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 :-) |
24 |
+ |
bool FullMCAnalysis=false; |
25 |
|
|
26 |
|
float luminosity=4980;//4653.74;//p3523.18;//2096.0;//3172.73;//2096.0;//1936;//751.0;//486.0;//468.0//336.;//pb^{-1} |
27 |
|
// float luminosity=3523.18;//2096.0;//3172.73;//2096.0;//1936;//751.0;//486.0;//468.0//336.;//pb^{-1} |
46 |
|
TCut leptoncut("(pt1>20&&pt2>20)"); |
47 |
|
TCut passtrig("(passed_triggers||!is_data)"&&leptoncut); |
48 |
|
|
49 |
< |
TCut openmasscut("mll>12"); // this is the mass cut used in the off peak analysis! |
50 |
< |
TCut openGenmasscut("genMll>12"); |
49 |
> |
TCut openmasscut("mll>20"); // this is the mass cut used in the off peak analysis! |
50 |
> |
TCut openGenmasscut("genMll>20"); |
51 |
|
TCut opensidebandcut("mll>500&&mll<100"); // this won't let anything thru ... |
52 |
|
|
53 |
|
TCut Restrmasscut("abs(mll-91.2)<20"); // this is the mass cut used in the off peak analysis! |
59 |
|
|
60 |
|
|
61 |
|
//TCut cutmass("mll>2"); |
62 |
< |
TCut basiccut("mll>2"&&leptoncut);//basically nothing. |
62 |
> |
// TCut basiccut("mll>2"&&leptoncut&&"(bTagProbTHighEff[0]>1.7)&&(bTagProbTHighEff[1]>1.7)");//basically nothing. |
63 |
> |
TCut basiccut("mll>2"&&leptoncut);//basically nothing |
64 |
|
TCut basicqualitycut("(pfJetGoodNum>=2&&pfJetGoodID[0]!=0)&&(pfJetGoodNum>=2&&pfJetGoodID[1]!=0)"&&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. |
65 |
|
//TCut jetqualitycut("(pfJetGoodNum>=2&&pfJetGoodID[0])&&(pfJetGoodNum>=2&&pfJetGoodID[1])");//now part of the basiccut |
66 |
|
//TCut jetqualitycut("mll>0"); |
67 |
|
|
68 |
|
TCut cutnJets("pfJetGoodNum>=3"&&basicqualitycut); |
69 |
< |
TCut cutnJetsJESdown("pfJetGoodNum25>=3"&&basicqualitycut); |
70 |
< |
TCut cutnJetsJESup("pfJetGoodNum35>=3"&&basicqualitycut); |
69 |
> |
TCut cutnJetsJESdown("pfJetGoodNumn1sigma>=3"&&basicqualitycut); |
70 |
> |
TCut cutnJetsJESup("pfJetGoodNump1sigma>=3"&&basicqualitycut); |
71 |
|
TCut cutOSOF("(id1!=id2)&&(ch1*ch2<0)"); |
72 |
|
TCut cutOSSF("(id1==id2)&&(ch1*ch2<0)"); |
73 |
|
TCut sidebandcut(Restrsidebandcut); |
82 |
|
int ScanXzones=15; // number of zones in x for (mSUGRA) scans |
83 |
|
int ScanYzones=15; // number of zones in y for (mSUGRA) scans |
84 |
|
|
85 |
< |
string mSUGRAxsFile="../../Plotting/Modules/external/msugra_m0_m12_10_0_1_NLO_1.0.txt"; |
85 |
> |
string mSUGRAxsFile="Plotting/Modules/external/msugra_m0_m12_10_0_1_NLO_1.0.txt"; |
86 |
> |
string SMSReferenceXSFile="Plotting/Modules/external/reference_xSec_SMS-new.root"; |
87 |
|
|
88 |
|
//scan types: |
89 |
|
int mSUGRA=1; |
121 |
|
int JESdown=1; |
122 |
|
int JESup=2; |
123 |
|
|
124 |
+ |
// Uncertainties ONPEAK |
125 |
+ |
float zjetsestimateuncertONPEAK=0.25; |
126 |
+ |
float emuncertONPEAK=0.25; |
127 |
+ |
float emsidebanduncertONPEAK=0.25; |
128 |
+ |
float eemmsidebanduncertONPEAK=0.25; |
129 |
+ |
|
130 |
+ |
// Uncertainties OFFPEAK (iJZB) |
131 |
+ |
float zjetsestimateuncertOFFPEAK=0.25; |
132 |
+ |
float emuncertOFFPEAK=0.1; |
133 |
+ |
float emsidebanduncertOFFPEAK=0.0; |
134 |
+ |
float eemmsidebanduncertOFFPEAK=0.0; |
135 |
+ |
|
136 |
+ |
|
137 |
|
//some refinement: nicer color gradient |
138 |
|
Double_t stops[5] = { 0.00, 0.34, 0.61, 0.84, 1.00 }; |
139 |
|
Double_t red[5] = { 0.00, 0.00, 0.87, 1.00, 0.51 }; |
159 |
|
bool computeJZBresponse=false; |
160 |
|
|
161 |
|
bool requireZ=true;//this is switched to "off" automatically when doing offpeak! |
162 |
+ |
|
163 |
+ |
string ScanSampleDirectory="SMS_T5zzh_newslots";//DileptonmSUGRAScan__pieces |
164 |
|
|
165 |
< |
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" |
165 |
> |
// 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" |
166 |
> |
// 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" |
167 |
|
|
168 |
|
string FilterEfficiencyFile = "/shome/buchmann/JellyfishCBAF/DistributedModelCalculations/Limits/FilterEfficiencyv3.root"; |
169 |
+ |
string CMSSW_dir = "/shome/buchmann/final_production_2011/CMSSW_4_2_8/src/"; |
170 |
|
|
171 |
|
ResultLibrary allresults; |
172 |
|
} |