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

Comparing UserCode/cbrown/AnalysisFramework/Plotting/various_studies.C (file contents):
Revision 1.7 by buchmann, Wed Jul 20 09:15:27 2011 UTC vs.
Revision 1.8 by buchmann, Wed Jul 20 14:34:53 2011 UTC

# Line 14 | Line 14 | to compile: just run make ... that shoul
14   #include "Modules/ActiveSamples.C"
15   #include "Modules/UpperLimitsWithShape.C"
16   #include "Modules/Setup.C"
17 + #include "Modules/external/cl95cms.c"
18   #include "Modules/Plotting_Functions.C" //also included for peak finding etc.
19   #include "Modules/StudyModule.C" //also included for peak finding etc.
20 + #include "Modules/ExperimentalModule.C"
21 + #include "Modules/ResultModule.C"
22 + #include "Modules/LimitCalculation.C"
23  
24   #include <TCut.h>
25   #include <TROOT.h>
# Line 66 | Line 70 | int main()
70    bool do_jzb_correction=false; //use this to find out the correction factor; originally this was done automatically but now you should do this separately and update Modules/Setup.C
71    bool do_find_sideband_definition=false;
72    bool do_pick_up_events=true; ///DONE
73 +  bool do_ttbar_with_shapes=false;
74 +  bool do_upper_limits=false;
75    
76    //**** part 1 : peak finding
77    float MCPeak=0,MCPeakError=0,DataPeak=0,DataPeakError=0,MCSigma=10,DataSigma=10;
# Line 97 | Line 103 | int main()
103    
104    if(study_sidebands) look_at_sidebands(mcjzb.str(),datajzb.str());
105    
106 +  vector<float>jzb_limit_bins;
107 +  jzb_limit_bins.push_back(50);jzb_limit_bins.push_back(75);jzb_limit_bins.push_back(100);
108 +  jzb_limit_bins.push_back(150);jzb_limit_bins.push_back(200);jzb_limit_bins.push_back(500);
109 +
110 +  
111 +  if(do_ttbar_with_shapes) prepare_ttbar_limits(mcjzb.str(),datajzb.str(),DataPeakError,MCPeakError,jzb_limit_bins);
112 +  
113 +  if(do_upper_limits) calculate_upper_limits(mcjzb.str(),datajzb.str());
114 +  
115    if(do_pick_up_events) {
116      
117      dout << "Observed: " << endl;

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines