1 |
– |
/**** |
2 |
– |
|
3 |
– |
Off peak status (RestrictToMassPeak) : |
4 |
– |
|
5 |
– |
x Necessary adaptations identified |
6 |
– |
x Started working on necessary adaptations |
7 |
– |
x Necessary adaptations implemented |
8 |
– |
x Necessary adaptations tested |
9 |
– |
|
10 |
– |
DONE! |
11 |
– |
|
12 |
– |
|
13 |
– |
****/ |
14 |
– |
|
1 |
|
#include <iostream> |
2 |
|
#include <vector> |
3 |
|
#include <sys/stat.h> |
39 |
|
|
40 |
|
/// NEW: Only indicate the very first file - please place scans in ntuples/SMS/ and adapt the name here and in SUSYScan.C (the files are loaded when needed avoiding high memory usage and startup times) |
41 |
|
/// If you have different scans (e.g. you have a T5zzl one as well) you can create an additional directory such as SMS_T5zzl, but you need to change that below and in SUSYScan.C as well (where the files are loaded) |
42 |
< |
scansample.AddSample(prefix+"buchmann/ntuples/"+PlottingSetup::ScanSampleDirectory+"/"+mSUGRAorSMS+"_clean_splitup_0_0.root",mSUGRAorSMS,1,1,false,true,0,kRed); // there is only one sample in the scan; we give each event weight "1" |
42 |
> |
// scansample.AddSample(prefix+"buchmann/ntuples/"+PlottingSetup::ScanSampleDirectory+"/"+mSUGRAorSMS+"_clean_splitup_0_0.root",mSUGRAorSMS,1,1,false,true,0,kRed); // there is only one sample in the scan; we give each event weight "1" |
43 |
> |
// scansample.AddSample("/scratch/buchmann/buchmann/AllSamplesStartingAt10GeV_rest5_mSUGRA/mSUGRA_dilepton_m0-220to3000_m12-100to1000_tanb-10andA0-0_7TeV-Pythia6Z/output_8.root","mSUGRA",1,1,false,true,0,kRed); // there is only one sample in the scan; we give each event weight "1" |
44 |
> |
// scansample.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_400__M12_400.root","mSUGRA",1,1,false,true,0,kRed); // there is only one sample in the scan; we give each event weight "1" |
45 |
> |
|
46 |
> |
|
47 |
> |
|
48 |
|
|
49 |
|
if(loadall) { |
50 |
|
for(int ix=0;ix<PlottingSetup::ScanXzones;ix++) { |
56 |
|
|
57 |
|
|
58 |
|
flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak |
59 |
< |
if(!PlottingSetup::RestrictToMassPeak) { |
59 |
> |
if(PlottingSetup::RestrictToMassPeak==false) { |
60 |
|
PlottingSetup::cutmass=PlottingSetup::openmasscut; |
61 |
|
// PlottingSetup::sidebandcut=TCut("mll<2&&mll>3");//impossible cut - because we don't want any sidebands! |
62 |
|
PlottingSetup::sidebandcut=TCut("SideBandsHaveBeenCalledForOffPeakAnalysis__ERROR==1");//this will trigger errors whenever sidebands are called - perfect for checking any remaining uses (checked, and no uses were found). |
63 |
|
PlottingSetup::genMassCut=PlottingSetup::openGenmasscut;//impossible cut - because we don't want any sidebands! |
64 |
+ |
PlottingSetup::requireZ=false; |
65 |
|
} |
66 |
+ |
dout << "Currently configuration: The mass cut is \"" << PlottingSetup::cutmass << "\" and RestrictToMassPeak is currently set to " << PlottingSetup::RestrictToMassPeak << endl; |
67 |
+ |
|
68 |
+ |
|
69 |
|
} |
70 |
|
|
71 |
|
|
72 |
< |
void define_samples(bool showList, samplecollection &allsamples, samplecollection &signalsamples, samplecollection &scansample, samplecollection &raresample, samplecollection &systsamples) |
72 |
> |
void define_samples(bool showList, samplecollection &allsamples, samplecollection &signalsamples, samplecollection &scansample, samplecollection &raresample, samplecollection &systsamples, samplecollection &qcdsamples) |
73 |
|
{ |
74 |
+ |
|
75 |
+ |
/* |
76 |
+ |
Good sources for XS: |
77 |
+ |
https://twiki.cern.ch/twiki/bin/view/CMS/ProductionSummer2011 |
78 |
+ |
e.g. |
79 |
+ |
https://twiki.cern.ch/twiki/bin/view/CMS/ProductionSummer2011#PYTHIA6 |
80 |
+ |
|
81 |
+ |
In general, all the xsec for datasets are in PREP, eg: |
82 |
+ |
http://cms.cern.ch/iCMS/prep/campaignmanagement?id=Summer11 |
83 |
+ |
|
84 |
+ |
or |
85 |
+ |
|
86 |
+ |
http://cms.cern.ch/iCMS/prep/requestmanagement?dsn=QCD_Pt-15to30_TuneZ2_7TeV_pythia6&campid=Summer11 |
87 |
+ |
*/ |
88 |
+ |
|
89 |
+ |
// source for 8 TeV: https://twiki.cern.ch/twiki/bin/viewauth/CMS/StandardModelCrossSectionsat8TeV |
90 |
+ |
|
91 |
|
flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak |
92 |
|
if(!PlottingSetup::RestrictToMassPeak) { |
93 |
|
PlottingSetup::cutmass=PlottingSetup::openmasscut; |
96 |
|
} |
97 |
|
|
98 |
|
float ZJetsCrossSection = 3048.0; //NNLO---------------------------------- |
99 |
< |
// float ZJetsCrossSectiontiny = 2659.0; |
100 |
< |
float ZJetsCrossSectionlow = 4998.0; |
101 |
< |
// float ZJetsCrossSectionM20 = 4998.0;//summer11 |
102 |
< |
// float ZJetsCrossSectionD6Tlow = 400.0;//summer11 |
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 |
> |
|
103 |
|
float TTbarCrossSection = 165.0;//(NLO) ---- 165.0; // approx. NNLO----- |
104 |
+ |
if(PlottingSetup::is2012) TTbarCrossSection=225.197; |
105 |
+ |
// at 8 TeV: 225.197 (https://twiki.cern.ch/twiki/bin/view/CMS/Zpt8TeV) |
106 |
|
float WJetsCrossSection = 31314.0;//NNLO-------3.131e4; //NNLO------------ |
107 |
|
float ZnunuCrossSection = 5760.0;//NNLO -------4.5e+3; //(LO);------------ |
108 |
|
//Single top x-sections from here: |
111 |
|
float SingleTopTCrossSection = 64.57; // NLO;-------------------------------- |
112 |
|
float SingleTopUCrossSection = 15.74; // NLO;---------------------------------- |
113 |
|
float VVJetsCrossSection = 4.8; // LO;------------------------------------- |
114 |
+ |
/* float LM0CrossSection = 38.93 *1.41 ; |
115 |
+ |
float LM1CrossSection = 4.888 *1.34 ; |
116 |
+ |
float LM2CrossSection =0.6027 *1.33 ; |
117 |
+ |
float LM3CrossSection = 3.438*1.40 ; // LO*k*/ |
118 |
|
float LM4CrossSection = 1.879*1.35; // k*LO |
119 |
< |
float LM8CrossSection = 0.7300*1.41 ; // k*LO |
119 |
> |
if(PlottingSetup::is2012) LM4CrossSection=2.88199; //(only LO, no k factor) |
120 |
> |
/* |
121 |
> |
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) |
124 |
> |
/* |
125 |
> |
float LM7CrossSection = 1.209*1.11; // k*LO*/ |
126 |
> |
float LM8CrossSection = 0.7300*1.41; // k*LO |
127 |
> |
/* float LM9CrossSection = 7.134*1.48; // k*LO |
128 |
> |
float LM10CrossSection = 0.04778*1.22 ; // k*LO |
129 |
> |
float LM11CrossSection = 0.8236 *1.35 ; // k*LO |
130 |
> |
float LM12CrossSection = 4.414 *1.34 ; // k*LO |
131 |
> |
float LM13CrossSection = 6.899*1.42 ; // k*LO*/ |
132 |
|
//the following numbers are from the MadGraphStandardModel210Summary as linked on the GeneratorProduction2011 page |
133 |
|
// float QCD50to100CrossSection=30000000; // not used |
134 |
< |
float QCD100to250CrossSection=7000000; |
135 |
< |
float QCD250to500CrossSection=171000; |
136 |
< |
float QCD500to1000CrossSection=5200; |
137 |
< |
float QCD1000toInfCrossSection=83.0; |
134 |
> |
//float QCD100to250CrossSection=7000000; |
135 |
> |
//float QCD250to500CrossSection=171000; |
136 |
> |
//float QCD500to1000CrossSection=5200; |
137 |
> |
//float QCD1000toInfCrossSection=83.0; |
138 |
|
|
139 |
|
// float totEventsQCD50to100=207418.0; // not used |
140 |
< |
long totEventsQCD100to250 = 638792; |
141 |
< |
long totEventsQCD250to500 = 344454; |
142 |
< |
long totEventsQCD500to1000 = (long)10e10; //unknown but now suppressed. (no event passes anyway) |
143 |
< |
long totEventsQCD1000toInf = 163185; |
140 |
> |
//long totEventsQCD100to250 = 638792; |
141 |
> |
//long totEventsQCD250to500 = 344454; |
142 |
> |
//long totEventsQCD500to1000 = (long)10e10; //unknown but now suppressed. (no event passes anyway) |
143 |
> |
//long totEventsQCD1000toInf = 163185; |
144 |
|
|
145 |
|
long totEventsZjets = 32005720;//summer11 |
146 |
< |
long totEventsZjetsEElow = 29169462; |
117 |
< |
long totEventsZjetsEEtiny = 9310870; |
118 |
< |
long totEventsZjetsMMlow = 29479180; |
119 |
< |
long totEventsZjetsMMtiny = 9758899; |
120 |
< |
long totEventsZjetsTTlow = 19937479; |
121 |
< |
long totEventsZjetsTTtiny = 2200000; |
146 |
> |
long totEventslowZjets = 10272470; |
147 |
|
|
148 |
< |
long totEventsZjetsD6T = 2543706;//summer11 |
149 |
< |
long totEventsZjetsD6Tlow = 181799;//summer11 |
148 |
> |
// long totEventsZjetsD6T = 2543706;//summer11 |
149 |
> |
// long totEventsZjetsD6Tlow = 181799;//summer11 |
150 |
|
|
151 |
|
//spring:2313911.0;//2743142.0; // DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola |
152 |
|
long totEventsTTbar = 3701947;//summer11 |
158 |
|
long totEventsSingleTopS = 493868;//489472.0; |
159 |
|
long totEventsSingleTopT = 475460;//477610.0; |
160 |
|
long totEventsSingleTopU = 489417;//477599.0; |
161 |
< |
long totEventsLM4 = 218380; |
162 |
< |
long totEventsLM8 = 180960; |
161 |
> |
/* long totEventsLM0 = 1234567;//bad guess, will be replaced. |
162 |
> |
long totEventsLM1 = 1234567;//bad guess, will be replaced. |
163 |
> |
long totEventsLM2 = 1234567;//bad guess, will be replaced.*/ |
164 |
> |
long totEventsLM3 = 220000;//bad guess, will be replaced. |
165 |
> |
/* long totEventsLM4 = 218380; |
166 |
> |
long totEventsLM5 = 220000; |
167 |
> |
long totEventsLM6 = 220000; |
168 |
> |
long totEventsLM7 = 217800;*/ |
169 |
> |
long totEventsLM8 = 217800;/* |
170 |
> |
long totEventsLM9 = 220000; |
171 |
> |
long totEventsLM10 = 1234567; |
172 |
> |
long totEventsLM11 = 216999; |
173 |
> |
long totEventsLM12 = 217400; |
174 |
> |
long totEventsLM13 = 1234567;*/ |
175 |
|
|
176 |
|
|
177 |
|
Int_t nice_blue = TColor::GetColor("#2E9AFE"); |
178 |
|
Int_t nice_green = TColor::GetColor("#81f781"); |
179 |
< |
Int_t nice_red = TColor::GetColor("#F78181"); |
179 |
> |
// Int_t nice_red = TColor::GetColor("#F78181"); |
180 |
|
Int_t nice_pink = TColor::GetColor("#F781BE"); |
181 |
|
Int_t nice_orange= TColor::GetColor("#F7BE81");//#FA9624 |
182 |
|
Int_t nice_black = TColor::GetColor("#2A0A1B");//black with a bit of purple :-) |
198 |
|
diboson_color = kGreen+2; |
199 |
|
lm_color = kViolet+7; |
200 |
|
|
201 |
< |
Int_t rare_color = nice_blue; |
201 |
> |
// Int_t rare_color = nice_blue; |
202 |
|
|
203 |
|
char hostname[1023]; |
204 |
|
gethostname(hostname,1023); |
206 |
|
string prefix="/shome/"; |
207 |
|
if((int)((string)hostname).find("falda")>-1) prefix="/scratch/"; |
208 |
|
|
209 |
< |
allsamples.AddSample(prefix+"buchmann/ntuples/Data/AllData_Complete_4653_ipb_PAPER.root","Data",1,1,true,false,0,kBlack); |
210 |
< |
|
211 |
< |
|
212 |
< |
if(!PlottingSetup::RestrictToMassPeak) { |
213 |
< |
allsamples.AddSample(prefix+"lbaeni/jzb/ntuples/MC/JZBplus/DYToEE_M-20_CT10_TuneZ2_7TeV-powheg-pythia__Summer11-PU_S4_START42_V11-v1__AODSIM.root","Z+Jets",totEventsZjetsEElow,ZJetsCrossSectionlow,false,false,7,dy_color); |
214 |
< |
// allsamples.AddSample(prefix+"lbaeni/jzb/ntuples/MC/JZBplus/DYToEE_M-10To20_CT10_TuneZ2_7TeV-powheg-pythia__Summer11-PU_S4_START42_V11-v1__AODSIM.root","lmZ+Jets",totEventsZjetsEEtiny,ZJetsCrossSectiontiny,false,false,9,nice_blue); |
215 |
< |
allsamples.AddSample(prefix+"lbaeni/jzb/ntuples/MC/JZBplus/DYToMuMu_M-20_CT10_TuneZ2_7TeV-powheg-pythia__Summer11-PU_S4_START42_V11-v1__AODSIM.root","Z+Jets",totEventsZjetsMMlow,ZJetsCrossSectionlow,false,false,7,dy_color); |
216 |
< |
// allsamples.AddSample(prefix+"lbaeni/jzb/ntuples/MC/JZBplus/DYToMuMu_M-10To20_CT10_TuneZ2_7TeV-powheg-pythia__Summer11-PU_S4_START42_V11-v1__AODSIM.root","lmZ+Jets",totEventsZjetsMMtiny,ZJetsCrossSectiontiny,false,false,9,nice_blue); |
217 |
< |
allsamples.AddSample(prefix+"lbaeni/jzb/ntuples/MC/JZBplus/DYToTauTau_M-20_CT10_TuneZ2_7TeV-powheg-pythia-tauola__Summer11-PU_S4_START42_V11-v1__AODSIM.root","Z+Jets",totEventsZjetsTTlow,ZJetsCrossSectionlow,false,false,7,dy_color); |
218 |
< |
// allsamples.AddSample(prefix+"lbaeni/jzb/ntuples/MC/JZBplus/DYToTauTau_M-10To20_TuneZ2_7TeV-pythia6-tauola__Summer11-PU_S3_START42_V11-v2__AODSIM.root","lmZ+Jets",totEventsZjetsTTtiny,ZJetsCrossSectiontiny,false,false,9,nice_blue); |
209 |
> |
if(!PlottingSetup::is2012) { |
210 |
> |
allsamples.AddSample(prefix+"buchmann/ntuples/Data/AllData_Complete_4653_ipb_PAPER.root","Data",1,1,true,false,0,kBlack); |
211 |
> |
allsamples.AddSample(prefix+"buchmann//ntuples/MC/NewSelectionForPaper/NewlyReweighted/DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola_Summer11___NTupleV020300.root","Z+Jets",totEventsZjets,ZJetsCrossSection,false,false,7,dy_color); |
212 |
> |
allsamples.AddSample(prefix+"buchmann/ntuples/MC/Partial_DYJetsToLL_M-10To50_TuneZ2_7TeV-madgraph__Fall11-PU_S6_START42_V14B-v1__AODSIM.root","Z+Jets",totEventslowZjets,LowZJetsCrossSection,false,false,7,dy_color); |
213 |
> |
allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/TTJets_TuneZ2_7TeV-madgraph-tauola_Summer11_3.root","t#bar{t}",totEventsTTbar,TTbarCrossSection,false,false,1,ttbar_color);//summer11 |
214 |
> |
allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/WJetsToLNu_TuneZ2_7TeV-madgraph-tauola_Summer11_2.root","W+Jets",totEventsWJets,WJetsCrossSection,false,false,3,wjets_color);//summer |
215 |
> |
allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/VVJetsTo4L_TuneD6T_7TeV-madgraph-tauola.root","WW/WZ/ZZ",totEventsVVJets,VVJetsCrossSection,false,false,4,diboson_color); |
216 |
> |
allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/TToBLNu_TuneZ2_s-channel_7TeV-madgraph_2.root","Single-top",totEventsSingleTopS,SingleTopSCrossSection,false,false,2,singletop_color); |
217 |
> |
allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/TToBLNu_TuneZ2_t-channel_7TeV-madgraph.root","Single-top",totEventsSingleTopT,SingleTopTCrossSection,false,false,2,singletop_color); |
218 |
> |
allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/TToBLNu_TuneZ2_tW-channel_7TeV-madgraph.root","Single-top",totEventsSingleTopU,SingleTopUCrossSection,false,false,2,singletop_color); |
219 |
> |
allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/ZinvisibleJets_7TeV-madgraph.root","Z nunu",totEventsZnunu,ZnunuCrossSection,false,false,4,diboson_color); |
220 |
> |
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); |
221 |
> |
allsamples.AddSample(prefix+"/buchmann/ntuples/MC/JZBplus/GVJets_7TeV-madgraph-Summer11-PU_S4.root","#gamma+VJets",1067879,56.64 ,false,false,6,kRed); |
222 |
> |
//allsamples.AddSample("/shome/buchmann/ntuples/MC/LM4_SUSY_sftsht_7TeV-pythia6__Summer11-PU_S4_START42_V11-v2_novtx_fixed3___WITHGENINFO.root","LM4",totEventsLM4,LM4CrossSection,false,true,8,lm_color); |
223 |
> |
allsamples.AddSample("/shome//lbaeni/jzb/ntuples/MC/JZBplus/SignalSamples/LM8_SUSY_sftsht_7TeV-pythia6.root","LM8",totEventsLM8,LM8CrossSection,false,true,8,lm_color); |
224 |
|
} else { |
225 |
< |
allsamples.AddSample(prefix+"buchmann//ntuples/MC/NewSelectionForPaper/NewlyReweighted/DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola_Summer11___NTupleV020300.root","Z+Jets",totEventsZjets,ZJetsCrossSection,false,false,7,dy_color); |
225 |
> |
|
226 |
> |
write_info(__FUNCTION__,"Setting up 2012 configuration - once you've stopped looking at 2011 data please make this the default!"); |
227 |
> |
cutWeight=TCut("weight3D"); |
228 |
> |
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. |
229 |
> |
PlottingSetup::cutnJets=TCut("pfJetGoodNum40>=3"&&PlottingSetup::basicqualitycut); |
230 |
> |
PlottingSetup::luminosity=PlottingSetup::luminosity2012; |
231 |
> |
PlottingSetup::lumiuncert=PlottingSetup::lumiuncert2012; |
232 |
> |
|
233 |
> |
// cJZB response correction (705 /pb) : |
234 |
> |
PlottingSetup::jzbvariabledata="(jzb[1]+0.0656534*pt)"; |
235 |
> |
PlottingSetup::jzbvariablemc="(jzb[1]+0.0322238*pt)"; |
236 |
> |
|
237 |
> |
allsamples.AddSample("/shome/buchmann/ntuples/Data/Data_920ipb_UpdatedTriggers_updatedJetID_v2.root","Data",1,1,true,false,0,kBlack); |
238 |
> |
|
239 |
> |
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); |
240 |
> |
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); |
241 |
> |
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); |
242 |
> |
// 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); |
243 |
> |
|
244 |
> |
/* if(PlottingSetup::DoBTag) { |
245 |
> |
PlottingSetup::jzbvariabledata="(jzb[1]+0.0924584*pt)"; |
246 |
> |
PlottingSetup::jzbvariablemc="(jzb[1]+0.0477569*pt)"; |
247 |
> |
}*/ |
248 |
> |
|
249 |
|
} |
185 |
– |
|
186 |
– |
allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/TTJets_TuneZ2_7TeV-madgraph-tauola_Summer11_3.root","t#bar{t}",totEventsTTbar,TTbarCrossSection,false,false,1,ttbar_color);//summer11 |
187 |
– |
allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/WJetsToLNu_TuneZ2_7TeV-madgraph-tauola_Summer11_2.root","W+Jets",totEventsWJets,WJetsCrossSection,false,false,3,wjets_color);//summer |
188 |
– |
allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/VVJetsTo4L_TuneD6T_7TeV-madgraph-tauola.root","WW/WZ/ZZ",totEventsVVJets,VVJetsCrossSection,false,false,4,diboson_color); |
189 |
– |
allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/TToBLNu_TuneZ2_s-channel_7TeV-madgraph_2.root","Single-top",totEventsSingleTopS,SingleTopSCrossSection,false,false,2,singletop_color); |
190 |
– |
allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/TToBLNu_TuneZ2_t-channel_7TeV-madgraph.root","Single-top",totEventsSingleTopT,SingleTopTCrossSection,false,false,2,singletop_color); |
191 |
– |
allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/TToBLNu_TuneZ2_tW-channel_7TeV-madgraph.root","Single-top",totEventsSingleTopU,SingleTopUCrossSection,false,false,2,singletop_color); |
192 |
– |
allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/ZinvisibleJets_7TeV-madgraph.root","Z nunu",totEventsZnunu,ZnunuCrossSection,false,false,4,diboson_color); |
193 |
– |
allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/LM4_SUSY_sftsht_7TeV-pythia6.root","LM4",totEventsLM4,LM4CrossSection,false,true,8,lm_color); |
194 |
– |
|
250 |
|
if(!PlottingSetup::RestrictToMassPeak) { |
251 |
|
// signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM0_SUSY_sftsht_7TeV-pythia6.root","LM0",219796,38.93*1.41 ,false,true,3,lm_color); |
252 |
|
// signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM1_SUSY_sftsht_7TeV-pythia6.root","LM1",218176,4.888*1.34,false,true,4,lm_color); |
258 |
|
signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM11_SUSY_sftsht_7TeV-pythia6.root","LM11",218703,0.8236*1.35,false,true,9,lm_color); |
259 |
|
signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM12_SUSY_sftsht_7TeV-pythia6.root","LM12",223491,4.414*1.34,false,true,10,lm_color); |
260 |
|
signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM13_SUSY_sftsht_7TeV-pythia6.root","LM13",173240,6.899*1.42,false,true,11,lm_color);*/ |
261 |
+ |
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); |
262 |
|
} |
263 |
< |
signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM4_SUSY_sftsht_7TeV-pythia6.root","LM4",totEventsLM4,LM4CrossSection,false,true,1,lm_color); |
264 |
< |
signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM8_SUSY_sftsht_7TeV-pythia6.root","LM8",totEventsLM8,LM8CrossSection,false,true,2,lm_color); |
263 |
> |
// signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM4_SUSY_sftsht_7TeV-pythia6.root","LM4",totEventsLM4,LM4CrossSection,false,true,1,lm_color); |
264 |
> |
// signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM8_SUSY_sftsht_7TeV-pythia6.root","LM8",totEventsLM8,LM8CrossSection,false,true,2,lm_color); |
265 |
> |
|
266 |
> |
/* signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_200__M12_400.root","mSUGRA_M0_200__M12_400",0,1,false,true,1,lm_color); |
267 |
> |
signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_200__M12_300.root","mSUGRA_M0_200__M12_300",0,1,false,true,2,lm_color); |
268 |
> |
signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_400__M12_400.root","mSUGRA_M0_400__M12_400",0,1,false,true,3,lm_color); |
269 |
> |
signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_400__M12_300.root","mSUGRA_M0_400__M12_300",0,1,false,true,4,lm_color); |
270 |
> |
signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_600__M12_400.root","mSUGRA_M0_600__M12_400",0,1,false,true,5,lm_color); |
271 |
> |
signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_600__M12_300.root","mSUGRA_M0_600__M12_300",0,1,false,true,6,lm_color); |
272 |
> |
signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_800__M12_400.root","mSUGRA_M0_800__M12_400",0,1,false,true,7,lm_color); |
273 |
> |
signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_800__M12_300.root","mSUGRA_M0_800__M12_300",0,1,false,true,8,lm_color); |
274 |
> |
signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_1000__M12_400.root","mSUGRA_M0_1000__M12_400",0,1,false,true,9,lm_color); |
275 |
> |
signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_1000__M12_300.root","mSUGRA_M0_1000__M12_300",0,1,false,true,10,lm_color); |
276 |
> |
signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_1500__M12_400.root","mSUGRA_M0_1500__M12_400",0,1,false,true,11,lm_color); |
277 |
> |
signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_1500__M12_300.root","mSUGRA_M0_1500__M12_300",0,1,false,true,12,lm_color);*/ |
278 |
> |
//signalsamples.AddSample("/shome//buchmann/ntuples/SMS_T1lh/SMS_MassGlu_400__MassLSP_200.root","SMS_MassGlu_400__MassLSP_200",0,1,false,true,12,lm_color); |
279 |
> |
//signalsamples.AddSample("/shome//buchmann/ntuples/SMS_T1lh/SMS_MassGlu_500__MassLSP_200.root","SMS_MassGlu_500__MassLSP_200",0,1,false,true,12,lm_color); |
280 |
> |
//signalsamples.AddSample("/shome//buchmann/ntuples/SMS_T1lh/SMS_MassGlu_600__MassLSP_200.root","SMS_MassGlu_600__MassLSP_200",0,1,false,true,12,lm_color); |
281 |
> |
//signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_2000__M12_200.root","M_{0}=2000,M_{1/2}=200",0,1,false,true,12,lm_color); |
282 |
> |
//signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_300__M12_800.root","M_{0}=300,M_{1/2}=800",0,1,false,true,12,lm_color); |
283 |
> |
|
284 |
> |
if(PlottingSetup::is2012) { |
285 |
> |
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); |
286 |
> |
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); |
287 |
> |
} |
288 |
> |
/*signalsamples.AddSample("/shome//buchmann/ntuples/SMS_T1lh/SMS_MassGlu_400__MassLSP_300.root","SMS_MassGlu_400__MassLSP_200",0,1,false,true,12,lm_color); |
289 |
> |
signalsamples.AddSample("/shome//buchmann/ntuples/SMS_T1lh/SMS_MassGlu_500__MassLSP_300.root","SMS_MassGlu_400__MassLSP_200",0,1,false,true,12,lm_color); |
290 |
> |
signalsamples.AddSample("/shome//buchmann/ntuples/SMS_T1lh/SMS_MassGlu_600__MassLSP_300.root","SMS_MassGlu_400__MassLSP_200",0,1,false,true,12,lm_color);*/ |
291 |
|
|
210 |
– |
// systsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/HighStatSamples/DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola_Fall11-PU_S6_START42_V14B-v1.root","Z+Jets",18368890,ZJetsCrossSection,false,false,7,dy_color); |
292 |
|
|
293 |
+ |
// systsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/HighStatSamples/DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola_Fall11-PU_S6_START42_V14B-v1.root","Z+Jets",18368890,ZJetsCrossSection,false,false,7,dy_color); |
294 |
|
|
295 |
< |
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); |
296 |
< |
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 |
295 |
> |
if(!PlottingSetup::is2012) { |
296 |
> |
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); |
297 |
> |
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 |
298 |
> |
} else { |
299 |
> |
systsamples.AddSample("/shome/pablom/tmp/CMSSW_5_2_5/src/DiLeptonAnalysis/NTupleProducer/macros/mc/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball-Summer12-PU_S7_START52_V9-v2.root","Z+Jets",0,ZJetsCrossSection,false,false,7,dy_color); |
300 |
> |
systsamples.AddSample("/shome/pablom/tmp/CMSSW_5_2_5/src/DiLeptonAnalysis/NTupleProducer/macros/mc/V03-03-00_TTJets_TuneZ2star_8TeV-madgraph-tauola-Summer12-PU_S7_START52_V9-v1.root","t#bar{t}",0,TTbarCrossSection,false,false,1,ttbar_color); |
301 |
> |
} |
302 |
|
|
303 |
|
// |
304 |
|
/* |
305 |
|
|
306 |
|
|
307 |
|
// We do not use the following (rare) samples; nonetheless we include them in a separate category to prove that they are negligible |
308 |
< |
write_warning(__FUNCTION__,"Rare samples not yet reprocessed with new selectioN");/* |
308 |
> |
write_warning(__FUNCTION__,"Rare samples not yet reprocessed with new selectioN"); |
309 |
|
raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTWWTo2Lplus2Nu_NTP.root","Rare SS",39791,4.76E-5,false,false,6,rare_color);//summer11 |
310 |
|
raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/WWminusTo2L2Nu_NTP.root","Rare SS",51063,0.001559,false,false,6,rare_color);//summer11 |
311 |
|
raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTWWTo2Lminus2Nu_NTP.root","Rare SS",39496 ,4.719E-5,false,false,6,rare_color);//summer11 |
318 |
|
raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTZTo2Lminus2Nu_NTP.root","Rare SS",36781,0.001946,false,false,6,rare_color);//summer11 |
319 |
|
*/ |
320 |
|
//scan files are only defined when necessary! |
321 |
+ |
|
322 |
+ |
// XS taken from: |
323 |
+ |
// https://twiki.cern.ch/twiki/bin/view/CMS/QCDGammaPlusJetQCD11005 |
324 |
|
|
325 |
+ |
if(!PlottingSetup::RestrictToMassPeak) { |
326 |
+ |
PlottingSetup::cutmass=PlottingSetup::openmasscut; |
327 |
+ |
PlottingSetup::sidebandcut=TCut("SideBandsHaveBeenCalledForOffPeakAnalysis__ERROR==1");//this will trigger errors whenever sidebands are called - perfect for checking any remaining uses (checked, and no uses were found). |
328 |
+ |
PlottingSetup::genMassCut=PlottingSetup::openGenmasscut; |
329 |
+ |
PlottingSetup::requireZ=false; |
330 |
+ |
} |
331 |
+ |
|
332 |
+ |
|
333 |
+ |
if(PlottingSetup::DoBTag) { |
334 |
+ |
PlottingSetup::cutnJets=PlottingSetup::cutnJets&&PlottingSetup::bTagRequirement; |
335 |
+ |
} |
336 |
+ |
|
337 |
+ |
|
338 |
+ |
|
339 |
|
if(showList) allsamples.ListSamples(); |
340 |
|
|
341 |
|
} |