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_2096";

  float luminosity=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");
  int data=1;
  int mc=0;
  int mcwithsignal=2;
  TCut passtrig("(passed_triggers||!is_data)");
  TCut cutmass("abs(mll-91.2)<20");
  //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; //guessed values for official file
//  float mLSPstart=25;float mLSPend=1200;float mLSPstep=25; //guessed values for official file
  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;

  
  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="bayesian";//what method to use to set limits
  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.
  bool ConsiderSignalContaminationForLimits=false; //whether or not to consider signal contamination when computing limits (standard:true)
  
  float JZBPeakPositionData=-999;
  float JZBPeakPositionMC=-999;
  float JZBPeakWidthData=-999;
  float JZBPeakWidthMC=-999;
  
}
Revision:	1.35
Committed:	Mon Sep 5 22:03:59 2011 UTC (13 years, 8 months ago) by buchmann
Content type:	text/plain
Branch:	MAIN
Changes since 1.34:	+2 -2 lines
Log Message:	Updated lumi result in acc. w/ CR's calcs
#	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.35	string directoryname="official_2096";
18	buchmann	1.34
19	buchmann	1.35	float luminosity=2096.0;//1936;//751.0;//486.0;//468.0//336.;//pb^{-1}
20	buchmann	1.34	float lumiuncert=0.045;// to be indicated in [0,1] range, e.g. for 4% write 0.04
21	buchmann	1.4	string jzbvariabledata="jzb[1]+0.06*pt";
22			string jzbvariablemc="jzb[1]+0.04*pt";
23	fronga	1.25	float jzbHigh = 350.; // Range for JZB plots
24	buchmann	1.34
25	buchmann	1.1	samplecollection allsamples("completesamplecollection");
26	buchmann	1.9	samplecollection signalsamples("signalsamplecollection");
27	buchmann	1.14	samplecollection scansample("scansamplecollection");
28	buchmann	1.1	int data=1;
29			int mc=0;
30			int mcwithsignal=2;
31			TCut passtrig("(passed_triggers\|\|!is_data)");
32	buchmann	1.20	TCut cutmass("abs(mll-91.2)<20");
33	buchmann	1.1	//TCut cutmass("mll>2");
34			TCut basiccut("mll>2");//basically nothing.
35	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.
36	buchmann	1.1	//TCut jetqualitycut("(pfJetGoodNum>=2&&pfJetGoodID[0])&&(pfJetGoodNum>=2&&pfJetGoodID[1])");//now part of the basiccut
37			//TCut jetqualitycut("mll>0");
38	buchmann	1.7
39	buchmann	1.20	TCut cutnJets("pfJetGoodNum>=3"&&basicqualitycut);
40	buchmann	1.7	TCut cutnJetsJESdown("pfJetGoodNum25>=3"&&basicqualitycut);
41			TCut cutnJetsJESup("pfJetGoodNum35>=3"&&basicqualitycut);
42	buchmann	1.1	TCut cutOSOF("(id1!=id2)&&(ch1*ch2<0)");
43			TCut cutOSSF("(id1==id2)&&(ch1*ch2<0)");
44	buchmann	1.22	TCut sidebandcut("(mll>55&&mll<70)\|\|(mll>112&&mll<160)");
45	fronga	1.25
46	buchmann	1.22	//TCut sidebandcut("(mll>61&&mll<70)\|\|(mll>112&&mll<190)");
47	buchmann	1.1	//TCut basiccut("(passed_triggers\|\|!is_data)");
48
49	buchmann	1.15	// SUSY scan parameters
50			// float mglustart=25;float mgluend=1200;float mglustep=25; //guessed values for official file
51			// float mLSPstart=25;float mLSPend=1200;float mLSPstep=25; //guessed values for official file
52	buchmann	1.16	float mglustart=25;float mgluend=1200;float mglustep=25;
53			float mLSPstart=25;float mLSPend=1200;float mLSPstep=25;
54	buchmann	1.31
55			float m0start=20; float m0end=2000; float m0step=20;
56			float m12start=20; float m12end=760; float m12step=20;
57	buchmann	1.15
58
59	buchmann	1.1	TCut essential(passtrig);//add here any cuts you ALWAYS want
60			int dogaus=0;
61			int doKM=1;
62			int dogaus2sigma=2;
63			int dogaus3sigma=3;
64			int Kostasmethod=-99;
65
66	buchmann	1.13	float fitresultconstdata=0;//this is the result when fitting in the 0-30 GeV range
67			float fitresultconstmc=0;//this is the result when fitting in the 0-30 GeV range
68	buchmann	1.1
69	buchmann	1.33	int method=-1;//Fitting method
70	buchmann	1.8
71			//now some style issues:
72			float DataMarkerSize=1.2;
73	buchmann	1.10
74			//here we save our number of predicted and observed events (with errors)
75			vector<float> Nobs;
76			vector<float> Npred;
77			vector<float> Nprederr;
78	buchmann	1.13
79			int noJES=0;
80			int JESdown=1;
81			int JESup=2;
82	buchmann	1.14
83	buchmann	1.17	//some refinement: nicer color gradient
84			Double_t stops[5] = { 0.00, 0.34, 0.61, 0.84, 1.00 };
85			Double_t red[5] = { 0.00, 0.00, 0.87, 1.00, 0.51 };
86			Double_t green[5] = { 0.00, 0.81, 1.00, 0.20, 0.00 };
87			Double_t blue[5] = { 0.51, 1.00, 0.12, 0.00, 0.00 };
88	buchmann	1.18	int fi=TColor::CreateGradientColorTable(5, stops, red, green,blue, 255);
89	buchmann	1.17
90	buchmann	1.32	// LIMITS
91	buchmann	1.33	int nlimittoys=1000; // how many toys for setting limits
92	buchmann	1.32	string limitmethod="bayesian";//what method to use to set limits
93	buchmann	1.24	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.
94	buchmann	1.30	bool ConsiderSignalContaminationForLimits=false; //whether or not to consider signal contamination when computing limits (standard:true)
95	buchmann	1.23
96	buchmann	1.28	float JZBPeakPositionData=-999;
97			float JZBPeakPositionMC=-999;
98			float JZBPeakWidthData=-999;
99			float JZBPeakWidthMC=-999;
100
101	buchmann	1.1	}