| 8 |
|
#include "SampleClass.C" |
| 9 |
|
#endif |
| 10 |
|
#define SetupLoaded |
| 11 |
+ |
#ifndef ResultLibraryClassLoaded |
| 12 |
+ |
#include "ResultLibraryClass.C" |
| 13 |
+ |
#endif |
| 14 |
+ |
|
| 15 |
|
|
| 16 |
|
using namespace std; |
| 17 |
|
|
| 18 |
|
|
| 19 |
|
namespace PlottingSetup { |
| 20 |
|
|
| 21 |
< |
// string directoryname="official_2096ipb___forPASv5_withPDF"; |
| 18 |
< |
string directoryname="official_3523ipb___forAN_prepaper_NewSelection"; |
| 21 |
> |
string directoryname="Results_With_4653ipb__Step1"; |
| 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 |
|
|
| 25 |
< |
float luminosity=3523.18;//2096.0;//3172.73;//2096.0;//1936;//751.0;//486.0;//468.0//336.;//pb^{-1} |
| 25 |
> |
float luminosity=4653.74;//p3523.18;//2096.0;//3172.73;//2096.0;//1936;//751.0;//486.0;//468.0//336.;//pb^{-1} |
| 26 |
|
// float luminosity=3523.18;//2096.0;//3172.73;//2096.0;//1936;//751.0;//486.0;//468.0//336.;//pb^{-1} |
| 27 |
|
// float luminosity=2096.0;//3172.73;//2096.0;//1936;//751.0;//486.0;//468.0//336.;//pb^{-1} |
| 28 |
|
float lumiuncert=0.045;// to be indicated in [0,1] range, e.g. for 4% write 0.04 |
| 33 |
|
|
| 34 |
|
samplecollection allsamples("completesamplecollection"); |
| 35 |
|
samplecollection signalsamples("signalsamplecollection"); |
| 36 |
+ |
samplecollection systsamples("systematicssamplecollection"); |
| 37 |
|
samplecollection scansample("scansamplecollection"); |
| 38 |
|
samplecollection raresample("raresamplecollection"); |
| 39 |
|
int data=1; |
| 57 |
|
TCut cutOSSF("(id1==id2)&&(ch1*ch2<0)"); |
| 58 |
|
TCut sidebandcut("(mll>55&&mll<70)||(mll>112&&mll<160)"); |
| 59 |
|
|
| 56 |
– |
//TCut sidebandcut("(mll>61&&mll<70)||(mll>112&&mll<190)"); |
| 57 |
– |
//TCut basiccut("(passed_triggers||!is_data)"); |
| 58 |
– |
|
| 60 |
|
// SUSY scan parameters |
| 61 |
|
float mglustart=25;float mgluend=1200;float mglustep=25; |
| 62 |
|
float mLSPstart=25;float mLSPend=1200;float mLSPstep=25; |
| 64 |
|
float m0start=20; float m0end=2000; float m0step=20; |
| 65 |
|
float m12start=20; float m12end=760; float m12step=20; |
| 66 |
|
|
| 67 |
< |
int ScanXzones=10; // number of zones in x for (mSUGRA) scans |
| 68 |
< |
int ScanYzones=10; // number of zones in y for (mSUGRA) scans |
| 67 |
> |
int ScanXzones=15; // number of zones in x for (mSUGRA) scans |
| 68 |
> |
int ScanYzones=15; // number of zones in y for (mSUGRA) scans |
| 69 |
> |
|
| 70 |
> |
//scan types: |
| 71 |
> |
int mSUGRA=1; |
| 72 |
> |
int SMS=2; |
| 73 |
> |
int GMSB=3; |
| 74 |
|
|
| 75 |
|
TCut essential(passtrig);//add here any cuts you ALWAYS want |
| 76 |
|
int dogaus=0; |
| 93 |
|
vector<float> Nprederr; |
| 94 |
|
|
| 95 |
|
//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 |
| 96 |
+ |
bool alwaysflip=false; |
| 97 |
|
vector<float> flippedNobs; |
| 98 |
|
vector<float> flippedNpred; |
| 99 |
|
vector<float> flippedNprederr; |
| 112 |
|
// LIMITS |
| 113 |
|
int nlimittoys=1000; // how many toys for setting limits |
| 114 |
|
string limitmethod="cls";//what method to use to set limits |
| 115 |
< |
int limitpatience=15; // for how many minutes should the limit calculation (for one configuration!) be allowed to run before being aborted? this only has an effect when when running on the grid. |
| 115 |
> |
int limitpatience=15; // for how many minutes should the limit calculation (for one configuration!) be allowed to run before being aborted? this only has an effect when when running on the t3 batch; the value for CRAB is on the next line |
| 116 |
> |
int limitpatienceCRAB=60; |
| 117 |
|
bool ConsiderSignalContaminationForLimits=true; //whether or not to consider signal contamination when computing limits (standard:true) |
| 118 |
|
int nuisancemodel=1; |
| 119 |
|
|
| 126 |
|
bool computeJZBefficiency=false; |
| 127 |
|
bool computeJZBresponse=false; |
| 128 |
|
|
| 129 |
+ |
bool requireZ=true; |
| 130 |
+ |
|
| 131 |
+ |
string ScanSampleDirectory="GMSB/"; // possibilities (atm) : SMS/, SMS_T5zzl/, GMSB/, mSUGRA/ note: this string needs to either contain "SMS", "GMSB", or "mSUGRA" |
| 132 |
|
|
| 133 |
< |
//watch out, the cbafbasedir string is in GeneralToolBox |
| 133 |
> |
ResultLibrary allresults; |
| 134 |
|
} |
