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

Comparing UserCode/cbrown/Development/Plotting/Modules/ActiveSamples.C (file contents):
Revision 1.1 by buchmann, Mon Jan 30 14:45:54 2012 UTC vs.
Revision 1.10 by buchmann, Wed May 23 15:13:58 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
1   #include <iostream>
2   #include <vector>
3   #include <sys/stat.h>
# Line 53 | Line 39 | void define_SMS_sample(bool showList, sa
39      
40        /// 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)
41        /// 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)
42 <        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"
42 > //      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"
43 > //    scansample.AddSample("/scratch/buchmann/buchmann/AllSamplesStartingAt10GeV_rest5_mSUGRA/mSUGRA_dilepton_m0-220to3000_m12-100to1000_tanb-10andA0-0_7TeV-Pythia6Z/output_8.root","mSUGRA",1,1,false,true,0,kRed); // there is only one sample in the scan; we give each event weight "1"
44 > //    scansample.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_400__M12_400.root","mSUGRA",1,1,false,true,0,kRed); // there is only one sample in the scan; we give each event weight "1"
45 >    
46 >    
47 >    
48  
49      if(loadall) {
50        for(int ix=0;ix<PlottingSetup::ScanXzones;ix++) {
# Line 65 | Line 56 | void define_SMS_sample(bool showList, sa
56          
57  
58      flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak
59 <    if(!PlottingSetup::RestrictToMassPeak) {
59 >    if(PlottingSetup::RestrictToMassPeak==false) {
60        PlottingSetup::cutmass=PlottingSetup::openmasscut;
61   //      PlottingSetup::sidebandcut=TCut("mll<2&&mll>3");//impossible cut - because we don't want any sidebands!
62        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).
63        PlottingSetup::genMassCut=PlottingSetup::openGenmasscut;//impossible cut - because we don't want any sidebands!
64 +      PlottingSetup::requireZ=false;
65      }
66 +    dout << "Currently configuration: The mass cut is \"" << PlottingSetup::cutmass << "\" and RestrictToMassPeak is currently set to " << PlottingSetup::RestrictToMassPeak << endl;
67 +    if(PlottingSetup::DoBTag) PlottingSetup::cutnJets=PlottingSetup::cutnJets&&PlottingSetup::bTagRequirement;
68   }
69  
70  
71 < void define_samples(bool showList, samplecollection &allsamples, samplecollection &signalsamples, samplecollection &scansample, samplecollection &raresample, samplecollection &systsamples)
71 > void define_samples(bool showList, samplecollection &allsamples, samplecollection &signalsamples, samplecollection &scansample, samplecollection &raresample, samplecollection &systsamples, samplecollection &qcdsamples)
72   {
73 +
74 + /*
75 + Good sources for XS:
76 + https://twiki.cern.ch/twiki/bin/view/CMS/ProductionSummer2011
77 + e.g.
78 + https://twiki.cern.ch/twiki/bin/view/CMS/ProductionSummer2011#PYTHIA6
79 +
80 + In general, all the xsec for datasets are in PREP, eg:
81 + http://cms.cern.ch/iCMS/prep/campaignmanagement?id=Summer11
82 +
83 + or
84 +
85 + http://cms.cern.ch/iCMS/prep/requestmanagement?dsn=QCD_Pt-15to30_TuneZ2_7TeV_pythia6&campid=Summer11
86 + */
87 +
88 + // source for 8 TeV: https://twiki.cern.ch/twiki/bin/viewauth/CMS/StandardModelCrossSectionsat8TeV
89 +
90      flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak
91      if(!PlottingSetup::RestrictToMassPeak) {
92        PlottingSetup::cutmass=PlottingSetup::openmasscut;
# Line 84 | Line 95 | void define_samples(bool showList, sampl
95      }
96    
97      float ZJetsCrossSection         = 3048.0; //NNLO----------------------------------
98 < //    float ZJetsCrossSectionM20      = 4998.0;//summer11
99 < //    float ZJetsCrossSectionD6Tlow   = 400.0;//summer11
98 >    if(PlottingSetup::is2012) ZJetsCrossSection=3503.71;
99 >    // at 8 TeV: 3503.71 (https://twiki.cern.ch/twiki/bin/view/CMS/Zpt8TeV)
100 >    float LowZJetsCrossSection      = 11908; //https://twiki.cern.ch/twiki/bin/view/CMS/TWikiTop2011DataMCTrig
101 >    
102      float TTbarCrossSection         = 165.0;//(NLO) ---- 165.0; // approx. NNLO-----
103 +    if(PlottingSetup::is2012) TTbarCrossSection=225.197;
104 +    // at 8 TeV: 225.197 (https://twiki.cern.ch/twiki/bin/view/CMS/Zpt8TeV)
105      float WJetsCrossSection         = 31314.0;//NNLO-------3.131e4; //NNLO------------
106      float ZnunuCrossSection         = 5760.0;//NNLO -------4.5e+3; //(LO);------------
107   //Single top x-sections from here:
# Line 95 | Line 110 | void define_samples(bool showList, sampl
110      float SingleTopTCrossSection     = 64.57; // NLO;--------------------------------
111      float SingleTopUCrossSection     = 15.74; // NLO;----------------------------------
112      float VVJetsCrossSection        = 4.8; // LO;-------------------------------------
113 <    float LM4CrossSection           = 1.879*1.35; // k*LO
114 <    float LM8CrossSection           = 0.7300*1.41 ; // k*LO  
113 > /*    float LM0CrossSection           =         38.93 *1.41 ;
114 >    float LM1CrossSection           = 4.888 *1.34 ;
115 >    float LM2CrossSection           =0.6027 *1.33 ;*/
116 >    float LM3CrossSection           = 3.438*1.40 ; // LO*k
117 > /*    float LM4CrossSection           = 1.879*1.35; // k*LO
118 >    float LM5CrossSection           = 0.473*1.34; // k*LO
119 >    float LM6CrossSection           = 0.3104*1.30; // k*LO
120 >    float LM7CrossSection           = 1.209*1.11; // k*LO
121 >    float LM8CrossSection           = 0.7300*1.41; // k*LO
122 >    float LM9CrossSection           = 7.134*1.48; // k*LO
123 >    float LM10CrossSection           = 0.04778*1.22 ; // k*LO
124 >    float LM11CrossSection           = 0.8236 *1.35 ; // k*LO
125 >    float LM12CrossSection           = 4.414 *1.34 ; // k*LO
126 >    float LM13CrossSection           = 6.899*1.42 ; // k*LO*/
127      //the following numbers are from the MadGraphStandardModel210Summary as linked on the GeneratorProduction2011 page    
128      //    float QCD50to100CrossSection=30000000; // not used
129 <    float QCD100to250CrossSection=7000000;
130 <    float QCD250to500CrossSection=171000;
131 <    float QCD500to1000CrossSection=5200;
132 <    float QCD1000toInfCrossSection=83.0;
129 >    //float QCD100to250CrossSection=7000000;
130 >    //float QCD250to500CrossSection=171000;
131 >    //float QCD500to1000CrossSection=5200;
132 >    //float QCD1000toInfCrossSection=83.0;
133      
134      //    float totEventsQCD50to100=207418.0; // not used
135 <    long totEventsQCD100to250   = 638792;
136 <    long totEventsQCD250to500   = 344454;
137 <    long totEventsQCD500to1000  = (long)10e10;  //unknown but now suppressed. (no event passes anyway)
138 <    long totEventsQCD1000toInf  = 163185;
135 >    //long totEventsQCD100to250   = 638792;
136 >    //long totEventsQCD250to500   = 344454;
137 >    //long totEventsQCD500to1000  = (long)10e10;  //unknown but now suppressed. (no event passes anyway)
138 >    //long totEventsQCD1000toInf  = 163185;
139      
140      long totEventsZjets         = 32005720;//summer11
141 <    long totEventsZjetsD6T      = 2543706;//summer11
142 <    long totEventsZjetsD6Tlow   = 181799;//summer11
143 < //    long totEventsZjetsEE       = 2738801;
144 < //    long totEventsZjetsMM       = 2148325;
118 < //    long totEventsZjetsTT       = 2032536;
141 >    long totEventslowZjets      = 10272470;
142 >
143 > //    long totEventsZjetsD6T      = 2543706;//summer11
144 > //    long totEventsZjetsD6Tlow   = 181799;//summer11
145  
146      //spring:2313911.0;//2743142.0; // DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola
147      long totEventsTTbar         =  3701947;//summer11
# Line 127 | Line 153 | void define_samples(bool showList, sampl
153      long totEventsSingleTopS    =   493868;//489472.0;
154      long totEventsSingleTopT    =   475460;//477610.0;
155      long totEventsSingleTopU    =   489417;//477599.0;
156 <    long totEventsLM4           =   218380;
157 <    long totEventsLM8           =   180960;
156 > /*    long totEventsLM0           =   1234567;//bad guess, will be replaced.
157 >    long totEventsLM1           =   1234567;//bad guess, will be replaced.
158 >    long totEventsLM2           =   1234567;//bad guess, will be replaced.*/
159 >    long totEventsLM3           =   220000;//bad guess, will be replaced.
160 > /*    long totEventsLM4           =   218380;
161 >    long totEventsLM5           =   220000;
162 >    long totEventsLM6           =   220000;
163 >    long totEventsLM7           =   217800;
164 >    long totEventsLM8           =   217800;
165 >    long totEventsLM9           =   220000;
166 >    long totEventsLM10          =   1234567;
167 >    long totEventsLM11          =   216999;
168 >    long totEventsLM12          =   217400;
169 >    long totEventsLM13          =   1234567;*/
170  
171      
172      Int_t nice_blue  = TColor::GetColor("#2E9AFE");
173      Int_t nice_green = TColor::GetColor("#81f781");
174 <    Int_t nice_red   = TColor::GetColor("#F78181");
174 > //    Int_t nice_red   = TColor::GetColor("#F78181");
175      Int_t nice_pink  = TColor::GetColor("#F781BE");
176      Int_t nice_orange= TColor::GetColor("#F7BE81");//#FA9624
177      Int_t nice_black = TColor::GetColor("#2A0A1B");//black with a bit of purple :-)    
# Line 155 | Line 193 | void define_samples(bool showList, sampl
193      diboson_color       = kGreen+2;
194      lm_color            = kViolet+7;
195  
196 <    Int_t rare_color            = nice_blue;
196 > //    Int_t rare_color          = nice_blue;
197    
198      char hostname[1023];
199      gethostname(hostname,1023);
# Line 163 | Line 201 | void define_samples(bool showList, sampl
201      string prefix="/shome/";
202      if((int)((string)hostname).find("falda")>-1) prefix="/scratch/";
203      
204 <    allsamples.AddSample(prefix+"buchmann/ntuples/Data/AllData_Complete_4653_ipb_PAPER.root","Data",1,1,true,false,0,kBlack);
205 <
206 <
207 <    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);
208 <    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
209 <    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
210 <    allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/VVJetsTo4L_TuneD6T_7TeV-madgraph-tauola.root","WW/WZ/ZZ",totEventsVVJets,VVJetsCrossSection,false,false,4,diboson_color);
211 <    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);
212 <    allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/TToBLNu_TuneZ2_t-channel_7TeV-madgraph.root","Single-top",totEventsSingleTopT,SingleTopTCrossSection,false,false,2,singletop_color);
213 <    allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/TToBLNu_TuneZ2_tW-channel_7TeV-madgraph.root","Single-top",totEventsSingleTopU,SingleTopUCrossSection,false,false,2,singletop_color);
214 <    allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/ZinvisibleJets_7TeV-madgraph.root","Z nunu",totEventsZnunu,ZnunuCrossSection,false,false,4,diboson_color);
215 <    allsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/LM4_SUSY_sftsht_7TeV-pythia6.root","LM4",totEventsLM4,LM4CrossSection,false,true,8,lm_color);
216 <
217 <
218 <
204 >    if(!PlottingSetup::is2012) {
205 >      allsamples.AddSample(prefix+"buchmann/ntuples/Data/AllData_Complete_4653_ipb_PAPER.root","Data",1,1,true,false,0,kBlack);
206 >      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);
207 >      allsamples.AddSample(prefix+"buchmann/ntuples/MC/Partial_DYJetsToLL_M-10To50_TuneZ2_7TeV-madgraph__Fall11-PU_S6_START42_V14B-v1__AODSIM.root","Z+Jets",totEventslowZjets,LowZJetsCrossSection,false,false,7,dy_color);
208 >      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
209 >      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
210 >      allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/VVJetsTo4L_TuneD6T_7TeV-madgraph-tauola.root","WW/WZ/ZZ",totEventsVVJets,VVJetsCrossSection,false,false,4,diboson_color);
211 >      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);
212 >      allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/TToBLNu_TuneZ2_t-channel_7TeV-madgraph.root","Single-top",totEventsSingleTopT,SingleTopTCrossSection,false,false,2,singletop_color);
213 >      allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/TToBLNu_TuneZ2_tW-channel_7TeV-madgraph.root","Single-top",totEventsSingleTopU,SingleTopUCrossSection,false,false,2,singletop_color);
214 >      allsamples.AddSample(prefix+"/buchmann/ntuples/MC/NewSelectionForPaper/CorrectFull2011PileUp/ZinvisibleJets_7TeV-madgraph.root","Z nunu",totEventsZnunu,ZnunuCrossSection,false,false,4,diboson_color);
215 >      allsamples.AddSample(prefix+"/buchmann/ntuples/MC/JZBplus/WGstarToLNu2E_TuneZ2_7TeV-madgraph-tauola__Summer11-PU_S4_START42_V11-v1__AODSIM.root","W+#gamma",266882,5.546 ,false,false,5,kGreen);
216 >      allsamples.AddSample(prefix+"/buchmann/ntuples/MC/JZBplus/GVJets_7TeV-madgraph-Summer11-PU_S4.root","#gamma+VJets",1067879,56.64 ,false,false,6,kRed);
217 >      //allsamples.AddSample("/shome/buchmann/ntuples/MC/LM4_SUSY_sftsht_7TeV-pythia6__Summer11-PU_S4_START42_V11-v2_novtx_fixed3___WITHGENINFO.root","LM4",totEventsLM4,LM4CrossSection,false,true,8,lm_color);
218 >      allsamples.AddSample("/shome//lbaeni/jzb/ntuples/MC/JZBplus/SignalSamples/LM8_SUSY_sftsht_7TeV-pythia6.root","LM3",totEventsLM3,LM3CrossSection,false,true,8,lm_color);
219 >    } else {
220 >      
221 >      write_info(__FUNCTION__,"Setting up 2012 configuration - once you've stopped looking at 2011 data please make this the default!");
222 >      cutWeight=TCut("weight3D");
223 >      basicqualitycut=TCut("(pfJetGoodNum40>=2&&pfJetGoodID[0]!=0)&&(pfJetGoodNum40>=2&&pfJetGoodID[1]!=0)"&&basiccut);//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.
224 >      cutnJets=TCut("pfJetGoodNum40>=3"&&basicqualitycut);
225 >      PlottingSetup::luminosity=PlottingSetup::luminosity2012;
226 >      PlottingSetup::lumiuncert=PlottingSetup::lumiuncert2012;
227 >      
228 >      // cJZB response correction (705 /pb) :
229 >      PlottingSetup::jzbvariabledata="(jzb[1]+0.0656534*pt)";
230 >      PlottingSetup::jzbvariablemc="(jzb[1]+0.0322238*pt)";
231 >
232 >      allsamples.AddSample("/shome/buchmann/ntuples/Data/DataProcessed705ipb_btagfix.root","Data",1,1,true,false,0,kBlack);
233 >      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/ETHchen/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball-Summer12-PU_S7_START52_V9-v2.root","Z+Jets",0,ZJetsCrossSection,false,false,7,dy_color);
234 >      allsamples.AddSample("/shome/buchmann/ntuples/MC8tev/ETHchen/TTJets_TuneZ2star_8TeV-madgraph-tauola-Summer12-PU_S7_START52_V9-v1.root","t#bar{t}",0,TTbarCrossSection,false,false,1,ttbar_color);
235 >    }
236      if(!PlottingSetup::RestrictToMassPeak) {
237   //      signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM0_SUSY_sftsht_7TeV-pythia6.root","LM0",219796,38.93*1.41 ,false,true,3,lm_color);
238   //      signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM1_SUSY_sftsht_7TeV-pythia6.root","LM1",218176,4.888*1.34,false,true,4,lm_color);
239   //      signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM2_SUSY_sftsht_7TeV-pythia6.root","LM2",206336,0.6027*1.33,false,true,5,lm_color);
240 <        signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM3_SUSY_sftsht_7TeV-pythia6.root","LM3",215120,3.438*1.40,false,true,6,lm_color);
240 > write_warning(__FUNCTION__,"CURRENTLY ALL SIGNAL SAMPLES HAVE BEEN DEACTIVATED (EXCEPT LM4 N LM8)");
241 > /*      signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM3_SUSY_sftsht_7TeV-pythia6.root","LM3",215120,3.438*1.40,false,true,6,lm_color);
242          signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM5_SUSY_sftsht_7TeV-pythia6.root","LM5",223992,0.473*1.34,false,true,7,lm_color);
243          signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM9_SUSY_sftsht_7TeV-pythia6.root","LM9",227808,7.134*1.48,false,true,8,lm_color);
244          signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM11_SUSY_sftsht_7TeV-pythia6.root","LM11",218703,0.8236*1.35,false,true,9,lm_color);
245          signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM12_SUSY_sftsht_7TeV-pythia6.root","LM12",223491,4.414*1.34,false,true,10,lm_color);
246 <        signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM13_SUSY_sftsht_7TeV-pythia6.root","LM13",173240,6.899*1.42,false,true,11,lm_color);
246 >        signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM13_SUSY_sftsht_7TeV-pythia6.root","LM13",173240,6.899*1.42,false,true,11,lm_color);*/
247      }
248 <    signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM4_SUSY_sftsht_7TeV-pythia6.root","LM4",totEventsLM4,LM4CrossSection,false,true,1,lm_color);
249 <    signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM8_SUSY_sftsht_7TeV-pythia6.root","LM8",totEventsLM8,LM8CrossSection,false,true,2,lm_color);
248 > //    signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM4_SUSY_sftsht_7TeV-pythia6.root","LM4",totEventsLM4,LM4CrossSection,false,true,1,lm_color);
249 > //    signalsamples.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/LM8_SUSY_sftsht_7TeV-pythia6.root","LM8",totEventsLM8,LM8CrossSection,false,true,2,lm_color);
250  
251 < //    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);
251 > /*      signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_200__M12_400.root","mSUGRA_M0_200__M12_400",0,1,false,true,1,lm_color);
252 >      signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_200__M12_300.root","mSUGRA_M0_200__M12_300",0,1,false,true,2,lm_color);
253 >      signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_400__M12_400.root","mSUGRA_M0_400__M12_400",0,1,false,true,3,lm_color);
254 >      signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_400__M12_300.root","mSUGRA_M0_400__M12_300",0,1,false,true,4,lm_color);
255 >      signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_600__M12_400.root","mSUGRA_M0_600__M12_400",0,1,false,true,5,lm_color);
256 >      signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_600__M12_300.root","mSUGRA_M0_600__M12_300",0,1,false,true,6,lm_color);
257 >      signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_800__M12_400.root","mSUGRA_M0_800__M12_400",0,1,false,true,7,lm_color);
258 >      signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_800__M12_300.root","mSUGRA_M0_800__M12_300",0,1,false,true,8,lm_color);
259 >      signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_1000__M12_400.root","mSUGRA_M0_1000__M12_400",0,1,false,true,9,lm_color);
260 >      signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_1000__M12_300.root","mSUGRA_M0_1000__M12_300",0,1,false,true,10,lm_color);
261 >      signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_1500__M12_400.root","mSUGRA_M0_1500__M12_400",0,1,false,true,11,lm_color);
262 >      signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_1500__M12_300.root","mSUGRA_M0_1500__M12_300",0,1,false,true,12,lm_color);*/
263 > //signalsamples.AddSample("/shome//buchmann/ntuples/SMS_T1lh/SMS_MassGlu_400__MassLSP_200.root","SMS_MassGlu_400__MassLSP_200",0,1,false,true,12,lm_color);
264 > //signalsamples.AddSample("/shome//buchmann/ntuples/SMS_T1lh/SMS_MassGlu_500__MassLSP_200.root","SMS_MassGlu_500__MassLSP_200",0,1,false,true,12,lm_color);
265 > //signalsamples.AddSample("/shome//buchmann/ntuples/SMS_T1lh/SMS_MassGlu_600__MassLSP_200.root","SMS_MassGlu_600__MassLSP_200",0,1,false,true,12,lm_color);
266 > signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_2000__M12_200.root","M_{0}=2000,M_{1/2}=200",0,1,false,true,12,lm_color);
267 > signalsamples.AddSample("/shome/lbaeni/jzb/DileptonmSUGRAScan__pieces/mSUGRA_M0_300__M12_800.root","M_{0}=300,M_{1/2}=800",0,1,false,true,12,lm_color);
268 >
269 >
270 > /*signalsamples.AddSample("/shome//buchmann/ntuples/SMS_T1lh/SMS_MassGlu_400__MassLSP_300.root","SMS_MassGlu_400__MassLSP_200",0,1,false,true,12,lm_color);
271 > signalsamples.AddSample("/shome//buchmann/ntuples/SMS_T1lh/SMS_MassGlu_500__MassLSP_300.root","SMS_MassGlu_400__MassLSP_200",0,1,false,true,12,lm_color);
272 > signalsamples.AddSample("/shome//buchmann/ntuples/SMS_T1lh/SMS_MassGlu_600__MassLSP_300.root","SMS_MassGlu_400__MassLSP_200",0,1,false,true,12,lm_color);*/
273  
274  
275 + //    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);
276 +
277      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);
278      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
279  
# Line 203 | Line 282 | void define_samples(bool showList, sampl
282      
283  
284      // We do not use the following (rare) samples; nonetheless we include them in a separate category to prove that they are negligible
285 <    write_warning(__FUNCTION__,"Rare samples not yet reprocessed with new selectioN");/*
285 >    write_warning(__FUNCTION__,"Rare samples not yet reprocessed with new selectioN");
286      raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTWWTo2Lplus2Nu_NTP.root","Rare SS",39791,4.76E-5,false,false,6,rare_color);//summer11
287      raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/WWminusTo2L2Nu_NTP.root","Rare SS",51063,0.001559,false,false,6,rare_color);//summer11
288      raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTWWTo2Lminus2Nu_NTP.root","Rare SS",39496 ,4.719E-5,false,false,6,rare_color);//summer11
# Line 216 | Line 295 | void define_samples(bool showList, sampl
295      raresample.AddSample(prefix+"buchmann/ntuples/MC/NewSelectionForPaper/TTZTo2Lminus2Nu_NTP.root","Rare SS",36781,0.001946,false,false,6,rare_color);//summer11
296      */
297      //scan files are only defined when necessary!
298 +    
299 + // XS taken from:
300 + //                https://twiki.cern.ch/twiki/bin/view/CMS/QCDGammaPlusJetQCD11005
301  
302 +    if(!PlottingSetup::RestrictToMassPeak) {
303 +      PlottingSetup::cutmass=PlottingSetup::openmasscut;
304 +      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).
305 +      PlottingSetup::genMassCut=PlottingSetup::openGenmasscut;
306 +      PlottingSetup::requireZ=false;
307 +    }
308 +    
309      if(showList) allsamples.ListSamples();
310  
311   }

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines