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> |
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; |
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; |