Source counts

[giuspugl]

First implementation of integrating radio source spectra.
given a flux cut and freq, it looks for the closest number count model from Lagache et al 2019 (which is an updated version of Tucci et al. 2011) and returns the level of poissonian noise level due to undetected sources .

[dpole]

Hey Peppe, thanks for the PR. From what I understand your routine computes only the autospectra, storing on disk the total contribution from all the galaxy populations at some reference frequencies. Is that correct?

[giuspugl]

Yes the estimation is only for autospectra.To estimate the level of Poissonian noise you need to integrate as a function of flux ,S^2 *n(s) with n encodes the differential number counts coming from the updated Tucci model, published in Lagache et al. 2019 . The files in data folder encode the number counts

[dpole]

Do you think that it is possible to come up with a recipe for the cross spectra? I presume you need number counts and an SED for each population. BTW, we already have tools for several types of scaling as well as tools for summing several independent populations.

[giuspugl]

That's a good point ! I think we can follow the prescriptions presented in Reichardt et al 2012 (eq. 18) . Do you want me to keep committing in this PR to include Xspectra too ?

[dpole]

Sure, you can keep committing to this PR.

[codecov-commenter]

Codecov Report

❌ Patch coverage is 18.18182% with 54 lines in your changes missing coverage. Please review.
✅ Project coverage is 55.49%. Comparing base (46a6151) to head (edc014a).

Files with missing lines	Patch %	Lines
fgspectra/frequency.py	18.18%	54 Missing ⚠️

❗ There is a different number of reports uploaded between BASE (46a6151) and HEAD (edc014a). Click for more details.

HEAD has 1 upload less than BASE

Flag	BASE (46a6151)	HEAD (edc014a)
	2	1

Additional details and impacted files

@@             Coverage Diff             @@
##           fluxcut      #34      +/-   ##
===========================================
- Coverage    63.60%   55.49%   -8.12%     
===========================================
  Files            4        4              
  Lines          294      355      +61     
===========================================
+ Hits           187      197      +10     
- Misses         107      158      +51     

☔ View full report in Codecov by Harness.
📢 Have feedback on the report? Share it here.

🚀 New features to boost your workflow:

• ❄️ Test Analytics: Detect flaky tests, report on failures, and find test suite problems.

[giuspugl]

This is a very first implementation of Cross Spectra. It relies on the definition of Cross spectra implemented in Reichardt et al 2012 (see eq.18 therein). A major difference is in the estimates of spectral indices, in this case. I used the spectral indices for Radio sources estimated from Planck Catalogue from 30 to 353 GHz.