search

zblz / naima

Derivation of non-thermal particle distributions through MCMC spectral fitting
http://naima.readthedocs.io
BSD 3-Clause "New" or "Revised" License
45 stars 54 forks source link

How to deal with the upper limit data? #153

Closed 55hyan closed 7 years ago

55hyan commented 7 years ago

How can we deal with the upper limit datas when we use naima for fitting? (Such as RX J1713.7-3946, there are several upper limit datas in the Fermi observation, when we do the fitting, how should we consider?).

zblz commented 7 years ago

naima deals with them by adding a certain value to the log-likelihood whenever the UL is violated by the model, otherwise they have no effect on the likelihood. This value depends on the confidence level of the UL.

I'm not sure what you question is. Do you mean whether you should include them in the spectrum used for the fitting? See the data format docs page for more info on how to do that.

55hyan commented 7 years ago

Yes, I don't know how to include the upper limits in the data format, but I have got it now. Thanks.

55hyan commented 7 years ago

I am also interested in the secondary electron spectrum from proton-proton interaction, did you have any update for this part?

zblz commented 7 years ago

Secondary spectrum is not yet computed by naima, there is an issue to track progress on that: #70.

I'll close this issue as the UL issue has been resolved.