ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/cbrown/Development/Plotting/Selective_Plot_Generator.C
Revision: 1.15
Committed: Tue May 29 11:26:38 2012 UTC (12 years, 11 months ago) by buchmann
Content type: text/plain
Branch: MAIN
Changes since 1.14: +8 -17 lines
Log Message: 
added shorthand for btagging (-b)

File Contents

# User Rev Content
1 buchmann 1.1 #include <iostream>
2     #include <vector>
3     #include <sys/stat.h>
4     #include <getopt.h>
5     #include <stdio.h>
6     #include <stdlib.h>
7     #include "Modules/GeneralToolBox.C"
8     #include "Modules/SampleClass.C"
9     #include "Modules/setTDRStyle.C"
10     #include "Modules/Setup.C"
11     #include "Modules/Poisson_Calculator.C"
12     #include "Modules/ActiveSamples.C"
13     #include "Modules/PeakFinder.C"
14     #include "Modules/UpperLimitsWithShape.C"
15     #include "Modules/Plotting_Functions.C"
16     #include "Modules/LimitCalculation.C"
17     #include "Modules/ResultModule.C"
18 buchmann 1.4 #include "Modules/CrossSectionReader.C"
19 buchmann 1.1 #include "Modules/Systematics.C"
20 buchmann 1.10 #include "Modules/SugarCoating.C"
21     #include "Modules/ExclusionPlot.C"
22 buchmann 1.1 #include "Modules/SUSYScan.C"
23    
24 buchmann 1.10
25    
26 buchmann 1.1 #include <TCut.h>
27     #include <TROOT.h>
28     #include <TCanvas.h>
29     #include <TMath.h>
30     #include <TColor.h>
31     #include <TPaveText.h>
32     #include <TRandom.h>
33     #include <TH1.h>
34     #include <TH2.h>
35     #include <TF1.h>
36     #include <TSQLResult.h>
37    
38     #ifndef Verbosity
39     #define Verbosity 0
40     #endif
41     #ifndef HUSH
42     #define HUSH 1
43     #endif
44    
45     using namespace PlottingSetup;
46    
47     void usage(int passed=0 ) {
48     std::cout << "USAGE : " << std::endl;
49     std::cout << "You can use different options when running this program : " << std::endl;
50     std::cout << std::endl;
51     std::cout << "\033[1;34m all\033[0m \t\t All processes of the standard workflow" << std::endl;
52     std::cout << "\033[1;34m peak\033[0m \t\t Peak finding" << std::endl;
53     std::cout << "\033[1;34m kin\033[0m \t\t Kinematic plots" << std::endl;
54     std::cout << "\033[1;34m pfkin\033[0m \t Kinematic plots for PF variables" << std::endl;
55     std::cout << "\033[1;34m region\033[0m \t Region comparison plots" << std::endl;
56     std::cout << "\033[1;34m lepton\033[0m \t Lepton comparison plots" << std::endl;
57     std::cout << "\033[1;34m jzb\033[0m \t\t JZB plots" << std::endl;
58     std::cout << "\033[1;34m pred\033[0m \t\t Prediction plots" << std::endl;
59     std::cout << "\033[1;34m ratio\033[0m \t Ratio plots (standalone ratio plot)" << std::endl;
60     std::cout << "\033[1;34m signalbg\033[0m \t Signal vs. background plots (shape comparison)" << std::endl;
61     std::cout << "\033[1;34m ttbar\033[0m \t TTbar sideband comparison" << std::endl;
62     std::cout << "\033[1;34m zjets\033[0m \t ZJets plots" << std::endl;
63 buchmann 1.7 std::cout << "\033[1;34m diboson\033[0m \t Diboson plots" << std::endl;
64 buchmann 1.1 std::cout << "\033[1;34m results\033[0m \t Compute results (JZB>X)" << std::endl;
65     std::cout << "\033[1;34m syst\033[0m \t\t Compute systematics" << std::endl;
66     std::cout << "\033[1;34m eff\033[0m \t\t Compute all (partial) efficiencies" << std::endl;
67     std::cout << "\033[1;34m shapes\033[0m \t Prepare ROOT files for computing limits using shapes" << std::endl;
68     std::cout << "\033[1;34m counting\033[0m \t Establish limits with counting experiment" << std::endl;
69 buchmann 1.14 std::cout << "\033[1;34m shapelimit\033[0m \t Establish limits with shapes" << std::endl;
70 buchmann 1.15 std::cout << "\033[1;34m qcd\033[0m \t\t Show that QCD is negligible" << endl;
71 buchmann 1.1 std::cout << "\033[1;34m response\033[0m \t Compute response correction" << std::endl;
72     std::cout << "\033[1;34m pickup\033[0m \t Pick up interesting events" << std::endl;
73     std::cout << "\033[1;34m yields\033[0m \t Compute MC yields for this lumi" << std::endl;
74     std::cout << "\033[1;34m sideassess\033[0m\t Assessment of the mass sidebands definition in the data" << std::endl;
75     std::cout << "\033[1;34m jes\033[0m \t\t Make Jet Energy Scale (JES) plot" << std::endl;
76 buchmann 1.15 std::cout << "\033[1;34m effcurve\033[0m \t\t Make JZB efficiency plot" << std::endl;
77 buchmann 1.1 std::cout << "\033[1;34m save\033[0m \t\t Save results in a template file (for reuse when establishing limits for scans)" << std::endl;
78     std::cout << "\033[1;34m png\033[0m \t\t Save all plots as pngs (also available: --pdf, --root, --eps, --C) standard: png&pdf" << std::endl;
79     std::cout << "\033[1;34m offpeak\033[0m \t Sets offpeak to true (i.e. carry out off-peak analysis!)" << std::endl;
80     std::cout << "\033[1;34m onpeak\033[0m \t Sets onpeak to true (i.e. carry out standard on-peak analysis!)" << std::endl;
81 buchmann 1.2 std::cout << "\033[1;34m jzbcuts\033[0m \t Sets the jzb cuts (please provide them in the format 50,100,150,200 and so on (comma separated)" << endl;
82     std::cout << "\033[1;34m ptsanity\033[0m \t Checks whether all samples have the correct pt range (additional canvas for you to check)" << endl;
83 buchmann 1.15 std::cout << "\033[1;34m mlls\033[0m \t Draw dilepton invariant mass for all signals" << endl;
84 buchmann 1.1 std::cout << std::endl;
85     std::cout << "\033[1;34m dir (d)\033[0m\t Directory where all plots will be saved" << std::endl;
86     std::cout << "\033[1;34m public \033[0m\t Switch to be activated when using the algorithm in public (no branding etc.)" << std::endl;
87     std::cout << std::endl;
88     std::cout << std::endl;
89     std::cout << "\033[1;34m met\033[0m \t\t Compute obs/pred for met cut vs jzb cut (in data, lm4 and lm8)" << std::endl;
90 buchmann 1.15 std::cout << "\033[1;34m metvsjzb\033[0m \t\t Produce MET:JZB plots" << std::endl;
91 buchmann 1.1 std::cout << std::endl;
92     std::cout << std::endl;
93     std::cout << "\033[1;34m paper\033[0m \t\t Switch paper mode on (no \" Preliminary \" on plots)" << std::endl;
94 buchmann 1.15 std::cout << "\033[1;34m approved\033[0m \t\t Only approved plots" << std::endl;
95 buchmann 1.1 exit(-1);
96     }
97    
98     int main (int argc, char ** argv)
99     {
100     int do_all = false;/// DONE
101     int do_peak_finding = false; /// DONE
102     int do_kinematic_variables = false; ///DONE
103     int do_kinematic_PF_variables = false; ///DONE
104     int do_region_comparison = false; /// Comparison between signal and control regions
105     int do_lepton_comparison = false; ///DONE
106     int do_jzb_plots = false; /// DONE
107     int do_pred = false; /// DONE
108     bool overlay_signal = false; /// Overlay LM4 signal to predictions
109     int do_ratio = false; /// DONE ------------------------------------------- still need to adapt the legend!
110     int do_signal_bg_comparison_plot = false; /// DONE
111     int do_ttbar_comparison = false; /// DONE
112     int do_jzb_efficiency_curve = false; /// DONE
113     int do_zjets_comparison = false; /// DONE
114 buchmann 1.7 int do_diboson_plots = false;
115 buchmann 1.1 // int calculate_pred_and_observed = false; /// DONE -- now rerouted to the new way of computing results
116     // int calculate_yields = false; ///DONE /// superfluous (replaced by get_new_results)
117     int get_new_results = false; ///DONE
118     bool dopoisson = false;//should we calculate stat err with poisson? ///DONE
119     bool verbose = false;//prints out a lot more information ... ///DONE
120     //--------------------------------------------
121     // Systematics and limits
122     int do_compute_systematics=false; /// DONE (might want to add something related to the peak?)
123     bool requireZ=true; // should we require a Z for MC efficiency?
124     int do_compute_efficiency=false; /// DONE (might want to add something related to the peak?)
125     int do_prepare_limits_using_shapes=false; /// ALMOST DONE **************************************
126     int do_compute_upper_limits_from_counting_experiment=false; /// DONE
127     bool doobserved=false; // want to get the observed limit as well?
128     int doquick=1; // this will cause the number of predicted/observed events to hurry up (0 = do mc and all results, 1 = all results but no mc, 2 = don't save results in library and don't compute mc stuff)
129 buchmann 1.14 int do_compute_upper_limits_from_shapes=false;
130 buchmann 1.1
131     //--------------------------------------------
132     // More not-so-standard stuff
133     int do_model_scan=false; // this can take days - use the grid version!
134     int do_zjet_ttbar_shapes=false; /// DONE
135     int do_test=false;/// DONE (just any test you wish)
136     int do_pick_up_events=false; ///DONE
137    
138     //--------------------------------------------
139     // To be run manually !
140     int do_response_correction=false; ///DONE
141     //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
142    
143     int do_sideband_assessment=false;
144 buchmann 1.2 int do_qcd=false;
145 buchmann 1.1
146     int do_save_template=false;
147     int do_mc_yields=false; // compute MC yields
148     int do_metjzb=false;
149    
150     int savepdf=true;
151     int saveC=true;
152     int saveRoot=true;
153     int savepng=true;
154     int saveeps=false;
155    
156 buchmann 1.3 int dooffpeak=true;
157 buchmann 1.1 int doonpeak=false;
158    
159     int do_JES=false;
160 buchmann 1.2 int do_ptsanity=false;
161 buchmann 1.3 int do_mlls=false;
162 buchmann 1.5 int do_met_vs_jzb=false;
163 buchmann 1.1
164     std::string directory="";
165     int option_iterator;
166     int option_counter=0;
167     bool moreoptions=true;
168    
169     string jzbcuts_string="";
170    
171     while(moreoptions) {
172     static struct option long_options[] =
173     {
174     /* These options set a flag. */
175     {"all", no_argument, &do_all, 1},
176     {"peak", no_argument, &do_peak_finding, 1},
177     {"kin", no_argument, &do_kinematic_variables, 1},
178     {"pfkin", no_argument, &do_kinematic_PF_variables, 1},
179     {"region", no_argument, &do_region_comparison, 1},
180     {"lepton", no_argument, &do_lepton_comparison, 1},
181     {"jzb", no_argument, &do_jzb_plots, 1},
182     {"pred", no_argument, &do_pred, 1},
183     {"ratio", no_argument, &do_ratio, 1},
184     {"signalbg", no_argument, &do_signal_bg_comparison_plot, 1},
185     {"ttbar", no_argument, &do_ttbar_comparison, 1},
186     {"zjets", no_argument, &do_zjets_comparison, 1},
187 buchmann 1.7 {"diboson", no_argument, &do_diboson_plots, 1},
188 buchmann 1.1 {"results", no_argument, &get_new_results, 1},
189     {"syst", no_argument, &do_compute_systematics, 1},
190     {"eff", no_argument, &do_compute_efficiency, 1},
191     {"shapes", no_argument, &do_prepare_limits_using_shapes, 1},
192     {"counting", no_argument, &do_compute_upper_limits_from_counting_experiment, 1},
193 buchmann 1.14 {"shapelimit",no_argument, &do_compute_upper_limits_from_shapes, 1},
194 buchmann 1.2 {"ptsanity", no_argument, &do_ptsanity, 1},
195 buchmann 1.1 {"response", no_argument, &do_response_correction, 1},
196     {"pickup", no_argument, &do_pick_up_events, 1},
197     {"yields", no_argument, &do_mc_yields, 1},
198     {"sideassess",no_argument, &do_sideband_assessment, 1},
199 buchmann 1.5 {"metvsjzb", no_argument, &do_met_vs_jzb, 1},
200 buchmann 1.2 {"qcd", no_argument, &do_qcd, 1},
201 buchmann 1.1 {"save", no_argument, &do_save_template,1},
202     {"png", no_argument, &savepng,1},
203     {"eps", no_argument, &saveeps,1},
204     {"pdf", no_argument, &savepdf,1},
205     {"root", no_argument, &saveRoot,1},
206     {"C", no_argument, &saveC,1},
207     {"offpeak", no_argument, &dooffpeak,1},
208     {"onpeak", no_argument, &doonpeak,1},
209     {"met", no_argument, &do_metjzb,1},
210     {"jes", no_argument, &do_JES,1},
211     {"effcurve", no_argument, &do_jzb_efficiency_curve,1},
212     {"public", no_argument, &PlottingSetup::publicmode,1},
213     {"paper", no_argument, &PlottingSetup::PaperMode,1},
214 buchmann 1.9 {"approved", no_argument, &PlottingSetup::Approved,1},
215 buchmann 1.3 {"mlls", no_argument, &do_mlls,1},
216 buchmann 1.1 /* The following options store values.*/
217     {"jzbcuts", required_argument, 0, 'j'},
218     {"dir", required_argument, 0, 'd'},
219     {0, 0, 0, 0}
220     };
221     int option_index = 0;
222 buchmann 1.15 option_iterator = getopt_long(argc, argv, "ad:j:",long_options, &option_index);
223 buchmann 1.1 if(option_iterator == -1) moreoptions=false;
224     else {
225     option_counter++;
226     switch (option_iterator)
227     {
228     case 0:
229     if (long_options[option_index].flag != 0)
230     break;
231     printf ("option %s", long_options[option_index].name);
232     if (optarg)
233     printf (" with arg %s", optarg);
234     printf ("\n");
235     break;
236     case 'a':
237     do_all=true;
238     break;
239     case 'd':
240     directory=(std::string)optarg;
241     std::cout<<"Option directory was passed with argument " << optarg << std::endl;
242     break;
243     case 'j':
244     jzbcuts_string=(std::string)optarg;
245     std::cout<<"JZB cuts were manually defined:" << optarg << std::endl;
246     break;
247     case '?':
248     usage(option_iterator);
249     break;
250     default:
251     usage(option_iterator);
252     }
253     }
254     }
255    
256     if(directory!="") PlottingSetup::directoryname=directory;
257     if(dooffpeak) PlottingSetup::RestrictToMassPeak=false;
258     if(doonpeak) PlottingSetup::RestrictToMassPeak=true;
259     if(option_counter==0) usage();
260    
261     ///----------------------------------- BELOW THIS LINE: NO MORE OPTIONS BUT ACTUAL FUNCTION CALLS! ---------------------------------------------------------
262     gROOT->SetStyle("Plain");
263     bool do_fat_line=false; // if you want to have HistLineWidth=1 and FuncWidth=1 as it was before instead of 2
264     setTDRStyle(do_fat_line);
265     gStyle->SetTextFont(42);
266     bool showList=true;
267     set_directory(directoryname);//Indicate the directory name where you'd like to save the output files in Setup.C
268     set_treename("events");//you can set the treename here to be used; options are "events" (for reco) for "PFevents" (for particle flow)
269     // define_samples(showList,allsamples);
270 buchmann 1.2 define_samples(showList,allsamples,signalsamples,scansample,raresample,systsamples,qcdsamples);
271 buchmann 1.1 setlumi(luminosity);
272     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!
273     stringstream resultsummary;
274    
275 buchmann 1.11 write_analysis_type(PlottingSetup::RestrictToMassPeak,PlottingSetup::DoBTag);
276 buchmann 1.1 do_png(savepng);
277     do_pdf(savepdf);
278     do_eps(saveeps);
279     do_C(saveC);
280     do_root(saveRoot);
281    
282    
283     global_ratio_binning.push_back(0);
284     global_ratio_binning.push_back(5);
285     global_ratio_binning.push_back(10);
286     global_ratio_binning.push_back(20);
287     global_ratio_binning.push_back(50);
288     global_ratio_binning.push_back(100);
289     global_ratio_binning.push_back(200);
290     global_ratio_binning.push_back(400);
291     //global_ratio_binning.push_back(500);
292    
293     //these are the JZB cuts defining our search regions
294     vector<float> jzb_cut; //starting where, please?
295     if(jzbcuts_string=="") {
296     jzb_cut.push_back(50);
297     // jzb_cut.push_back(75);
298     jzb_cut.push_back(100);
299     // jzb_cut.push_back(125);
300     jzb_cut.push_back(150);
301     // jzb_cut.push_back(175);
302     jzb_cut.push_back(200);
303     // jzb_cut.push_back(225);
304     jzb_cut.push_back(250);
305     // jzb_cut.push_back(275);
306     // jzb_cut.push_back(300);
307     } else {
308     manually_set_jzb_cuts(jzb_cut,jzbcuts_string);
309     if(jzb_cut.size()==0) {
310     write_error(__FUNCTION__,"There are no jzb cuts!");
311     return -1;
312     }
313     }
314    
315     todo();
316 buchmann 1.2
317 buchmann 1.1 //**** part 1 : peak finding
318     float MCPeak=0,MCPeakError=0,DataPeak=0,DataPeakError=0,MCSigma=10,DataSigma=10;
319     method=Kostasmethod;//Kostasmethod;//dogaus3sigma;// options: dogaus,doKM,dogaus2sigma,dogaus3sigma
320 buchmann 1.8 //write_warning(__FUNCTION__,"NOT DOING ANY PEAK FINDING ATM");
321 buchmann 1.14 if(do_peak_finding||PlottingSetup::Approved||do_mlls||do_zjets_comparison||do_sideband_assessment||do_metjzb||do_zjet_ttbar_shapes||do_diboson_plots||do_region_comparison||do_jzb_plots||do_lepton_comparison||do_signal_bg_comparison_plot||do_pred||do_ratio||do_mc_yields||do_save_template||do_compute_upper_limits_from_counting_experiment||do_compute_upper_limits_from_shapes||do_compute_systematics||do_model_scan||do_prepare_limits_using_shapes||do_all||get_new_results||do_ttbar_comparison||do_kinematic_variables||do_kinematic_PF_variables) find_peaks(MCPeak,MCPeakError, DataPeak, DataPeakError,MCSigma,DataSigma,resultsummary);
322 buchmann 1.1
323     stringstream datajzb;
324     if(DataPeak>0) datajzb<<"("<<jzbvariabledata<<"-"<<TMath::Abs(DataPeak)<<")";
325     else datajzb<<"("<<jzbvariabledata<<"+"<<TMath::Abs(DataPeak)<<")";
326     stringstream mcjzb;
327     if(MCPeak>0) mcjzb<<"("<<jzbvariablemc<<"-"<<TMath::Abs(MCPeak)<<")";
328     else mcjzb<<"("<<jzbvariablemc<<"+"<<TMath::Abs(MCPeak)<<")";
329    
330     dout << "With peak correction, we get : " << endl;
331     dout << " Data : " << datajzb.str() << endl;
332     dout << " MC : " << mcjzb.str() << endl;
333    
334    
335    
336     //**** part 2 : kinematic plots
337 buchmann 1.9 if(do_kinematic_variables||(do_all&&!PlottingSetup::Approved)) do_kinematic_plots(mcjzb.str(),datajzb.str());
338 buchmann 1.1 if(do_kinematic_PF_variables) do_kinematic_PF_plots(mcjzb.str(),datajzb.str());
339    
340     //**** part 3: Leptonic comparison (ee vs mm, eemm vs emu)
341 buchmann 1.9 if(do_lepton_comparison||(do_all&&!PlottingSetup::Approved)) lepton_comparison_plots();
342 buchmann 1.1
343    
344     //**** part 3b: comparison between control regions (SFZPJZBPOS,SFZPJZBNEG,...)
345 buchmann 1.9 if (do_region_comparison||(do_all&&!PlottingSetup::Approved)) region_comparison_plots(mcjzb.str(),datajzb.str(),jzb_cut);
346 buchmann 1.1
347     //**** part 4: JZB plots (OSOF, OSSF)
348     if(do_jzb_plots||do_all) jzb_plots(mcjzb.str(),datajzb.str(),global_ratio_binning);
349    
350     //**** part 5 : Prediction plots
351     if(do_pred||do_all) do_prediction_plots(mcjzb.str(),datajzb.str(),DataSigma,MCSigma,overlay_signal);
352    
353     //**** part 6: Ratio plots
354 buchmann 1.9 if(do_ratio||(do_all&&!PlottingSetup::Approved)) do_ratio_plots(mcjzb.str(),datajzb.str(),global_ratio_binning);
355 buchmann 1.1
356     //**** part 7: Some systematics plots
357 buchmann 1.9 if(do_signal_bg_comparison_plot||(do_all&&!PlottingSetup::Approved)) signal_bg_comparison();
358     if(do_ttbar_comparison||(do_all&&!PlottingSetup::Approved)) ttbar_sidebands_comparison(mcjzb.str(),global_ratio_binning);
359 buchmann 1.13 if(do_zjets_comparison||(do_all&&!PlottingSetup::Approved)) zjets_prediction_comparison(mcjzb.str());
360 buchmann 1.7
361     if(do_diboson_plots) diboson_plots(mcjzb.str(),datajzb.str(),global_ratio_binning);
362 buchmann 1.1
363 buchmann 1.9 if(do_JES||(do_all&&!PlottingSetup::Approved)) make_JES_plot();
364 buchmann 1.1
365 buchmann 1.9 if (do_jzb_efficiency_curve) plot_jzb_sel_eff(mcjzb.str(),signalsamples,jzb_cut);
366 buchmann 1.1
367     //**** part 8: observed and predicted!
368     if(do_all||get_new_results) doquick=0;
369     if(do_save_template) doquick=2;
370     if(do_all||do_save_template||do_compute_upper_limits_from_counting_experiment||do_model_scan||get_new_results) get_result(mcjzb.str(),datajzb.str(),DataPeakError,MCPeakError,jzb_cut,verbose,dopoisson,doquick);
371    
372     vector<float>jzb_shape_limit_bins;
373 buchmann 1.14 jzb_shape_limit_bins.push_back(50);
374     jzb_shape_limit_bins.push_back(65.05);
375     jzb_shape_limit_bins.push_back(73.4);
376     jzb_shape_limit_bins.push_back(88.7);
377     jzb_shape_limit_bins.push_back(115.25);
378 buchmann 1.12 jzb_shape_limit_bins.push_back(200);
379     // jzb_shape_limit_bins.push_back(250);
380 buchmann 1.1 if(do_prepare_limits_using_shapes) prepare_limits(mcjzb.str(),datajzb.str(),DataPeakError,MCPeakError,jzb_shape_limit_bins);
381    
382     //------------------------ end of standard functions
383    
384     // if(calculate_yields||do_all) calculate_all_yields(mcjzb.str(),jzb_cut); // now outdated - get_results now replaces this.
385    
386    
387     vector<vector<float> > all_systematics;
388     vector<vector<float> > all_efficiency;
389    
390     if(do_compute_systematics||do_compute_upper_limits_from_counting_experiment||do_all) all_systematics=compute_systematics(mcjzb.str(),MCPeakError,alwaysflip,datajzb.str(),signalsamples,jzb_cut,requireZ);
391     // if(do_compute_efficiency) compute_efficiency(mcjzb.str(),allsamples, 50, -200, 300);
392     if(do_compute_efficiency) write_warning(__FUNCTION__,"efficiency computation deactivated");
393    
394     if(do_compute_upper_limits_from_counting_experiment||do_all) compute_upper_limits_from_counting_experiment(all_systematics,jzb_cut,mcjzb.str(),doobserved,alwaysflip);
395    
396     //------------------------ end of analysis parts - below: special applications
397    
398 buchmann 1.5 if(do_model_scan) scan_SUSY_parameter_space(mcjzb.str(),datajzb.str(),jzb_cut,requireZ,MCPeakError,DataPeakError);
399 buchmann 1.1
400     if(do_zjet_ttbar_shapes) draw_ttbar_and_zjets_shape(mcjzb.str(),datajzb.str());
401    
402    
403     //now do predicted and observed.
404     //to that end, we need to consider systematics!
405     // systematics come from:
406     // Poisson error on pred. contributions (easy)
407     // +20%/-40% from MC closure test
408     // from peak position --> whenever we want to investigate systematics, peak position must have been set!
409    
410     //jzb response correction
411     if(do_response_correction) {
412     find_correction_factors(jzbvariabledata,jzbvariablemc);
413     dout << "Please update Modules/Setup.C to reflect the following values (round them first ... )" <<endl;
414     dout << "Corrected JZB variable definition " << endl << " Data: " << jzbvariabledata << " and MC: " << jzbvariablemc << endl;
415     dout << endl;
416     dout << "If you're feeling lazy, copy & paste this : " << endl;
417     dout << " string jzbvariabledata=\""<<jzbvariabledata<<"\";"<<endl;
418     dout << " string jzbvariablemc=\"" <<jzbvariablemc<<"\";"<<endl;
419     }
420    
421     if(do_pick_up_events) {
422    
423     dout << "Observed: " << endl;
424     pick_up_events((const char*)(cutmass&&cutOSSF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)>100)"));
425     dout << "Predicted (JZB<-100) OSSF" << endl;
426     pick_up_events((const char*)(cutmass&&cutOSSF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)<-100)"));
427     dout << "Predicted (emu, JZB>100) OSOF" << endl;
428     pick_up_events((const char*)(cutmass&&cutOSOF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)>100)"));
429     dout << "Predicted (emu, JZB<-100) OSOF" << endl;
430     pick_up_events((const char*)(cutmass&&cutOSOF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)<-100)"));
431     dout << "Predicted (SB SF, JZB>100) OSSF" << endl;
432     pick_up_events((const char*)(sidebandcut&&cutOSSF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)>100)"));
433     dout << "Predicted (SB SF, JZB<-100) OSSF" << endl;
434     pick_up_events((const char*)(sidebandcut&&cutOSSF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)<-100)"));
435     dout << "Predicted (SB emu, JZB>100) OSOF" << endl;
436     pick_up_events((const char*)(sidebandcut&&cutOSOF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)>100)"));
437     dout << "Predicted (SB emu, JZB<-100) OSOF" << endl;
438     pick_up_events((const char*)(cutOSOF&&cutnJets&&basiccut&&sidebandcut&&"((jzb[1]+0.06*pt-2.84727)<-100)"));
439     }
440    
441 buchmann 1.14 if(do_save_template||do_compute_upper_limits_from_shapes) {
442     write_error(__FUNCTION__,"Currently commented out shape storage");
443 buchmann 1.15 /* save_template(mcjzb.str(),datajzb.str(),jzb_cut,MCPeakError,DataPeakError,jzb_shape_limit_bins);
444     PrepareDataShapes(mcjzb.str(), datajzb.str(), jzb_shape_limit_bins, DataPeakError);*/
445 buchmann 1.6 }
446 buchmann 1.14
447     if(do_compute_upper_limits_from_shapes) compute_upper_limits_from_shapes(mcjzb.str(),datajzb.str(),jzb_cut,requireZ, MCPeakError, DataPeakError);
448 buchmann 1.1
449     if(do_sideband_assessment) sideband_assessment(datajzb.str(),30.0,50.0);
450    
451     if(do_mc_yields) compute_MC_yields(mcjzb.str(),jzb_cut);
452    
453     if(do_metjzb) met_jzb_cut(datajzb.str(),mcjzb.str(),jzb_cut);
454    
455 buchmann 1.2 if(do_qcd||do_all) qcd_plots(datajzb.str(),mcjzb.str(),jzb_cut);
456    
457 buchmann 1.3 if(do_mlls) do_mlls_plot(mcjzb.str());
458    
459 buchmann 1.2 if(do_ptsanity) check_ptsanity();
460 buchmann 1.3
461 buchmann 1.5 if(do_met_vs_jzb) met_vs_jzb_plots();
462 buchmann 1.1 if(do_test) test();
463    
464     return 0;
465     }
466    
467