20 |
|
|
21 |
|
string directoryname="Jellyfish_tests"; |
22 |
|
|
23 |
< |
bool RestrictToMassPeak=false; //if you want to switch between offpeak ("false") and onpeak ("true") analysis please use this switch; the masscut below will be adapted automatically when adding samples :-) |
23 |
> |
bool RestrictToMassPeak=true; //if you want to switch between offpeak ("false") and onpeak ("true") analysis please use this switch; the masscut below will be adapted automatically when adding samples :-) |
24 |
> |
bool FullMCAnalysis=false; |
25 |
> |
bool DoBTag=false; |
26 |
|
|
27 |
+ |
// the 2012 switch is in GeneralToolBox |
28 |
+ |
|
29 |
|
float luminosity=4980;//4653.74;//p3523.18;//2096.0;//3172.73;//2096.0;//1936;//751.0;//486.0;//468.0//336.;//pb^{-1} |
30 |
|
// float luminosity=3523.18;//2096.0;//3172.73;//2096.0;//1936;//751.0;//486.0;//468.0//336.;//pb^{-1} |
31 |
|
// float luminosity=2096.0;//3172.73;//2096.0;//1936;//751.0;//486.0;//468.0//336.;//pb^{-1} |
32 |
|
float lumiuncert=0.022;//0.045;// to be indicated in [0,1] range, e.g. for 4% write 0.04 |
33 |
+ |
|
34 |
+ |
float luminosity2012=5100; |
35 |
+ |
float lumiuncert2012=0.022; |
36 |
|
|
37 |
|
// string jzbvariabledata="jzb[1]+0.06*pt"; |
38 |
|
// string jzbvariablemc="jzb[1]+0.04*pt"; |
39 |
< |
string jzbvariabledata="(jzb[1]+0.053*pt)"; |
40 |
< |
string jzbvariablemc="(jzb[1]+0.043*pt)"; |
41 |
< |
float jzbHigh = 400.; // Range for JZB plots |
39 |
> |
string jzbvariabledata="(jzb[1]+0.062*pt)"; |
40 |
> |
string jzbvariablemc="(jzb[1]+0.036*pt)"; |
41 |
> |
float jzbHigh = 450.; // Range for JZB plots |
42 |
> |
float iMllLow = 20.; // Range for Edge Fitting (and mll plot) for iJZB |
43 |
> |
float iMllHigh = 400.; // Range for Edge Fitting (and mll plot) for iJZB |
44 |
|
|
45 |
|
samplecollection allsamples("completesamplecollection"); |
46 |
|
samplecollection qcdsamples("QCDcollection"); |
51 |
|
int data=1; |
52 |
|
int mc=0; |
53 |
|
int mcwithsignal=2; |
54 |
< |
TCut leptoncut("(pt1>20&&pt2>20)"); |
54 |
> |
TCut leptoncut("(abs(eta1)<1.4 && abs(eta2)<1.4 && pt1>20&&pt2>10)"); |
55 |
|
TCut passtrig("(passed_triggers||!is_data)"&&leptoncut); |
56 |
|
|
57 |
|
TCut openmasscut("mll>20"); // this is the mass cut used in the off peak analysis! |
65 |
|
TCut cutmass(Restrmasscut); |
66 |
|
TCut genMassCut(RestrGenmasscut); |
67 |
|
|
68 |
+ |
TCut bTagRequirement("bTagProbCSVBP[0]>0.679"); |
69 |
|
|
70 |
|
//TCut cutmass("mll>2"); |
71 |
|
// TCut basiccut("mll>2"&&leptoncut&&"(bTagProbTHighEff[0]>1.7)&&(bTagProbTHighEff[1]>1.7)");//basically nothing. |
75 |
|
//TCut jetqualitycut("mll>0"); |
76 |
|
|
77 |
|
TCut cutnJets("pfJetGoodNum>=3"&&basicqualitycut); |
78 |
< |
TCut cutnJetsJESdown("pfJetGoodNumm1sigma>=3"&&basicqualitycut); |
78 |
> |
TCut cutnJetsJESdown("pfJetGoodNumn1sigma>=3"&&basicqualitycut); |
79 |
|
TCut cutnJetsJESup("pfJetGoodNump1sigma>=3"&&basicqualitycut); |
80 |
|
TCut cutOSOF("(id1!=id2)&&(ch1*ch2<0)"); |
81 |
|
TCut cutOSSF("(id1==id2)&&(ch1*ch2<0)"); |
91 |
|
int ScanXzones=15; // number of zones in x for (mSUGRA) scans |
92 |
|
int ScanYzones=15; // number of zones in y for (mSUGRA) scans |
93 |
|
|
94 |
< |
string mSUGRAxsFile="../../Plotting/Modules/external/msugra_m0_m12_10_0_1_NLO_1.0.txt"; |
94 |
> |
string mSUGRAxsFile="Plotting/Modules/external/msugra_m0_m12_10_0_1_NLO_1.0.txt"; |
95 |
> |
string SMSReferenceXSFile="Plotting/Modules/external/reference_xSec_SMS-new.root"; |
96 |
|
|
97 |
|
//scan types: |
98 |
|
int mSUGRA=1; |
130 |
|
int JESdown=1; |
131 |
|
int JESup=2; |
132 |
|
|
133 |
+ |
// Uncertainties ONPEAK |
134 |
+ |
float zjetsestimateuncertONPEAK=0.25; |
135 |
+ |
float emuncertONPEAK=0.25; |
136 |
+ |
float emsidebanduncertONPEAK=0.25; |
137 |
+ |
float eemmsidebanduncertONPEAK=0.25; |
138 |
+ |
|
139 |
+ |
// Uncertainties OFFPEAK (iJZB) |
140 |
+ |
float zjetsestimateuncertOFFPEAK=0.25; |
141 |
+ |
float emuncertOFFPEAK=0.1; |
142 |
+ |
float emsidebanduncertOFFPEAK=0.0; |
143 |
+ |
float eemmsidebanduncertOFFPEAK=0.0; |
144 |
+ |
|
145 |
+ |
|
146 |
|
//some refinement: nicer color gradient |
147 |
|
Double_t stops[5] = { 0.00, 0.34, 0.61, 0.84, 1.00 }; |
148 |
|
Double_t red[5] = { 0.00, 0.00, 0.87, 1.00, 0.51 }; |
168 |
|
bool computeJZBresponse=false; |
169 |
|
|
170 |
|
bool requireZ=true;//this is switched to "off" automatically when doing offpeak! |
171 |
+ |
|
172 |
+ |
// string ScanSampleDirectory="SMS_T5zzh_newslots";//DileptonmSUGRAScan__pieces |
173 |
|
|
174 |
< |
// string ScanSampleDirectory="DileptonScan__pieces"; // possibilities (atm) : SMS_T5zz/, SMS_T5zzl/, SMS_T5zzh/, GMSB/, mSUGRA/ note: this string needs to either contain "SMS", "GMSB", or "mSUGRA" |
175 |
< |
string ScanSampleDirectory="SMS_T1lh"; // possibilities (atm) : SMS_T5zz/, SMS_T5zzl/, SMS_T5zzh/, GMSB/, mSUGRA/ note: this string needs to either contain "SMS", "GMSB", or "mSUGRA" |
174 |
> |
string ScanSampleDirectory="DileptonmSUGRAScan__pieces"; // possibilities (atm) : SMS_T5zz/, SMS_T5zzl/, SMS_T5zzh/, GMSB/, mSUGRA/ note: this string needs to either contain "SMS", "GMSB", or "mSUGRA" |
175 |
> |
// string ScanSampleDirectory="SMS_T1lh"; // possibilities (atm) : SMS_T5zz/, SMS_T5zzl/, SMS_T5zzh/, GMSB/, mSUGRA/ note: this string needs to either contain "SMS", "GMSB", or "mSUGRA" |
176 |
|
|
177 |
|
string FilterEfficiencyFile = "/shome/buchmann/JellyfishCBAF/DistributedModelCalculations/Limits/FilterEfficiencyv3.root"; |
178 |
|
string CMSSW_dir = "/shome/buchmann/final_production_2011/CMSSW_4_2_8/src/"; |