ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/cbrown/AnalysisFramework/Plotting/Modules/Setup.C

Comparing UserCode/cbrown/AnalysisFramework/Plotting/Modules/Setup.C (file contents):
Revision 1.17 by buchmann, Tue Jul 26 13:08:54 2011 UTC vs.
Revision 1.59 by buchmann, Wed Nov 30 10:33:22 2011 UTC

#	Line 1 | Line 1
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	<	float luminosity=1078;//751.0;//486.0;//468.0//336.;//pb^{-1}
21	<	float lumiuncert=0.04;// to be indicated in [0,1] range, e.g. for 4% write 0.04
20	>
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=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
29	>
30		string jzbvariabledata="jzb[1]+0.06*pt";
31		string jzbvariablemc="jzb[1]+0.04*pt";
32	<
32	>	float jzbHigh = 350.; // Range for JZB plots
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;
40		int mc=0;
41		int mcwithsignal=2;
42		TCut passtrig("(passed_triggers||!is_data)");
43		TCut cutmass("abs(mll-91.2)<20");
44	+	TCut genMassCut("abs(genMll-91.2)<20");
45	+	TCut openmasscut("mll>40"); // this is the mass cut used in the off peak analysis!
46	+	TCut openGenmasscut("genMll>40");
47		//TCut cutmass("mll>2");
48		TCut basiccut("mll>2");//basically nothing.
49		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.
#	Line 36 | Line 56 | namespace PlottingSetup {
56		TCut cutOSOF("(id1!=id2)&&(ch1*ch2<0)");
57		TCut cutOSSF("(id1==id2)&&(ch1*ch2<0)");
58		TCut sidebandcut("(mll>55&&mll<70)||(mll>112&&mll<160)");
59	+
60	+	TCut Weight("weight"); // this is the weight that is used everywhere
61	+
62		//TCut sidebandcut("(mll>61&&mll<70)||(mll>112&&mll<190)");
63		//TCut basiccut("(passed_triggers||!is_data)");
64
65		// SUSY scan parameters
43	–	//  float mglustart=25;float mgluend=1200;float mglustep=25; //guessed values for official file
44	–	//  float mLSPstart=25;float mLSPend=1200;float mLSPstep=25; //guessed values for official file
66		float mglustart=25;float mgluend=1200;float mglustep=25;
67		float mLSPstart=25;float mLSPend=1200;float mLSPstep=25;
47	–
68
69	+	float m0start=20;  float m0end=2000;  float m0step=20;
70	+	float m12start=20; float m12end=760;  float m12step=20;
71	+
72	+	int ScanXzones=15; // number of zones in x for (mSUGRA) scans
73	+	int ScanYzones=15; // number of zones in y for (mSUGRA) scans
74	+
75	+	//scan types:
76	+	int mSUGRA=1;
77	+	int SMS=2;
78	+	int GMSB=3;
79	+
80		TCut essential(passtrig);//add here any cuts you ALWAYS want
81		int dogaus=0;
82		int doKM=1;
#	Line 56 | Line 87 | namespace PlottingSetup {
87		float fitresultconstdata=0;//this is the result when fitting in the 0-30 GeV range
88		float fitresultconstmc=0;//this is the result when fitting in the 0-30 GeV range
89
90	<	int method=-1;
90	>	int method=-1;//Fitting method
91
92		//now some style issues:
93		float DataMarkerSize=1.2;
#	Line 66 | Line 97 | namespace PlottingSetup {
97		vector<float> Npred;
98		vector<float> Nprederr;
99
100	+	//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
101	+	bool alwaysflip=false;
102	+	vector<float> flippedNobs;
103	+	vector<float> flippedNpred;
104	+	vector<float> flippedNprederr;
105	+
106		int noJES=0;
107		int JESdown=1;
108		int JESup=2;
#	Line 75 | Line 112 | namespace PlottingSetup {
112		Double_t red[5]   = { 0.00, 0.00, 0.87, 1.00, 0.51 };
113		Double_t green[5] = { 0.00, 0.81, 1.00, 0.20, 0.00 };
114		Double_t blue[5]  = { 0.51, 1.00, 0.12, 0.00, 0.00 };
115	<	TColor::CreateGradientColorTable(5, stops, red, green,blue, 255);
115	>	int fi=TColor::CreateGradientColorTable(5, stops, red, green,blue, 255);
116
117	+	// LIMITS
118	+	int nlimittoys=1000; // how many toys for setting limits
119	+	string limitmethod="cls";//what method to use to set limits
120	+	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
121	+	int limitpatienceCRAB=60;
122	+	bool ConsiderSignalContaminationForLimits=true; //whether or not to consider signal contamination when computing limits (standard:true)
123	+	int nuisancemodel=1;
124	+
125	+	float JZBPeakPositionData=-999;
126	+	float JZBPeakPositionMC=-999;
127	+	float JZBPeakWidthData=-999;
128	+	float JZBPeakWidthMC=-999;
129	+
130	+	// two possible future systematics that ATM only take up CPU time in SUSY scans
131	+	bool computeJZBefficiency=false;
132	+	bool computeJZBresponse=false;
133	+
134	+	bool requireZ=true;
135	+
136	+	string ScanSampleDirectory="GMSB/"; // possibilities (atm) : SMS/, SMS_T5zzl/, GMSB/, mSUGRA/     note: this string needs to either contain "SMS", "GMSB", or "mSUGRA"
137	+
138	+	ResultLibrary allresults;
139		}

Diff Legend

–	Removed lines
+	Added lines
<	Changed lines
>	Changed lines