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

Comparing UserCode/cbrown/Development/Plotting/Modules/Setup.C (file contents):
Revision 1.36 by buchmann, Thu Oct 11 08:57:44 2012 UTC vs.
Revision 1.43 by buchmann, Fri Jun 28 15:02:21 2013 UTC

# Line 22 | Line 22 | namespace PlottingSetup {
22    
23    bool RestrictToMassPeak=true; //if you want to switch between offpeak ("false") and onpeak ("true") analysis please use this switch; the masscut below will be adapted automatically when adding samples :-)
24    
25 <  bool UseSidebandsForcJZB=true;
25 >  bool UseSidebandsForcJZB=false;
26    bool FullMCAnalysis=false;
27    bool DoBTag=false;
28    
# Line 33 | Line 33 | namespace PlottingSetup {
33   //  float luminosity=2096.0;//3172.73;//2096.0;//1936;//751.0;//486.0;//468.0//336.;//pb^{-1}
34    float lumiuncert=0.022;//0.045;// to be indicated in [0,1] range, e.g. for 4% write 0.04
35    
36 <  float luminosity2012=9200;
36 >  //float luminosity2012=19490;//9200;
37 >  float luminosity2012=9200;//9200;
38    float lumiuncert2012=0.022;
39  
40   //  string jzbvariabledata="jzb[1]+0.06*pt";
# Line 42 | Line 43 | namespace PlottingSetup {
43    string jzbvariablemc="(jzb[1]+0.036*pt)";
44    float jzbHigh = 450.; // Range for JZB plots
45    float iMllLow = 20.; // Range for Edge Fitting (and mll plot) for iJZB
46 <  float iMllHigh = 400.; // Range for Edge Fitting (and mll plot) for iJZB
46 >  float iMllHigh = 300.; // Range for Edge Fitting (and mll plot) for iJZB
47  
48    samplecollection allsamples("completesamplecollection");
49    samplecollection qcdsamples("QCDcollection");
# Line 51 | Line 52 | namespace PlottingSetup {
52    samplecollection scansample("scansamplecollection");
53    samplecollection raresample("raresamplecollection");
54    samplecollection comparesamples("comparesamplecollection");
55 +  samplecollection SingleLeptonData("SingleLeptonData");
56    int data=1;
57    int mc=0;
58    int mcwithsignal=2;
59 <  TCut leptoncut("(abs(eta1)<2.4 && abs(eta2)<2.4 && pt1>20 && pt2>20 )");
60 < //   TCut passtrig("(passed_filter && (!is_data||  (passed_triggers && ( (id1==id2)*(id1==0)*passed_ee + (id1==id2)*(id1==1)*passed_mm + (id1!=id2)*passed_em) )   ))"&&leptoncut);
61 <  TCut passtrig("((passed_triggers *  (  (id1==id2)*(id1==0)*(trigger_bit&1) + (id1==id2)*(id1==1)*(trigger_bit&2) + (id1!=id2)*(trigger_bit&4|trigger_bit&8)   )   ||!is_data))"&&leptoncut);
59 >   TCut leptoncut("(abs(eta1)<1.4 && abs(eta2)<1.4 && pt1>20 && pt2>20)");
60 >  TCut passtrig("((passed_triggers *  (  (id1==id2)*(id1==0)*(trigger_bit&1) + (id1==id2)*(id1==1)*(trigger_bit&2) + (id1!=id2)*(trigger_bit&4|trigger_bit&8)   )  ||!is_data))"&&leptoncut);
61 > //   TCut passtrig("((passed_triggers ))"&&leptoncut);
62 > //  TCut passtrig(leptoncut);
63    
64    TCut openmasscut("mll>20"); // this is the mass cut used in the off peak analysis!
65    TCut openGenmasscut("genMll>15");
# Line 75 | Line 78 | namespace PlottingSetup {
78    //TCut cutmass("mll>2");
79   //  TCut basiccut("mll>2"&&leptoncut&&"(bTagProbTHighEff[0]>1.7)&&(bTagProbTHighEff[1]>1.7)");//basically nothing.
80    TCut basiccut("mll>2"&&leptoncut);//basically nothing
81 <  TCut basicqualitycut("(pfJetGoodNum>=2&&pfJetGoodID[0]!=0)&&(pfJetGoodNum>=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.
81 >  TCut basicqualitycut("(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.
82    
83    TCut cutnJets("pfJetGoodNum40>=3"&&basicqualitycut);
84    TCut cutnJetsJESdown("pfJetGoodNumn1sigma>=3"&&basicqualitycut);
# Line 83 | Line 86 | namespace PlottingSetup {
86    TCut cutOSOF("(id1!=id2)&&(ch1*ch2<0)");
87    TCut cutOSSF("(id1==id2)&&(ch1*ch2<0)");
88    TCut sidebandcut(Restrsidebandcut);
89 +  
90 +  TCut TriLepSF("tri_id2==tri_id3");
91 +  TCut TriLepOF("tri_id2!=tri_id3");
92 +
93  
94    // SUSY scan parameters
95    float mglustart=25;float mgluend=1200;float mglustep=25;
# Line 135 | Line 142 | namespace PlottingSetup {
142    
143    // Uncertainties ONPEAK
144    float zjetsestimateuncertONPEAK=0.2;
145 <  float emuncertONPEAK=0.1;
146 <  float emsidebanduncertONPEAK=0.1;
147 <  float eemmsidebanduncertONPEAK=0.1;
145 >  float emuncertONPEAK=0.10;
146 >  float emsidebanduncertONPEAK=0.10;
147 >  float eemmsidebanduncertONPEAK=0.10;
148  
149    // Uncertainties OFFPEAK (iJZB)
150 <  float zjetsestimateuncertOFFPEAK=0.25;
151 <  float emuncertOFFPEAK=0.1;
150 >  float zjetsestimateuncertOFFPEAK=0.2;
151 >  float emuncertOFFPEAK=0.10;
152    float emsidebanduncertOFFPEAK=0.0;
153    float eemmsidebanduncertOFFPEAK=0.0;
154  
# Line 181 | Line 188 | namespace PlottingSetup {
188    string CMSSW_dir = "/shome/buchmann/final_production_2011/CMSSW_4_2_8/src/";
189  
190    ResultLibrary allresults;
191 +  
192 +  bool IsZbAnalysis=false;
193 +  bool IsWZAnalysis=false;
194 +  bool IsJZBAnalysis=true;
195 +  
196 +  int do_CleanCache=false;
197   }

Diff Legend

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