1 |
|
#include <iostream> |
2 |
+ |
#include <fstream> |
3 |
|
#include <TCut.h> |
4 |
+ |
#include <TColor.h> |
5 |
+ |
#include <TStyle.h> |
6 |
+ |
|
7 |
|
#ifndef SampleClassLoaded |
8 |
|
#include "SampleClass.C" |
9 |
|
#endif |
14 |
|
|
15 |
|
namespace PlottingSetup { |
16 |
|
float luminosity=1078;//751.0;//486.0;//468.0//336.;//pb^{-1} |
17 |
+ |
float lumiuncert=0.04;// to be indicated in [0,1] range, e.g. for 4% write 0.04 |
18 |
|
string jzbvariabledata="jzb[1]+0.06*pt"; |
19 |
|
string jzbvariablemc="jzb[1]+0.04*pt"; |
20 |
+ |
float jzbHigh = 350.; // Range for JZB plots |
21 |
|
|
22 |
|
samplecollection allsamples("completesamplecollection"); |
23 |
+ |
samplecollection signalsamples("signalsamplecollection"); |
24 |
+ |
samplecollection scansample("scansamplecollection"); |
25 |
|
int data=1; |
26 |
|
int mc=0; |
27 |
|
int mcwithsignal=2; |
39 |
|
TCut cutOSOF("(id1!=id2)&&(ch1*ch2<0)"); |
40 |
|
TCut cutOSSF("(id1==id2)&&(ch1*ch2<0)"); |
41 |
|
TCut sidebandcut("(mll>55&&mll<70)||(mll>112&&mll<160)"); |
42 |
+ |
|
43 |
|
//TCut sidebandcut("(mll>61&&mll<70)||(mll>112&&mll<190)"); |
44 |
|
//TCut basiccut("(passed_triggers||!is_data)"); |
45 |
|
|
46 |
+ |
// SUSY scan parameters |
47 |
+ |
// float mglustart=25;float mgluend=1200;float mglustep=25; //guessed values for official file |
48 |
+ |
// float mLSPstart=25;float mLSPend=1200;float mLSPstep=25; //guessed values for official file |
49 |
+ |
float mglustart=25;float mgluend=1200;float mglustep=25; |
50 |
+ |
float mLSPstart=25;float mLSPend=1200;float mLSPstep=25; |
51 |
+ |
|
52 |
+ |
|
53 |
|
TCut essential(passtrig);//add here any cuts you ALWAYS want |
54 |
|
int dogaus=0; |
55 |
|
int doKM=1; |
57 |
|
int dogaus3sigma=3; |
58 |
|
int Kostasmethod=-99; |
59 |
|
|
60 |
< |
float fitresultconstdata=0;//this is the result when fitting in the 0-30 GeV range! |
61 |
< |
float fitresultconstmc=0;//this is the result when fitting in the 0-30 GeV range! |
60 |
> |
float fitresultconstdata=0;//this is the result when fitting in the 0-30 GeV range |
61 |
> |
float fitresultconstmc=0;//this is the result when fitting in the 0-30 GeV range |
62 |
|
|
63 |
|
int method=-1; |
64 |
+ |
|
65 |
+ |
//now some style issues: |
66 |
+ |
float DataMarkerSize=1.2; |
67 |
+ |
|
68 |
+ |
//here we save our number of predicted and observed events (with errors) |
69 |
+ |
vector<float> Nobs; |
70 |
+ |
vector<float> Npred; |
71 |
+ |
vector<float> Nprederr; |
72 |
+ |
|
73 |
+ |
int noJES=0; |
74 |
+ |
int JESdown=1; |
75 |
+ |
int JESup=2; |
76 |
+ |
|
77 |
+ |
//some refinement: nicer color gradient |
78 |
+ |
Double_t stops[5] = { 0.00, 0.34, 0.61, 0.84, 1.00 }; |
79 |
+ |
Double_t red[5] = { 0.00, 0.00, 0.87, 1.00, 0.51 }; |
80 |
+ |
Double_t green[5] = { 0.00, 0.81, 1.00, 0.20, 0.00 }; |
81 |
+ |
Double_t blue[5] = { 0.51, 1.00, 0.12, 0.00, 0.00 }; |
82 |
+ |
int fi=TColor::CreateGradientColorTable(5, stops, red, green,blue, 255); |
83 |
+ |
|
84 |
+ |
int limitpatience=50; // 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. |
85 |
+ |
|
86 |
|
} |