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

Comparing UserCode/cbrown/Development/Plotting/Modules/ActiveSamples.C (file contents):
Revision 1.20 by buchmann, Tue Jun 19 16:26:46 2012 UTC vs.
Revision 1.25 by buchmann, Tue Jul 31 16:12:27 2012 UTC

# Line 96 | Line 96 | http://cms.cern.ch/iCMS/prep/requestmana
96      }
97    
98      float ZJetsCrossSection         = 3048.0; //NNLO----------------------------------
99 <    if(PlottingSetup::is2012) ZJetsCrossSection=3503.71;
99 >    if (PlottingSetup::is2012) ZJetsCrossSection=3503.71;
100      // at 8 TeV: 3503.71 (https://twiki.cern.ch/twiki/bin/view/CMS/Zpt8TeV)
101      float LowZJetsCrossSection      = 11908; //https://twiki.cern.ch/twiki/bin/view/CMS/TWikiTop2011DataMCTrig
102 +    if (PlottingSetup::is2012) LowZJetsCrossSection = 11050*0.069*1.15; // LO xs * filter eff * k-factor (same as ZJets)
103      
104      float TTbarCrossSection         = 165.0;//(NLO) ---- 165.0; // approx. NNLO-----
105      if(PlottingSetup::is2012) TTbarCrossSection=225.197;
# Line 111 | Line 112 | http://cms.cern.ch/iCMS/prep/requestmana
112      float SingleTopTCrossSection     = 64.57; // NLO;--------------------------------
113      float SingleTopUCrossSection     = 15.74; // NLO;----------------------------------
114      float VVJetsCrossSection        = 4.8; // LO;-------------------------------------
115 < /*    float LM0CrossSection           =         38.93 *1.41 ;
116 <    float LM1CrossSection           = 4.888 *1.34 ;
116 <    float LM2CrossSection           =0.6027 *1.33 ;
117 <    float LM3CrossSection           = 3.438*1.40 ; // LO*k*/
115 >    float LM3CrossSection           = 0.0; // 2011: don't know the x-section
116 >    if(PlottingSetup::is2012) LM3CrossSection=5.21678; //(only LO, no k factor)
117      float LM4CrossSection           = 1.879*1.35; // k*LO
118      if(PlottingSetup::is2012) LM4CrossSection=2.88199; //(only LO, no k factor)
120 –    /*
121 –    float LM5CrossSection           = 0.473*1.34; // k*LO*/
119      float LM6CrossSection           = 0.3104*1.30; // k*LO
120      if(PlottingSetup::is2012) LM6CrossSection=0.501773;//(only LO, no k factor)
124 –    /*
125 –    float LM7CrossSection           = 1.209*1.11; // k*LO*/
121      float LM8CrossSection           = 0.7300*1.41; // k*LO
122      if(PlottingSetup::is2012) LM8CrossSection=1.14;
128 – /*    float LM9CrossSection           = 7.134*1.48; // k*LO
129 –    float LM10CrossSection           = 0.04778*1.22 ; // k*LO
130 –    float LM11CrossSection           = 0.8236 *1.35 ; // k*LO
131 –    float LM12CrossSection           = 4.414 *1.34 ; // k*LO
132 –    float LM13CrossSection           = 6.899*1.42 ; // k*LO*/
133 –    //the following numbers are from the MadGraphStandardModel210Summary as linked on the GeneratorProduction2011 page    
134 –    //    float QCD50to100CrossSection=30000000; // not used
135 –    //float QCD100to250CrossSection=7000000;
136 –    //float QCD250to500CrossSection=171000;
137 –    //float QCD500to1000CrossSection=5200;
138 –    //float QCD1000toInfCrossSection=83.0;
139 –    
140 –    //    float totEventsQCD50to100=207418.0; // not used
141 –    //long totEventsQCD100to250   = 638792;
142 –    //long totEventsQCD250to500   = 344454;
143 –    //long totEventsQCD500to1000  = (long)10e10;  //unknown but now suppressed. (no event passes anyway)
144 –    //long totEventsQCD1000toInf  = 163185;
123      
124      long totEventsZjets         = 32005720;//summer11
125      long totEventslowZjets      = 10272470;
126  
149 – //    long totEventsZjetsD6T      = 2543706;//summer11
150 – //    long totEventsZjetsD6Tlow   = 181799;//summer11
151 –
152 –    //spring:2313911.0;//2743142.0; // DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola
127      long totEventsTTbar         =  3701947;//summer11
154 –    //1161621.0;//Spring//1144028.0; // PabloV13/TTJets_TuneZ2_7TeV-madgraph-tauola
128      long totEventsWJets         = 56674902;//summer11
156 –    //15010237.0;//spring//14350756.0; /// WJetsToLNu_TuneZ2_7TeV-madgraph-tauola.root
129      long totEventsZnunu         =  2106977;//2167964.0; //ZinvisibleJets_7TeV-madgraph.root
130      long totEventsVVJets        =   959076;//509072.0;
131      long totEventsSingleTopS    =   493868;//489472.0;
132      long totEventsSingleTopT    =   475460;//477610.0;
133      long totEventsSingleTopU    =   489417;//477599.0;
162 – /*    long totEventsLM0           =   1234567;//bad guess, will be replaced.
163 –    long totEventsLM1           =   1234567;//bad guess, will be replaced.
164 –    long totEventsLM2           =   1234567;//bad guess, will be replaced.*/
165 –    long totEventsLM3           =   220000;//bad guess, will be replaced.
166 – /*    long totEventsLM4           =   218380;
167 –    long totEventsLM5           =   220000;
168 –    long totEventsLM6           =   220000;
169 –    long totEventsLM7           =   217800;*/
170 –    long totEventsLM8           =   217800;/*
171 –    long totEventsLM9           =   220000;
172 –    long totEventsLM10          =   1234567;
173 –    long totEventsLM11          =   216999;
174 –    long totEventsLM12          =   217400;
175 –    long totEventsLM13          =   1234567;*/
176 –
134      
135      Int_t nice_blue  = TColor::GetColor("#2E9AFE");
136      Int_t nice_green = TColor::GetColor("#81f781");
# Line 225 | Line 182 | http://cms.cern.ch/iCMS/prep/requestmana
182        
183   //      allsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_200__M12_440.root","LM8",10000,0.132,false,true,8,lm_color);
184        
185 +      
186 +      PlottingSetup::emuncertOFFPEAK=0.1;
187 +      write_info(__FUNCTION__,"Have set emu uncertainty to 10% which is the 2011 value - for 2012 this is 25% at the moment due to the low stats ttbar sample");
188      } else {
189        
190        write_info(__FUNCTION__,"Setting up 2012 configuration - once you've stopped looking at 2011 data please make this the default!");
191        PlottingSetup::basicqualitycut=TCut("(pfJetGoodNum40>=2&&pfJetGoodID[0]!=0)&&(pfJetGoodNum40>=2&&pfJetGoodID[1]!=0)"&&PlottingSetup::basiccut);//don't use this for the "essential cut", because we want to plot nJets as well as mll in the inclusive case; we thus use it as an addition nJets cut.
192        PlottingSetup::cutnJets=TCut("pfJetGoodNum40>=3"&&PlottingSetup::basicqualitycut);
193 +
194 +      // For Kostas' control region
195 + //       PlottingSetup::basicqualitycut=TCut("(pfJetGoodNum40<=2&&pfJetGoodID[0]!=0)&&(pfJetGoodNum40>=2&&pfJetGoodID[1]!=0)"&&PlottingSetup::basiccut);//don't use this for the "essential cut", because we want to plot nJets as well as mll in the inclusive case; we thus use it as an addition nJets cut.
196 + //       PlottingSetup::cutnJets=TCut("pfJetGoodNum40<3"&&PlottingSetup::basicqualitycut);
197 +
198        PlottingSetup::luminosity=PlottingSetup::luminosity2012;
199        PlottingSetup::lumiuncert=PlottingSetup::lumiuncert2012;
200        
201 <      // cJZB response correction (705 /pb) :
202 <      PlottingSetup::jzbvariabledata="(jzb[1]+0.0656534*pt)";
203 <      PlottingSetup::jzbvariablemc="(jzb[1]+0.0322238*pt)";
201 >      // Values for 5.1 /fb (2012, cJZB)
202 >      PlottingSetup::jzbvariabledata="(jzb[1]+0.064814*pt)";
203 >      PlottingSetup::jzbvariablemc="(jzb[1]+0.036854*pt)";
204        
205 +        
206 +        
207        if(!PlottingSetup::RestrictToMassPeak) {
208 <        PlottingSetup::jzbvariabledata="(((id1==0&&id1==id2)*(jzb[1]+0.0698233*pt))+((id1==1&&id1==id2)*(jzb[1]+0.0572628*pt))+((id1!=id2)*(jzb[1]+0.0629109*pt)))";
209 <        PlottingSetup::jzbvariablemc="(((id1==0&&id1==id2)*(jzb[1]+0.0409392*pt))+((id1==1&&id1==id2)*(jzb[1]+0.0327419*pt))+((id1!=id2)*(jzb[1]+0.0361354*pt)))";
208 >        // Values for 5.1 /fb (2012, iJZB)
209 >        PlottingSetup::jzbvariabledata="(jzb[1]+0.059979*pt)";
210 >        PlottingSetup::jzbvariablemc="(jzb[1]+0.034665*pt)";
211  
212        }
213 +      
214 + //       //  Note: Using same response for bJZB as for cJZB
215 + //       if(PlottingSetup::DoBTag) {
216 + //         // Values for 5.1 /fb (2012, bJZB)
217 + //      PlottingSetup::jzbvariabledata="(jzb[1]+0.092*pt)";
218 + //         PlottingSetup::jzbvariablemc="(jzb[1]+0.040*pt)";
219 + //       }
220 +      
221          
222 +      SingleTopSCrossSection=5.55;//Top mass: m(top)=173GeV, Default scales: muF=muR=mt, PDF: MSTW2008 NNLO
223 +      SingleTopTCrossSection=87.1;//Top mass: m(top)=173GeV, Default scales: muF=muR=mt, PDF: MSTW2008 NNLO
224 +      SingleTopUCrossSection=22.2;//Top mass: m(top)=173GeV, Default scales: muF=muR=mt, PDF: MSTW2008 NNLO
225 +      
226 +      allsamples.AddSample("/shome/buchmann/ntuples/Data/Data_5100ipb.root","Data",1,1,true,false,0,kBlack);
227 +
228 +      allsamples.AddSample("/shome/fronga/ntuples/8TeV/DYJetsToLLNoTau_M-50_TuneZ2Star_8TeV-madgraph-tarball-Summer12-PU_S7_START52_V9-v2.root","Z+Jets",-30346766,ZJetsCrossSection,false,false,8,dy_color);
229 +      allsamples.AddSample("/shome/fronga/ntuples/8TeV/DYJetsToLLNoTau_M-10To50filter_8TeV-madgraph-Summer12-PU_S7_START52_V9-v1.root","Z+Jets",-6955271,LowZJetsCrossSection,false,false,8,dy_color);
230 +      allsamples.AddSample("/shome/fronga/ntuples/8TeV/DYJetsToTauTau_M-50_TuneZ2Star_8TeV-madgraph-tarball-Summer12-PU_S7_START52_V9-v2.root","Z+Jets (#tau#tau)",-30346766,ZJetsCrossSection,false,false,7,kRed-4);
231 +      allsamples.AddSample("/shome/fronga/ntuples/8TeV/DYJetsToTauTau_M-10To50filter_8TeV-madgraph-Summer12-PU_S7_START52_V9-v1.root","Z+Jets (#tau#tau)",-6955271,LowZJetsCrossSection,false,false,7,kRed-4);
232 +
233 + //       allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball-Summer12-PU_S7_START52_V9-v2.root","Z+Jets",0,ZJetsCrossSection,false,false,7,dy_color);
234 + //       allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/DYJetsToLL_M-10To50filter_8TeV-madgraph-Summer12-PU_S7_START52_V9-v1.root","Z+Jets",0,LowZJetsCrossSection,false,false,7,dy_color);
235  
236 < //      allsamples.AddSample("/shome/buchmann/ntuples/Data/Data_920ipb_UpdatedTriggers_updatedJetID_v2.root","Data",1,1,true,false,0,kBlack);
237 <      allsamples.AddSample("/shome/buchmann/ntuples/Data/Data_2420ipb_correctJSON.root","Data",1,1,true,false,0,kBlack);
236 >      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/TTJets_TuneZ2star_8TeV-madgraph-tauola-Summer12-PU_S7_START52_V5-v1.root","t#bar{t}",0,TTbarCrossSection,false,false,1,ttbar_color);//summer11
237 >      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/WJetsToLNu_TuneZ2Star_8TeV-madgraph-tarball-Summer12-PU_S7_START52_V9-v1.root","W+Jets",0,36257.2,false,false,3,wjets_color);//summer
238 >
239 >      // SINGLE-TOP
240 >      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/T_s-channel_TuneZ2star_8TeV-powheg-tauola-Summer12-PU_S7_START52_V9-v1.root","Single-top",0,1.76,false,false,2,singletop_color);
241 >      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/Tbar_s-channel_TuneZ2star_8TeV-powheg-tauola-Summer12-PU_S7_START52_V9-v1.root","Single-top",0,3.79,false,false,2,singletop_color);
242 >      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/Tbar_t-channel_TuneZ2star_8TeV-powheg-tauola-Summer12-PU_S7_START52_V9-v1.root","Single-top",0,56.4,false,false,2,singletop_color);
243 >      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/Tbar_tW-channel-DR_TuneZ2star_8TeV-powheg-tauola-Summer12-PU_S7_START52_V9-v1.root","Single-top",0,SingleTopUCrossSection,false,false,2,singletop_color);
244        
245 <      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball-Summer12-PU_S7_START52_V9-v2.root","Z+Jets",0,ZJetsCrossSection,false,false,7,dy_color);
246 <      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/TTJets_TuneZ2star_8TeV-madgraph-tauola-Summer12-PU_S7_START52_V9-v1.root","t#bar{t}",0,TTbarCrossSection,false,false,1,ttbar_color);
247 < //      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/LM4_leo____leo-PYTHIA6_SUSY_LM4_sftsht_8TeV_GEN_FASTSIM-2012_Startup_inTimeOnly-V3-84.root","LM4",0,LM4CrossSection,false,true,1,lm_color);
248 < //        allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/leo-PYTHIA6_SUSY_LM6_sftsht_8TeV_GEN_FASTSIM-2012_Startup_inTimeOnly.root","LM6",0,LM6CrossSection,false,true,1,lm_color);
249 <        
250 < /*      if(PlottingSetup::DoBTag) {
251 <        PlottingSetup::jzbvariabledata="(jzb[1]+0.0924584*pt)";
252 <        PlottingSetup::jzbvariablemc="(jzb[1]+0.0477569*pt)";
253 <      }*/
245 >      // DIBOSONS: LO with a k-factor
246 >      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/ZZJetsTo4L_TuneZ2star_8TeV-madgraph-tauola-Summer12-PU_S7_START52_V9-v3.root","Dibosons",0,0.181,false,false,4,diboson_color);
247 >      allsamples.AddSample("/shome/fronga/ntuples/8TeV/ZZJetsTo2L2Nu_TuneZ2star_8TeV-madgraph-tauola-Summer12-PU_S7_START52_V9-v3.root","Dibosons",0,0.3,false,false,4,diboson_color);
248 >      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/WZJetsTo3LNu_TuneZ2_8TeV-madgraph-tauola-Summer12-PU_S7_START52_V9-v2.root","Dibosons",0,1.06,false,false,4,diboson_color);
249 >      allsamples.AddSample("/shome/fronga/ntuples/8TeV/WWJetsTo2L2Nu_TuneZ2star_8TeV-madgraph-tauola-Summer12-PU_S7_START52_V9-v1.root","Dibosons",0,5.0,false,false,4,diboson_color);
250 >
251 >      //ttV
252 >      //allsamples.AddSample("/shome/fronga/ntuples/8TeV/TTZJets_8TeV-madgraph_v2-Summer12-PU_S7_START52_V9-v1.root","ttV",0,0.172,false,false,5,kGreen);
253 >      //allsamples.AddSample("/shome/fronga/ntuples/8TeV/TTGJets_8TeV-madgraph-Summer12-PU_S7_START52_V9-v1.root",   "ttV",0,1.444,false,false,5,kGreen);
254 >      //allsamples.AddSample("/shome/fronga/ntuples/8TeV/TTWJets_8TeV-madgraph-Summer12-PU_S7_START52_V9-v1.root",   "ttV",0,0.2149,false,false,5,kGreen);
255 >      
256 >      
257 >      //  allsamples.AddSample(prefix+"/buchmann/ntuples/MC/JZBplus/WGstarToLNu2E_TuneZ2_7TeV-madgraph-tauola__Summer11-PU_S4_START42_V11-v1__AODSIM.root","W+#gamma",266882,5.546 ,false,false,5,kGreen);
258 >      
259 >      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/leo-PYTHIA6_SUSY_LM6_sftsht_8TeV_GEN_FASTSIM-2012_Startup_inTimeOnly.root","LM6",0,LM4CrossSection,false,true,1,lm_color);
260          
261      }
261 –    if(!PlottingSetup::RestrictToMassPeak) {
262 –        if(!PlottingSetup::is2012) signalsamples.AddSample("/shome/buchmann//ntuples/MC/NewSelectionForPaper/NewlyReweighted/DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola_Summer11___NTupleV020300.root","LM8",totEventsLM8,LM8CrossSection,false,true,8,lm_color);
263 –    }
264 –
262  
263      if(PlottingSetup::is2012) {
264 < //        signalsamples.AddSample("/shome/buchmann/ntuples/MC8tev/LM4_leo____leo-PYTHIA6_SUSY_LM4_sftsht_8TeV_GEN_FASTSIM-2012_Startup_inTimeOnly-V3-84.root","LM4",0,LM4CrossSection,false,true,1,lm_color);
265 < //        signalsamples.AddSample("/shome/buchmann/ntuples/MC8tev/leo-PYTHIA6_SUSY_LM6_sftsht_8TeV_GEN_FASTSIM-2012_Startup_inTimeOnly.root","LM6",0,LM6CrossSection,false,true,2,lm_color);
264 >        signalsamples.AddSample("/shome/buchmann/ntuples/MC8tev/trash/LM4_leo____leo-PYTHIA6_SUSY_LM4_sftsht_8TeV_GEN_FASTSIM-2012_Startup_inTimeOnly-V3-84.root","LM4",0,LM4CrossSection,false,true,1,lm_color);
265 >        signalsamples.AddSample("/shome/buchmann/ntuples/MC8tev/leo-PYTHIA6_SUSY_LM6_sftsht_8TeV_GEN_FASTSIM-2012_Startup_inTimeOnly.root","LM6",0,LM6CrossSection,false,true,2,lm_color);
266 >        signalsamples.AddSample("/shome/buchmann/ntuples/MC8tev/PYTHIA6_SUSY_LM3_sftsht_8TeV_GEN_FASTSIM-2012_Startup_inTimeOnly-V3-leo_AllWeightOne.root","LM3",0,LM3CrossSection,false,true,3,lm_color);
267   //        signalsamples.AddSample("/shome/buchmann/ntuples/MC8tev/leo-PYTHIA6_SUSY_LM8_sftsht_8TeV_GEN_FASTSIM-2012_Startup_inTimeOnly.root","LM8",0,LM6CrossSection,false,true,3,lm_color);
268 <        cutWeight=TCut("PUweight");
268 >        cutWeight=TCut("(weight*(weight<1000)*(is_data+(!is_data)*((id1==id2)*0.95+(id1!=id2)*0.94)))");
269 >        //FR cutWeight=TCut("(weight*(weight<1000)*(is_data+(!is_data)*((id1==id2&&id1==0)*0.95+(id1==id2&&id1==1)*0.88+(id1!=id2)*0.92)))");
270 >        //cutWeight=TCut("(weight*(weight<1000)*((id1==id2&&id1==0)*0.95+(id1==id2&&id1==1)*0.88+(id1!=id2)*0.92))");
271      }
272  
273  
# Line 277 | Line 277 | http://cms.cern.ch/iCMS/prep/requestmana
277        systsamples.AddSample(prefix+"buchmann//ntuples/MC/NewSelectionForPaper/NewlyReweighted/DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola_Summer11___NTupleV020300.root","Z+Jets",32005720,ZJetsCrossSection,false,false,7,dy_color);
278        systsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/HighStatSamples/TTJets_Fall11-PU_S6_START42_V14B.root","t#bar{t}",58370118,TTbarCrossSection,false,false,1,ttbar_color);//summer11
279      } else {
280 +      
281        systsamples.AddSample("/shome/buchmann/ntuples/MC8tev/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball-Summer12-PU_S7_START52_V9-v2.root","Z+Jets",0,ZJetsCrossSection,false,false,7,dy_color);
282 <      systsamples.AddSample("/shome/buchmann/ntuples/MC8tev/TTJets_TuneZ2star_8TeV-madgraph-tauola-Summer12-PU_S7_START52_V9-v1.root","t#bar{t}",0,TTbarCrossSection,false,false,1,ttbar_color);
282 >      systsamples.AddSample("/shome/buchmann/ntuples/MC8tev/DYJetsToLL_M-10To50filter_8TeV-madgraph-Summer12-PU_S7_START52_V9-v1.root","Z+Jets",0,LowZJetsCrossSection,false,false,7,dy_color);
283 >      systsamples.AddSample("/shome/buchmann/ntuples/MC8tev/TTJets_TuneZ2star_8TeV-madgraph-tauola-Summer12-PU_S7_START52_V5-v1.root","t#bar{t}",0,TTbarCrossSection,false,false,1,ttbar_color);//summer11
284 >      
285      }
286  
287  
# Line 299 | Line 302 | http://cms.cern.ch/iCMS/prep/requestmana
302  
303      
304      
305 <    if(showList) allsamples.ListSamples();
305 >    if(showList) {
306 >      allsamples.ListSamples();
307 >      signalsamples.ListSamples();
308 >      systsamples.ListSamples();
309 >    }
310  
311   }

Diff Legend

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