| 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"); |
| 66 |
|
float m0start=200; float m0end=3000; float m0step=20; |
| 67 |
|
float m12start=100; float m12end=1000; float m12step=20; |
| 68 |
|
|
| 69 |
< |
int ScanXzones=15; // number of zones in x for (mSUGRA) scans |
| 70 |
< |
int ScanYzones=15; // number of zones in y for (mSUGRA) scans |
| 69 |
> |
int ScanXzones=10; // number of zones in x for (mSUGRA) scans |
| 70 |
> |
int ScanYzones=10; // number of zones in y for (mSUGRA) scans |
| 71 |
|
|
| 72 |
|
//scan types: |
| 73 |
|
int mSUGRA=1; |
| 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; |
