bla

analysis/plot/CreateHisto.C

@@ -21,9 +21,6 @@ void histo(double hours = 168, int binWidth = 100)   // 168 hours = 7 days
 
     double sf, norm;    
     
-    // Veto cut
-    TCut veto("M1<=1 && E1<140+0.21*E0"");
-
     // BACKGROUND    
     TChain * backgrChain = new TChain("data");    
     chain->Add(datadir+"*0192*");    
@@ -37,6 +34,7 @@ void histo(double hours = 168, int binWidth = 100)   // 168 hours = 7 days
     TH1F * hBackgr = new TH1F("hBackgr", "hBackgr", bins, emin, emax);
     TH1F * hBackgrVeto = new TH1F("hBackgrVeto", "hBackgrVeto", bins, emin, emax);
     backgrChain -> Draw("EDepSignal >> hBackgr", "Z0 != 10");
+    TCut veto("M1 <= 1 && E1 < 140 + 0.21*E0");
     backgrChain -> Draw("EDepSignal >> hBackgrVeto", "Z0 != 10" && veto);
     sf = hours/backgrHours;
     hBackgr->Scale(sf);