24 |
|
{ |
25 |
|
char hostname[1023]; |
26 |
|
gethostname(hostname,1023); |
27 |
+ |
|
28 |
+ |
SanityChecks(); |
29 |
|
|
30 |
|
/* |
31 |
|
|
71 |
|
} |
72 |
|
|
73 |
|
|
74 |
< |
void define_samples(bool showList, samplecollection &allsamples, samplecollection &signalsamples, samplecollection &scansample, samplecollection &raresample, samplecollection &systsamples, samplecollection &qcdsamples) |
74 |
> |
void define_samples(bool showList, samplecollection &allsamples, samplecollection &signalsamples, samplecollection &scansample, samplecollection &raresample, samplecollection &systsamples, samplecollection &qcdsamples, samplecollection &comparesamples) |
75 |
|
{ |
76 |
< |
|
77 |
< |
/* |
78 |
< |
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 |
< |
|
76 |
> |
SanityChecks(); |
77 |
> |
|
78 |
> |
|
79 |
|
flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak |
80 |
|
if(!PlottingSetup::RestrictToMassPeak) { |
81 |
|
PlottingSetup::cutmass=PlottingSetup::openmasscut; |
82 |
|
PlottingSetup::sidebandcut=TCut("mll<2&&mll>3");//impossible cut - because we don't want any sidebands! |
83 |
|
PlottingSetup::genMassCut=PlottingSetup::openGenmasscut;//impossible cut - because we don't want any sidebands! |
84 |
|
} |
85 |
< |
|
86 |
< |
float ZJetsCrossSection = 3048.0; //NNLO---------------------------------- |
87 |
< |
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: |
109 |
< |
//https://twiki.cern.ch/twiki/bin/view/CMS/SingleTopSigma |
110 |
< |
float SingleTopSCrossSection = 4.63; // NLO;---------------------------------- |
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 |
< |
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; |
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; |
144 |
< |
|
145 |
< |
long totEventsZjets = 32005720;//summer11 |
146 |
< |
long totEventslowZjets = 10272470; |
147 |
< |
|
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 |
153 |
< |
//1161621.0;//Spring//1144028.0; // PabloV13/TTJets_TuneZ2_7TeV-madgraph-tauola |
154 |
< |
long totEventsWJets = 56674902;//summer11 |
155 |
< |
//15010237.0;//spring//14350756.0; /// WJetsToLNu_TuneZ2_7TeV-madgraph-tauola.root |
156 |
< |
long totEventsZnunu = 2106977;//2167964.0; //ZinvisibleJets_7TeV-madgraph.root |
157 |
< |
long totEventsVVJets = 959076;//509072.0; |
158 |
< |
long totEventsSingleTopS = 493868;//489472.0; |
159 |
< |
long totEventsSingleTopT = 475460;//477610.0; |
160 |
< |
long totEventsSingleTopU = 489417;//477599.0; |
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"); |
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 :-) |
183 |
< |
|
184 |
< |
|
185 |
< |
Int_t dy_color = nice_blue; |
186 |
< |
Int_t wjets_color = kGray; |
187 |
< |
Int_t ttbar_color = nice_green; |
188 |
< |
Int_t singletop_color = nice_orange; |
189 |
< |
Int_t qcd_color = nice_black; |
190 |
< |
Int_t diboson_color = nice_pink; |
191 |
< |
Int_t lm_color = kRed; |
192 |
< |
|
193 |
< |
dy_color = kYellow; |
194 |
< |
wjets_color = kGray; |
195 |
< |
ttbar_color = kMagenta+2; |
196 |
< |
singletop_color = kBlue; |
197 |
< |
qcd_color = kPink; |
198 |
< |
diboson_color = kGreen+2; |
199 |
< |
lm_color = kViolet+7; |
200 |
< |
|
201 |
< |
// Int_t rare_color = nice_blue; |
202 |
< |
|
203 |
< |
char hostname[1023]; |
204 |
< |
gethostname(hostname,1023); |
85 |
> |
if(PlottingSetup::RestrictToMassPeak&&!PlottingSetup::UseSidebandsForcJZB) { |
86 |
> |
PlottingSetup::sidebandcut=TCut("mll<2&&mll>3"); |
87 |
> |
} |
88 |
|
|
206 |
– |
string prefix="/shome/"; |
207 |
– |
if((int)((string)hostname).find("falda")>-1) prefix="/scratch/"; |
89 |
|
|
90 |
< |
if(!PlottingSetup::is2012) { |
91 |
< |
allsamples.AddSample(prefix+"buchmann/ntuples/Data/AllData_Complete_4653_ipb_PAPER.root","Data",1,1,true,false,0,kBlack); |
92 |
< |
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); |
93 |
< |
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); |
94 |
< |
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 |
95 |
< |
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 |
96 |
< |
allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/VVJetsTo4L_TuneD6T_7TeV-madgraph-tauola.root","WW/WZ/ZZ",totEventsVVJets,VVJetsCrossSection,false,false,4,diboson_color); |
97 |
< |
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); |
98 |
< |
allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/TToBLNu_TuneZ2_t-channel_7TeV-madgraph.root","Single-top",totEventsSingleTopT,SingleTopTCrossSection,false,false,2,singletop_color); |
99 |
< |
allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/TToBLNu_TuneZ2_tW-channel_7TeV-madgraph.root","Single-top",totEventsSingleTopU,SingleTopUCrossSection,false,false,2,singletop_color); |
100 |
< |
allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/ZinvisibleJets_7TeV-madgraph.root","Z nunu",totEventsZnunu,ZnunuCrossSection,false,false,4,diboson_color); |
101 |
< |
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); |
102 |
< |
allsamples.AddSample(prefix+"/buchmann/ntuples/MC/JZBplus/GVJets_7TeV-madgraph-Summer11-PU_S4.root","#gamma+VJets",1067879,56.64 ,false,false,6,kRed); |
103 |
< |
//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); |
104 |
< |
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 { |
90 |
> |
|
91 |
> |
|
92 |
> |
|
93 |
> |
bool SplitByJetFlavor=false; |
94 |
> |
bool SplitByLeptonFlavor=true; |
95 |
> |
bool SplitDibosons=false; |
96 |
> |
bool SplitOthers=true; // splits others into "Dibosons","W+Jets", and "Single-top" |
97 |
> |
|
98 |
> |
if(PlottingSetup::IsJZBAnalysis) { |
99 |
> |
//our main analysis |
100 |
> |
SplitByJetFlavor=false; |
101 |
> |
SplitByLeptonFlavor=true; |
102 |
> |
SplitDibosons=false; |
103 |
> |
SplitOthers=false; |
104 |
> |
} |
105 |
|
|
106 |
+ |
if(PlottingSetup::IsWZAnalysis) { |
107 |
+ |
//WZ cross check for main analysis |
108 |
+ |
SplitByJetFlavor=false; |
109 |
+ |
SplitByLeptonFlavor=false; |
110 |
+ |
SplitDibosons=true; |
111 |
+ |
SplitOthers=true; |
112 |
+ |
} |
113 |
+ |
if(PlottingSetup::IsZbAnalysis) { |
114 |
+ |
//b-jet residual correction calculation |
115 |
+ |
SplitByJetFlavor=true; |
116 |
+ |
SplitByLeptonFlavor=false; |
117 |
+ |
SplitDibosons=false; |
118 |
+ |
SplitOthers=true; |
119 |
+ |
} |
120 |
+ |
|
121 |
+ |
if(!PlottingSetup::is2012) { |
122 |
+ |
// 2011 Data and MC |
123 |
+ |
PlottingSetup::emuncertOFFPEAK=0.1; |
124 |
+ |
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"); |
125 |
+ |
LoadSamplesFromJSON("Samples2011.samples",SplitByJetFlavor,SplitByLeptonFlavor,SplitDibosons,SplitOthers); |
126 |
+ |
} else { |
127 |
+ |
// 2012 Data and MC |
128 |
|
write_info(__FUNCTION__,"Setting up 2012 configuration - once you've stopped looking at 2011 data please make this the default!"); |
129 |
< |
cutWeight=TCut("weight3D"); |
130 |
< |
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. |
129 |
> |
// 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. |
130 |
> |
PlottingSetup::basicqualitycut=TCut("(pfJetGoodNum40>=2&&pfJetGoodID[0]!=0)&&(pfJetGoodNum40>=2&&pfJetGoodID[1]!=0)"&&PlottingSetup::basiccut); |
131 |
|
PlottingSetup::cutnJets=TCut("pfJetGoodNum40>=3"&&PlottingSetup::basicqualitycut); |
132 |
+ |
|
133 |
|
PlottingSetup::luminosity=PlottingSetup::luminosity2012; |
134 |
|
PlottingSetup::lumiuncert=PlottingSetup::lumiuncert2012; |
135 |
|
|
136 |
< |
// cJZB response correction (705 /pb) : |
137 |
< |
PlottingSetup::jzbvariabledata="(jzb[1]+0.0656534*pt)"; |
138 |
< |
PlottingSetup::jzbvariablemc="(jzb[1]+0.0322238*pt)"; |
139 |
< |
|
140 |
< |
allsamples.AddSample("/shome/buchmann/ntuples/Data/Data_920ipb_UpdatedTriggers_updatedJetID_v2.root","Data",1,1,true,false,0,kBlack); |
136 |
> |
// Values for 5.1 /fb (2012, cJZB) |
137 |
> |
PlottingSetup::jzbvariabledata="(jzb[1]+0.061016*pt)"; |
138 |
> |
PlottingSetup::jzbvariablemc="(jzb[1]+0.035264*pt)"; |
139 |
> |
|
140 |
> |
// Values for 5.1 /fb (2012, iJZB) |
141 |
> |
if(!PlottingSetup::RestrictToMassPeak) { |
142 |
> |
PlottingSetup::jzbvariabledata="(jzb[1]+0.061016*pt)"; |
143 |
> |
PlottingSetup::jzbvariablemc="(jzb[1]+0.035264*pt)"; |
144 |
> |
} |
145 |
> |
|
146 |
|
|
147 |
< |
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); |
148 |
< |
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); |
149 |
< |
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 |
< |
}*/ |
147 |
> |
if(PlottingSetup::RestrictToMassPeak) { |
148 |
> |
PlottingSetup::passtrig=TCut("pt1>20&&pt2>20"&&PlottingSetup::passtrig); |
149 |
> |
} |
150 |
|
|
151 |
+ |
if ( !PlottingSetup::is53reco ) { |
152 |
+ |
//52 reco samples |
153 |
+ |
LoadSamplesFromJSON("Samples52.samples",SplitByJetFlavor,SplitByLeptonFlavor,SplitDibosons,SplitOthers); |
154 |
+ |
} else { |
155 |
+ |
if(PlottingSetup::IsZbAnalysis) LoadSamplesFromJSON("Samples53_20ifb.samples",SplitByJetFlavor,SplitByLeptonFlavor,SplitDibosons,SplitOthers); |
156 |
+ |
else LoadSamplesFromJSON("Samples53.samples",SplitByJetFlavor,SplitByLeptonFlavor,SplitDibosons,SplitOthers); |
157 |
+ |
} |
158 |
|
} |
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); |
253 |
– |
// signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM2_SUSY_sftsht_7TeV-pythia6.root","LM2",206336,0.6027*1.33,false,true,5,lm_color); |
254 |
– |
write_warning(__FUNCTION__,"CURRENTLY ALL SIGNAL SAMPLES HAVE BEEN DEACTIVATED (EXCEPT LM4 N LM8)"); |
255 |
– |
/* signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM3_SUSY_sftsht_7TeV-pythia6.root","LM3",215120,3.438*1.40,false,true,6,lm_color); |
256 |
– |
signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM5_SUSY_sftsht_7TeV-pythia6.root","LM5",223992,0.473*1.34,false,true,7,lm_color); |
257 |
– |
signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM9_SUSY_sftsht_7TeV-pythia6.root","LM9",227808,7.134*1.48,false,true,8,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); |
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); |
159 |
|
|
160 |
|
if(PlottingSetup::is2012) { |
161 |
< |
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); |
162 |
< |
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); |
161 |
> |
cutWeight=TCut("(weight*(weight<1000)*(is_data+(!is_data)*((id1==id2)*0.95+(id1!=id2)*0.94)))"); |
162 |
> |
// 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)))"); |
163 |
|
} |
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 |
– |
|
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); |
164 |
|
|
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"); |
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 |
312 |
– |
raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTWTo2Lplus2Nu_NTP.root","Rare SS",48764,0.006841,false,false,6,rare_color);//summer11 |
313 |
– |
raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/WWplusTo2L2Nu_NTP.root","Rare SS",48244 ,0.004216,false,false,6,rare_color);//summer11 |
314 |
– |
raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTZTo2Lplus2Nu_NTP.root","Rare SS",36740,0.002024,false,false,6,rare_color);//summer11 |
315 |
– |
raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/WWWTo2Lplus2Nu_NTP.root","Rare SS",43322,0.008957,false,false,6,rare_color);//summer11 |
316 |
– |
raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTWTo2Lminus2Nu_NTP.root","Rare SS",50435,0.002705,false,false,6,rare_color);//summer11 |
317 |
– |
raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/WWWTo2Lminus2Nu_NTP.root","Rare SS",17147,0.004109,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 |
– |
|
165 |
|
// XS taken from: |
166 |
|
// https://twiki.cern.ch/twiki/bin/view/CMS/QCDGammaPlusJetQCD11005 |
167 |
|
|
171 |
|
PlottingSetup::genMassCut=PlottingSetup::openGenmasscut; |
172 |
|
PlottingSetup::requireZ=false; |
173 |
|
} |
174 |
< |
|
174 |
> |
if(PlottingSetup::RestrictToMassPeak&&!PlottingSetup::UseSidebandsForcJZB) { |
175 |
> |
PlottingSetup::sidebandcut=TCut("SideBandsHaveBeenCalledForOnPeakAnalysisWithoutSidebands__ERROR==1"); |
176 |
> |
PlottingSetup::requireZ=true; |
177 |
> |
} |
178 |
|
|
179 |
|
if(PlottingSetup::DoBTag) { |
180 |
|
PlottingSetup::cutnJets=PlottingSetup::cutnJets&&PlottingSetup::bTagRequirement; |
181 |
|
} |
336 |
– |
|
182 |
|
|
183 |
< |
|
184 |
< |
if(showList) allsamples.ListSamples(); |
185 |
< |
|
183 |
> |
if(showList) { |
184 |
> |
allsamples.ListSamples(); |
185 |
> |
signalsamples.ListSamples(); |
186 |
> |
systsamples.ListSamples(); |
187 |
> |
} |
188 |
|
} |