ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/cbrown/Development/DistributedModelCalculations/Systematics/SystematicsWorkerScript.C
Revision: 1.2
Committed: Wed Apr 11 08:55:56 2012 UTC (13 years ago) by buchmann
Content type: text/plain
Branch: MAIN
Changes since 1.1: +2 -1 lines
Log Message:
reading previous config a bit earlier for complete setup

File Contents

# Content
1 #include <iostream>
2 #include <vector>
3 #include <sys/stat.h>
4 #include "../../Plotting/Modules/GeneralToolBox.C"
5 #include "../../Plotting/Modules/SampleClass.C"
6 #include "../../Plotting/Modules/Setup.C"
7 #include "../../Plotting/Modules/PeakFinder.C"
8 #include "../../Plotting/Modules/Poisson_Calculator.C"
9 #include "../../Plotting/Modules/setTDRStyle.C"
10 #include "../../Plotting/Modules/ActiveSamples.C"
11 #include "../../Plotting/Modules/UpperLimitsWithShape.C"
12 #include "../../Plotting/Modules/Plotting_Functions.C"
13 #include "../../Plotting/Modules/LimitCalculation.C"
14 #include "../../Plotting/Modules/ResultModule.C"
15 #include "../../Plotting/Modules/CrossSectionReader.C"
16 #include "../../Plotting/Modules/Systematics.C"
17 #include "../../Plotting/Modules/SUSYScan.C"
18 #include "../last_configuration.C"
19
20 #include <TCut.h>
21 #include <TROOT.h>
22 #include <TCanvas.h>
23 #include <TMath.h>
24 #include <TColor.h>
25 #include <TPaveText.h>
26 #include <TRandom.h>
27 #include <TH1.h>
28 #include <TH2.h>
29 #include <TF1.h>
30 #include <TSQLResult.h>
31
32 #ifndef Verbosity
33 #define Verbosity 0
34 #endif
35 #ifndef HUSH
36 #define HUSH 1
37 #endif
38
39 using namespace PlottingSetup;
40
41 int main(int narg, char *args[])
42 {
43 if(narg<3||narg>3) {
44 cout << "USAGE : " << endl;
45 cout << args[0] << " Njobs thisjob" << endl;
46 cout << " where Njobs is the number of jobs, and thisjob correspond to this jobs' job number" << endl;
47 return -1;
48 }
49
50 float njobs=atoi(args[1]);
51 float jobnumber=atoi(args[2]);
52
53 cout << "Number of jobs: " << njobs << endl;
54 cout << "This is job number : " << jobnumber << endl;
55
56 read_config();
57
58 gROOT->SetStyle("Plain");
59 setTDRStyle();
60 gStyle->SetTextFont(42);
61 bool showList=false;
62 set_directory(PlottingSetup::directoryname);
63 set_treename("events");//you can set the treename here to be used; options are "events" (for reco) for "PFevents" (for particle flow)
64 define_SMS_sample(showList,allsamples,signalsamples,scansample);
65 setlumi(luminosity);
66 do_png(true);
67 do_pdf(false);
68 do_eps(false);
69 do_C(false);
70 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!
71
72
73 bool doquick=true;
74 bool systonly=true;
75 scan_SUSY_parameter_space(mcjzb,datajzb,jzb_cuts,PlottingSetup::requireZ,MCPeakError,DataPeakError, njobs, jobnumber,systonly);
76
77 return 0;
78 }
79