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.21 by fronga, Thu Jun 21 11:37:29 2012 UTC vs.
Revision 1.26 by buchmann, Tue Aug 14 12:14:59 2012 UTC

# Line 24 | Line 24 | void define_SMS_sample(bool showList, sa
24   {
25      char hostname[1023];
26      gethostname(hostname,1023);
27 +    
28 +    SanityChecks();
29  
30      /*
31      
# Line 71 | Line 73 | void define_SMS_sample(bool showList, sa
73  
74   void define_samples(bool showList, samplecollection &allsamples, samplecollection &signalsamples, samplecollection &scansample, samplecollection &raresample, samplecollection &systsamples, samplecollection &qcdsamples)
75   {
76 +  SanityChecks();
77  
78   /*
79   Good sources for XS:
# Line 112 | Line 115 | http://cms.cern.ch/iCMS/prep/requestmana
115      float SingleTopTCrossSection     = 64.57; // NLO;--------------------------------
116      float SingleTopUCrossSection     = 15.74; // NLO;----------------------------------
117      float VVJetsCrossSection        = 4.8; // LO;-------------------------------------
118 < /*    float LM0CrossSection           =         38.93 *1.41 ;
119 <    float LM1CrossSection           = 4.888 *1.34 ;
117 <    float LM2CrossSection           =0.6027 *1.33 ;
118 <    float LM3CrossSection           = 3.438*1.40 ; // LO*k*/
118 >    float LM3CrossSection           = 0.0; // 2011: don't know the x-section
119 >    if(PlottingSetup::is2012) LM3CrossSection=5.21678; //(only LO, no k factor)
120      float LM4CrossSection           = 1.879*1.35; // k*LO
121      if(PlottingSetup::is2012) LM4CrossSection=2.88199; //(only LO, no k factor)
121 –    /*
122 –    float LM5CrossSection           = 0.473*1.34; // k*LO*/
122      float LM6CrossSection           = 0.3104*1.30; // k*LO
123      if(PlottingSetup::is2012) LM6CrossSection=0.501773;//(only LO, no k factor)
125 –    /*
126 –    float LM7CrossSection           = 1.209*1.11; // k*LO*/
124      float LM8CrossSection           = 0.7300*1.41; // k*LO
125      if(PlottingSetup::is2012) LM8CrossSection=1.14;
129 – /*    float LM9CrossSection           = 7.134*1.48; // k*LO
130 –    float LM10CrossSection           = 0.04778*1.22 ; // k*LO
131 –    float LM11CrossSection           = 0.8236 *1.35 ; // k*LO
132 –    float LM12CrossSection           = 4.414 *1.34 ; // k*LO
133 –    float LM13CrossSection           = 6.899*1.42 ; // k*LO*/
134 –    //the following numbers are from the MadGraphStandardModel210Summary as linked on the GeneratorProduction2011 page    
135 –    //    float QCD50to100CrossSection=30000000; // not used
136 –    //float QCD100to250CrossSection=7000000;
137 –    //float QCD250to500CrossSection=171000;
138 –    //float QCD500to1000CrossSection=5200;
139 –    //float QCD1000toInfCrossSection=83.0;
140 –    
141 –    //    float totEventsQCD50to100=207418.0; // not used
142 –    //long totEventsQCD100to250   = 638792;
143 –    //long totEventsQCD250to500   = 344454;
144 –    //long totEventsQCD500to1000  = (long)10e10;  //unknown but now suppressed. (no event passes anyway)
145 –    //long totEventsQCD1000toInf  = 163185;
126      
127      long totEventsZjets         = 32005720;//summer11
128      long totEventslowZjets      = 10272470;
129  
150 – //    long totEventsZjetsD6T      = 2543706;//summer11
151 – //    long totEventsZjetsD6Tlow   = 181799;//summer11
152 –
153 –    //spring:2313911.0;//2743142.0; // DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola
130      long totEventsTTbar         =  3701947;//summer11
155 –    //1161621.0;//Spring//1144028.0; // PabloV13/TTJets_TuneZ2_7TeV-madgraph-tauola
131      long totEventsWJets         = 56674902;//summer11
157 –    //15010237.0;//spring//14350756.0; /// WJetsToLNu_TuneZ2_7TeV-madgraph-tauola.root
132      long totEventsZnunu         =  2106977;//2167964.0; //ZinvisibleJets_7TeV-madgraph.root
133      long totEventsVVJets        =   959076;//509072.0;
134      long totEventsSingleTopS    =   493868;//489472.0;
135      long totEventsSingleTopT    =   475460;//477610.0;
136      long totEventsSingleTopU    =   489417;//477599.0;
163 – /*    long totEventsLM0           =   1234567;//bad guess, will be replaced.
164 –    long totEventsLM1           =   1234567;//bad guess, will be replaced.
165 –    long totEventsLM2           =   1234567;//bad guess, will be replaced.*/
166 –    long totEventsLM3           =   220000;//bad guess, will be replaced.
167 – /*    long totEventsLM4           =   218380;
168 –    long totEventsLM5           =   220000;
169 –    long totEventsLM6           =   220000;
170 –    long totEventsLM7           =   217800;*/
171 –    long totEventsLM8           =   217800;/*
172 –    long totEventsLM9           =   220000;
173 –    long totEventsLM10          =   1234567;
174 –    long totEventsLM11          =   216999;
175 –    long totEventsLM12          =   217400;
176 –    long totEventsLM13          =   1234567;*/
177 –
137      
138      Int_t nice_blue  = TColor::GetColor("#2E9AFE");
139      Int_t nice_green = TColor::GetColor("#81f781");
# Line 226 | Line 185 | http://cms.cern.ch/iCMS/prep/requestmana
185        
186   //      allsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_200__M12_440.root","LM8",10000,0.132,false,true,8,lm_color);
187        
188 +      
189 +      PlottingSetup::emuncertOFFPEAK=0.1;
190 +      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");
191      } else {
192        
193        write_info(__FUNCTION__,"Setting up 2012 configuration - once you've stopped looking at 2011 data please make this the default!");
194        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.
195        PlottingSetup::cutnJets=TCut("pfJetGoodNum40>=3"&&PlottingSetup::basicqualitycut);
196 +
197 +      // For Kostas' control region
198 + //       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.
199 + //       PlottingSetup::cutnJets=TCut("pfJetGoodNum40<3"&&PlottingSetup::basicqualitycut);
200 +
201        PlottingSetup::luminosity=PlottingSetup::luminosity2012;
202        PlottingSetup::lumiuncert=PlottingSetup::lumiuncert2012;
203        
204 <      // cJZB response correction (705 /pb) :
205 <      PlottingSetup::jzbvariabledata="(jzb[1]+0.0656534*pt)";
206 <      PlottingSetup::jzbvariablemc="(jzb[1]+0.0322238*pt)";
204 >      // Values for 5.1 /fb (2012, cJZB)
205 >      PlottingSetup::jzbvariabledata="(jzb[1]+0.064814*pt)";
206 >      PlottingSetup::jzbvariablemc="(jzb[1]+0.036854*pt)";
207        
208 +        
209 +        
210        if(!PlottingSetup::RestrictToMassPeak) {
211 <        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)))";
212 <        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)))";
211 >        // Values for 5.1 /fb (2012, iJZB)
212 >        PlottingSetup::jzbvariabledata="(jzb[1]+0.059979*pt)";
213 >        PlottingSetup::jzbvariablemc="(jzb[1]+0.034665*pt)";
214  
215        }
216 +      
217 + //       //  Note: Using same response for bJZB as for cJZB
218 + //       if(PlottingSetup::DoBTag) {
219 + //         // Values for 5.1 /fb (2012, bJZB)
220 + //      PlottingSetup::jzbvariabledata="(jzb[1]+0.092*pt)";
221 + //         PlottingSetup::jzbvariablemc="(jzb[1]+0.040*pt)";
222 + //       }
223 +      
224          
225 <
226 <      allsamples.AddSample("/shome/buchmann/ntuples/Data/Data_3553ifb.root","Data",1,1,true,false,0,kBlack);
225 >      SingleTopSCrossSection=5.55;//Top mass: m(top)=173GeV, Default scales: muF=muR=mt, PDF: MSTW2008 NNLO
226 >      SingleTopTCrossSection=87.1;//Top mass: m(top)=173GeV, Default scales: muF=muR=mt, PDF: MSTW2008 NNLO
227 >      SingleTopUCrossSection=22.2;//Top mass: m(top)=173GeV, Default scales: muF=muR=mt, PDF: MSTW2008 NNLO
228 >      
229 > //       allsamples.AddSample("/shome/buchmann/ntuples/Data/METdataset_5100ipb.root","Data",1,1,true,false,0,kBlack);
230 >      allsamples.AddSample("/shome/buchmann/ntuples/Data/Data_5100ipb_ZbReady.root","Data",1,1,true,false,0,kBlack);
231        
250 –      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball-Summer12-PU_S7_START52_V9-v2.root","Z+Jets",30346766,ZJetsCrossSection,false,false,7,dy_color);
251 –      allsamples.AddSample("/scratch/fronga/ntuples/V030300v3/DYJetsToLL_M-10To50filter_8TeV-madgraph-Summer12-PU_S7_START52_V9-v1.root","Z+Jets",7123299,LowZJetsCrossSection,false,false,7,dy_color);      
252 –      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/TTJets_TuneZ2star_8TeV-madgraph-tauola-Summer12-PU_S7_START52_V9-v1.root","t#bar{t}",6719290,TTbarCrossSection,false,false,1,ttbar_color);
232  
233 < //      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);
234 < //        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);
235 <        
236 < /*      if(PlottingSetup::DoBTag) {
237 <        PlottingSetup::jzbvariabledata="(jzb[1]+0.0924584*pt)";
238 <        PlottingSetup::jzbvariablemc="(jzb[1]+0.0477569*pt)";
239 <      }*/
233 >      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);
234 >      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);
235 >      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);
236 >      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);
237 >
238 > //       allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/OLD/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball-Summer12-PU_S7_START52_V9-v2.root","Z+Jets",0,ZJetsCrossSection,false,false,7,dy_color);
239 > //       allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/OLD/DYJetsToLL_M-10To50filter_8TeV-madgraph-Summer12-PU_S7_START52_V9-v1.root","Z+Jets",0,LowZJetsCrossSection,false,false,7,dy_color);
240 >
241 >      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/OLD/TTJets_TuneZ2star_8TeV-madgraph-tauola-Summer12-PU_S7_START52_V5-v1.root","t#bar{t}",0,TTbarCrossSection,false,false,1,ttbar_color);//summer11
242 >      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/OLD/WJetsToLNu_TuneZ2Star_8TeV-madgraph-tarball-Summer12-PU_S7_START52_V9-v1.root","W+Jets",0,36257.2,false,false,3,wjets_color);//summer
243 >
244 >      // SINGLE-TOP
245 >      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/OLD/T_s-channel_TuneZ2star_8TeV-powheg-tauola-Summer12-PU_S7_START52_V9-v1.root","Single-top",0,1.76,false,false,2,singletop_color);
246 >      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/OLD/Tbar_s-channel_TuneZ2star_8TeV-powheg-tauola-Summer12-PU_S7_START52_V9-v1.root","Single-top",0,3.79,false,false,2,singletop_color);
247 >      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/OLD/Tbar_t-channel_TuneZ2star_8TeV-powheg-tauola-Summer12-PU_S7_START52_V9-v1.root","Single-top",0,56.4,false,false,2,singletop_color);
248 >      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/OLD/Tbar_tW-channel-DR_TuneZ2star_8TeV-powheg-tauola-Summer12-PU_S7_START52_V9-v1.root","Single-top",0,SingleTopUCrossSection,false,false,2,singletop_color);
249 >      
250 >      // DIBOSONS: LO with a k-factor
251 >      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/OLD/ZZJetsTo4L_TuneZ2star_8TeV-madgraph-tauola-Summer12-PU_S7_START52_V9-v3.root","Dibosons",0,0.181,false,false,4,diboson_color);
252 >      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);
253 >      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/OLD/WZJetsTo3LNu_TuneZ2_8TeV-madgraph-tauola-Summer12-PU_S7_START52_V9-v2.root","Dibosons",0,1.06,false,false,4,diboson_color);
254 >      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);
255 >
256 >      //ttV
257 >      //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);
258 >      //allsamples.AddSample("/shome/fronga/ntuples/8TeV/TTGJets_8TeV-madgraph-Summer12-PU_S7_START52_V9-v1.root",   "ttV",0,1.444,false,false,5,kGreen);
259 >      //allsamples.AddSample("/shome/fronga/ntuples/8TeV/TTWJets_8TeV-madgraph-Summer12-PU_S7_START52_V9-v1.root",   "ttV",0,0.2149,false,false,5,kGreen);
260 >      
261 >      
262 >      //  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);
263 >      
264 >      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/OLD/leo-PYTHIA6_SUSY_LM6_sftsht_8TeV_GEN_FASTSIM-2012_Startup_inTimeOnly.root","LM6",0,LM4CrossSection,false,true,1,lm_color);
265          
266      }
263 –    if(!PlottingSetup::RestrictToMassPeak) {
264 –        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);
265 –    }
266 –
267  
268      if(PlottingSetup::is2012) {
269 < //        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);
270 < //        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);
271 < //        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);
272 <      //cutWeight=TCut("PUweight");
269 >        signalsamples.AddSample("/shome/buchmann/ntuples/MC8tev/OLD/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);
270 >        signalsamples.AddSample("/shome/buchmann/ntuples/MC8tev/OLD/leo-PYTHIA6_SUSY_LM6_sftsht_8TeV_GEN_FASTSIM-2012_Startup_inTimeOnly.root","LM6",0,LM6CrossSection,false,true,2,lm_color);
271 >        signalsamples.AddSample("/shome/buchmann/ntuples/MC8tev/OLD/PYTHIA6_SUSY_LM3_sftsht_8TeV_GEN_FASTSIM-2012_Startup_inTimeOnly-V3-leo_AllWeightOne.root","LM3",0,LM3CrossSection,false,true,3,lm_color);
272 > //        signalsamples.AddSample("/shome/buchmann/ntuples/MC8tev/OLD/leo-PYTHIA6_SUSY_LM8_sftsht_8TeV_GEN_FASTSIM-2012_Startup_inTimeOnly.root","LM8",0,LM6CrossSection,false,true,3,lm_color);
273 >        cutWeight=TCut("(weight*(weight<1000)*(is_data+(!is_data)*((id1==id2)*0.95+(id1!=id2)*0.94)))");
274 >        //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)))");
275 >        //cutWeight=TCut("(weight*(weight<1000)*((id1==id2&&id1==0)*0.95+(id1==id2&&id1==1)*0.88+(id1!=id2)*0.92))");
276      }
277  
278  
# Line 279 | Line 282 | http://cms.cern.ch/iCMS/prep/requestmana
282        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);
283        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
284      } else {
285 <      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);
286 <      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);
285 >      
286 >      systsamples.AddSample("/shome/buchmann/ntuples/MC8tev/OLD/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball-Summer12-PU_S7_START52_V9-v2.root","Z+Jets",0,ZJetsCrossSection,false,false,7,dy_color);
287 >      systsamples.AddSample("/shome/buchmann/ntuples/MC8tev/OLD/DYJetsToLL_M-10To50filter_8TeV-madgraph-Summer12-PU_S7_START52_V9-v1.root","Z+Jets",0,LowZJetsCrossSection,false,false,7,dy_color);
288 >      systsamples.AddSample("/shome/buchmann/ntuples/MC8tev/OLD/TTJets_TuneZ2star_8TeV-madgraph-tauola-Summer12-PU_S7_START52_V5-v1.root","t#bar{t}",0,TTbarCrossSection,false,false,1,ttbar_color);//summer11
289 >      
290      }
291  
292  
# Line 301 | Line 307 | http://cms.cern.ch/iCMS/prep/requestmana
307  
308      
309      
310 <    if(showList) allsamples.ListSamples();
310 >    if(showList) {
311 >      allsamples.ListSamples();
312 >      signalsamples.ListSamples();
313 >      systsamples.ListSamples();
314 >    }
315  
316   }

Diff Legend

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