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

using namespace std;


namespace PlottingSetup {

//  string directoryname="official_2096ipb___forPASv5_withPDF";
  string directoryname="official_3523ipb___forAN_prepaper_NewSelection";
  
  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=3523.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.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";
  float jzbHigh = 350.; // Range for JZB plots

  samplecollection allsamples("completesamplecollection");
  samplecollection signalsamples("signalsamplecollection");
  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)");

  //TCut sidebandcut("(mll>61&&mll<70)||(mll>112&&mll<190)");
  //TCut basiccut("(passed_triggers||!is_data)");
  
  // 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=10; // number of zones in x for (mSUGRA) scans
  int ScanYzones=10; // number of zones in y for (mSUGRA) scans

  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;
  
  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 grid.
  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;

  //watch out, the cbafbasedir string is in GeneralToolBox
}
Revision:	1.46
Committed:	Mon Oct 24 15:05:37 2011 UTC (13 years, 6 months ago) by buchmann
Content type:	text/plain
Branch:	MAIN
Changes since 1.45:	+3 -3 lines
Log Message:	Upgrade from HoneyPot to IceCreamBowl: merged in offpeak stuff, different warning color for frederic, saving to rootfile, only 3 attempts when computing limits and much, much more
#	User	Rev	Content
1	buchmann	1.1	#include <iostream>
2	buchmann	1.11	#include <fstream>
3	buchmann	1.1	#include <TCut.h>
4	buchmann	1.18	#include <TColor.h>
5			#include <TStyle.h>
6	buchmann	1.11
7	buchmann	1.1	#ifndef SampleClassLoaded
8			#include "SampleClass.C"
9			#endif
10	buchmann	1.5	#define SetupLoaded
11	buchmann	1.1
12			using namespace std;
13
14
15			namespace PlottingSetup {
16	buchmann	1.34
17	buchmann	1.45	// string directoryname="official_2096ipb___forPASv5_withPDF";
18			string directoryname="official_3523ipb___forAN_prepaper_NewSelection";
19	buchmann	1.39
20	buchmann	1.43	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 :-)
21	buchmann	1.39
22	buchmann	1.45	float luminosity=3523.18;//2096.0;//3172.73;//2096.0;//1936;//751.0;//486.0;//468.0//336.;//pb^{-1}
23	buchmann	1.44	// float luminosity=3523.18;//2096.0;//3172.73;//2096.0;//1936;//751.0;//486.0;//468.0//336.;//pb^{-1}
24	buchmann	1.45	// float luminosity=2096.0;//3172.73;//2096.0;//1936;//751.0;//486.0;//468.0//336.;//pb^{-1}
25	buchmann	1.34	float lumiuncert=0.045;// to be indicated in [0,1] range, e.g. for 4% write 0.04
26	buchmann	1.36
27	buchmann	1.4	string jzbvariabledata="jzb[1]+0.06*pt";
28			string jzbvariablemc="jzb[1]+0.04*pt";
29	fronga	1.25	float jzbHigh = 350.; // Range for JZB plots
30	buchmann	1.34
31	buchmann	1.1	samplecollection allsamples("completesamplecollection");
32	buchmann	1.9	samplecollection signalsamples("signalsamplecollection");
33	buchmann	1.14	samplecollection scansample("scansamplecollection");
34	buchmann	1.41	samplecollection raresample("raresamplecollection");
35	buchmann	1.1	int data=1;
36			int mc=0;
37			int mcwithsignal=2;
38			TCut passtrig("(passed_triggers\|\|!is_data)");
39	buchmann	1.20	TCut cutmass("abs(mll-91.2)<20");
40	buchmann	1.44	TCut genMassCut("abs(genMll-91.2)<20");
41	buchmann	1.46	TCut openmasscut("mll>40"); // this is the mass cut used in the off peak analysis!
42			TCut openGenmasscut("genMll>40");
43	buchmann	1.1	//TCut cutmass("mll>2");
44			TCut basiccut("mll>2");//basically nothing.
45	buchmann	1.7	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.
46	buchmann	1.1	//TCut jetqualitycut("(pfJetGoodNum>=2&&pfJetGoodID[0])&&(pfJetGoodNum>=2&&pfJetGoodID[1])");//now part of the basiccut
47			//TCut jetqualitycut("mll>0");
48	buchmann	1.7
49	buchmann	1.20	TCut cutnJets("pfJetGoodNum>=3"&&basicqualitycut);
50	buchmann	1.7	TCut cutnJetsJESdown("pfJetGoodNum25>=3"&&basicqualitycut);
51			TCut cutnJetsJESup("pfJetGoodNum35>=3"&&basicqualitycut);
52	buchmann	1.1	TCut cutOSOF("(id1!=id2)&&(ch1*ch2<0)");
53			TCut cutOSSF("(id1==id2)&&(ch1*ch2<0)");
54	buchmann	1.22	TCut sidebandcut("(mll>55&&mll<70)\|\|(mll>112&&mll<160)");
55	fronga	1.25
56	buchmann	1.22	//TCut sidebandcut("(mll>61&&mll<70)\|\|(mll>112&&mll<190)");
57	buchmann	1.1	//TCut basiccut("(passed_triggers\|\|!is_data)");
58
59	buchmann	1.15	// SUSY scan parameters
60	buchmann	1.16	float mglustart=25;float mgluend=1200;float mglustep=25;
61			float mLSPstart=25;float mLSPend=1200;float mLSPstep=25;
62	buchmann	1.31
63			float m0start=20; float m0end=2000; float m0step=20;
64			float m12start=20; float m12end=760; float m12step=20;
65	buchmann	1.42
66			int ScanXzones=10; // number of zones in x for (mSUGRA) scans
67			int ScanYzones=10; // number of zones in y for (mSUGRA) scans
68	buchmann	1.15
69	buchmann	1.1	TCut essential(passtrig);//add here any cuts you ALWAYS want
70			int dogaus=0;
71			int doKM=1;
72			int dogaus2sigma=2;
73			int dogaus3sigma=3;
74			int Kostasmethod=-99;
75
76	buchmann	1.13	float fitresultconstdata=0;//this is the result when fitting in the 0-30 GeV range
77			float fitresultconstmc=0;//this is the result when fitting in the 0-30 GeV range
78	buchmann	1.1
79	buchmann	1.33	int method=-1;//Fitting method
80	buchmann	1.8
81			//now some style issues:
82			float DataMarkerSize=1.2;
83	buchmann	1.10
84			//here we save our number of predicted and observed events (with errors)
85			vector<float> Nobs;
86			vector<float> Npred;
87			vector<float> Nprederr;
88	buchmann	1.13
89			int noJES=0;
90			int JESdown=1;
91			int JESup=2;
92	buchmann	1.14
93	buchmann	1.17	//some refinement: nicer color gradient
94			Double_t stops[5] = { 0.00, 0.34, 0.61, 0.84, 1.00 };
95			Double_t red[5] = { 0.00, 0.00, 0.87, 1.00, 0.51 };
96			Double_t green[5] = { 0.00, 0.81, 1.00, 0.20, 0.00 };
97			Double_t blue[5] = { 0.51, 1.00, 0.12, 0.00, 0.00 };
98	buchmann	1.18	int fi=TColor::CreateGradientColorTable(5, stops, red, green,blue, 255);
99	buchmann	1.17
100	buchmann	1.32	// LIMITS
101	buchmann	1.33	int nlimittoys=1000; // how many toys for setting limits
102	buchmann	1.40	string limitmethod="cls";//what method to use to set limits
103	buchmann	1.46	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.
104	buchmann	1.37	bool ConsiderSignalContaminationForLimits=true; //whether or not to consider signal contamination when computing limits (standard:true)
105	buchmann	1.36	int nuisancemodel=1;
106	buchmann	1.23
107	buchmann	1.28	float JZBPeakPositionData=-999;
108			float JZBPeakPositionMC=-999;
109			float JZBPeakWidthData=-999;
110			float JZBPeakWidthMC=-999;
111	buchmann	1.36
112	buchmann	1.40	// two possible future systematics that ATM only take up CPU time in SUSY scans
113			bool computeJZBefficiency=false;
114			bool computeJZBresponse=false;
115
116	buchmann	1.36	//watch out, the cbafbasedir string is in GeneralToolBox
117	buchmann	1.1	}