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.22 by buchmann, Mon Jul 9 12:37:11 2012 UTC

# Line 112 | 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 ;
117    float LM2CrossSection           =0.6027 *1.33 ;
118    float LM3CrossSection           = 3.438*1.40 ; // LO*k*/
115      float LM4CrossSection           = 1.879*1.35; // k*LO
116      if(PlottingSetup::is2012) LM4CrossSection=2.88199; //(only LO, no k factor)
121    /*
122    float LM5CrossSection           = 0.473*1.34; // k*LO*/
117      float LM6CrossSection           = 0.3104*1.30; // k*LO
118      if(PlottingSetup::is2012) LM6CrossSection=0.501773;//(only LO, no k factor)
125    /*
126    float LM7CrossSection           = 1.209*1.11; // k*LO*/
119      float LM8CrossSection           = 0.7300*1.41; // k*LO
120      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;
121      
122      long totEventsZjets         = 32005720;//summer11
123      long totEventslowZjets      = 10272470;
124  
150 //    long totEventsZjetsD6T      = 2543706;//summer11
151 //    long totEventsZjetsD6Tlow   = 181799;//summer11
152
153    //spring:2313911.0;//2743142.0; // DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola
125      long totEventsTTbar         =  3701947;//summer11
155    //1161621.0;//Spring//1144028.0; // PabloV13/TTJets_TuneZ2_7TeV-madgraph-tauola
126      long totEventsWJets         = 56674902;//summer11
157    //15010237.0;//spring//14350756.0; /// WJetsToLNu_TuneZ2_7TeV-madgraph-tauola.root
127      long totEventsZnunu         =  2106977;//2167964.0; //ZinvisibleJets_7TeV-madgraph.root
128      long totEventsVVJets        =   959076;//509072.0;
129      long totEventsSingleTopS    =   493868;//489472.0;
130      long totEventsSingleTopT    =   475460;//477610.0;
131      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
132      
133      Int_t nice_blue  = TColor::GetColor("#2E9AFE");
134      Int_t nice_green = TColor::GetColor("#81f781");
# Line 226 | Line 180 | http://cms.cern.ch/iCMS/prep/requestmana
180        
181   //      allsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_200__M12_440.root","LM8",10000,0.132,false,true,8,lm_color);
182        
183 +      
184 +      PlottingSetup::emuncertOFFPEAK=0.1;
185 +      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");
186      } else {
187        
188        write_info(__FUNCTION__,"Setting up 2012 configuration - once you've stopped looking at 2011 data please make this the default!");
# Line 234 | Line 191 | http://cms.cern.ch/iCMS/prep/requestmana
191        PlottingSetup::luminosity=PlottingSetup::luminosity2012;
192        PlottingSetup::lumiuncert=PlottingSetup::lumiuncert2012;
193        
194 <      // cJZB response correction (705 /pb) :
195 <      PlottingSetup::jzbvariabledata="(jzb[1]+0.0656534*pt)";
196 <      PlottingSetup::jzbvariablemc="(jzb[1]+0.0322238*pt)";
194 >      // Values for 5.1 /fb (2012, cJZB)
195 >      PlottingSetup::jzbvariabledata="(jzb[1]+0.064814*pt)";
196 >      PlottingSetup::jzbvariablemc="(jzb[1]+0.036854*pt)";
197        
198 +        
199 +        
200        if(!PlottingSetup::RestrictToMassPeak) {
201 <        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)))";
202 <        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)))";
201 >        // Values for 5.1 /fb (2012, iJZB)
202 >        PlottingSetup::jzbvariabledata="(jzb[1]+0.059979*pt)";
203 >        PlottingSetup::jzbvariablemc="(jzb[1]+0.034665*pt)";
204  
205        }
206 +      
207 + //  Note: Using same response for bJZB as for cJZB
208 + /*      if(PlottingSetup::DoBTag) {
209 +        // Values for 5.1 /fb (2012, bJZB)
210 +        PlottingSetup::jzbvariabledata="(jzb[1]+0.082732*pt)";
211 +        PlottingSetup::jzbvariablemc="(jzb[1]+0.075*pt)";
212 +      }*/
213 +      
214          
215 + SingleTopSCrossSection=5.55;//Top mass: m(top)=173GeV, Default scales: muF=muR=mt, PDF: MSTW2008 NNLO
216 + SingleTopTCrossSection=87.1;//Top mass: m(top)=173GeV, Default scales: muF=muR=mt, PDF: MSTW2008 NNLO
217 + SingleTopUCrossSection=22.2;//Top mass: m(top)=173GeV, Default scales: muF=muR=mt, PDF: MSTW2008 NNLO
218  
248      allsamples.AddSample("/shome/buchmann/ntuples/Data/Data_3553ifb.root","Data",1,1,true,false,0,kBlack);
219        
220 <      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);
221 <      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);      
222 <      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);
220 >      
221 >      
222 >      allsamples.AddSample("/shome/buchmann/ntuples/Data/Data_5100ipb.root","Data",1,1,true,false,0,kBlack);
223 >      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);
224 >      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);
225 >      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
226 >      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
227 > //       allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/VVJetsTo4L_TuneD6T_7TeV-madgraph-tauola.root","WW/WZ/ZZ",totEventsVVJets,VVJetsCrossSection,false,false,4,diboson_color);
228 >      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);
229 >      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);
230 >      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);
231 >      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);
232 >      
233 >      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);
234 >      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);
235 >      
236 >      
237 > //       allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/ZinvisibleJets_7TeV-madgraph.root","Z nunu",totEventsZnunu,ZnunuCrossSection,false,false,4,diboson_color);
238 > //       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);
239 >      
240 >      
241  
242 < //      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);
255 < //        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);
256 <        
257 < /*      if(PlottingSetup::DoBTag) {
258 <        PlottingSetup::jzbvariabledata="(jzb[1]+0.0924584*pt)";
259 <        PlottingSetup::jzbvariablemc="(jzb[1]+0.0477569*pt)";
260 <      }*/
242 >      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);
243          
244      }
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
245  
246      if(PlottingSetup::is2012) {
247 < //        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);
248 < //        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);
247 >        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);
248 >        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);
249   //        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);
250 <      //cutWeight=TCut("PUweight");
250 >        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)))");
251 > //        cutWeight=TCut("(weight*(weight<1000)*((id1==id2&&id1==0)*0.95+(id1==id2&&id1==1)*0.88+(id1!=id2)*0.92))");
252      }
253  
254  
# Line 279 | Line 258 | http://cms.cern.ch/iCMS/prep/requestmana
258        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);
259        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
260      } else {
261 +      
262        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);
263 <      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);
263 >      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);
264 >      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
265 >      
266      }
267  
268  

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines