| 9 |
|
#include "Modules/GeneralToolBox.C" |
| 10 |
|
#include "Modules/SampleClass.C" |
| 11 |
|
#include "Modules/PeakFinder.C" |
| 12 |
– |
#include "Modules/HistoDB.C" |
| 12 |
|
#include "Modules/Poisson_Calculator.C" |
| 13 |
|
#include "Modules/setTDRStyle.C" |
| 14 |
|
#include "Modules/ActiveSamples.C" |
| 79 |
|
if(MCPeak>0) mcjzb<<"("<<jzbvariablemc<<"-"<<TMath::Abs(MCPeak)<<")"; |
| 80 |
|
else mcjzb<<"("<<jzbvariablemc<<"+"<<TMath::Abs(MCPeak)<<")"; |
| 81 |
|
|
| 82 |
< |
cout << "With peak correction, we get : " << endl; |
| 83 |
< |
cout << " Data : " << datajzb.str() << endl; |
| 84 |
< |
cout << " MC : " << datajzb.str() << endl; |
| 82 |
> |
dout << "With peak correction, we get : " << endl; |
| 83 |
> |
dout << " Data : " << datajzb.str() << endl; |
| 84 |
> |
dout << " MC : " << datajzb.str() << endl; |
| 85 |
|
|
| 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; |
| 90 |
> |
dout << "Please update Modules/Setup.C to reflect the following values (round them first ... )" <<endl; |
| 91 |
> |
dout << "Corrected JZB variable definition " << endl << " Data: " << jzbvariabledata << " and MC: " << jzbvariablemc << endl; |
| 92 |
> |
dout << endl; |
| 93 |
> |
dout << "If you're feeling lazy, copy & paste this : " << endl; |
| 94 |
> |
dout << " string jzbvariabledata=\""<<jzbvariabledata<<"\";"<<endl; |
| 95 |
> |
dout << " string jzbvariablemc=\"" <<jzbvariablemc<<"\";"<<endl; |
| 96 |
|
} |
| 97 |
|
|
| 98 |
|
if(study_sidebands) look_at_sidebands(mcjzb.str(),datajzb.str()); |
| 99 |
|
|
| 100 |
|
if(do_pick_up_events) { |
| 101 |
|
|
| 102 |
< |
cout << "Observed: " << endl; |
| 102 |
> |
dout << "Observed: " << endl; |
| 103 |
|
pick_up_events((const char*)(cutmass&&cutOSSF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)>100)")); |
| 104 |
< |
cout << "Predicted (JZB<-100) OSSF" << endl; |
| 104 |
> |
dout << "Predicted (JZB<-100) OSSF" << endl; |
| 105 |
|
pick_up_events((const char*)(cutmass&&cutOSSF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)<-100)")); |
| 106 |
< |
cout << "Predicted (emu, JZB>100) OSOF" << endl; |
| 106 |
> |
dout << "Predicted (emu, JZB>100) OSOF" << endl; |
| 107 |
|
pick_up_events((const char*)(cutmass&&cutOSOF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)>100)")); |
| 108 |
< |
cout << "Predicted (emu, JZB<-100) OSOF" << endl; |
| 108 |
> |
dout << "Predicted (emu, JZB<-100) OSOF" << endl; |
| 109 |
|
pick_up_events((const char*)(cutmass&&cutOSOF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)<-100)")); |
| 110 |
< |
cout << "Predicted (SB emu, JZB>100) OSOF" << endl; |
| 110 |
> |
dout << "Predicted (SB emu, JZB>100) OSOF" << endl; |
| 111 |
|
pick_up_events((const char*)(sidebandcut&&cutOSOF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)>100)")); |
| 112 |
< |
cout << "Predicted (SB emu, JZB<-100) OSOF" << endl; |
| 112 |
> |
dout << "Predicted (SB emu, JZB<-100) OSOF" << endl; |
| 113 |
|
pick_up_events((const char*)(cutOSOF&&cutnJets&&basiccut&&sidebandcut&&"((jzb[1]+0.06*pt-2.84727)<-100)")); |
| 114 |
< |
cout << "Predicted (SB SF, JZB>100) OSSF" << endl; |
| 114 |
> |
dout << "Predicted (SB SF, JZB>100) OSSF" << endl; |
| 115 |
|
pick_up_events((const char*)(sidebandcut&&cutOSSF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)>100)")); |
| 116 |
< |
cout << "Predicted (SB SF, JZB<-100) OSSF" << endl; |
| 116 |
> |
dout << "Predicted (SB SF, JZB<-100) OSSF" << endl; |
| 117 |
|
pick_up_events((const char*)(sidebandcut&&cutOSSF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)<-100)")); |
| 118 |
|
} |
| 119 |
|
|
