ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

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

Comparing UserCode/cbrown/AnalysisFramework/Plotting/Modules/ActiveSamples.C (file contents):
Revision 1.35 by buchmann, Wed Sep 28 16:35:29 2011 UTC vs.
Revision 1.56 by fronga, Wed Jan 18 14:14:29 2012 UTC

#	Line 1 | Line 1
1	+	/****
2	+
3	+	Off peak status (RestrictToMassPeak) :
4	+
5	+	x  Necessary adaptations identified
6	+	x  Started working on necessary adaptations
7	+	x  Necessary adaptations implemented
8	+	x  Necessary adaptations tested
9	+
10	+	DONE!
11	+
12	+
13	+	****/
14	+
15		#include <iostream>
16		#include <vector>
17		#include <sys/stat.h>
#	Line 20 | Line 34
34		#define HUSH 1
35		#endif
36
37	<	void define_SMS_sample(bool showList, samplecollection &allsamples, samplecollection &signalsamples, samplecollection &scansample)
37	>	void define_SMS_sample(bool showList, samplecollection &allsamples, samplecollection &signalsamples, samplecollection &scansample, bool loadall=false)
38		{
39		char hostname[1023];
40		gethostname(hostname,1023);
41
42		/*
43
44	<	CRUCIAL: The NAME needs to contain either "SMS" or "mSUGRA" for the SUSY scan algorithm (to be able to distinguish the two)
44	>	CRUCIAL: The NAME needs to contain either "SMS" or "mSUGRA" for the SUSY scan algorithm (to be able to distinguish the two) -- this is now done automatically based on PlottingSetup::ScanSampleDirectory
45
46		*/
47	<
48	<	if((int)((string)hostname).find("falda")>-1) {
49	<	scansample.AddSample("/scratch/buchmann/mSUGRA3.root","SMS_T6Z",1,1,false,true,1,kRed); // there is only one sample in the scan; we give each event weight "1"
50	<	} else {
51	<	//      scansample.AddSample("/shome/pablom/SimplifiedModels.root","SMS",1,1,false,true,1,kRed); // there is only one sample in the scan; we give each event weight "1"
52	<	//      scansample.AddSample("/shome/pablom/SimplifiedModels.root","SMS",1,1,false,true,1,kRed); // there is only one sample in the scan; we give each event weight "1"
53	<	scansample.AddSample("/shome/buchmann/Complete_mSUGRA_Scan_cleaned.root","mSUGRA",1,1,false,true,1,kRed); // there is only one sample in the scan; we give each event weight "1"
54	<	for(int ix=0;ix<PlottingSetup::ScanXzones;ix++) {
55	<	for(int iy=0;iy<PlottingSetup::ScanYzones;iy++) {
56	<	scansample.AddSample("/shome/buchmann/ntuples/mSUGRA/mSUGRA_clean_splitup_"+any2string(ix)+"_"+any2string(iy)+".root","mSUGRA",1,1,false,true,PlottingSetup::ScanYzones*ix+iy,kRed); // there is only one sample in the scan; we give each event weight "1"
57	<	}
47	>
48	>	string prefix="/shome/";
49	>	if((int)((string)hostname).find("falda")>-1) prefix="/scratch/";
50	>	string mSUGRAorSMS="mSUGRA";
51	>	if((int)((string)PlottingSetup::ScanSampleDirectory).find("SMS")>-1) mSUGRAorSMS="SMS";
52	>	if((int)((string)PlottingSetup::ScanSampleDirectory).find("GMSB")>-1) mSUGRAorSMS="GMSB";
53	>
54	>	/// NEW: Only indicate the very first file - please place scans in ntuples/SMS/ and adapt the name here and in SUSYScan.C (the files are loaded when needed avoiding high memory usage and startup times)
55	>	/// If you have different scans (e.g. you have a T5zzl one as well) you can create an additional directory such as SMS_T5zzl, but you need to change that below and in SUSYScan.C as well (where the files are loaded)
56	>	scansample.AddSample(prefix+"buchmann/ntuples/"+PlottingSetup::ScanSampleDirectory+"/"+mSUGRAorSMS+"_clean_splitup_0_0.root",mSUGRAorSMS,1,1,false,true,0,kRed); // there is only one sample in the scan; we give each event weight "1"
57	>
58	>	if(loadall) {
59	>	for(int ix=0;ix<PlottingSetup::ScanXzones;ix++) {
60	>	for(int iy=1;iy<PlottingSetup::ScanYzones;iy++) {//note: the "zero-zero" sample is already loaded above!
61	>	scansample.AddSample((prefix+"buchmann/ntuples/"+PlottingSetup::ScanSampleDirectory+"/"+mSUGRAorSMS+"_clean_splitup_"+any2string(ix)+"_"+any2string(iy)+".root"),mSUGRAorSMS,1,1,false,true,ix*PlottingSetup::ScanXzones+iy,kRed); // there is only one sample in the scan; we give each event weight "1"
62		}
63	+	}
64	+	}
65	+
66	+
67	+	flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak
68	+	if(!PlottingSetup::RestrictToMassPeak) {
69	+	PlottingSetup::cutmass=PlottingSetup::openmasscut;
70	+	//      PlottingSetup::sidebandcut=TCut("mll<2&&mll>3");//impossible cut - because we don't want any sidebands!
71	+	PlottingSetup::sidebandcut=TCut("SideBandsHaveBeenCalledForOffPeakAnalysis__ERROR==1");//this will trigger errors whenever sidebands are called - perfect for checking any remaining uses (checked, and no uses were found).
72	+	PlottingSetup::genMassCut=PlottingSetup::openGenmasscut;//impossible cut - because we don't want any sidebands!
73		}
74		}
75
76	<	void define_samples(bool showList, samplecollection &allsamples, samplecollection &signalsamples, samplecollection &scansample, samplecollection &raresample)
76	>
77	>	void define_samples(bool showList, samplecollection &allsamples, samplecollection &signalsamples, samplecollection &scansample, samplecollection &raresample, samplecollection &systsamples)
78		{
79	+	flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak
80	+	if(!PlottingSetup::RestrictToMassPeak) {
81	+	PlottingSetup::cutmass=PlottingSetup::openmasscut;
82	+	PlottingSetup::sidebandcut=TCut("mll<2&&mll>3");//impossible cut - because we don't want any sidebands!
83	+	PlottingSetup::genMassCut=PlottingSetup::openGenmasscut;//impossible cut - because we don't want any sidebands!
84	+	}
85
86		float ZJetsCrossSection         = 3048.0; //NNLO----------------------------------
87	+	//    float ZJetsCrossSectionM20      = 4998.0;//summer11
88	+	//    float ZJetsCrossSectionD6Tlow   = 400.0;//summer11
89		float TTbarCrossSection         = 165.0;//(NLO) ---- 165.0; // approx. NNLO-----
90		float WJetsCrossSection         = 31314.0;//NNLO-------3.131e4; //NNLO------------
91		float ZnunuCrossSection         = 5760.0;//NNLO -------4.5e+3; //(LO);------------
92	<	float SingleTopSCrossSection     = 1.49; // NLO;----------------------------------
93	<	float SingleTopTCrossSection     = 20.540; // NLO;--------------------------------
94	<	float SingleTopUCrossSection     = 10.6; // NLO;----------------------------------
92	>	//Single top x-sections from here:
93	>	//https://twiki.cern.ch/twiki/bin/view/CMS/SingleTopSigma
94	>	float SingleTopSCrossSection     = 4.63; // NLO;----------------------------------
95	>	float SingleTopTCrossSection     = 64.57; // NLO;--------------------------------
96	>	float SingleTopUCrossSection     = 15.74; // NLO;----------------------------------
97		float VVJetsCrossSection        = 4.8; // LO;-------------------------------------
98	<	float LM4CrossSection           = 2.537; // k*LO
98	>	float LM4CrossSection           = 1.879*1.35; // k*LO
99	>	float LM8CrossSection           = 0.7300*1.41 ; // k*LO
100		//the following numbers are from the MadGraphStandardModel210Summary as linked on the GeneratorProduction2011 page
101		//    float QCD50to100CrossSection=30000000; // not used
102		float QCD100to250CrossSection=7000000;
#	Line 67 | Line 107 | void define_samples(bool showList, sampl
107		//    float totEventsQCD50to100=207418.0; // not used
108		long totEventsQCD100to250   = 638792;
109		long totEventsQCD250to500   = 344454;
110	<	long totEventsQCD500to1000  = 10e10;  //unknown but now suppressed. (no event passes anyway)
110	>	long totEventsQCD500to1000  = (long)10e10;  //unknown but now suppressed. (no event passes anyway)
111		long totEventsQCD1000toInf  = 163185;
112
113		long totEventsZjets         = 32005720;//summer11
114	+	long totEventsZjetsD6T      = 2543706;//summer11
115	+	long totEventsZjetsD6Tlow   = 181799;//summer11
116	+	//    long totEventsZjetsEE       = 2738801;
117	+	//    long totEventsZjetsMM       = 2148325;
118	+	//    long totEventsZjetsTT       = 2032536;
119	+
120		//spring:2313911.0;//2743142.0; // DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola
121	<	long totEventsTTbar         =  3701872;//summer11
121	>	long totEventsTTbar         =  3701947;//summer11
122		//1161621.0;//Spring//1144028.0; // PabloV13/TTJets_TuneZ2_7TeV-madgraph-tauola
123	<	long totEventsWJets         = 53739703;//summer11
123	>	long totEventsWJets         = 56674902;//summer11
124		//15010237.0;//spring//14350756.0; /// WJetsToLNu_TuneZ2_7TeV-madgraph-tauola.root
125		long totEventsZnunu         =  2106977;//2167964.0; //ZinvisibleJets_7TeV-madgraph.root
126		long totEventsVVJets        =   959076;//509072.0;
127		long totEventsSingleTopS    =   493868;//489472.0;
128		long totEventsSingleTopT    =   475460;//477610.0;
129		long totEventsSingleTopU    =   489417;//477599.0;
130	+	long totEventsLM4           =   218380;
131	+	long totEventsLM8           =   180960;
132
133
134		Int_t nice_blue  = TColor::GetColor("#2E9AFE");
#	Line 104 | Line 152 | void define_samples(bool showList, sampl
152		ttbar_color         = kMagenta+2;
153		singletop_color     = kBlue;
154		qcd_color           = kPink;
155	<	diboson_color       = kGreen;
155	>	diboson_color       = kGreen+2;
156		lm_color            = kViolet+7;
157
158		Int_t rare_color            = nice_blue;
#	Line 112 | Line 160 | void define_samples(bool showList, sampl
160		char hostname[1023];
161		gethostname(hostname,1023);
162
163	<	if((int)((string)hostname).find("falda")>-1) {
164	<	//not running locally on my computer
165	<	//      allsamples.AddSample("/scratch/buchmann/ntuples2011/data/AllData_Jul6_1078pb_boolean.root","Data",1,1,true,false,0,kBlack); // this is the 1078 /pb version
166	<	//      allsamples.AddSample("/scratch/fronga/v1.47/AllData_1.4fb_stripped.root","Data",1,1,true,false,0,kBlack);
119	<	//      allsamples.AddSample("/scratch/buchmann/Data_Aug26__1936.root","Data",1,1,true,false,0,kBlack);
120	<	allsamples.AddSample("/scratch/buchmann/ntuples2011/data/AllData_Aug26__2048_ipb__NewFormat_Complete_official.root","Data",1,1,true,false,0,kBlack);
121	<	//      allsamples.AddSample("/scratch/buchmann/ntuples2011/data/Data_v6_test.root","Data",1,1,true,false,0,kBlack);
122	<	//      allsamples.AddSample("/scratch/fronga/JSON_160404-163869_v1.29/AllCertified191_v1.29.root","Data",1,1,true,false,0,kBlack); // this is the version used for the 2011 PAS (update), i.e. SUS-11-012
123	<
124	<	allsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/TTJets_TuneZ2_7TeV-madgraph-tauola_Summer11_2.root","TTbar",totEventsTTbar,TTbarCrossSection,false,false,1,ttbar_color);//summer11
125	<	allsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/WJetsToLNu_TuneZ2_7TeV-madgraph-tauola.root","WJets",totEventsWJets,WJetsCrossSection,false,false,3,wjets_color);//summer
126	<	allsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/VVJetsTo4L_TuneD6T_7TeV-madgraph-tauola.root","DiBosons",totEventsVVJets,VVJetsCrossSection,false,false,4,diboson_color);
127	<	allsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/TToBLNu_TuneZ2_s-channel_7TeV-madgraph_2.root","SingleTop",totEventsSingleTopS,SingleTopSCrossSection,false,false,2,singletop_color);
128	<	allsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/TToBLNu_TuneZ2_t-channel_7TeV-madgraph.root","SingleTop",totEventsSingleTopT,SingleTopTCrossSection,false,false,2,singletop_color);
129	<	allsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/TToBLNu_TuneZ2_tW-channel_7TeV-madgraph.root","SingleTop",totEventsSingleTopU,SingleTopUCrossSection,false,false,2,singletop_color);
130	<	allsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/ZinvisibleJets_7TeV-madgraph.root","Z nunu",totEventsZnunu,ZnunuCrossSection,false,false,4,diboson_color);
131	<	allsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola_Summer11.root","ZJets",totEventsZjets,ZJetsCrossSection,false,false,6,dy_color);//summer11
132	<	allsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM4_SUSY_sftsht_7TeV-pythia6.root","LM4",218536,1.879,false,true,7,lm_color);
133	<
134	<	signalsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM0_SUSY_sftsht_7TeV-pythia6.root","LM0",219796,38.93 ,false,true,3,lm_color);
135	<	signalsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM1_SUSY_sftsht_7TeV-pythia6.root","LM1",218176,4.888,false,true,4,lm_color);
136	<	signalsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM2_SUSY_sftsht_7TeV-pythia6.root","LM2",206336,0.6027,false,true,5,lm_color);
137	<	signalsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM3_SUSY_sftsht_7TeV-pythia6.root","LM3",215120,3.438,false,true,6,lm_color);
138	<	signalsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM4_SUSY_sftsht_7TeV-pythia6.root","LM4",218536,1.879,false,true,1,lm_color);
139	<	signalsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM5_SUSY_sftsht_7TeV-pythia6.root","LM5",223992,0.473,false,true,7,lm_color);
140	<	signalsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM8_SUSY_sftsht_7TeV-pythia6.root","LM8",180960,0.73,false,true,2,lm_color);
141	<	signalsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM9_SUSY_sftsht_7TeV-pythia6.root","LM9",227808,7.134,false,true,8,lm_color);
142	<	signalsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM11_SUSY_sftsht_7TeV-pythia6.root","LM11",218703,0.8236,false,true,9,lm_color);
143	<	signalsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM12_SUSY_sftsht_7TeV-pythia6.root","LM12",223491,4.414,false,true,10,lm_color);
144	<	signalsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM13_SUSY_sftsht_7TeV-pythia6.root","LM13",173240,6.899,false,true,11,lm_color);
145	<
146	<	} else {
147	<
148	<	allsamples.AddSample("/shome/buchmann/ntuples/Data/AllData_2096pb_v158.root","Data",1,1,true,false,0,kBlack);
149	<	//      allsamples.AddSample("/shome/buchmann/ntuples/Data/Latest_assemblage.root","Data",1,1,true,false,0,kBlack);
150	<
151	<	allsamples.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/TTJets_TuneZ2_7TeV-madgraph-tauola_Summer11_2.root","TTbar",totEventsTTbar,TTbarCrossSection,false,false,1,ttbar_color);//summer11
152	<
153	<	//write_warning(__FUNCTION__,"TEMPORARILY USING POWHEG TT to 2l 2nu 2b sample!");
154	<	//allsamples.AddSample("/scratch/buchmann/ntuples/NewPOWHEG_sample/TTTo2L2Nu2B_7TeV-powheg-pythia6__Summer11-PU_S4_START42_V11-v1___previously_TTJets.root","TTbar",10339374,16.53,false,false,1,ttbar_color);//summer11
155	<	//xs from https://twiki.cern.ch/twiki/bin/view/CMS/CrossSectionDetails
156	<
157	<	allsamples.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/WJetsToLNu_TuneZ2_7TeV-madgraph-tauola_Summer11.root","WJets",totEventsWJets,WJetsCrossSection,false,false,3,wjets_color);//summer
158	<	allsamples.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/VVJetsTo4L_TuneD6T_7TeV-madgraph-tauola.root","DiBosons",totEventsVVJets,VVJetsCrossSection,false,false,4,diboson_color);
159	<	allsamples.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/TToBLNu_TuneZ2_s-channel_7TeV-madgraph_2.root","SingleTop",totEventsSingleTopS,SingleTopSCrossSection,false,false,2,singletop_color);
160	<	allsamples.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/TToBLNu_TuneZ2_t-channel_7TeV-madgraph.root","SingleTop",totEventsSingleTopT,SingleTopTCrossSection,false,false,2,singletop_color);
161	<	allsamples.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/TToBLNu_TuneZ2_tW-channel_7TeV-madgraph.root","SingleTop",totEventsSingleTopU,SingleTopUCrossSection,false,false,2,singletop_color);
162	<	allsamples.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/ZinvisibleJets_7TeV-madgraph.root","Z nunu",totEventsZnunu,ZnunuCrossSection,false,false,4,diboson_color);
163	<
164	<
165	<
166	<
167	<	allsamples.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola_Summer11.root","ZJets",totEventsZjets,ZJetsCrossSection,false,false,7,dy_color);//summer11
168	<	//      allsamples.AddSample("/scratch/buchmann/Skimmed_DY.root","ZJets",totEventsZjets,ZJetsCrossSection,false,false,6,dy_color);//summer11
169	<
170	<
171	<	allsamples.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/LM4_SUSY_sftsht_7TeV-pythia6.root","LM4",218536,1.879,false,true,8,lm_color);
172	<
173	<
174	<	/*      signalsamples.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/LM0_SUSY_sftsht_7TeV-pythia6.root","LM0",219796,38.93 ,false,true,3,lm_color);
175	<	signalsamples.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/LM1_SUSY_sftsht_7TeV-pythia6.root","LM1",218176,4.888,false,true,4,lm_color);
176	<	signalsamples.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/LM2_SUSY_sftsht_7TeV-pythia6.root","LM2",206336,0.6027,false,true,5,lm_color);
177	<	signalsamples.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/LM3_SUSY_sftsht_7TeV-pythia6.root","LM3",215120,3.438,false,true,6,lm_color);*/
178	<	signalsamples.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/LM4_SUSY_sftsht_7TeV-pythia6.root","LM4",218536,1.879,false,true,1,lm_color);/*
179	<	signalsamples.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/LM5_SUSY_sftsht_7TeV-pythia6.root","LM5",223992,0.473,false,true,7,lm_color);*/
180	<	signalsamples.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/LM8_SUSY_sftsht_7TeV-pythia6.root","LM8",180960,0.73,false,true,2,lm_color);/*
181	<	signalsamples.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/LM9_SUSY_sftsht_7TeV-pythia6.root","LM9",227808,7.134,false,true,8,lm_color);
182	<	signalsamples.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/LM11_SUSY_sftsht_7TeV-pythia6.root","LM11",218703,0.8236,false,true,9,lm_color);
183	<	signalsamples.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/LM12_SUSY_sftsht_7TeV-pythia6.root","LM12",223491,4.414,false,true,10,lm_color);
184	<	signalsamples.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/LM13_SUSY_sftsht_7TeV-pythia6.root","LM13",173240,6.899,false,true,11,lm_color);*/
185	<
186	<
187	<	// We do not use the following (rare) samples; nonetheless we include them in a separate category to prove that they are negligible
188	<	raresample.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/TTWWTo2Lplus2Nu_NTP.root","Rare SS",39791,4.76E-5,false,false,6,rare_color);//summer11
189	<	raresample.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/WWminusTo2L2Nu_NTP.root","Rare SS",51063,0.001559,false,false,6,rare_color);//summer11
190	<	raresample.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/TTWWTo2Lminus2Nu_NTP.root","Rare SS",39496 ,4.719E-5,false,false,6,rare_color);//summer11
191	<	raresample.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/TTWTo2Lplus2Nu_NTP.root","Rare SS",48764,0.006841,false,false,6,rare_color);//summer11
192	<	raresample.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/WWplusTo2L2Nu_NTP.root","Rare SS",48244 ,0.004216,false,false,6,rare_color);//summer11
193	<	raresample.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/TTZTo2Lplus2Nu_NTP.root","Rare SS",36740,0.002024,false,false,6,rare_color);//summer11
194	<	raresample.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/WWWTo2Lplus2Nu_NTP.root","Rare SS",43322,0.008957,false,false,6,rare_color);//summer11
195	<	raresample.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/TTWTo2Lminus2Nu_NTP.root","Rare SS",50435,0.002705,false,false,6,rare_color);//summer11
196	<	raresample.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/WWWTo2Lminus2Nu_NTP.root","Rare SS",17147,0.004109,false,false,6,rare_color);//summer11
197	<	raresample.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/TTZTo2Lminus2Nu_NTP.root","Rare SS",36781,0.001946,false,false,6,rare_color);//summer11
198	<
163	>	string prefix="/shome/";
164	>	if((int)((string)hostname).find("falda")>-1) prefix="/scratch/";
165	>
166	>	allsamples.AddSample(prefix+"buchmann/ntuples/Data/AllData_Complete_4653_ipb_PAPER.root","Data",1,1,true,false,0,kBlack);
167
168
169	+	allsamples.AddSample(prefix+"buchmann//ntuples/MC/NewSelectionForPaper/NewlyReweighted/DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola_Summer11___NTupleV020300.root","Z+Jets",totEventsZjets,ZJetsCrossSection,false,false,7,dy_color);
170	+	allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/TTJets_TuneZ2_7TeV-madgraph-tauola_Summer11_3.root","t#bar{t}",totEventsTTbar,TTbarCrossSection,false,false,1,ttbar_color);//summer11
171	+	allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/WJetsToLNu_TuneZ2_7TeV-madgraph-tauola_Summer11_2.root","W+Jets",totEventsWJets,WJetsCrossSection,false,false,3,wjets_color);//summer
172	+	allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/VVJetsTo4L_TuneD6T_7TeV-madgraph-tauola.root","WW/WZ/ZZ",totEventsVVJets,VVJetsCrossSection,false,false,4,diboson_color);
173	+	allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/TToBLNu_TuneZ2_s-channel_7TeV-madgraph_2.root","Single-top",totEventsSingleTopS,SingleTopSCrossSection,false,false,2,singletop_color);
174	+	allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/TToBLNu_TuneZ2_t-channel_7TeV-madgraph.root","Single-top",totEventsSingleTopT,SingleTopTCrossSection,false,false,2,singletop_color);
175	+	allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/TToBLNu_TuneZ2_tW-channel_7TeV-madgraph.root","Single-top",totEventsSingleTopU,SingleTopUCrossSection,false,false,2,singletop_color);
176	+	allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/ZinvisibleJets_7TeV-madgraph.root","Z nunu",totEventsZnunu,ZnunuCrossSection,false,false,4,diboson_color);
177	+	allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/LM4_SUSY_sftsht_7TeV-pythia6.root","LM4",totEventsLM4,LM4CrossSection,false,true,8,lm_color);
178
202	–	}
203	–	/**/
204	–	//    signalsamples.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/WithTwoSix/LM4_SUSY_sftsht_7TeV-pythia6.root","REALLY OLD LM4!",218536,1.879,false,true,1,lm_color);
205	–	//    signalsamples.AddSample("/shome/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/MCSummer_PURW_172802_LP__JetID3p0/WithTwoSix/LM8_SUSY_sftsht_7TeV-pythia6.root","REALLY OLD LM8!",218536,0.73,false,true,1,lm_color);
206	–
207	–	//    scansample.AddSample("/scratch/fronga/SMS/T6z_GlGlToQQZQQZ_38xFall10.root","T6Z",1,1,false,true,1,lm_color); // there is only one sample in the scan; we give each event weight "1"
208	–	//scansample.AddSample("/scratch/buchmann/allfinalevent_NTupleProducer_42X_mc.root","T6Z",1,1,false,true,1,lm_color); // there is only one sample in the scan; we give each event weight "1"
209	–	//    scansample.AddSample("/scratch/buchmann/SUSYScanTestJZB.root","T6Z",1,1,false,true,1,lm_color); // there is only one sample in the scan; we give each event weight "1"
179
180
181	<	define_SMS_sample(showList, allsamples, signalsamples, scansample);
181	>	if(!PlottingSetup::RestrictToMassPeak) {
182	>	//      signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM0_SUSY_sftsht_7TeV-pythia6.root","LM0",219796,38.93*1.41 ,false,true,3,lm_color);
183	>	//      signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM1_SUSY_sftsht_7TeV-pythia6.root","LM1",218176,4.888*1.34,false,true,4,lm_color);
184	>	//      signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM2_SUSY_sftsht_7TeV-pythia6.root","LM2",206336,0.6027*1.33,false,true,5,lm_color);
185	>	signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM3_SUSY_sftsht_7TeV-pythia6.root","LM3",215120,3.438*1.40,false,true,6,lm_color);
186	>	signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM5_SUSY_sftsht_7TeV-pythia6.root","LM5",223992,0.473*1.34,false,true,7,lm_color);
187	>	signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM9_SUSY_sftsht_7TeV-pythia6.root","LM9",227808,7.134*1.48,false,true,8,lm_color);
188	>	signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM11_SUSY_sftsht_7TeV-pythia6.root","LM11",218703,0.8236*1.35,false,true,9,lm_color);
189	>	signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM12_SUSY_sftsht_7TeV-pythia6.root","LM12",223491,4.414*1.34,false,true,10,lm_color);
190	>	signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM13_SUSY_sftsht_7TeV-pythia6.root","LM13",173240,6.899*1.42,false,true,11,lm_color);
191	>	}
192	>	signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM4_SUSY_sftsht_7TeV-pythia6.root","LM4",totEventsLM4,LM4CrossSection,false,true,1,lm_color);
193	>	signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM8_SUSY_sftsht_7TeV-pythia6.root","LM8",totEventsLM8,LM8CrossSection,false,true,2,lm_color);
194
195	<	if(showList) allsamples.ListSamples();
195	>	//    systsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/HighStatSamples/DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola_Fall11-PU_S6_START42_V14B-v1.root","Z+Jets",18368890,ZJetsCrossSection,false,false,7,dy_color);
196
216	–	}
197
198	+	systsamples.AddSample(prefix+"buchmann//ntuples/MC/NewSelectionForPaper/NewlyReweighted/DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola_Summer11___NTupleV020300.root","Z+Jets",32005720,ZJetsCrossSection,false,false,7,dy_color);
199	+	systsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/HighStatSamples/TTJets_Fall11-PU_S6_START42_V14B.root","t#bar{t}",58370118,TTbarCrossSection,false,false,1,ttbar_color);//summer11
200
201	<	void define_spring_samples(bool showList, samplecollection &allsamples)
202	<	{
203	<
204	<	float ZJetsCrossSection         = 3048.0; //NNLO----------------------------------
205	<	float TTbarCrossSection         = 165.0;//(NLO) ---- 165.0; // approx. NNLO-----
206	<	float WJetsCrossSection         = 31314.0;//NNLO-------3.131e4; //NNLO------------
207	<	float ZnunuCrossSection         = 5760.0;//NNLO -------4.5e+3; //(LO);------------
208	<	float SingleTopSCrossSection     = 1.49; // NLO;----------------------------------
209	<	float SingleTopTCrossSection     = 20.540; // NLO;--------------------------------
210	<	float SingleTopUCrossSection     = 10.6; // NLO;----------------------------------
211	<	float VVJetsCrossSection        = 4.8; // LO;-------------------------------------
212	<	float LM4CrossSection           = 2.537; // k*LO
213	<	//the following numbers are from the MadGraphStandardModel210Summary as linked on the GeneratorProduction2011 page
214	<	//    float QCD50to100CrossSection=30000000; // not used
215	<	float QCD100to250CrossSection=7000000;
216	<	float QCD250to500CrossSection=171000;
217	<	float QCD500to1000CrossSection=5200;
218	<	float QCD1000toInfCrossSection=83.0;
237	<
238	<	//    float totEventsQCD50to100=207418.0; // not used
239	<	float totEventsQCD100to250=638792.0;
240	<	float totEventsQCD250to500=344454.0;
241	<	float totEventsQCD500to1000=10e10;//unknown but now suppressed. (no event passes anyway)
242	<	float totEventsQCD1000toInf=163185.0;
243	<
244	<	float totEventsZjets         = 32005720.0;//summer11
245	<	//spring:2313911.0;//2743142.0; // DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola
246	<	float totEventsTTbar         = 3701947.0;//summer11
247	<	//1161621.0;//Spring//1144028.0; // PabloV13/TTJets_TuneZ2_7TeV-madgraph-tauola
248	<	float totEventsWJets         = 53739703.0;//summer11
249	<	//15010237.0;//spring//14350756.0; /// WJetsToLNu_TuneZ2_7TeV-madgraph-tauola.root
250	<	float totEventsZnunu         = 2106977.0;//2167964.0; //ZinvisibleJets_7TeV-madgraph.root
251	<	float totEventsVVJets        = 959076.0;//509072.0;
252	<	float totEventsSingleTopS    = 493868.0;//489472.0;
253	<	float totEventsSingleTopT    = 475460.0;//477610.0;
254	<	float totEventsSingleTopU    = 489417.0;//477599.0;
201	>	//
202	>	/*
203	>
204	>
205	>	// We do not use the following (rare) samples; nonetheless we include them in a separate category to prove that they are negligible
206	>	write_warning(__FUNCTION__,"Rare samples not yet reprocessed with new selectioN");/*
207	>	raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTWWTo2Lplus2Nu_NTP.root","Rare SS",39791,4.76E-5,false,false,6,rare_color);//summer11
208	>	raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/WWminusTo2L2Nu_NTP.root","Rare SS",51063,0.001559,false,false,6,rare_color);//summer11
209	>	raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTWWTo2Lminus2Nu_NTP.root","Rare SS",39496 ,4.719E-5,false,false,6,rare_color);//summer11
210	>	raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTWTo2Lplus2Nu_NTP.root","Rare SS",48764,0.006841,false,false,6,rare_color);//summer11
211	>	raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/WWplusTo2L2Nu_NTP.root","Rare SS",48244 ,0.004216,false,false,6,rare_color);//summer11
212	>	raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTZTo2Lplus2Nu_NTP.root","Rare SS",36740,0.002024,false,false,6,rare_color);//summer11
213	>	raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/WWWTo2Lplus2Nu_NTP.root","Rare SS",43322,0.008957,false,false,6,rare_color);//summer11
214	>	raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTWTo2Lminus2Nu_NTP.root","Rare SS",50435,0.002705,false,false,6,rare_color);//summer11
215	>	raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/WWWTo2Lminus2Nu_NTP.root","Rare SS",17147,0.004109,false,false,6,rare_color);//summer11
216	>	raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTZTo2Lminus2Nu_NTP.root","Rare SS",36781,0.001946,false,false,6,rare_color);//summer11
217	>	*/
218	>	//scan files are only defined when necessary!
219
256	–
257	–	Int_t nice_blue  = TColor::GetColor("#2E9AFE");
258	–	Int_t nice_green = TColor::GetColor("#81f781");
259	–	Int_t nice_red   = TColor::GetColor("#F78181");
260	–	Int_t nice_pink  = TColor::GetColor("#F781BE");
261	–	Int_t nice_orange= TColor::GetColor("#F7BE81");//#FA9624
262	–	Int_t nice_black = TColor::GetColor("#2A0A1B");//black with a bit of purple :-)
263	–
264	–
265	–	Int_t dy_color        = nice_blue;
266	–	Int_t wjets_color     = kGray;
267	–	Int_t ttbar_color     = nice_green;
268	–	Int_t singletop_color = nice_orange;
269	–	Int_t qcd_color       = nice_black;
270	–	Int_t diboson_color   = nice_pink;
271	–	Int_t lm_color        = kRed;
272	–
273	–	dy_color            = kYellow;
274	–	wjets_color         = kGray;
275	–	ttbar_color         = kMagenta+2;
276	–	singletop_color     = kBlue;
277	–	qcd_color           = kPink;
278	–	diboson_color       = kGreen;
279	–	lm_color            = kViolet+7;
280	–
281	–	//    allsamples.AddSample("/scratch/buchmann/AllData_Jun10___DCS_648pb_MoreTriggers4_DCS.root","Data",1,1,true,false,0,kBlack);
282	–	//    allsamples.AddSample("/scratch/buchmann/AllData_Jun10___486pb_MoreTriggers3.root","Data",1,1,true,false,0,kBlack);
283	–	//allsamples.AddSample("/scratch/buchmann/AllData_Jun17_160404-166861_real_complete_2.root","Data",1,1,true,false,0,kBlack);
284	–	allsamples.AddSample("/scratch/buchmann/ntuples2011/data/AllData_Jul6_1078pb.root","Data",1,1,true,false,0,kBlack);
285	–
286	–	//    allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/TTJets_TuneZ2_7TeV-madgraph-tauola.root","t#bar{t}",totEventsTTbar,TTbarCrossSection,false,false,1,ttbar_color);//spring11
287	–	allsamples.AddSample("/scratch/fronga/MC_v1.29/TTJets_TuneZ2_7TeV-madgraph-tauola_Summer11_3.root","t#bar{t}",totEventsTTbar,TTbarCrossSection,false,false,1,ttbar_color);//summer11
288	–	//    allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/WJetsToLNu_TuneZ2_7TeV-madgraph-tauola.root","WJets",totEventsWJets,WJetsCrossSection,false,false,3,wjets_color); //spring11
289	–	allsamples.AddSample("/scratch/fronga/MC_v1.29/WJetsToLNu_TuneZ2_7TeV-madgraph-tauola.root","WJets",totEventsWJets,WJetsCrossSection,false,false,3,wjets_color);//summer
290	–	allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/VVJetsTo4L_TuneD6T_7TeV-madgraph-tauola.root","Dibosons",totEventsVVJets,VVJetsCrossSection,false,false,4,diboson_color);
291	–	allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/TToBLNu_TuneZ2_s-channel_7TeV-madgraph_2.root","SingleTop",totEventsSingleTopS,SingleTopSCrossSection,false,false,2,singletop_color);
292	–	allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/TToBLNu_TuneZ2_t-channel_7TeV-madgraph.root","SingleTop",totEventsSingleTopT,SingleTopTCrossSection,false,false,2,singletop_color);
293	–	allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/TToBLNu_TuneZ2_tW-channel_7TeV-madgraph.root","SingleTop",totEventsSingleTopU,SingleTopUCrossSection,false,false,2,singletop_color);
294	–	allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/ZinvisibleJets_7TeV-madgraph.root","Z nunu",totEventsZnunu,ZnunuCrossSection,false,false,4,diboson_color);
295	–	//    allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/QCD_TuneD6T_HT-100To250_7TeV-madgraph.root","QCD",totEventsQCD100to250,QCD100to250CrossSection,false,false,5,qcd_color);
296	–	//    allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/QCD_TuneD6T_HT-250To500_7TeV-madgraph.root","QCD",totEventsQCD250to500,QCD250to500CrossSection,false,false,5,qcd_color);
297	–	//    allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/QCD_TuneD6T_HT-500To1000_7TeV-madgraph.root","QCD",totEventsQCD500to1000,QCD500to1000CrossSection,false,false,5,qcd_color);
298	–	//    allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/QCD_TuneD6T_HT-1000_7TeV-madgraph.root","QCD",totEventsQCD1000toInf,QCD1000toInfCrossSection,false,false,5,qcd_color);
299	–	//    allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola.root","ZJets",totEventsZjets,ZJetsCrossSection,false,false,6,dy_color);//spring11
300	–	allsamples.AddSample("/scratch/fronga/MC_v1.29/DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola_Summer11.root","ZJets",totEventsZjets,ZJetsCrossSection,false,false,6,dy_color);//summer11
301	–	allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/LM4_SUSY_sftsht_7TeV-pythia6.root","LM4",218536,1.879,false,true,7,lm_color);
302	–
220		if(showList) allsamples.ListSamples();
221
222		}
306	–

Diff Legend

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