Dear experts,
I have a question about the calculation about upperlimit by using Fedlman method. I want to know that why if I choose the oberserve and background value are same and the value is more than 40, the upperlimits given are same. And for 40, the upperlimit value is 0 .
Hi,
The class is designed with some default parameters to be used for small values of n and n-background, where it is typically used.
In your case, for example for nobs=40
and nbkg=40
, you need to set larger values of the scan muMax
. The default is 50.
You need to set a value of muMax
much larger than nobs
, at least larger than nobs+nbkg
. You might need to increase the muStep
(default is 0.005) if it takes too long to compute.
TFeldmanCousins fc;
fc.SetCL(0.9);
fc.SetMuMax(100);
fc.CalculateUpperLimit(40,40);
will give you a value around 11.5
Cheers
Lorenzo
Dear moneta,
Thank you for your suggestion. So if the values of background and obs range from 1 to 100, can I choose a maximum value like MuMax (200) ? For small values of background and obs like nobs +nbkg is 10, if I setMuMax is 100, is it ok?
Hi,
I think muMax should be at least large as nobs+nbkg
. You should also note that for large values the Poisson is like a normal distribution and you can just use the limit computed from a normal distribution.
For example for a 90% CL ,
double upper_limit = nobs - b + sqrt(nobs) * ROOT::Math::normal_quantile_c(0.05,1);
Lorenzo
Oh, I see. That’s very useful! Thank you for your help!
This topic was automatically closed 14 days after the last reply. New replies are no longer allowed.