Changeset 19910 for trunk/FACT++

Timestamp:

12/15/19 18:54:47 (5 years ago)

Author:

tbretz

Message:

Forgot to adapt the upper limit calculation to the fact that the BG is already the background scaled to only one wobble position.

File:

• trunk/FACT++/src/spectrum.cc (modified) (2 diffs)

trunk/FACT++/src/spectrum.cc

@@ -735 +735 @@
         cout << "Histograms will be written to " << out << ".hist.root\n";
     }
+    if (verbose>0)
+        cout << "\nCalculating upper limits for a confidence interval of " << confidence << endl;
 #endif
 
@@ -2088 +2090 @@
 
                 if (feldman)
-                    feldman_ul[bin] = fc.CalculateUpperLimit(dat_sig, dat_bg/5)/width/areatime/eff;
-
-                rolke.SetPoissonBkgKnownEff(dat_sig,  dat_bg, 5, eff);
+                    feldman_ul[bin] = fc.CalculateUpperLimit(dat_sig, dat_bg)/width/areatime/eff;
+
+                rolke.SetPoissonBkgKnownEff(dat_sig,  dat_bg*5, 5, eff);
                 rolke_ll[bin]  = rolke.GetLowerLimit()/width/areatime;
                 rolke_ul[bin]  = rolke.GetUpperLimit()/width/areatime;
 
-                rolke.SetPoissonBkgKnownEff(dat_isig, dat_ibg, 5, ieff);
+                rolke.SetPoissonBkgKnownEff(dat_isig, dat_ibg*5, 5, ieff);
                 rolke_int[bin] = rolke.GetUpperLimit()/areatime;
 #endif