- Timestamp:
- 12/15/19 18:54:47 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/FACT++/src/spectrum.cc
r19908 r19910 735 735 cout << "Histograms will be written to " << out << ".hist.root\n"; 736 736 } 737 if (verbose>0) 738 cout << "\nCalculating upper limits for a confidence interval of " << confidence << endl; 737 739 #endif 738 740 … … 2088 2090 2089 2091 if (feldman) 2090 feldman_ul[bin] = fc.CalculateUpperLimit(dat_sig, dat_bg /5)/width/areatime/eff;2091 2092 rolke.SetPoissonBkgKnownEff(dat_sig, dat_bg , 5, eff);2092 feldman_ul[bin] = fc.CalculateUpperLimit(dat_sig, dat_bg)/width/areatime/eff; 2093 2094 rolke.SetPoissonBkgKnownEff(dat_sig, dat_bg*5, 5, eff); 2093 2095 rolke_ll[bin] = rolke.GetLowerLimit()/width/areatime; 2094 2096 rolke_ul[bin] = rolke.GetUpperLimit()/width/areatime; 2095 2097 2096 rolke.SetPoissonBkgKnownEff(dat_isig, dat_ibg , 5, ieff);2098 rolke.SetPoissonBkgKnownEff(dat_isig, dat_ibg*5, 5, ieff); 2097 2099 rolke_int[bin] = rolke.GetUpperLimit()/areatime; 2098 2100 #endif
Note:
See TracChangeset
for help on using the changeset viewer.