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;
  
  //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 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.48
Committed:	Mon Nov 7 15:23:24 2011 UTC (13 years, 6 months ago) by buchmann
Content type:	text/plain
Branch:	MAIN
Changes since 1.47:	+1 -0 lines
Log Message:	Added option to use flipping
#	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	buchmann	1.47	//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
90	buchmann	1.48	bool alwaysflip=false;
91	buchmann	1.47	vector<float> flippedNobs;
92			vector<float> flippedNpred;
93			vector<float> flippedNprederr;
94
95	buchmann	1.13	int noJES=0;
96			int JESdown=1;
97			int JESup=2;
98	buchmann	1.14
99	buchmann	1.17	//some refinement: nicer color gradient
100			Double_t stops[5] = { 0.00, 0.34, 0.61, 0.84, 1.00 };
101			Double_t red[5] = { 0.00, 0.00, 0.87, 1.00, 0.51 };
102			Double_t green[5] = { 0.00, 0.81, 1.00, 0.20, 0.00 };
103			Double_t blue[5] = { 0.51, 1.00, 0.12, 0.00, 0.00 };
104	buchmann	1.18	int fi=TColor::CreateGradientColorTable(5, stops, red, green,blue, 255);
105	buchmann	1.17
106	buchmann	1.32	// LIMITS
107	buchmann	1.33	int nlimittoys=1000; // how many toys for setting limits
108	buchmann	1.40	string limitmethod="cls";//what method to use to set limits
109	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.
110	buchmann	1.37	bool ConsiderSignalContaminationForLimits=true; //whether or not to consider signal contamination when computing limits (standard:true)
111	buchmann	1.36	int nuisancemodel=1;
112	buchmann	1.23
113	buchmann	1.28	float JZBPeakPositionData=-999;
114			float JZBPeakPositionMC=-999;
115			float JZBPeakWidthData=-999;
116			float JZBPeakWidthMC=-999;
117	buchmann	1.36
118	buchmann	1.40	// two possible future systematics that ATM only take up CPU time in SUSY scans
119			bool computeJZBefficiency=false;
120			bool computeJZBresponse=false;
121
122	buchmann	1.47
123	buchmann	1.36	//watch out, the cbafbasedir string is in GeneralToolBox
124	buchmann	1.1	}