ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/cbrown/AnalysisFramework/Plotting/Modules/Systematics.C
(Generate patch)

Comparing UserCode/cbrown/AnalysisFramework/Plotting/Modules/Systematics.C (file contents):
Revision 1.2 by buchmann, Tue Jul 19 07:04:41 2011 UTC vs.
Revision 1.3 by buchmann, Tue Jul 19 09:28:30 2011 UTC

# Line 114 | Line 114 | float pileup(TTree *events, string infor
114          jzbMax = myJzbMax;
115          
116          // Acceptance cuts
117 <        TCut kbase("abs(genMllSel-91.2)<20&&pfJetGoodNum>2&&genZPtSel>0&&abs(mll-91.2)<20&&((id1+1)*(id2+1)*ch1*ch2)!=-2");
117 >        TCut kbase("abs(genMllSel-91.2)<20&&genNjets>2&&genZPtSel>0&&abs(mll-91.2)<20&&((id1+1)*(id2+1)*ch1*ch2)!=-2");
118          
119          TH1F* hLM4 = plotEff(events,kbase,informalname);
120          hLM4->SetMinimum(0.);
# Line 141 | Line 141 | void MCefficiency(TTree *events,float &r
141          
142          char jzbSelStr[256]; sprintf(jzbSelStr,"%f",jzbSel);
143          // All acceptance cuts at gen. level
144 <        TCut kbase("abs(genMllSel-91.2)<20&&pfJetGoodNum>2&&genZPt>0&&genJZB"+geq_or_leq()+TString(jzbSelStr)+"&&genId1==-genId2");
144 >        TCut kbase("abs(genMllSel-91.2)<20&&genNjets>2&&genZPt>0&&genJZB"+geq_or_leq()+TString(jzbSelStr)+"&&genId1==-genId2");
145          // Corresponding reco. cuts
146          TCut ksel("abs(mll-91.2)<20&&id1==id2&&"+TString(mcjzb)+geq_or_leq()+TString(jzbSelStr));
147          
# Line 156 | Line 156 | void MCefficiency(TTree *events,float &r
156   }
157  
158   float JZBefficiency(TTree *events, string informalname) {
159 <        TCut kbase("abs(genMllSel-91.2)<20&&pfJetGoodNum>2&&genZPt>0&&abs(mll-91.2)<20&&((id1+1)*(id2+1)*ch1*ch2)!=-2");
159 >        TCut kbase("abs(genMllSel-91.2)<20&&genNjets>2&&genZPt>0&&abs(mll-91.2)<20&&((id1+1)*(id2+1)*ch1*ch2)!=-2");
160          TH1F* hLM4 = plotEff(events,kbase,informalname);
161          Int_t bin = hLM4->FindBin(jzbSel); // To get the error
162          std::cout << "  Efficiency at JZB==" << jzbSel  << std::endl;
# Line 204 | Line 204 | void JZBjetScale(TTree *events, float &j
204   // Effect of energy scale on JZB efficiency
205   void doJZBscale(TTree *events, float &down, float &up, float &syst, float systematic, string informalname) {
206          
207 <        TCut kbase("abs(genMllSel-91.2)<20&&genZPt>0&&pfJetGoodNum>2");
207 >        TCut kbase("abs(genMllSel-91.2)<20&&genZPt>0&&genNjets>2");
208          TCut ksel("abs(mll-91.2)<20&&((id1+1)*(id2+1)*ch1*ch2)!=-2");
209          
210          nBins =    50;
# Line 229 | Line 229 | void doJZBscale(TTree *events, float &do
229   void JZBresponse(TTree *events, bool isMET = kFALSE, Float_t myJzbMax = 200., Int_t nPeriods = 9 ) {
230          
231          jzbMin = 20;
232 <        TCut kbase("abs(genMllSel-91.2)<20&&genZPtSel>0&&pfJetGoodNum>2");
232 >        TCut kbase("abs(genMllSel-91.2)<20&&genZPtSel>0&&genNjets>2");
233          TCut ksel("abs(mll-91.2)<20&&((id1+1)*(id2+1)*ch1*ch2)!=-2");
234          
235          TProfile* hJzbResp = new TProfile("hJzbResp","JZB response  ; JZB true (GeV/c); JZB reco. / JZB true",

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines