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.12 by buchmann, Thu Mar 29 19:07:05 2012 UTC vs.
Revision 1.17 by buchmann, Fri May 4 12:04:30 2012 UTC

# Line 21 | Line 21 | namespace PlottingSetup {
21    string directoryname="Jellyfish_tests";
22    
23    bool RestrictToMassPeak=false; //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 +  bool FullMCAnalysis=false;
25    
26    float luminosity=4980;//4653.74;//p3523.18;//2096.0;//3172.73;//2096.0;//1936;//751.0;//486.0;//468.0//336.;//pb^{-1}
27   //  float luminosity=3523.18;//2096.0;//3172.73;//2096.0;//1936;//751.0;//486.0;//468.0//336.;//pb^{-1}
# Line 45 | Line 46 | namespace PlottingSetup {
46    TCut leptoncut("(pt1>20&&pt2>20)");
47    TCut passtrig("(passed_triggers||!is_data)"&&leptoncut);
48    
49 <  TCut openmasscut("mll>12"); // this is the mass cut used in the off peak analysis!
50 <  TCut openGenmasscut("genMll>12");
49 >  TCut openmasscut("mll>20"); // this is the mass cut used in the off peak analysis!
50 >  TCut openGenmasscut("genMll>20");
51    TCut opensidebandcut("mll>500&&mll<100"); // this won't let anything thru ...
52    
53    TCut Restrmasscut("abs(mll-91.2)<20"); // this is the mass cut used in the off peak analysis!
# Line 58 | Line 59 | namespace PlottingSetup {
59    
60    
61    //TCut cutmass("mll>2");
62 <  TCut basiccut("mll>2"&&leptoncut);//basically nothing.
62 > //  TCut basiccut("mll>2"&&leptoncut&&"(bTagProbTHighEff[0]>1.7)&&(bTagProbTHighEff[1]>1.7)");//basically nothing.
63 >  TCut basiccut("mll>2"&&leptoncut);//basically nothing
64    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.
65    //TCut jetqualitycut("(pfJetGoodNum>=2&&pfJetGoodID[0])&&(pfJetGoodNum>=2&&pfJetGoodID[1])");//now part of the basiccut
66    //TCut jetqualitycut("mll>0");
67    
68    TCut cutnJets("pfJetGoodNum>=3"&&basicqualitycut);
69 <  TCut cutnJetsJESdown("pfJetGoodNum25>=3"&&basicqualitycut);
70 <  TCut cutnJetsJESup("pfJetGoodNum35>=3"&&basicqualitycut);
69 >  TCut cutnJetsJESdown("pfJetGoodNumn1sigma>=3"&&basicqualitycut);
70 >  TCut cutnJetsJESup("pfJetGoodNump1sigma>=3"&&basicqualitycut);
71    TCut cutOSOF("(id1!=id2)&&(ch1*ch2<0)");
72    TCut cutOSSF("(id1==id2)&&(ch1*ch2<0)");
73    TCut sidebandcut(Restrsidebandcut);
# Line 80 | Line 82 | namespace PlottingSetup {
82    int ScanXzones=15; // number of zones in x for (mSUGRA) scans
83    int ScanYzones=15; // number of zones in y for (mSUGRA) scans
84    
85 <  string mSUGRAxsFile="../../Plotting/Modules/external/msugra_m0_m12_10_0_1_NLO_1.0.txt";
85 >  string mSUGRAxsFile="Plotting/Modules/external/msugra_m0_m12_10_0_1_NLO_1.0.txt";
86 >  string SMSReferenceXSFile="Plotting/Modules/external/reference_xSec_SMS-new.root";
87  
88    //scan types:
89    int mSUGRA=1;
# Line 118 | Line 121 | namespace PlottingSetup {
121    int JESdown=1;
122    int JESup=2;
123    
124 +  // Uncertainties ONPEAK
125 +  float zjetsestimateuncertONPEAK=0.25;
126 +  float emuncertONPEAK=0.25;
127 +  float emsidebanduncertONPEAK=0.25;
128 +  float eemmsidebanduncertONPEAK=0.25;
129 +
130 +  // Uncertainties OFFPEAK (iJZB)
131 +  float zjetsestimateuncertOFFPEAK=0.25;
132 +  float emuncertOFFPEAK=0.1;
133 +  float emsidebanduncertOFFPEAK=0.0;
134 +  float eemmsidebanduncertOFFPEAK=0.0;
135 +
136 +  
137    //some refinement: nicer color gradient
138    Double_t stops[5] = { 0.00, 0.34, 0.61, 0.84, 1.00 };
139    Double_t red[5]   = { 0.00, 0.00, 0.87, 1.00, 0.51 };
# Line 143 | Line 159 | namespace PlottingSetup {
159    bool computeJZBresponse=false;
160  
161    bool requireZ=true;//this is switched to "off" automatically when doing offpeak!
162 +  
163 +  string ScanSampleDirectory="SMS_T5zzh_newslots";//DileptonmSUGRAScan__pieces
164  
165 <  string ScanSampleDirectory="DileptonScan__pieces"; // possibilities (atm) : SMS_T5zz/, SMS_T5zzl/, SMS_T5zzh/, GMSB/, mSUGRA/     note: this string needs to either contain "SMS", "GMSB", or "mSUGRA"
165 > //  string ScanSampleDirectory="DileptonScan__pieces"; // possibilities (atm) : SMS_T5zz/, SMS_T5zzl/, SMS_T5zzh/, GMSB/, mSUGRA/     note: this string needs to either contain "SMS", "GMSB", or "mSUGRA"
166 > //  string ScanSampleDirectory="SMS_T1lh"; // possibilities (atm) : SMS_T5zz/, SMS_T5zzl/, SMS_T5zzh/, GMSB/, mSUGRA/     note: this string needs to either contain "SMS", "GMSB", or "mSUGRA"
167  
168    string FilterEfficiencyFile = "/shome/buchmann/JellyfishCBAF/DistributedModelCalculations/Limits/FilterEfficiencyv3.root";
169 +  string CMSSW_dir = "/shome/buchmann/final_production_2011/CMSSW_4_2_8/src/";
170  
171    ResultLibrary allresults;
172   }

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines