18 |
|
|
19 |
|
namespace PlottingSetup { |
20 |
|
|
21 |
< |
string directoryname="Results_With_4653ipb__Step8___AllInOne___PAPERMODE"; |
21 |
> |
string directoryname="Jellyfish_tests"; |
22 |
|
|
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 |
|
|
32 |
|
float jzbHigh = 400.; // Range for JZB plots |
33 |
|
|
34 |
|
samplecollection allsamples("completesamplecollection"); |
35 |
+ |
samplecollection qcdsamples("QCDcollection"); |
36 |
|
samplecollection signalsamples("signalsamplecollection"); |
37 |
|
samplecollection systsamples("systematicssamplecollection"); |
38 |
|
samplecollection scansample("scansamplecollection"); |
41 |
|
int mc=0; |
42 |
|
int mcwithsignal=2; |
43 |
|
TCut passtrig("(passed_triggers||!is_data)"); |
44 |
+ |
TCut leptoncut("(pt1>20&&pt2>20)"); |
45 |
|
TCut cutmass("abs(mll-91.2)<20"); |
46 |
|
TCut genMassCut("abs(genMll-91.2)<20"); |
47 |
|
TCut openmasscut("mll>40"); // this is the mass cut used in the off peak analysis! |
48 |
|
TCut openGenmasscut("genMll>40"); |
49 |
|
//TCut cutmass("mll>2"); |
50 |
< |
TCut basiccut("mll>2");//basically nothing. |
50 |
> |
TCut basiccut("mll>2"&&leptoncut);//basically nothing. |
51 |
|
TCut basicqualitycut("(pfJetGoodNum>=2&&pfJetGoodID[0]!=0)&&(pfJetGoodNum>=2&&pfJetGoodID[1]!=0)");//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. |
52 |
|
//TCut jetqualitycut("(pfJetGoodNum>=2&&pfJetGoodID[0])&&(pfJetGoodNum>=2&&pfJetGoodID[1])");//now part of the basiccut |
53 |
|
//TCut jetqualitycut("mll>0"); |
96 |
|
|
97 |
|
//here we save our "flipped" number of predicted and observed events (with errors) -- this means that we consider JZB<-X as observed and construct the corresponding prediction |
98 |
|
bool alwaysflip=false; |
99 |
+ |
bool allowflipping=false; |
100 |
|
vector<float> flippedNobs; |
101 |
|
vector<float> flippedNpred; |
102 |
|
vector<float> flippedNprederr; |