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 |
< |
if (PlottingSetup::is2012) LowZJetsCrossSection = 11050*0.069*1.15; // LO xs * filter eff * k-factor (same as ZJets) |
103 |
< |
|
104 |
< |
float TTbarCrossSection = 165.0;//(NLO) ---- 165.0; // approx. NNLO----- |
105 |
< |
if(PlottingSetup::is2012) TTbarCrossSection=225.197; |
106 |
< |
// at 8 TeV: 225.197 (https://twiki.cern.ch/twiki/bin/view/CMS/Zpt8TeV) |
107 |
< |
float WJetsCrossSection = 31314.0;//NNLO-------3.131e4; //NNLO------------ |
108 |
< |
float ZnunuCrossSection = 5760.0;//NNLO -------4.5e+3; //(LO);------------ |
109 |
< |
//Single top x-sections from here: |
110 |
< |
//https://twiki.cern.ch/twiki/bin/view/CMS/SingleTopSigma |
111 |
< |
float SingleTopSCrossSection = 4.63; // NLO;---------------------------------- |
112 |
< |
float SingleTopTCrossSection = 64.57; // NLO;-------------------------------- |
113 |
< |
float SingleTopUCrossSection = 15.74; // NLO;---------------------------------- |
114 |
< |
float VVJetsCrossSection = 4.8; // LO;------------------------------------- |
115 |
< |
float LM3CrossSection = 0.0; // 2011: don't know the x-section |
116 |
< |
if(PlottingSetup::is2012) LM3CrossSection=5.21678; //(only LO, no k factor) |
117 |
< |
float LM4CrossSection = 1.879*1.35; // k*LO |
118 |
< |
if(PlottingSetup::is2012) LM4CrossSection=2.88199; //(only LO, no k factor) |
119 |
< |
float LM6CrossSection = 0.3104*1.30; // k*LO |
120 |
< |
if(PlottingSetup::is2012) LM6CrossSection=0.501773;//(only LO, no k factor) |
121 |
< |
float LM8CrossSection = 0.7300*1.41; // k*LO |
122 |
< |
if(PlottingSetup::is2012) LM8CrossSection=1.14; |
123 |
< |
|
124 |
< |
long totEventsZjets = 32005720;//summer11 |
125 |
< |
long totEventslowZjets = 10272470; |
126 |
< |
|
127 |
< |
long totEventsTTbar = 3701947;//summer11 |
128 |
< |
long totEventsWJets = 56674902;//summer11 |
129 |
< |
long totEventsZnunu = 2106977;//2167964.0; //ZinvisibleJets_7TeV-madgraph.root |
130 |
< |
long totEventsVVJets = 959076;//509072.0; |
131 |
< |
long totEventsSingleTopS = 493868;//489472.0; |
132 |
< |
long totEventsSingleTopT = 475460;//477610.0; |
133 |
< |
long totEventsSingleTopU = 489417;//477599.0; |
134 |
< |
|
135 |
< |
Int_t nice_blue = TColor::GetColor("#2E9AFE"); |
136 |
< |
Int_t nice_green = TColor::GetColor("#81f781"); |
137 |
< |
// Int_t nice_red = TColor::GetColor("#F78181"); |
138 |
< |
Int_t nice_pink = TColor::GetColor("#F781BE"); |
139 |
< |
Int_t nice_orange= TColor::GetColor("#F7BE81");//#FA9624 |
140 |
< |
Int_t nice_black = TColor::GetColor("#2A0A1B");//black with a bit of purple :-) |
141 |
< |
|
142 |
< |
|
143 |
< |
Int_t dy_color = nice_blue; |
144 |
< |
Int_t wjets_color = kGray; |
145 |
< |
Int_t ttbar_color = nice_green; |
146 |
< |
Int_t singletop_color = nice_orange; |
147 |
< |
Int_t qcd_color = nice_black; |
148 |
< |
Int_t diboson_color = nice_pink; |
149 |
< |
Int_t lm_color = kRed; |
150 |
< |
|
151 |
< |
dy_color = kYellow; |
152 |
< |
wjets_color = kGray; |
153 |
< |
ttbar_color = kMagenta+2; |
154 |
< |
singletop_color = kBlue; |
155 |
< |
qcd_color = kPink; |
156 |
< |
diboson_color = kGreen+2; |
157 |
< |
lm_color = kViolet+7; |
158 |
< |
|
159 |
< |
// Int_t rare_color = nice_blue; |
160 |
< |
|
161 |
< |
char hostname[1023]; |
162 |
< |
gethostname(hostname,1023); |
85 |
> |
if(PlottingSetup::RestrictToMassPeak&&!PlottingSetup::UseSidebandsForcJZB) { |
86 |
> |
PlottingSetup::sidebandcut=TCut("mll<2&&mll>3"); |
87 |
> |
} |
88 |
|
|
164 |
– |
string prefix="/shome/"; |
165 |
– |
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); |
182 |
< |
|
183 |
< |
// allsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_200__M12_440.root","LM8",10000,0.132,false,true,8,lm_color); |
184 |
< |
|
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 |
< |
|
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 |
|
|
194 |
– |
// For Kostas' control region |
195 |
– |
// 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. |
196 |
– |
// PlottingSetup::cutnJets=TCut("pfJetGoodNum40<3"&&PlottingSetup::basicqualitycut); |
197 |
– |
|
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 |
< |
|
206 |
< |
|
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)"; |
210 |
< |
PlottingSetup::jzbvariablemc="(jzb[1]+0.034665*pt)"; |
211 |
< |
|
142 |
> |
PlottingSetup::jzbvariabledata="(jzb[1]+0.061016*pt)"; |
143 |
> |
PlottingSetup::jzbvariablemc="(jzb[1]+0.035264*pt)"; |
144 |
|
} |
145 |
|
|
214 |
– |
// // Note: Using same response for bJZB as for cJZB |
215 |
– |
// if(PlottingSetup::DoBTag) { |
216 |
– |
// // Values for 5.1 /fb (2012, bJZB) |
217 |
– |
// PlottingSetup::jzbvariabledata="(jzb[1]+0.092*pt)"; |
218 |
– |
// PlottingSetup::jzbvariablemc="(jzb[1]+0.040*pt)"; |
219 |
– |
// } |
146 |
|
|
147 |
+ |
if(PlottingSetup::RestrictToMassPeak) { |
148 |
+ |
PlottingSetup::passtrig=TCut("pt1>20&&pt2>20"&&PlottingSetup::passtrig); |
149 |
+ |
} |
150 |
|
|
151 |
< |
SingleTopSCrossSection=5.55;//Top mass: m(top)=173GeV, Default scales: muF=muR=mt, PDF: MSTW2008 NNLO |
152 |
< |
SingleTopTCrossSection=87.1;//Top mass: m(top)=173GeV, Default scales: muF=muR=mt, PDF: MSTW2008 NNLO |
153 |
< |
SingleTopUCrossSection=22.2;//Top mass: m(top)=173GeV, Default scales: muF=muR=mt, PDF: MSTW2008 NNLO |
154 |
< |
|
155 |
< |
|
156 |
< |
|
157 |
< |
|
229 |
< |
allsamples.AddSample("/shome/buchmann/ntuples/Data/Data_5100ipb.root","Data",1,1,true,false,0,kBlack); |
230 |
< |
|
231 |
< |
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,7,dy_color); |
232 |
< |
allsamples.AddSample("/shome/fronga/ntuples/8TeV/DYJetsToLLNoTau_M-10To50filter_8TeV-madgraph-Summer12-PU_S7_START52_V9-v1.root","Z+Jets",-6955271,LowZJetsCrossSection,false,false,7,dy_color); |
233 |
< |
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,8,kRed-4); |
234 |
< |
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,8,kRed-4); |
235 |
< |
|
236 |
< |
// 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); |
237 |
< |
// 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); |
238 |
< |
|
239 |
< |
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 |
240 |
< |
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 |
241 |
< |
// allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/VVJetsTo4L_TuneD6T_7TeV-madgraph-tauola.root","WW/WZ/ZZ",totEventsVVJets,VVJetsCrossSection,false,false,4,diboson_color); |
242 |
< |
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); |
243 |
< |
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); |
244 |
< |
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); |
245 |
< |
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); |
246 |
< |
|
247 |
< |
// DIBOSONS: LO with a k-factor |
248 |
< |
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); |
249 |
< |
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); |
250 |
< |
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); |
251 |
< |
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); |
252 |
< |
|
253 |
< |
//ttZ: irrelevant allsamples.AddSample("/shome/fronga/ntuples/8TeV/TTZJets_8TeV-madgraph_v2-Summer12-PU_S7_START52_V9-v1.root","ttZ",0,0.172,false,false,5,kOrange-2); |
254 |
< |
|
255 |
< |
|
256 |
< |
// allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/ZinvisibleJets_7TeV-madgraph.root","Z nunu",totEventsZnunu,ZnunuCrossSection,false,false,4,diboson_color); |
257 |
< |
// 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); |
258 |
< |
|
259 |
< |
|
260 |
< |
|
261 |
< |
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); |
262 |
< |
|
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 |
|
} |
159 |
|
|
160 |
|
if(PlottingSetup::is2012) { |
161 |
< |
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); |
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); |
268 |
< |
signalsamples.AddSample("/shome/fronga/ntuples/8TeV/PYTHIA6_SUSY_LM3_sftsht_8TeV_GEN_FASTSIM-2012_Startup_inTimeOnly-V3-leo.root","LM3",0,LM3CrossSection,false,true,3,lm_color); |
269 |
< |
// 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); |
270 |
< |
cutWeight=TCut("(weight*(weight<1000)*(is_data+(!is_data)*((id1==id2)*0.95+(id1!=id2)*0.94)))"); |
271 |
< |
//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)))"); |
272 |
< |
//cutWeight=TCut("(weight*(weight<1000)*((id1==id2&&id1==0)*0.95+(id1==id2&&id1==1)*0.88+(id1!=id2)*0.92))"); |
273 |
< |
} |
274 |
< |
|
275 |
< |
|
276 |
< |
// 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); |
277 |
< |
|
278 |
< |
if(!PlottingSetup::is2012) { |
279 |
< |
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); |
280 |
< |
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 |
281 |
< |
} else { |
282 |
< |
|
283 |
< |
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); |
284 |
< |
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); |
285 |
< |
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 |
286 |
< |
|
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 |
|
} |
164 |
|
|
289 |
– |
|
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 |
|
} |
304 |
– |
|
305 |
– |
|
182 |
|
|
183 |
|
if(showList) { |
184 |
|
allsamples.ListSamples(); |
185 |
|
signalsamples.ListSamples(); |
186 |
|
systsamples.ListSamples(); |
187 |
|
} |
312 |
– |
|
188 |
|
} |