Model for the power from radio sources as a function of the flux cut

[dpole]

We have received a request for this model, but details are not clear.
To potential users, please express any need that could be not obvious (or feel free to mark down also the obvious ones).
To radio sources experts, please provide information about the functional form (and, of course, the parameters) that we should use in the implementation
Thanks!

[ajvanengelen]

I have a routine that does this for the Tucci et al 2010 C2EX model, which is the model used by the ACT team for this purpose. I am happy to write it into the fgspectra library. It would be helpful for someone who is already working with this package to write a function prototype (i.e. just the definition line for the routine) so that it meshes with the rest of the software in the package,

[dpole]

Great, thanks. We are finalizing the API in these days. Once we are done, I'll prepare a prototype for you.
In the meanwhile people should feel free to provide further inputs.

[ajvanengelen]

Sounds great! Like I said in an earlier email, what I have boils down to maybe three or four lines of code, so the work will really be in just setting up the structure to put it into the right place.

[erminiacalabrese]

Hi both, sorry for being late in responding. People are requesting the ability of predicting radio power given a specific flux cut. Alex you have code for it and in terms of interface with fgspectra this should simply predict the amplitude of the radio sources which would enter as one of the parameters of the full foreground model.

[dpole]

Hi, we have converged on the API (and also the internal structure). We can now prepare a template for this component. In order to do that, we need to know something about the model. In particular, it is important to understand if you can factorize the multipole dependence and the frequency dependence.
@ajvanengelen or someone else, can you give us some information in this respect?

[erminiacalabrese]

Yes, these sources are just shot noise, so a flat Cl (ell(ell+1)/2piCl in Dl) and scale in frequency with a beta=-0.5 .

[dpole]

So the only variable is the Cl value as a function of the flux cut?

[erminiacalabrese]

The variable is the amplitude of Cl at ell=3000 (which can be predicted for a given flux cut) but varies as free amplitude for the CMB likelihood. The frequency index can be fixed for now.

[ajvanengelen]

Hi Erminia, For SO, will the point source mask be the same at all frequencies, and will it be determined from the sources found at a single frequency? Alex

…

On 2019-06-13 11:37 a.m., erminiacalabrese wrote: The variable is the amplitude of Cl at ell=3000 (which can be predicted for a given flux cut) but varies as free amplitude for the CMB likelihood. The frequency index can be fixed for now. — You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub <#21?email_source=notifications&email_token=ADYBOHFUP3XAQTP42NYFYFLP2JSS5A5CNFSM4HVHJM2KYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGODXUDENA#issuecomment-501756468>, or mute the thread <https://github.com/notifications/unsubscribe-auth/ADYBOHF7F5MIB6OVZQBYGFLP2JSS5ANCNFSM4HVHJM2A>.

[erminiacalabrese]

I don't know yet. We could think of implementing some automatic calculation going via flux cut but for now I'd keep it simple and give people the option of doing predictions but also varying an amplitude.

[dpole]

In this case we don't need any special tool. We just have to combine a power-law SED with a white cl. We only need the latter, I can take care of it.

[ajvanengelen]

Yes, you just have to integrate over the number of counts at each flux level to get the amplitude of the white C_l. Alex

…

On 2019-06-13 11:51 a.m., Davide Poletti wrote: In this case we don't need any special tool. We just have to combine a power-law SED with a white cl. We only need the latter, I can take care of it. — You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub <#21?email_source=notifications&email_token=ADYBOHG3K3M2LMU4YUXQKODP2JUJDA5CNFSM4HVHJM2KYY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGODXUESPY#issuecomment-501762367>, or mute the thread <https://github.com/notifications/unsubscribe-auth/ADYBOHGAY4OEOOOW3LIGRHLP2JUJDANCNFSM4HVHJM2A>.

[erminiacalabrese]

Exactly, and that can then be used for dusty sources and radio sources separately.

[dpole]

Ok, so I think we should implement two separate models. One that takes the amplitude of the white spectrum as a free parameter and the another that takes the flux cut as free parameter. Right?
I can easily take care of the former and prepare the skeleton of the latter

[jcolinhill]

Just to check, this should be a power-law in flux units, right? (c.f. Sec. 2.5, Eq. 13 of Dunkley et al 2013) It seems like the power-law implementation in fgspectra currently assumes a power-law in RJ temperature units.