| 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> |
| 24 |
|
{ |
| 25 |
|
char hostname[1023]; |
| 26 |
|
gethostname(hostname,1023); |
| 27 |
+ |
|
| 28 |
+ |
SanityChecks(); |
| 29 |
|
|
| 30 |
|
/* |
| 31 |
|
|
| 41 |
|
|
| 42 |
|
/// 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) |
| 43 |
|
/// 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) |
| 44 |
< |
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" |
| 44 |
> |
// 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" |
| 45 |
> |
// 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" |
| 46 |
> |
// 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" |
| 47 |
> |
|
| 48 |
> |
|
| 49 |
> |
|
| 50 |
|
|
| 51 |
|
if(loadall) { |
| 52 |
|
for(int ix=0;ix<PlottingSetup::ScanXzones;ix++) { |
| 58 |
|
|
| 59 |
|
|
| 60 |
|
flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak |
| 61 |
< |
if(!PlottingSetup::RestrictToMassPeak) { |
| 61 |
> |
if(PlottingSetup::RestrictToMassPeak==false) { |
| 62 |
|
PlottingSetup::cutmass=PlottingSetup::openmasscut; |
| 63 |
|
// PlottingSetup::sidebandcut=TCut("mll<2&&mll>3");//impossible cut - because we don't want any sidebands! |
| 64 |
|
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). |
| 65 |
|
PlottingSetup::genMassCut=PlottingSetup::openGenmasscut;//impossible cut - because we don't want any sidebands! |
| 66 |
+ |
PlottingSetup::requireZ=false; |
| 67 |
|
} |
| 68 |
+ |
dout << "Currently configuration: The mass cut is \"" << PlottingSetup::cutmass << "\" and RestrictToMassPeak is currently set to " << PlottingSetup::RestrictToMassPeak << endl; |
| 69 |
+ |
|
| 70 |
+ |
|
| 71 |
|
} |
| 72 |
|
|
| 73 |
|
|
| 74 |
< |
void define_samples(bool showList, samplecollection &allsamples, samplecollection &signalsamples, samplecollection &scansample, samplecollection &raresample, samplecollection &systsamples) |
| 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 |
|
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 |
+ |
if(PlottingSetup::RestrictToMassPeak&&!PlottingSetup::UseSidebandsForcJZB) { |
| 86 |
+ |
PlottingSetup::sidebandcut=TCut("mll<2&&mll>3"); |
| 87 |
+ |
} |
| 88 |
+ |
|
| 89 |
+ |
|
| 90 |
+ |
|
| 91 |
+ |
|
| 92 |
|
|
| 93 |
< |
float ZJetsCrossSection = 3048.0; //NNLO---------------------------------- |
| 94 |
< |
// float ZJetsCrossSectionM20 = 4998.0;//summer11 |
| 95 |
< |
// float ZJetsCrossSectionD6Tlow = 400.0;//summer11 |
| 96 |
< |
float TTbarCrossSection = 165.0;//(NLO) ---- 165.0; // approx. NNLO----- |
| 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 |
> |
// 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 |
> |
// 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 |
> |
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 |
> |
LoadSamplesFromJSON("Samples53.samples",SplitByJetFlavor,SplitByLeptonFlavor,SplitDibosons,SplitOthers); |
| 156 |
> |
} |
| 157 |
> |
} |
| 158 |
> |
|
| 159 |
> |
if(PlottingSetup::is2012) { |
| 160 |
> |
cutWeight=TCut("(weight*(weight<1000)*(is_data+(!is_data)*((id1==id2)*0.95+(id1!=id2)*0.94)))"); |
| 161 |
> |
} |
| 162 |
> |
|
| 163 |
> |
// XS taken from: |
| 164 |
> |
// https://twiki.cern.ch/twiki/bin/view/CMS/QCDGammaPlusJetQCD11005 |
| 165 |
> |
|
| 166 |
> |
if(!PlottingSetup::RestrictToMassPeak) { |
| 167 |
> |
PlottingSetup::cutmass=PlottingSetup::openmasscut; |
| 168 |
> |
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). |
| 169 |
> |
PlottingSetup::genMassCut=PlottingSetup::openGenmasscut; |
| 170 |
> |
PlottingSetup::requireZ=false; |
| 171 |
> |
} |
| 172 |
> |
if(PlottingSetup::RestrictToMassPeak&&!PlottingSetup::UseSidebandsForcJZB) { |
| 173 |
> |
PlottingSetup::sidebandcut=TCut("SideBandsHaveBeenCalledForOnPeakAnalysisWithoutSidebands__ERROR==1"); |
| 174 |
> |
PlottingSetup::requireZ=true; |
| 175 |
> |
} |
| 176 |
> |
|
| 177 |
> |
if(PlottingSetup::DoBTag) { |
| 178 |
> |
PlottingSetup::cutnJets=PlottingSetup::cutnJets&&PlottingSetup::bTagRequirement; |
| 179 |
> |
} |
| 180 |
> |
|
| 181 |
> |
|
| 182 |
> |
|
| 183 |
> |
if(showList) { |
| 184 |
> |
allsamples.ListSamples(); |
| 185 |
> |
signalsamples.ListSamples(); |
| 186 |
> |
systsamples.ListSamples(); |
| 187 |
> |
} |
| 188 |
> |
|
| 189 |
> |
} |
| 190 |
> |
|
| 191 |
> |
|
| 192 |
> |
|
| 193 |
> |
void define_Zb_samples(bool showList, samplecollection &allsamples, samplecollection &signalsamples, samplecollection &scansample, samplecollection &raresample, samplecollection &systsamples, samplecollection &qcdsamples) |
| 194 |
> |
{ |
| 195 |
> |
SanityChecks(); |
| 196 |
> |
|
| 197 |
> |
float ZJetsCrossSection = 3503.71; |
| 198 |
> |
float LowZJetsCrossSection = 11050*0.069*1.15; // LO xs * filter eff * k-factor (same as ZJets) |
| 199 |
> |
|
| 200 |
> |
float TTbarCrossSection = 225.197; |
| 201 |
|
float WJetsCrossSection = 31314.0;//NNLO-------3.131e4; //NNLO------------ |
| 202 |
|
float ZnunuCrossSection = 5760.0;//NNLO -------4.5e+3; //(LO);------------ |
| 92 |
– |
//Single top x-sections from here: |
| 93 |
– |
//https://twiki.cern.ch/twiki/bin/view/CMS/SingleTopSigma |
| 203 |
|
float SingleTopSCrossSection = 4.63; // NLO;---------------------------------- |
| 204 |
|
float SingleTopTCrossSection = 64.57; // NLO;-------------------------------- |
| 205 |
|
float SingleTopUCrossSection = 15.74; // NLO;---------------------------------- |
| 206 |
|
float VVJetsCrossSection = 4.8; // LO;------------------------------------- |
| 98 |
– |
float LM4CrossSection = 1.879*1.35; // k*LO |
| 99 |
– |
float LM8CrossSection = 0.7300*1.41 ; // k*LO |
| 100 |
– |
//the following numbers are from the MadGraphStandardModel210Summary as linked on the GeneratorProduction2011 page |
| 101 |
– |
// float QCD50to100CrossSection=30000000; // not used |
| 102 |
– |
float QCD100to250CrossSection=7000000; |
| 103 |
– |
float QCD250to500CrossSection=171000; |
| 104 |
– |
float QCD500to1000CrossSection=5200; |
| 105 |
– |
float QCD1000toInfCrossSection=83.0; |
| 106 |
– |
|
| 107 |
– |
// float totEventsQCD50to100=207418.0; // not used |
| 108 |
– |
long totEventsQCD100to250 = 638792; |
| 109 |
– |
long totEventsQCD250to500 = 344454; |
| 110 |
– |
long totEventsQCD500to1000 = (long)10e10; //unknown but now suppressed. (no event passes anyway) |
| 111 |
– |
long totEventsQCD1000toInf = 163185; |
| 207 |
|
|
| 208 |
|
long totEventsZjets = 32005720;//summer11 |
| 209 |
< |
long totEventsZjetsD6T = 2543706;//summer11 |
| 115 |
< |
long totEventsZjetsD6Tlow = 181799;//summer11 |
| 116 |
< |
// long totEventsZjetsEE = 2738801; |
| 117 |
< |
// long totEventsZjetsMM = 2148325; |
| 118 |
< |
// long totEventsZjetsTT = 2032536; |
| 209 |
> |
long totEventslowZjets = 10272470; |
| 210 |
|
|
| 120 |
– |
//spring:2313911.0;//2743142.0; // DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola |
| 211 |
|
long totEventsTTbar = 3701947;//summer11 |
| 122 |
– |
//1161621.0;//Spring//1144028.0; // PabloV13/TTJets_TuneZ2_7TeV-madgraph-tauola |
| 212 |
|
long totEventsWJets = 56674902;//summer11 |
| 124 |
– |
//15010237.0;//spring//14350756.0; /// WJetsToLNu_TuneZ2_7TeV-madgraph-tauola.root |
| 213 |
|
long totEventsZnunu = 2106977;//2167964.0; //ZinvisibleJets_7TeV-madgraph.root |
| 214 |
|
long totEventsVVJets = 959076;//509072.0; |
| 215 |
|
long totEventsSingleTopS = 493868;//489472.0; |
| 216 |
|
long totEventsSingleTopT = 475460;//477610.0; |
| 217 |
|
long totEventsSingleTopU = 489417;//477599.0; |
| 130 |
– |
long totEventsLM4 = 218380; |
| 131 |
– |
long totEventsLM8 = 180960; |
| 132 |
– |
|
| 218 |
|
|
| 219 |
|
Int_t nice_blue = TColor::GetColor("#2E9AFE"); |
| 220 |
|
Int_t nice_green = TColor::GetColor("#81f781"); |
| 136 |
– |
Int_t nice_red = TColor::GetColor("#F78181"); |
| 221 |
|
Int_t nice_pink = TColor::GetColor("#F781BE"); |
| 222 |
|
Int_t nice_orange= TColor::GetColor("#F7BE81");//#FA9624 |
| 223 |
|
Int_t nice_black = TColor::GetColor("#2A0A1B");//black with a bit of purple :-) |
| 224 |
|
|
| 141 |
– |
|
| 225 |
|
Int_t dy_color = nice_blue; |
| 226 |
|
Int_t wjets_color = kGray; |
| 227 |
|
Int_t ttbar_color = nice_green; |
| 238 |
|
diboson_color = kGreen+2; |
| 239 |
|
lm_color = kViolet+7; |
| 240 |
|
|
| 158 |
– |
Int_t rare_color = nice_blue; |
| 159 |
– |
|
| 241 |
|
char hostname[1023]; |
| 242 |
|
gethostname(hostname,1023); |
| 243 |
|
|
| 244 |
|
string prefix="/shome/"; |
| 245 |
|
if((int)((string)hostname).find("falda")>-1) prefix="/scratch/"; |
| 246 |
|
|
| 247 |
< |
allsamples.AddSample(prefix+"buchmann/ntuples/Data/AllData_Complete_4653_ipb_PAPER.root","Data",1,1,true,false,0,kBlack); |
| 248 |
< |
|
| 249 |
< |
|
| 250 |
< |
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); |
| 251 |
< |
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 |
| 252 |
< |
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 |
| 253 |
< |
allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/VVJetsTo4L_TuneD6T_7TeV-madgraph-tauola.root","WW/WZ/ZZ",totEventsVVJets,VVJetsCrossSection,false,false,4,diboson_color); |
| 254 |
< |
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); |
| 174 |
< |
allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/TToBLNu_TuneZ2_t-channel_7TeV-madgraph.root","Single-top",totEventsSingleTopT,SingleTopTCrossSection,false,false,2,singletop_color); |
| 175 |
< |
allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/TToBLNu_TuneZ2_tW-channel_7TeV-madgraph.root","Single-top",totEventsSingleTopU,SingleTopUCrossSection,false,false,2,singletop_color); |
| 176 |
< |
allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/ZinvisibleJets_7TeV-madgraph.root","Z nunu",totEventsZnunu,ZnunuCrossSection,false,false,4,diboson_color); |
| 177 |
< |
allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/LM4_SUSY_sftsht_7TeV-pythia6.root","LM4",totEventsLM4,LM4CrossSection,false,true,8,lm_color); |
| 178 |
< |
|
| 179 |
< |
|
| 180 |
< |
|
| 181 |
< |
if(!PlottingSetup::RestrictToMassPeak) { |
| 182 |
< |
// signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM0_SUSY_sftsht_7TeV-pythia6.root","LM0",219796,38.93*1.41 ,false,true,3,lm_color); |
| 183 |
< |
// signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM1_SUSY_sftsht_7TeV-pythia6.root","LM1",218176,4.888*1.34,false,true,4,lm_color); |
| 184 |
< |
// signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM2_SUSY_sftsht_7TeV-pythia6.root","LM2",206336,0.6027*1.33,false,true,5,lm_color); |
| 185 |
< |
signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM3_SUSY_sftsht_7TeV-pythia6.root","LM3",215120,3.438*1.40,false,true,6,lm_color); |
| 186 |
< |
signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM5_SUSY_sftsht_7TeV-pythia6.root","LM5",223992,0.473*1.34,false,true,7,lm_color); |
| 187 |
< |
signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM9_SUSY_sftsht_7TeV-pythia6.root","LM9",227808,7.134*1.48,false,true,8,lm_color); |
| 188 |
< |
signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM11_SUSY_sftsht_7TeV-pythia6.root","LM11",218703,0.8236*1.35,false,true,9,lm_color); |
| 189 |
< |
signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM12_SUSY_sftsht_7TeV-pythia6.root","LM12",223491,4.414*1.34,false,true,10,lm_color); |
| 190 |
< |
signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM13_SUSY_sftsht_7TeV-pythia6.root","LM13",173240,6.899*1.42,false,true,11,lm_color); |
| 191 |
< |
} |
| 192 |
< |
signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM4_SUSY_sftsht_7TeV-pythia6.root","LM4",totEventsLM4,LM4CrossSection,false,true,1,lm_color); |
| 193 |
< |
signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM8_SUSY_sftsht_7TeV-pythia6.root","LM8",totEventsLM8,LM8CrossSection,false,true,2,lm_color); |
| 194 |
< |
|
| 195 |
< |
// 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); |
| 196 |
< |
|
| 197 |
< |
|
| 198 |
< |
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); |
| 199 |
< |
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 |
| 200 |
< |
|
| 201 |
< |
// |
| 202 |
< |
/* |
| 203 |
< |
|
| 204 |
< |
|
| 205 |
< |
// We do not use the following (rare) samples; nonetheless we include them in a separate category to prove that they are negligible |
| 206 |
< |
write_warning(__FUNCTION__,"Rare samples not yet reprocessed with new selectioN");/* |
| 207 |
< |
raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTWWTo2Lplus2Nu_NTP.root","Rare SS",39791,4.76E-5,false,false,6,rare_color);//summer11 |
| 208 |
< |
raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/WWminusTo2L2Nu_NTP.root","Rare SS",51063,0.001559,false,false,6,rare_color);//summer11 |
| 209 |
< |
raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTWWTo2Lminus2Nu_NTP.root","Rare SS",39496 ,4.719E-5,false,false,6,rare_color);//summer11 |
| 210 |
< |
raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTWTo2Lplus2Nu_NTP.root","Rare SS",48764,0.006841,false,false,6,rare_color);//summer11 |
| 211 |
< |
raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/WWplusTo2L2Nu_NTP.root","Rare SS",48244 ,0.004216,false,false,6,rare_color);//summer11 |
| 212 |
< |
raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTZTo2Lplus2Nu_NTP.root","Rare SS",36740,0.002024,false,false,6,rare_color);//summer11 |
| 213 |
< |
raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/WWWTo2Lplus2Nu_NTP.root","Rare SS",43322,0.008957,false,false,6,rare_color);//summer11 |
| 214 |
< |
raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTWTo2Lminus2Nu_NTP.root","Rare SS",50435,0.002705,false,false,6,rare_color);//summer11 |
| 215 |
< |
raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/WWWTo2Lminus2Nu_NTP.root","Rare SS",17147,0.004109,false,false,6,rare_color);//summer11 |
| 216 |
< |
raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTZTo2Lminus2Nu_NTP.root","Rare SS",36781,0.001946,false,false,6,rare_color);//summer11 |
| 217 |
< |
*/ |
| 218 |
< |
//scan files are only defined when necessary! |
| 247 |
> |
PlottingSetup::luminosity=PlottingSetup::luminosity2012; |
| 248 |
> |
PlottingSetup::lumiuncert=PlottingSetup::lumiuncert2012; |
| 249 |
> |
|
| 250 |
> |
SingleTopSCrossSection=5.55;//Top mass: m(top)=173GeV, Default scales: muF=muR=mt, PDF: MSTW2008 NNLO |
| 251 |
> |
SingleTopTCrossSection=87.1;//Top mass: m(top)=173GeV, Default scales: muF=muR=mt, PDF: MSTW2008 NNLO |
| 252 |
> |
SingleTopUCrossSection=22.2;//Top mass: m(top)=173GeV, Default scales: muF=muR=mt, PDF: MSTW2008 NNLO |
| 253 |
> |
|
| 254 |
> |
allsamples.AddSample("/shome/buchmann/ntuples/Data/Data_5100ipb_ZbReady.root","Data",1,1,true,false,0,kBlack); |
| 255 |
|
|
| 256 |
< |
if(showList) allsamples.ListSamples(); |
| 256 |
> |
allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/MCSamples5p1/SplitUp/Zbb_Full_DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball-Summer12-PU_S7_START52_V9-v2.root","Z+b",0,ZJetsCrossSection,false,false,8,dy_color); |
| 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 |
> |
|
| 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 |
> |
|
| 263 |
> |
// SINGLE-TOP |
| 264 |
> |
allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/MCSamples5p1/T_s-channel_TuneZ2star_8TeV-powheg-tauola-Summer12-PU_S7_START52_V9-v1.root","Single-top",0,1.76,false,false,2,singletop_color); |
| 265 |
> |
allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/MCSamples5p1/Tbar_s-channel_TuneZ2star_8TeV-powheg-tauola-Summer12-PU_S7_START52_V9-v1.root","Single-top",0,3.79,false,false,2,singletop_color); |
| 266 |
> |
allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/MCSamples5p1/Tbar_t-channel_TuneZ2star_8TeV-powheg-tauola-Summer12-PU_S7_START52_V9-v1.root","Single-top",0,56.4,false,false,2,singletop_color); |
| 267 |
> |
allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/MCSamples5p1/Tbar_tW-channel-DR_TuneZ2star_8TeV-powheg-tauola-Summer12-PU_S7_START52_V9-v1.root","Single-top",0,SingleTopUCrossSection,false,false,2,singletop_color); |
| 268 |
> |
|
| 269 |
> |
// DIBOSONS: LO with a k-factor |
| 270 |
> |
allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/MCSamples5p1/ZZJetsTo4L_TuneZ2star_8TeV-madgraph-tauola-Summer12-PU_S7_START52_V9-v3.root","Dibosons",0,0.181,false,false,4,diboson_color); |
| 271 |
> |
|
| 272 |
> |
allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/MCSamples5p1/WZJetsTo3LNu_TuneZ2_8TeV-madgraph-tauola-Summer12-PU_S7_START52_V9-v2.root","Dibosons",0,1.06,false,false,4,diboson_color); |
| 273 |
|
|
| 274 |
+ |
// cutWeight=TCut("(weight*Efficiencyweightonly*(weight<1000)*(is_data+(!is_data)*((id1==id2&&id1==0)*0.95+(id1==id2&&id1==1)*0.88+(id1!=id2)*0.92)))"); |
| 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); |
| 279 |
+ |
write_error(__FUNCTION__,"Trigger efficiency off!!!!"); |
| 280 |
+ |
if(showList) { |
| 281 |
+ |
allsamples.ListSamples(); |
| 282 |
+ |
} |
| 283 |
|
} |
