Plotting/Modules/Setup.C

#include <iostream>
#include <fstream>
#include <TCut.h>
#include <TColor.h>
#include <TStyle.h>

#ifndef SampleClassLoaded
#include "SampleClass.C"
#endif
#define SetupLoaded
#ifndef ResultLibraryClassLoaded
#include "ResultLibraryClass.C"
#endif


using namespace std;


namespace PlottingSetup {

  string directoryname="Results_With_4980ipb__Step8___AllInOne___PAPERMODE2";
  
  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 :-)
  
  float luminosity=4980.0;//4653.74;//p3523.18;//2096.0;//3172.73;//2096.0;//1936;//751.0;//486.0;//468.0//336.;//pb^{-1}
//  float luminosity=4950.;//4653.74;//p3523.18;//2096.0;//3172.73;//2096.0;//1936;//751.0;//486.0;//468.0//336.;//pb^{-1}
//  float luminosity=3523.18;//2096.0;//3172.73;//2096.0;//1936;//751.0;//486.0;//468.0//336.;//pb^{-1}
//  float luminosity=2096.0;//3172.73;//2096.0;//1936;//751.0;//486.0;//468.0//336.;//pb^{-1}
  float lumiuncert=0.022;//0.022;// to be indicated in [0,1] range, e.g. for 4% write 0.04
//  float lumiuncert=0.045;// to be indicated in [0,1] range, e.g. for 4% write 0.04

  string jzbvariabledata="jzb[1]+0.06*pt";
  //string jzbvariablemc="jzb[1]+0.04*pt-2.44029"; // For efficiency curve
  string jzbvariablemc="jzb[1]+0.04*pt";
  float jzbHigh = 400.; // Range for JZB plots

  samplecollection allsamples("completesamplecollection");
  samplecollection signalsamples("signalsamplecollection");
  samplecollection systsamples("systematicssamplecollection");
  samplecollection scansample("scansamplecollection");
  samplecollection raresample("raresamplecollection");
  int data=1;
  int mc=0;
  int mcwithsignal=2;
  TCut passtrig("(passed_triggers||!is_data)");
  TCut cutmass("abs(mll-91.2)<20");
  TCut genMassCut("abs(genMll-91.2)<20");
  TCut openmasscut("mll>40"); // this is the mass cut used in the off peak analysis!
  TCut openGenmasscut("genMll>40");
  //TCut cutmass("mll>2");
  TCut basiccut("mll>2");//basically nothing.
  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.
  //TCut jetqualitycut("(pfJetGoodNum>=2&&pfJetGoodID[0])&&(pfJetGoodNum>=2&&pfJetGoodID[1])");//now part of the basiccut
  //TCut jetqualitycut("mll>0");
  
  TCut cutnJets("pfJetGoodNum>=3"&&basicqualitycut);
  TCut cutnJetsJESdown("pfJetGoodNum25>=3"&&basicqualitycut);
  TCut cutnJetsJESup("pfJetGoodNum35>=3"&&basicqualitycut);
  TCut cutOSOF("(id1!=id2)&&(ch1*ch2<0)");
  TCut cutOSSF("(id1==id2)&&(ch1*ch2<0)");
  TCut sidebandcut("(mll>55&&mll<70)||(mll>112&&mll<160)");

  // SUSY scan parameters
  float mglustart=25;float mgluend=1200;float mglustep=25;
  float mLSPstart=25;float mLSPend=1200;float mLSPstep=25;
  
  float m0start=20;  float m0end=2000;  float m0step=20;
  float m12start=20; float m12end=760;  float m12step=20;
  
  int ScanXzones=15; // number of zones in x for (mSUGRA) scans
  int ScanYzones=15; // number of zones in y for (mSUGRA) scans

  //scan types:
  int mSUGRA=1;
  int SMS=2;
  int GMSB=3;

  TCut essential(passtrig);//add here any cuts you ALWAYS want
  int dogaus=0;
  int doKM=1;
  int dogaus2sigma=2;
  int dogaus3sigma=3;
  int Kostasmethod=-99;
  
  float fitresultconstdata=0;//this is the result when fitting in the 0-30 GeV range
  float fitresultconstmc=0;//this is the result when fitting in the 0-30 GeV range
  
  int method=-1;//Fitting method
  
  //now some style issues:
  float DataMarkerSize=1.2;
  
  //here we save our number of predicted and observed events (with errors)
  vector<float> Nobs;
  vector<float> Npred;
  vector<float> Nprederr;
  
  //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
  bool alwaysflip=false;
  vector<float> flippedNobs;
  vector<float> flippedNpred;
  vector<float> flippedNprederr;
  
  int noJES=0;
  int JESdown=1;
  int JESup=2;
  
  //some refinement: nicer color gradient
  Double_t stops[5] = { 0.00, 0.34, 0.61, 0.84, 1.00 };
  Double_t red[5]   = { 0.00, 0.00, 0.87, 1.00, 0.51 };
  Double_t green[5] = { 0.00, 0.81, 1.00, 0.20, 0.00 };
  Double_t blue[5]  = { 0.51, 1.00, 0.12, 0.00, 0.00 };
  int fi=TColor::CreateGradientColorTable(5, stops, red, green,blue, 255);
  
  // LIMITS
  int nlimittoys=1000; // how many toys for setting limits
  string limitmethod="cls";//what method to use to set limits
  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
  int limitpatienceCRAB=60;
  bool ConsiderSignalContaminationForLimits=true; //whether or not to consider signal contamination when computing limits (standard:true)
  int nuisancemodel=1;
  
  float JZBPeakPositionData=-999;
  float JZBPeakPositionMC=-999;
  float JZBPeakWidthData=-999;
  float JZBPeakWidthMC=-999;

  // two possible future systematics that ATM only take up CPU time in SUSY scans
  bool computeJZBefficiency=false;
  bool computeJZBresponse=false;

  bool requireZ=true;

  string ScanSampleDirectory="SMS_T5zzl/"; // possibilities (atm) : SMS_T5zz/, SMS_T5zzl/, SMS_T5zzh/, GMSB/, mSUGRA/     note: this string needs to either contain "SMS", "GMSB", or "mSUGRA"

  ResultLibrary allresults;
}
Revision:	1.67
Committed:	Wed Apr 18 12:23:42 2012 UTC (13 years ago) by fronga
Content type:	text/plain
Branch:	MAIN
CVS Tags:	cbaf_4_98ifb_paper, HEAD
Changes since 1.66:	+2 -1 lines
Log Message:	Update for paper
#	Content
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
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="Results_With_4980ipb__Step8___AllInOne___PAPERMODE2";
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=4980.0;//4653.74;//p3523.18;//2096.0;//3172.73;//2096.0;//1936;//751.0;//486.0;//468.0//336.;//pb^{-1}
26	// float luminosity=4950.;//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}
28	// float luminosity=2096.0;//3172.73;//2096.0;//1936;//751.0;//486.0;//468.0//336.;//pb^{-1}
29	float lumiuncert=0.022;//0.022;// to be indicated in [0,1] range, e.g. for 4% write 0.04
30	// float lumiuncert=0.045;// to be indicated in [0,1] range, e.g. for 4% write 0.04
31
32	string jzbvariabledata="jzb[1]+0.06*pt";
33	//string jzbvariablemc="jzb[1]+0.04*pt-2.44029"; // For efficiency curve
34	string jzbvariablemc="jzb[1]+0.04*pt";
35	float jzbHigh = 400.; // Range for JZB plots
36
37	samplecollection allsamples("completesamplecollection");
38	samplecollection signalsamples("signalsamplecollection");
39	samplecollection systsamples("systematicssamplecollection");
40	samplecollection scansample("scansamplecollection");
41	samplecollection raresample("raresamplecollection");
42	int data=1;
43	int mc=0;
44	int mcwithsignal=2;
45	TCut passtrig("(passed_triggers\|\|!is_data)");
46	TCut cutmass("abs(mll-91.2)<20");
47	TCut genMassCut("abs(genMll-91.2)<20");
48	TCut openmasscut("mll>40"); // this is the mass cut used in the off peak analysis!
49	TCut openGenmasscut("genMll>40");
50	//TCut cutmass("mll>2");
51	TCut basiccut("mll>2");//basically nothing.
52	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.
53	//TCut jetqualitycut("(pfJetGoodNum>=2&&pfJetGoodID[0])&&(pfJetGoodNum>=2&&pfJetGoodID[1])");//now part of the basiccut
54	//TCut jetqualitycut("mll>0");
55
56	TCut cutnJets("pfJetGoodNum>=3"&&basicqualitycut);
57	TCut cutnJetsJESdown("pfJetGoodNum25>=3"&&basicqualitycut);
58	TCut cutnJetsJESup("pfJetGoodNum35>=3"&&basicqualitycut);
59	TCut cutOSOF("(id1!=id2)&&(ch1*ch2<0)");
60	TCut cutOSSF("(id1==id2)&&(ch1*ch2<0)");
61	TCut sidebandcut("(mll>55&&mll<70)\|\|(mll>112&&mll<160)");
62
63	// SUSY scan parameters
64	float mglustart=25;float mgluend=1200;float mglustep=25;
65	float mLSPstart=25;float mLSPend=1200;float mLSPstep=25;
66
67	float m0start=20; float m0end=2000; float m0step=20;
68	float m12start=20; float m12end=760; float m12step=20;
69
70	int ScanXzones=15; // number of zones in x for (mSUGRA) scans
71	int ScanYzones=15; // number of zones in y for (mSUGRA) scans
72
73	//scan types:
74	int mSUGRA=1;
75	int SMS=2;
76	int GMSB=3;
77
78	TCut essential(passtrig);//add here any cuts you ALWAYS want
79	int dogaus=0;
80	int doKM=1;
81	int dogaus2sigma=2;
82	int dogaus3sigma=3;
83	int Kostasmethod=-99;
84
85	float fitresultconstdata=0;//this is the result when fitting in the 0-30 GeV range
86	float fitresultconstmc=0;//this is the result when fitting in the 0-30 GeV range
87
88	int method=-1;//Fitting method
89
90	//now some style issues:
91	float DataMarkerSize=1.2;
92
93	//here we save our number of predicted and observed events (with errors)
94	vector<float> Nobs;
95	vector<float> Npred;
96	vector<float> Nprederr;
97
98	//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
99	bool alwaysflip=false;
100	vector<float> flippedNobs;
101	vector<float> flippedNpred;
102	vector<float> flippedNprederr;
103
104	int noJES=0;
105	int JESdown=1;
106	int JESup=2;
107
108	//some refinement: nicer color gradient
109	Double_t stops[5] = { 0.00, 0.34, 0.61, 0.84, 1.00 };
110	Double_t red[5] = { 0.00, 0.00, 0.87, 1.00, 0.51 };
111	Double_t green[5] = { 0.00, 0.81, 1.00, 0.20, 0.00 };
112	Double_t blue[5] = { 0.51, 1.00, 0.12, 0.00, 0.00 };
113	int fi=TColor::CreateGradientColorTable(5, stops, red, green,blue, 255);
114
115	// LIMITS
116	int nlimittoys=1000; // how many toys for setting limits
117	string limitmethod="cls";//what method to use to set limits
118	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
119	int limitpatienceCRAB=60;
120	bool ConsiderSignalContaminationForLimits=true; //whether or not to consider signal contamination when computing limits (standard:true)
121	int nuisancemodel=1;
122
123	float JZBPeakPositionData=-999;
124	float JZBPeakPositionMC=-999;
125	float JZBPeakWidthData=-999;
126	float JZBPeakWidthMC=-999;
127
128	// two possible future systematics that ATM only take up CPU time in SUSY scans
129	bool computeJZBefficiency=false;
130	bool computeJZBresponse=false;
131
132	bool requireZ=true;
133
134	string ScanSampleDirectory="SMS_T5zzl/"; // possibilities (atm) : SMS_T5zz/, SMS_T5zzl/, SMS_T5zzh/, GMSB/, mSUGRA/ note: this string needs to either contain "SMS", "GMSB", or "mSUGRA"
135
136	ResultLibrary allresults;
137	}