cositools / cosipy

The COSI high-level data analysis tools
Apache License 2.0
3 stars 16 forks source link

Revisit units of extended source fit tutorial #194

Open israelmcmc opened 6 days ago

israelmcmc commented 6 days ago

https://github.com/cositools/cosipy/blob/develop/docs/tutorials/spectral_fits/extended_source_fit/diffuse_511_spectral_fit.ipynb

The tutorial says that the units of the K normalization are differential: K = [0.00046, 0.0011, 0.0027, 4.5e-3]/u.cm/u.cm/u.s/u.keV

However, the code normalizes the spectrum based on the integral flux (assuming a 1 keV bin size) https://github.com/cositools/cosipy/blob/65706f739ea90e0d61103599908401e8d8f923d7/cosipy/threeml/custom_functions.py#L177

In addition, on #191 we're changing the definition of K to be unitless. For now, we changed the file in order to get the same numerical output, but we still need to fix this.