ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/cbrown/AnalysisFramework/Plotting/Modules/ActiveSamples.C
(Generate patch)

Comparing UserCode/cbrown/AnalysisFramework/Plotting/Modules/ActiveSamples.C (file contents):
Revision 1.38 by buchmann, Thu Sep 29 13:59:44 2011 UTC vs.
Revision 1.46 by buchmann, Tue Nov 8 09:20:00 2011 UTC

# Line 9 | Line 9 | x  Necessary adaptations tested
9  
10   DONE!
11  
12 +
13   ****/
14  
15   #include <iostream>
# Line 33 | Line 34 | DONE!
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);
# Line 43 | Line 44 | void define_SMS_sample(bool showList, sa
44      CRUCIAL: The NAME needs to contain either "SMS" or "mSUGRA" for the SUSY scan algorithm (to be able to distinguish the two)
45      
46      */
47 +
48 +    string prefix="/shome/";
49 +    if((int)((string)hostname).find("falda")>-1) prefix="/scratch/";
50      
47    if((int)((string)hostname).find("falda")>-1) {
48        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"
49    } else {
51        /// 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)
52 < //      scansample.AddSample("/shome/buchmann/ntuples/mSUGRA/mSUGRA_clean_splitup_0_0.root","mSUGRA",1,1,false,true,0,kRed); // there is only one sample in the scan; we give each event weight "1"
53 <        scansample.AddSample("/shome/buchmann/ntuples/SMS/SMS_clean_splitup_0_0.root","SMS",1,1,false,true,0,kRed); // there is only one sample in the scan; we give each event weight "1"
52 >      /// 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)
53 > //      scansample.AddSample(prefix+"buchmann/ntuples/mSUGRA/mSUGRA_clean_splitup_0_0.root","mSUGRA",1,1,false,true,0,kRed); // there is only one sample in the scan; we give each event weight "1"
54 >        scansample.AddSample(prefix+"buchmann/ntuples/SMS_T5zzl/SMS_clean_splitup_0_0.root","SMS",1,1,false,true,0,kRed); // there is only one sample in the scan; we give each event weight "1"
55 >
56 >    if(loadall) {
57 >      for(int ix=0;ix<PlottingSetup::ScanXzones;ix++) {
58 >        for(int iy=1;iy<PlottingSetup::ScanYzones;iy++) {//note: the "zero-zero" sample is already loaded above!
59 >          scansample.AddSample((prefix+"buchmann/ntuples/SMS/SMS_clean_splitup_"+any2string(ix)+"_"+any2string(iy)+".root"),"SMS",1,1,false,true,ix*PlottingSetup::ScanXzones+iy,kRed); // there is only one sample in the scan; we give each event weight "1"
60 >        }
61 >      }
62      }
63 +        
64  
65 <
65 >    flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak
66      if(!PlottingSetup::RestrictToMassPeak) {
67 <      PlottingSetup::cutmass=TCut("mll>55");
68 <      PlottingSetup::sidebandcut=TCut("mll<2&&mll>3");//impossible cut - because we don't want any sidebands!
67 >      PlottingSetup::cutmass=PlottingSetup::openmasscut;
68 > //      PlottingSetup::sidebandcut=TCut("mll<2&&mll>3");//impossible cut - because we don't want any sidebands!
69 >      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).
70 >      PlottingSetup::genMassCut=PlottingSetup::openGenmasscut;//impossible cut - because we don't want any sidebands!
71      }
72   }
73  
74 +
75   void define_samples(bool showList, samplecollection &allsamples, samplecollection &signalsamples, samplecollection &scansample, samplecollection &raresample)
76   {
77 <  write_warning(__FUNCTION__,"Gotten here");
77 >    flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak
78      if(!PlottingSetup::RestrictToMassPeak) {
79        PlottingSetup::cutmass=PlottingSetup::openmasscut;
80        PlottingSetup::sidebandcut=TCut("mll<2&&mll>3");//impossible cut - because we don't want any sidebands!
81 +      PlottingSetup::genMassCut=PlottingSetup::openGenmasscut;//impossible cut - because we don't want any sidebands!
82      }
83    
84      float ZJetsCrossSection         = 3048.0; //NNLO----------------------------------
85 + //    float ZJetsCrossSectionM20      = 4998.0;//summer11
86 + //    float ZJetsCrossSectionD6Tlow   = 400.0;//summer11
87      float TTbarCrossSection         = 165.0;//(NLO) ---- 165.0; // approx. NNLO-----
88      float WJetsCrossSection         = 31314.0;//NNLO-------3.131e4; //NNLO------------
89      float ZnunuCrossSection         = 5760.0;//NNLO -------4.5e+3; //(LO);------------
90 <    float SingleTopSCrossSection     = 1.49; // NLO;----------------------------------
91 <    float SingleTopTCrossSection     = 20.540; // NLO;--------------------------------
92 <    float SingleTopUCrossSection     = 10.6; // NLO;----------------------------------
90 > //Single top x-sections from here:
91 > //https://twiki.cern.ch/twiki/bin/view/CMS/SingleTopSigma
92 >    float SingleTopSCrossSection     = 4.63; // NLO;----------------------------------
93 >    float SingleTopTCrossSection     = 64.57; // NLO;--------------------------------
94 >    float SingleTopUCrossSection     = 15.74; // NLO;----------------------------------
95      float VVJetsCrossSection        = 4.8; // LO;-------------------------------------
96 <    float LM4CrossSection           = 2.537; // k*LO
96 >    float LM4CrossSection           = 1.879*1.35; // k*LO
97 >    float LM8CrossSection           = 0.7300*1.41 ; // k*LO  
98      //the following numbers are from the MadGraphStandardModel210Summary as linked on the GeneratorProduction2011 page    
99      //    float QCD50to100CrossSection=30000000; // not used
100      float QCD100to250CrossSection=7000000;
# Line 90 | Line 109 | void define_samples(bool showList, sampl
109      long totEventsQCD1000toInf  = 163185;
110      
111      long totEventsZjets         = 32005720;//summer11
112 +    long totEventsZjetsD6T      = 2543706;//summer11
113 +    long totEventsZjetsD6Tlow   = 181799;//summer11
114 + //    long totEventsZjetsEE       = 2738801;
115 + //    long totEventsZjetsMM       = 2148325;
116 + //    long totEventsZjetsTT       = 2032536;
117 +
118      //spring:2313911.0;//2743142.0; // DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola
119      long totEventsTTbar         =  3701872;//summer11
120      //1161621.0;//Spring//1144028.0; // PabloV13/TTJets_TuneZ2_7TeV-madgraph-tauola
121 <    long totEventsWJets         = 53739703;//summer11
121 >    long totEventsWJets         = 56674902;//summer11
122      //15010237.0;//spring//14350756.0; /// WJetsToLNu_TuneZ2_7TeV-madgraph-tauola.root
123      long totEventsZnunu         =  2106977;//2167964.0; //ZinvisibleJets_7TeV-madgraph.root
124      long totEventsVVJets        =   959076;//509072.0;
125      long totEventsSingleTopS    =   493868;//489472.0;
126      long totEventsSingleTopT    =   475460;//477610.0;
127      long totEventsSingleTopU    =   489417;//477599.0;
128 +    long totEventsLM4           =   218380;
129 +    long totEventsLM8           =   180960;
130  
131      
132      Int_t nice_blue  = TColor::GetColor("#2E9AFE");
# Line 131 | Line 158 | void define_samples(bool showList, sampl
158      char hostname[1023];
159      gethostname(hostname,1023);
160      
161 <    if((int)((string)hostname).find("falda")>-1) {
162 <      //not running locally on my computer
136 < //      allsamples.AddSample("/scratch/buchmann/ntuples2011/data/AllData_Jul6_1078pb_boolean.root","Data",1,1,true,false,0,kBlack); // this is the 1078 /pb version
137 < //      allsamples.AddSample("/scratch/fronga/v1.47/AllData_1.4fb_stripped.root","Data",1,1,true,false,0,kBlack);
138 < //      allsamples.AddSample("/scratch/buchmann/Data_Aug26__1936.root","Data",1,1,true,false,0,kBlack);
139 <        allsamples.AddSample("/scratch/buchmann/ntuples2011/data/AllData_Aug26__2048_ipb__NewFormat_Complete_official.root","Data",1,1,true,false,0,kBlack);
140 < //      allsamples.AddSample("/scratch/buchmann/ntuples2011/data/Data_v6_test.root","Data",1,1,true,false,0,kBlack);
141 < //      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
142 <      
143 <      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
144 <      allsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/WJetsToLNu_TuneZ2_7TeV-madgraph-tauola.root","WJets",totEventsWJets,WJetsCrossSection,false,false,3,wjets_color);//summer
145 <      allsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/VVJetsTo4L_TuneD6T_7TeV-madgraph-tauola.root","DiBosons",totEventsVVJets,VVJetsCrossSection,false,false,4,diboson_color);
146 <      allsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/TToBLNu_TuneZ2_s-channel_7TeV-madgraph_2.root","SingleTop",totEventsSingleTopS,SingleTopSCrossSection,false,false,2,singletop_color);
147 <      allsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/TToBLNu_TuneZ2_t-channel_7TeV-madgraph.root","SingleTop",totEventsSingleTopT,SingleTopTCrossSection,false,false,2,singletop_color);
148 <      allsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/TToBLNu_TuneZ2_tW-channel_7TeV-madgraph.root","SingleTop",totEventsSingleTopU,SingleTopUCrossSection,false,false,2,singletop_color);
149 <      allsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/ZinvisibleJets_7TeV-madgraph.root","Z nunu",totEventsZnunu,ZnunuCrossSection,false,false,4,diboson_color);
150 <      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
151 <      allsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM4_SUSY_sftsht_7TeV-pythia6.root","LM4",218536,1.879,false,true,7,lm_color);
152 <
153 <      signalsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM0_SUSY_sftsht_7TeV-pythia6.root","LM0",219796,38.93 ,false,true,3,lm_color);
154 <      signalsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM1_SUSY_sftsht_7TeV-pythia6.root","LM1",218176,4.888,false,true,4,lm_color);
155 <      signalsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM2_SUSY_sftsht_7TeV-pythia6.root","LM2",206336,0.6027,false,true,5,lm_color);
156 <      signalsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM3_SUSY_sftsht_7TeV-pythia6.root","LM3",215120,3.438,false,true,6,lm_color);
157 <      signalsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM4_SUSY_sftsht_7TeV-pythia6.root","LM4",218536,1.879,false,true,1,lm_color);
158 <      signalsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM5_SUSY_sftsht_7TeV-pythia6.root","LM5",223992,0.473,false,true,7,lm_color);
159 <      signalsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM8_SUSY_sftsht_7TeV-pythia6.root","LM8",180960,0.73,false,true,2,lm_color);
160 <      signalsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM9_SUSY_sftsht_7TeV-pythia6.root","LM9",227808,7.134,false,true,8,lm_color);
161 <      signalsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM11_SUSY_sftsht_7TeV-pythia6.root","LM11",218703,0.8236,false,true,9,lm_color);
162 <      signalsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM12_SUSY_sftsht_7TeV-pythia6.root","LM12",223491,4.414,false,true,10,lm_color);
163 <      signalsamples.AddSample("/scratch/buchmann/ntuples/MC/MixedSpringSummerMCepsPU_v2/LM13_SUSY_sftsht_7TeV-pythia6.root","LM13",173240,6.899,false,true,11,lm_color);
164 <
165 <    } else {
166 <
167 <      allsamples.AddSample("/shome/buchmann/ntuples/Data/AllData_2096pb_v158.root","Data",1,1,true,false,0,kBlack);
168 < //      allsamples.AddSample("/shome/buchmann/ntuples/Data/Latest_assemblage.root","Data",1,1,true,false,0,kBlack);
169 <
170 <      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
171 <
172 < //write_warning(__FUNCTION__,"TEMPORARILY USING POWHEG TT to 2l 2nu 2b sample!");
173 < //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
174 < //xs from https://twiki.cern.ch/twiki/bin/view/CMS/CrossSectionDetails
175 <
176 <      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
177 <      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);
178 <      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);
179 <      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);
180 <      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);
181 <      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);
182 <
183 <
184 <
185 <
186 <      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
187 < //      allsamples.AddSample("/scratch/buchmann/Skimmed_DY.root","ZJets",totEventsZjets,ZJetsCrossSection,false,false,6,dy_color);//summer11
188 <
189 <
190 <      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);
191 <
192 <
193 < /*      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);
194 <      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);
195 <      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);
196 <      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);*/
197 <      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);/*
198 <      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);*/
199 <      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);/*
200 <      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);
201 <      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);
202 <      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);
203 <      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);*/
204 <
205 <
206 <      // We do not use the following (rare) samples; nonetheless we include them in a separate category to prove that they are negligible
207 <      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
208 <      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
209 <      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
210 <      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
211 <      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
212 <      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
213 <      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
214 <      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
215 <      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
216 <      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
217 <
218 <
219 <
220 <
221 <    }
222 <    /**/
223 < //    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);
224 < //    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);
161 >    string prefix="/shome/";
162 >    if((int)((string)hostname).find("falda")>-1) prefix="/scratch/";
163      
164 < //    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"
165 <    //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"
166 < //    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"
167 <
164 > //    allsamples.AddSample(prefix+"buchmann/ntuples/Data/AllData_2809pb.root","Data",1,1,true,false,0,kBlack);
165 > //    if(PlottingSetup::RestrictToMassPeak) {
166 >        allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola_Summer11.root","ZJets",totEventsZjets,ZJetsCrossSection,false,false,7,dy_color);//summer11
167 > /*    } else {
168 >        allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/DYToEE_M-20_TuneZ2_7TeV-pythia6-tauola__Summer11-PU_S3_START42_V11-v2__AODSIM.root","ZJets",totEventsZjetsEE,ZJetsCrossSectionM20/3.0,false,false,7,dy_color);//summer11
169 >        allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/DYToMuMu_M-20_TuneZ2_7TeV-pythia6-tauola__Summer11-PU_S3_START42_V11-v2__AODSIM.root","ZJets",totEventsZjetsMM,ZJetsCrossSectionM20/3.0,false,false,7,dy_color);//summer11
170 >        allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/DYToTauTau_M-20_TuneZ2_7TeV-pythia6-tauola__Summer11-PU_S3_START42_V11-v2__AODSIM.root","ZJets",totEventsZjetsTT,ZJetsCrossSectionM20/3.0,false,false,7,dy_color);//summer11
171 >    }*/
172 >    allsamples.AddSample(prefix+"buchmann/ntuples/Data/NewSelectionForPaper/AllData_3p5_newselection__updated.root ","Data",1,1,true,false,0,kBlack);
173 >    allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTJets_TuneZ2_7TeV-madgraph-tauola_Summer11_2.root","TTbar",totEventsTTbar,TTbarCrossSection,false,false,1,ttbar_color);//summer11
174 >    allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/WJetsToLNu_TuneZ2_7TeV-madgraph-tauola_Summer11_2.root","WJets",totEventsWJets,WJetsCrossSection,false,false,3,wjets_color);//summer
175 >    allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/VVJetsTo4L_TuneD6T_7TeV-madgraph-tauola.root","DiBosons",totEventsVVJets,VVJetsCrossSection,false,false,4,diboson_color);
176 >    allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TToBLNu_TuneZ2_s-channel_7TeV-madgraph_2.root","SingleTop",totEventsSingleTopS,SingleTopSCrossSection,false,false,2,singletop_color);
177 >    allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TToBLNu_TuneZ2_t-channel_7TeV-madgraph.root","SingleTop",totEventsSingleTopT,SingleTopTCrossSection,false,false,2,singletop_color);
178 >    allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TToBLNu_TuneZ2_tW-channel_7TeV-madgraph.root","SingleTop",totEventsSingleTopU,SingleTopUCrossSection,false,false,2,singletop_color);
179 >    allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/ZinvisibleJets_7TeV-madgraph.root","Z nunu",totEventsZnunu,ZnunuCrossSection,false,false,4,diboson_color);
180 >    allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM4_SUSY_sftsht_7TeV-pythia6.root","LM4",totEventsLM4,LM4CrossSection,false,true,8,lm_color);
181  
231    define_SMS_sample(showList, allsamples, signalsamples, scansample);
232
233    if(showList) allsamples.ListSamples();
234
235 }
182  
183  
184 < void define_spring_samples(bool showList, samplecollection &allsamples)
185 < {
186 <  
187 <    float ZJetsCrossSection         = 3048.0; //NNLO----------------------------------
188 <    float TTbarCrossSection         = 165.0;//(NLO) ---- 165.0; // approx. NNLO-----
189 <    float WJetsCrossSection         = 31314.0;//NNLO-------3.131e4; //NNLO------------
190 <    float ZnunuCrossSection         = 5760.0;//NNLO -------4.5e+3; //(LO);------------
191 <    float SingleTopSCrossSection     = 1.49; // NLO;----------------------------------
192 <    float SingleTopTCrossSection     = 20.540; // NLO;--------------------------------
193 <    float SingleTopUCrossSection     = 10.6; // NLO;----------------------------------
194 <    float VVJetsCrossSection        = 4.8; // LO;-------------------------------------
195 <    float LM4CrossSection           = 2.537; // k*LO
196 <    //the following numbers are from the MadGraphStandardModel210Summary as linked on the GeneratorProduction2011 page    
197 <    //    float QCD50to100CrossSection=30000000; // not used
198 <    float QCD100to250CrossSection=7000000;
199 <    float QCD250to500CrossSection=171000;
200 <    float QCD500to1000CrossSection=5200;
201 <    float QCD1000toInfCrossSection=83.0;
202 <    
203 <    //    float totEventsQCD50to100=207418.0; // not used
204 <    float totEventsQCD100to250=638792.0;
205 <    float totEventsQCD250to500=344454.0;
206 <    float totEventsQCD500to1000=10e10;//unknown but now suppressed. (no event passes anyway)
207 <    float totEventsQCD1000toInf=163185.0;
208 <    
209 <    float totEventsZjets         = 32005720.0;//summer11
210 <    //spring:2313911.0;//2743142.0; // DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola
211 <    float totEventsTTbar         = 3701947.0;//summer11
212 <    //1161621.0;//Spring//1144028.0; // PabloV13/TTJets_TuneZ2_7TeV-madgraph-tauola
213 <    float totEventsWJets         = 53739703.0;//summer11
268 <    //15010237.0;//spring//14350756.0; /// WJetsToLNu_TuneZ2_7TeV-madgraph-tauola.root
269 <    float totEventsZnunu         = 2106977.0;//2167964.0; //ZinvisibleJets_7TeV-madgraph.root
270 <    float totEventsVVJets        = 959076.0;//509072.0;
271 <    float totEventsSingleTopS    = 493868.0;//489472.0;
272 <    float totEventsSingleTopT    = 475460.0;//477610.0;
273 <    float totEventsSingleTopU    = 489417.0;//477599.0;
184 >    if(!PlottingSetup::RestrictToMassPeak) {
185 > //      signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM0_SUSY_sftsht_7TeV-pythia6.root","LM0",219796,38.93*1.41 ,false,true,3,lm_color);
186 > //      signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM1_SUSY_sftsht_7TeV-pythia6.root","LM1",218176,4.888*1.34,false,true,4,lm_color);
187 > //      signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM2_SUSY_sftsht_7TeV-pythia6.root","LM2",206336,0.6027*1.33,false,true,5,lm_color);
188 >        signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM3_SUSY_sftsht_7TeV-pythia6.root","LM3",215120,3.438*1.40,false,true,6,lm_color);
189 >        signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM5_SUSY_sftsht_7TeV-pythia6.root","LM5",223992,0.473*1.34,false,true,7,lm_color);
190 >        signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM9_SUSY_sftsht_7TeV-pythia6.root","LM9",227808,7.134*1.48,false,true,8,lm_color);
191 >        signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM11_SUSY_sftsht_7TeV-pythia6.root","LM11",218703,0.8236*1.35,false,true,9,lm_color);
192 >        signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM12_SUSY_sftsht_7TeV-pythia6.root","LM12",223491,4.414*1.34,false,true,10,lm_color);
193 >        signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM13_SUSY_sftsht_7TeV-pythia6.root","LM13",173240,6.899*1.42,false,true,11,lm_color);
194 >    }
195 >    signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM4_SUSY_sftsht_7TeV-pythia6.root","LM4",totEventsLM4,LM4CrossSection,false,true,1,lm_color);
196 >    signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM8_SUSY_sftsht_7TeV-pythia6.root","LM8",totEventsLM8,LM8CrossSection,false,true,2,lm_color);
197 > /*
198 >    
199 >
200 >    // We do not use the following (rare) samples; nonetheless we include them in a separate category to prove that they are negligible
201 >    write_warning(__FUNCTION__,"Rare samples not yet reprocessed with new selectioN");/*
202 >    raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTWWTo2Lplus2Nu_NTP.root","Rare SS",39791,4.76E-5,false,false,6,rare_color);//summer11
203 >    raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/WWminusTo2L2Nu_NTP.root","Rare SS",51063,0.001559,false,false,6,rare_color);//summer11
204 >    raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTWWTo2Lminus2Nu_NTP.root","Rare SS",39496 ,4.719E-5,false,false,6,rare_color);//summer11
205 >    raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTWTo2Lplus2Nu_NTP.root","Rare SS",48764,0.006841,false,false,6,rare_color);//summer11
206 >    raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/WWplusTo2L2Nu_NTP.root","Rare SS",48244 ,0.004216,false,false,6,rare_color);//summer11
207 >    raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTZTo2Lplus2Nu_NTP.root","Rare SS",36740,0.002024,false,false,6,rare_color);//summer11
208 >    raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/WWWTo2Lplus2Nu_NTP.root","Rare SS",43322,0.008957,false,false,6,rare_color);//summer11
209 >    raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTWTo2Lminus2Nu_NTP.root","Rare SS",50435,0.002705,false,false,6,rare_color);//summer11
210 >    raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/WWWTo2Lminus2Nu_NTP.root","Rare SS",17147,0.004109,false,false,6,rare_color);//summer11
211 >    raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTZTo2Lminus2Nu_NTP.root","Rare SS",36781,0.001946,false,false,6,rare_color);//summer11
212 >    */
213 >    //scan files are only defined when necessary!
214  
275    
276    Int_t nice_blue  = TColor::GetColor("#2E9AFE");
277    Int_t nice_green = TColor::GetColor("#81f781");
278    Int_t nice_red   = TColor::GetColor("#F78181");
279    Int_t nice_pink  = TColor::GetColor("#F781BE");
280    Int_t nice_orange= TColor::GetColor("#F7BE81");//#FA9624
281    Int_t nice_black = TColor::GetColor("#2A0A1B");//black with a bit of purple :-)    
282    
283    
284    Int_t dy_color        = nice_blue;
285    Int_t wjets_color     = kGray;
286    Int_t ttbar_color     = nice_green;
287    Int_t singletop_color = nice_orange;
288    Int_t qcd_color       = nice_black;
289    Int_t diboson_color   = nice_pink;
290    Int_t lm_color        = kRed;
291    
292    dy_color            = kYellow;
293    wjets_color         = kGray;
294    ttbar_color         = kMagenta+2;
295    singletop_color     = kBlue;
296    qcd_color           = kPink;
297    diboson_color       = kGreen;
298    lm_color            = kViolet+7;
299  
300 //    allsamples.AddSample("/scratch/buchmann/AllData_Jun10___DCS_648pb_MoreTriggers4_DCS.root","Data",1,1,true,false,0,kBlack);
301 //    allsamples.AddSample("/scratch/buchmann/AllData_Jun10___486pb_MoreTriggers3.root","Data",1,1,true,false,0,kBlack);
302    //allsamples.AddSample("/scratch/buchmann/AllData_Jun17_160404-166861_real_complete_2.root","Data",1,1,true,false,0,kBlack);
303    allsamples.AddSample("/scratch/buchmann/ntuples2011/data/AllData_Jul6_1078pb.root","Data",1,1,true,false,0,kBlack);
304    
305 //    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
306    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
307 //    allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/WJetsToLNu_TuneZ2_7TeV-madgraph-tauola.root","WJets",totEventsWJets,WJetsCrossSection,false,false,3,wjets_color); //spring11
308    allsamples.AddSample("/scratch/fronga/MC_v1.29/WJetsToLNu_TuneZ2_7TeV-madgraph-tauola.root","WJets",totEventsWJets,WJetsCrossSection,false,false,3,wjets_color);//summer
309    allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/VVJetsTo4L_TuneD6T_7TeV-madgraph-tauola.root","Dibosons",totEventsVVJets,VVJetsCrossSection,false,false,4,diboson_color);
310    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);
311    allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/TToBLNu_TuneZ2_t-channel_7TeV-madgraph.root","SingleTop",totEventsSingleTopT,SingleTopTCrossSection,false,false,2,singletop_color);
312    allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/TToBLNu_TuneZ2_tW-channel_7TeV-madgraph.root","SingleTop",totEventsSingleTopU,SingleTopUCrossSection,false,false,2,singletop_color);
313    allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/ZinvisibleJets_7TeV-madgraph.root","Z nunu",totEventsZnunu,ZnunuCrossSection,false,false,4,diboson_color);
314 //    allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/QCD_TuneD6T_HT-100To250_7TeV-madgraph.root","QCD",totEventsQCD100to250,QCD100to250CrossSection,false,false,5,qcd_color);
315 //    allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/QCD_TuneD6T_HT-250To500_7TeV-madgraph.root","QCD",totEventsQCD250to500,QCD250to500CrossSection,false,false,5,qcd_color);
316 //    allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/QCD_TuneD6T_HT-500To1000_7TeV-madgraph.root","QCD",totEventsQCD500to1000,QCD500to1000CrossSection,false,false,5,qcd_color);
317 //    allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/QCD_TuneD6T_HT-1000_7TeV-madgraph.root","QCD",totEventsQCD1000toInf,QCD1000toInfCrossSection,false,false,5,qcd_color);
318 //    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
319    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
320    allsamples.AddSample("/scratch/buchmann/MCSpring2011PU__inclPFvarss__PURW_160404_166861/LM4_SUSY_sftsht_7TeV-pythia6.root","LM4",218536,1.879,false,true,7,lm_color);
321    
215      if(showList) allsamples.ListSamples();
216  
217 < }
325 <
217 > }

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines