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 |
|
SanityChecks(); |
77 |
< |
|
78 |
< |
/* |
79 |
< |
Good sources for XS: |
80 |
< |
https://twiki.cern.ch/twiki/bin/view/CMS/ProductionSummer2011 |
81 |
< |
e.g. |
82 |
< |
https://twiki.cern.ch/twiki/bin/view/CMS/ProductionSummer2011#PYTHIA6 |
83 |
< |
|
84 |
< |
In general, all the xsec for datasets are in PREP, eg: |
85 |
< |
http://cms.cern.ch/iCMS/prep/campaignmanagement?id=Summer11 |
86 |
< |
|
87 |
< |
or |
88 |
< |
|
89 |
< |
http://cms.cern.ch/iCMS/prep/requestmanagement?dsn=QCD_Pt-15to30_TuneZ2_7TeV_pythia6&campid=Summer11 |
90 |
< |
*/ |
91 |
< |
|
92 |
< |
// source for 8 TeV: https://twiki.cern.ch/twiki/bin/viewauth/CMS/StandardModelCrossSectionsat8TeV |
93 |
< |
|
77 |
> |
|
78 |
> |
|
79 |
|
flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak |
80 |
|
if(!PlottingSetup::RestrictToMassPeak) { |
81 |
|
PlottingSetup::cutmass=PlottingSetup::openmasscut; |
85 |
|
if(PlottingSetup::RestrictToMassPeak&&!PlottingSetup::UseSidebandsForcJZB) { |
86 |
|
PlottingSetup::sidebandcut=TCut("mll<2&&mll>3"); |
87 |
|
} |
103 |
– |
|
104 |
– |
|
105 |
– |
float ZJetsCrossSection = 3048.0; //NNLO---------------------------------- |
106 |
– |
if (PlottingSetup::is2012) ZJetsCrossSection=3503.71; |
107 |
– |
// at 8 TeV: 3503.71 (https://twiki.cern.ch/twiki/bin/view/CMS/Zpt8TeV) |
108 |
– |
float LowZJetsCrossSection = 11908; //https://twiki.cern.ch/twiki/bin/view/CMS/TWikiTop2011DataMCTrig |
109 |
– |
if (PlottingSetup::is2012) LowZJetsCrossSection = 11050*0.069*1.15; // LO xs * filter eff * k-factor (same as ZJets) |
110 |
– |
|
111 |
– |
float TTbarCrossSection = 165.0;//(NLO) ---- 165.0; // approx. NNLO----- |
112 |
– |
if(PlottingSetup::is2012) TTbarCrossSection=225.197; |
113 |
– |
// at 8 TeV: 225.197 (https://twiki.cern.ch/twiki/bin/view/CMS/Zpt8TeV) |
114 |
– |
float WJetsCrossSection = 31314.0;//NNLO-------3.131e4; //NNLO------------ |
115 |
– |
float ZnunuCrossSection = 5760.0;//NNLO -------4.5e+3; //(LO);------------ |
116 |
– |
//Single top x-sections from here: |
117 |
– |
//https://twiki.cern.ch/twiki/bin/view/CMS/SingleTopSigma |
118 |
– |
float SingleTopSCrossSection = 4.63; // NLO;---------------------------------- |
119 |
– |
float SingleTopTCrossSection = 64.57; // NLO;-------------------------------- |
120 |
– |
float SingleTopUCrossSection = 15.74; // NLO;---------------------------------- |
121 |
– |
float VVJetsCrossSection = 4.8; // LO;------------------------------------- |
122 |
– |
float LM3CrossSection = 0.0; // 2011: don't know the x-section |
123 |
– |
if(PlottingSetup::is2012) LM3CrossSection=5.21678; //(only LO, no k factor) |
124 |
– |
float LM4CrossSection = 1.879*1.35; // k*LO |
125 |
– |
if(PlottingSetup::is2012) LM4CrossSection=2.88199; //(only LO, no k factor) |
126 |
– |
float LM6CrossSection = 0.3104*1.30; // k*LO |
127 |
– |
if(PlottingSetup::is2012) LM6CrossSection=0.501773;//(only LO, no k factor) |
128 |
– |
float LM8CrossSection = 0.7300*1.41; // k*LO |
129 |
– |
if(PlottingSetup::is2012) LM8CrossSection=1.14; |
130 |
– |
|
131 |
– |
long totEventsZjets = 32005720;//summer11 |
132 |
– |
long totEventslowZjets = 10272470; |
133 |
– |
|
134 |
– |
long totEventsTTbar = 3701947;//summer11 |
135 |
– |
long totEventsWJets = 56674902;//summer11 |
136 |
– |
long totEventsZnunu = 2106977;//2167964.0; //ZinvisibleJets_7TeV-madgraph.root |
137 |
– |
long totEventsVVJets = 959076;//509072.0; |
138 |
– |
long totEventsSingleTopS = 493868;//489472.0; |
139 |
– |
long totEventsSingleTopT = 475460;//477610.0; |
140 |
– |
long totEventsSingleTopU = 489417;//477599.0; |
88 |
|
|
142 |
– |
Int_t nice_blue = TColor::GetColor("#2E9AFE"); |
143 |
– |
Int_t nice_green = TColor::GetColor("#81f781"); |
144 |
– |
// Int_t nice_red = TColor::GetColor("#F78181"); |
145 |
– |
Int_t nice_pink = TColor::GetColor("#F781BE"); |
146 |
– |
Int_t nice_orange= TColor::GetColor("#F7BE81");//#FA9624 |
147 |
– |
Int_t nice_black = TColor::GetColor("#2A0A1B");//black with a bit of purple :-) |
89 |
|
|
90 |
|
|
150 |
– |
Int_t dy_color = nice_blue; |
151 |
– |
Int_t wjets_color = kGray; |
152 |
– |
Int_t ttbar_color = nice_green; |
153 |
– |
Int_t singletop_color = nice_orange; |
154 |
– |
Int_t qcd_color = nice_black; |
155 |
– |
Int_t diboson_color = nice_pink; |
156 |
– |
Int_t lm_color = kRed; |
157 |
– |
|
158 |
– |
dy_color = kYellow; |
159 |
– |
wjets_color = kGray; |
160 |
– |
ttbar_color = kMagenta+2; |
161 |
– |
singletop_color = kBlue; |
162 |
– |
qcd_color = kPink; |
163 |
– |
diboson_color = kGreen+2; |
164 |
– |
lm_color = kViolet+7; |
165 |
– |
|
166 |
– |
// Int_t rare_color = nice_blue; |
91 |
|
|
92 |
< |
char hostname[1023]; |
93 |
< |
gethostname(hostname,1023); |
94 |
< |
|
95 |
< |
string prefix="/shome/"; |
96 |
< |
if((int)((string)hostname).find("falda")>-1) prefix="/scratch/"; |
97 |
< |
|
98 |
< |
if(!PlottingSetup::is2012) { |
99 |
< |
allsamples.AddSample(prefix+"buchmann/ntuples/Data/AllData_Complete_4653_ipb_PAPER.root","Data",1,1,true,false,0,kBlack); |
100 |
< |
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); |
101 |
< |
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); |
102 |
< |
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 |
103 |
< |
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 |
104 |
< |
allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/VVJetsTo4L_TuneD6T_7TeV-madgraph-tauola.root","WW/WZ/ZZ",totEventsVVJets,VVJetsCrossSection,false,false,4,diboson_color); |
181 |
< |
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); |
182 |
< |
allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/TToBLNu_TuneZ2_t-channel_7TeV-madgraph.root","Single-top",totEventsSingleTopT,SingleTopTCrossSection,false,false,2,singletop_color); |
183 |
< |
allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/TToBLNu_TuneZ2_tW-channel_7TeV-madgraph.root","Single-top",totEventsSingleTopU,SingleTopUCrossSection,false,false,2,singletop_color); |
184 |
< |
allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/ZinvisibleJets_7TeV-madgraph.root","Z nunu",totEventsZnunu,ZnunuCrossSection,false,false,4,diboson_color); |
185 |
< |
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); |
186 |
< |
// allsamples.AddSample(prefix+"/buchmann/ntuples/MC/JZBplus/GVJets_7TeV-madgraph-Summer11-PU_S4.root","#gamma+VJets",1067879,56.64 ,false,false,6,kRed); |
187 |
< |
//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); |
188 |
< |
// allsamples.AddSample("/shome//lbaeni/jzb/ntuples/MC/JZBplus/SignalSamples/LM8_SUSY_sftsht_7TeV-pythia6.root","LM8",totEventsLM8,LM8CrossSection,false,true,8,lm_color); |
189 |
< |
|
190 |
< |
// allsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_200__M12_440.root","LM8",10000,0.132,false,true,8,lm_color); |
191 |
< |
|
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 |
< |
|
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 |
< |
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 |
|
|
201 |
– |
// For Kostas' control region |
202 |
– |
// 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. |
203 |
– |
// PlottingSetup::cutnJets=TCut("pfJetGoodNum40<3"&&PlottingSetup::basicqualitycut); |
204 |
– |
|
133 |
|
PlottingSetup::luminosity=PlottingSetup::luminosity2012; |
134 |
|
PlottingSetup::lumiuncert=PlottingSetup::lumiuncert2012; |
135 |
|
|
136 |
|
// Values for 5.1 /fb (2012, cJZB) |
137 |
< |
PlottingSetup::jzbvariabledata="(jzb[1]+0.064814*pt)"; |
138 |
< |
PlottingSetup::jzbvariablemc="(jzb[1]+0.036854*pt)"; |
137 |
> |
PlottingSetup::jzbvariabledata="(jzb[1]+0.061016*pt)"; |
138 |
> |
PlottingSetup::jzbvariablemc="(jzb[1]+0.035264*pt)"; |
139 |
|
|
140 |
< |
|
213 |
< |
|
140 |
> |
// Values for 5.1 /fb (2012, iJZB) |
141 |
|
if(!PlottingSetup::RestrictToMassPeak) { |
142 |
< |
// Values for 5.1 /fb (2012, iJZB) |
143 |
< |
PlottingSetup::jzbvariabledata="(jzb[1]+0.059979*pt)"; |
217 |
< |
PlottingSetup::jzbvariablemc="(jzb[1]+0.034665*pt)"; |
218 |
< |
|
142 |
> |
PlottingSetup::jzbvariabledata="(jzb[1]+0.061016*pt)"; |
143 |
> |
PlottingSetup::jzbvariablemc="(jzb[1]+0.035264*pt)"; |
144 |
|
} |
145 |
|
|
146 |
|
|
147 |
|
if(PlottingSetup::RestrictToMassPeak) { |
148 |
|
PlottingSetup::passtrig=TCut("pt1>20&&pt2>20"&&PlottingSetup::passtrig); |
149 |
|
} |
225 |
– |
// // Note: Using same response for bJZB as for cJZB |
226 |
– |
// if(PlottingSetup::DoBTag) { |
227 |
– |
// // Values for 5.1 /fb (2012, bJZB) |
228 |
– |
// PlottingSetup::jzbvariabledata="(jzb[1]+0.092*pt)"; |
229 |
– |
// PlottingSetup::jzbvariablemc="(jzb[1]+0.040*pt)"; |
230 |
– |
// } |
231 |
– |
|
232 |
– |
|
233 |
– |
SingleTopSCrossSection=5.55;//Top mass: m(top)=173GeV, Default scales: muF=muR=mt, PDF: MSTW2008 NNLO |
234 |
– |
SingleTopTCrossSection=87.1;//Top mass: m(top)=173GeV, Default scales: muF=muR=mt, PDF: MSTW2008 NNLO |
235 |
– |
SingleTopUCrossSection=22.2;//Top mass: m(top)=173GeV, Default scales: muF=muR=mt, PDF: MSTW2008 NNLO |
236 |
– |
|
237 |
– |
// allsamples.AddSample("/shome/buchmann/ntuples/Data/METdataset_5100ipb.root","Data",1,1,true,false,0,kBlack); |
238 |
– |
allsamples.AddSample("/shome/buchmann/ntuples/Data/Data_5100ipb_ZbReady.root","Data",1,1,true,false,0,kBlack); |
239 |
– |
// allsamples.AddSample("/scratch/buchmann/2012ABrereco.root","Data",1,1,true,false,0,kBlack); |
240 |
– |
|
241 |
– |
|
242 |
– |
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); |
243 |
– |
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); |
244 |
– |
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); |
245 |
– |
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); |
246 |
– |
|
247 |
– |
// 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); |
248 |
– |
// 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); |
249 |
– |
|
250 |
– |
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 |
251 |
– |
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 |
252 |
– |
|
253 |
– |
// SINGLE-TOP |
254 |
– |
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); |
255 |
– |
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); |
256 |
– |
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); |
257 |
– |
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); |
258 |
– |
|
259 |
– |
// DIBOSONS: LO with a k-factor |
260 |
– |
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); |
261 |
– |
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); |
262 |
– |
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); |
263 |
– |
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); |
264 |
– |
|
265 |
– |
//ttV |
266 |
– |
//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); |
267 |
– |
//allsamples.AddSample("/shome/fronga/ntuples/8TeV/TTGJets_8TeV-madgraph-Summer12-PU_S7_START52_V9-v1.root", "ttV",0,1.444,false,false,5,kGreen); |
268 |
– |
//allsamples.AddSample("/shome/fronga/ntuples/8TeV/TTWJets_8TeV-madgraph-Summer12-PU_S7_START52_V9-v1.root", "ttV",0,0.2149,false,false,5,kGreen); |
269 |
– |
|
270 |
– |
|
271 |
– |
// 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); |
272 |
– |
|
273 |
– |
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); |
150 |
|
|
151 |
+ |
if ( !PlottingSetup::is53reco ) { |
152 |
+ |
//52 reco samples |
153 |
+ |
LoadSamplesFromJSON("Samples52.samples",SplitByJetFlavor,SplitByLeptonFlavor,SplitDibosons,SplitOthers); |
154 |
+ |
} else { |
155 |
+ |
LoadSamplesFromJSON("Samples53.samples",SplitByJetFlavor,SplitByLeptonFlavor,SplitDibosons,SplitOthers); |
156 |
+ |
} |
157 |
|
} |
158 |
|
|
159 |
|
if(PlottingSetup::is2012) { |
160 |
< |
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); |
279 |
< |
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); |
280 |
< |
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); |
281 |
< |
// 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); |
282 |
< |
cutWeight=TCut("(weight*(weight<1000)*(is_data+(!is_data)*((id1==id2)*0.95+(id1!=id2)*0.94)))"); |
283 |
< |
//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)))"); |
284 |
< |
//cutWeight=TCut("(weight*(weight<1000)*((id1==id2&&id1==0)*0.95+(id1==id2&&id1==1)*0.88+(id1!=id2)*0.92))"); |
285 |
< |
} |
286 |
< |
|
287 |
< |
|
288 |
< |
// 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); |
289 |
< |
|
290 |
< |
if(!PlottingSetup::is2012) { |
291 |
< |
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); |
292 |
< |
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 |
293 |
< |
} else { |
294 |
< |
|
295 |
< |
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); |
296 |
< |
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); |
297 |
< |
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 |
298 |
< |
|
160 |
> |
cutWeight=TCut("(weight*(weight<1000)*(is_data+(!is_data)*((id1==id2)*0.95+(id1!=id2)*0.94)))"); |
161 |
|
} |
162 |
|
|
301 |
– |
|
163 |
|
// XS taken from: |
164 |
|
// https://twiki.cern.ch/twiki/bin/view/CMS/QCDGammaPlusJetQCD11005 |
165 |
|
|
257 |
|
allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/MCSamples5p1/SplitUp/Zcc_Full_DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball-Summer12-PU_S7_START52_V9-v2.root","Z+c",0,ZJetsCrossSection,false,false,9,lm_color); |
258 |
|
allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/MCSamples5p1/SplitUp/Zll_Full_DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball-Summer12-PU_S7_START52_V9-v2.root","Z+l",0,ZJetsCrossSection,false,false,10,qcd_color); |
259 |
|
|
399 |
– |
|
260 |
|
allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/MCSamples5p1/TTJets_TuneZ2star_8TeV-madgraph-tauola-Summer12-PU_S7_START52_V9-v1.root","t#bar{t}",0,TTbarCrossSection,false,false,1,ttbar_color); |
261 |
|
allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/MCSamples5p1/WJetsToLNu_TuneZ2Star_8TeV-madgraph-tarball-Summer12-PU_S7_START52_V9-v1.root","W+Jets",0,36257.2,false,false,3,wjets_color); |
262 |
|
|
275 |
|
TCut bTagWeight=TCut("((EventFlavor!=5&&EventFlavor!=4)+(EventFlavor==5||EventFlavor==4)*(((Zb30_bTagProbCSVBP[0]>0.15&&Zb30_bTagProbCSVBP[0]<0.95)*(-0.113472343605*Zb30_bTagProbCSVBP[0] + 1.04926963159)) + ((Zb30_bTagProbCSVBP[0]>1.2&&Zb30_bTagProbCSVBP[0]<20.0)*(0.00066342906116*Zb30_bTagProbCSVBP[0] + 0.93334455507)) + (Zb30_bTagProbCSVBP[0]<0.15||Zb30_bTagProbCSVBP[0]>20) + (Zb30_bTagProbCSVBP[0]>0.95&&Zb30_bTagProbCSVBP[0]<1.2)))"); |
276 |
|
|
277 |
|
|
278 |
< |
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)))"*bTagWeight); |
278 |
> |
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)))"*bTagWeight); |
279 |
> |
write_error(__FUNCTION__,"Trigger efficiency off!!!!"); |
280 |
|
if(showList) { |
281 |
|
allsamples.ListSamples(); |
282 |
|
} |