4 |
|
#include <getopt.h> |
5 |
|
#include <stdio.h> |
6 |
|
#include <stdlib.h> |
7 |
+ |
|
8 |
|
#include "Modules/GeneralToolBox.C" |
9 |
|
#include "Modules/SampleClass.C" |
10 |
|
#include "Modules/setTDRStyle.C" |
11 |
< |
#include "Modules/Setup.C" |
11 |
< |
//#include "Modules/Poisson_Calculator.C" |
11 |
> |
#include "Modules/ZbSetup.C" |
12 |
|
#include "Modules/JSON/JSONSampleLoader.C" |
13 |
|
#include "Modules/ActiveSamples.C" |
14 |
– |
//#include "Modules/PeakFinder.C" |
15 |
– |
//#include "Modules/UpperLimitsWithShape.C" |
16 |
– |
//#include "Modules/Plotting_Functions.C" |
17 |
– |
//#include "Modules/LimitCalculation.C" |
18 |
– |
//#include "Modules/ResultModule.C" |
19 |
– |
//#include "Modules/CrossSectionReader.C" |
20 |
– |
//#include "Modules/Systematics.C" |
21 |
– |
//#include "Modules/SugarCoating.C" |
22 |
– |
//#include "Modules/ExclusionPlot.C" |
23 |
– |
//#include "Modules/SUSYScan.C" |
24 |
– |
|
14 |
|
#include "Modules/ValueClass.C" |
15 |
|
#include "Modules/ZbTools.C" |
16 |
|
|
133 |
|
PlottingSetup::RestrictToMassPeak=true; |
134 |
|
|
135 |
|
//define_Zb_samples(showList,allsamples,signalsamples,scansample,raresample,systsamples,qcdsamples); |
136 |
+ |
IsZbAnalysis=true; |
137 |
+ |
IsWZAnalysis=false; |
138 |
+ |
IsJZBAnalysis=false; |
139 |
|
write_warning(__FUNCTION__,"Need to adapt sample collection!"); |
140 |
|
define_samples(showList,allsamples,signalsamples,scansample,raresample,systsamples,qcdsamples,comparesamples); |
141 |
|
setlumi(luminosity); |
149 |
|
do_root(saveRoot); |
150 |
|
|
151 |
|
|
152 |
< |
|
153 |
< |
|
154 |
< |
IsZbAnalysis=true; |
163 |
< |
IsWZAnalysis=false; |
164 |
< |
IsJZBAnalysis=false; |
165 |
< |
|
166 |
< |
|
167 |
< |
|
168 |
< |
write_warning(__FUNCTION__,"Everything missing so far goes here: "); |
169 |
< |
cout << " Still need to carry out cross - checks" << endl; |
170 |
< |
cout << " Still need to define errors in final plot " << endl; |
171 |
< |
|
172 |
< |
|
152 |
> |
write_warning(__FUNCTION__,"Using preliminary weights!"); |
153 |
> |
//cutWeight=TCut("(weight*(weight<1000)*(is_data+(!is_data)*((id1==id2)*0.95+(id1!=id2)*0.94)))"); |
154 |
> |
cutWeight=TCut("(weight*(weight<1000)*(is_data+(!is_data)*((id1==id2&&id1==0)*0.95+(id1==id2&&id1==1)*0.88+(id1!=id2)*0.92)))"); |
155 |
|
do_basic_ZB_analysis(do_Full); |
156 |
|
return 0; |
157 |
|
} |