| 15 |
|
#include "Modules/ActiveSamples.C" |
| 16 |
|
#include "Modules/UpperLimitsWithShape.C" |
| 17 |
|
#include "Modules/Setup.C" |
| 18 |
– |
#include "Modules/Var_Studies.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 |
|
|
| 48 |
|
gStyle->SetTextFont(42); |
| 49 |
|
bool showList=false; |
| 50 |
|
set_treename("events");//you can set the treename here to be used; options are "events" (for reco) for "PFevents" (for particle flow) |
| 51 |
< |
define_samples(showList,allsamples); |
| 51 |
> |
define_samples(showList,allsamples,signalsamples); |
| 52 |
|
setlumi(luminosity); |
| 53 |
|
do_png(true); |
| 54 |
|
do_pdf(false); |
| 55 |
|
do_eps(false); |
| 56 |
|
do_C(false); |
| 57 |
< |
set_directory("comparison__corrected__fixed");//here you can give it a name, for instance "1234pb" or whatever you like :-) |
| 57 |
> |
set_directory("summer_mc");//here you can give it a name, for instance "1234pb" or whatever you like :-) |
| 58 |
|
setessentialcut(essential);//this sets the essential cut; this one is used in the draw command so it is AUTOMATICALLY applied everywhere. IMPORTANT: Do NOT store weights here! |
| 59 |
|
stringstream resultsummary; |
| 60 |
|
|
| 61 |
|
//what to do : |
| 62 |
|
bool do_all=false;/// DONE |
| 63 |
|
bool do_peak_finding=false; /// DONE |
| 65 |
– |
bool do_kinematic_variables=false; ///DONE |
| 66 |
– |
bool do_kinematic_PF_variables=false; ///DONE |
| 67 |
– |
bool do_lepton_comparison=false; ///DONE |
| 68 |
– |
bool do_jzb_plots=false; /// DONE |
| 69 |
– |
bool do_pred=false; /// DONE |
| 70 |
– |
bool do_ratio=false; /// DONE |
| 71 |
– |
bool do_signal_bg_comparison_plot=false; /// DONE |
| 64 |
|
bool calculate_pred_and_observed=false; /// DONE |
| 65 |
< |
bool calculate_yields=false; |
| 74 |
< |
|
| 75 |
< |
//-------------------------------------------- |
| 76 |
< |
// not so standard functions |
| 77 |
< |
bool do_upper_limits=false; |
| 78 |
< |
bool do_model_scan=false; |
| 79 |
< |
bool do_zjet_ttbar_shapes=false; /// DONE |
| 65 |
> |
bool study_sidebands=true; /// FOR REAL |
| 66 |
|
bool do_test=false;/// DONE (just any test you wish) |
| 81 |
– |
bool study_sidebands=false; |
| 82 |
– |
bool do_pick_up_events=true; |
| 83 |
– |
|
| 84 |
– |
//-------------------------------------------- |
| 67 |
|
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 |
| 68 |
+ |
bool do_find_sideband_definition=false; |
| 69 |
|
|
| 70 |
|
//**** part 1 : peak finding |
| 71 |
|
float MCPeak=0,MCPeakError=0,DataPeak=0,DataPeakError=0,MCSigma=10,DataSigma=10; |
| 83 |
|
cout << " Data : " << datajzb.str() << endl; |
| 84 |
|
cout << " MC : " << datajzb.str() << endl; |
| 85 |
|
|
| 86 |
< |
//jzb response correction |
| 87 |
< |
if(do_jzb_correction) find_correction_factors(jzbvariabledata,jzbvariablemc); |
| 88 |
< |
cout << "Please update Modules/Setup.C to reflect the following values (round them first ... )" <<endl; |
| 89 |
< |
cout << "Corrected JZB variable definition " << endl << " Data: " << jzbvariabledata << " and MC: " << jzbvariablemc << endl; |
| 90 |
< |
cout << endl; |
| 91 |
< |
cout << "If you're feeling lazy, copy & paste this : " << endl; |
| 92 |
< |
cout << " string jzbvariabledata=\""<<jzbvariabledata<<"\";"<<endl; |
| 93 |
< |
cout << " string jzbvariablemc=\"" <<jzbvariablemc<<"\";"<<endl; |
| 86 |
> |
if(do_find_sideband_definition) find_sideband_definition(); |
| 87 |
> |
|
| 88 |
> |
if(do_jzb_correction) { |
| 89 |
> |
find_correction_factors(jzbvariabledata,jzbvariablemc); |
| 90 |
> |
cout << "Please update Modules/Setup.C to reflect the following values (round them first ... )" <<endl; |
| 91 |
> |
cout << "Corrected JZB variable definition " << endl << " Data: " << jzbvariabledata << " and MC: " << jzbvariablemc << endl; |
| 92 |
> |
cout << endl; |
| 93 |
> |
cout << "If you're feeling lazy, copy & paste this : " << endl; |
| 94 |
> |
cout << " string jzbvariabledata=\""<<jzbvariabledata<<"\";"<<endl; |
| 95 |
> |
cout << " string jzbvariablemc=\"" <<jzbvariablemc<<"\";"<<endl; |
| 96 |
> |
} |
| 97 |
|
|
| 98 |
|
if(study_sidebands) look_at_sidebands(mcjzb.str(),datajzb.str()); |
| 99 |
|
|
