1 |
buchmann |
1.1 |
#include <iostream>
|
2 |
|
|
#include <vector>
|
3 |
|
|
#ifndef SampleClassLoaded
|
4 |
|
|
#include "SampleClass.C"
|
5 |
|
|
#endif
|
6 |
|
|
#define SetupLoaded
|
7 |
|
|
#ifndef ResultLibraryClassLoaded
|
8 |
|
|
#include "ResultLibraryClass.C"
|
9 |
|
|
#endif
|
10 |
|
|
|
11 |
|
|
using namespace std;
|
12 |
|
|
|
13 |
|
|
namespace PlottingSetup {
|
14 |
|
|
string datajzb="datajzb_ERROR";
|
15 |
|
|
string mcjzb="mcjzb_ERROR";
|
16 |
|
|
vector<float>jzb_cuts;
|
17 |
|
|
float MCPeakError=-999;
|
18 |
buchmann |
1.2 |
float DataPeakError=-999;
|
19 |
buchmann |
1.1 |
}
|
20 |
|
|
|
21 |
|
|
void read_config() {
|
22 |
buchmann |
1.2 |
datajzb="((jzb[1]+0.053*pt)-2.95688)";
|
23 |
|
|
mcjzb="((jzb[1]+0.043*pt)-3.05937)";
|
24 |
buchmann |
1.1 |
|
25 |
|
|
|
26 |
|
|
|
27 |
buchmann |
1.2 |
MCPeakError=0.159709;
|
28 |
|
|
DataPeakError=0.105771;
|
29 |
buchmann |
1.1 |
|
30 |
|
|
jzb_cuts.push_back(50); // JZB cut at 50
|
31 |
|
|
jzb_cuts.push_back(100); // JZB cut at 100
|
32 |
|
|
jzb_cuts.push_back(150); // JZB cut at 150
|
33 |
|
|
jzb_cuts.push_back(200); // JZB cut at 200
|
34 |
|
|
jzb_cuts.push_back(250); // JZB cut at 250
|
35 |
|
|
|
36 |
|
|
|
37 |
buchmann |
1.2 |
Nobs.push_back(1056); // JZB cut at 50
|
38 |
|
|
Nobs.push_back(272); // JZB cut at 100
|
39 |
|
|
Nobs.push_back(57); // JZB cut at 150
|
40 |
|
|
Nobs.push_back(17); // JZB cut at 200
|
41 |
|
|
Nobs.push_back(6); // JZB cut at 250
|
42 |
|
|
Npred.push_back(1157); // JZB cut at 50
|
43 |
|
|
Npred.push_back(283); // JZB cut at 100
|
44 |
|
|
Npred.push_back(86); // JZB cut at 150
|
45 |
|
|
Npred.push_back(33); // JZB cut at 200
|
46 |
|
|
Npred.push_back(11); // JZB cut at 250
|
47 |
|
|
Nprederr.push_back(134.014); // JZB cut at 50
|
48 |
|
|
Nprederr.push_back(32.8491); // JZB cut at 100
|
49 |
|
|
Nprederr.push_back(12.6686); // JZB cut at 150
|
50 |
|
|
Nprederr.push_back(6.58046); // JZB cut at 200
|
51 |
|
|
Nprederr.push_back(3.49428); // JZB cut at 250
|
52 |
|
|
|
53 |
|
|
|
54 |
|
|
flippedNobs.push_back(306); // JZB cut at 50
|
55 |
|
|
flippedNobs.push_back(21); // JZB cut at 100
|
56 |
|
|
flippedNobs.push_back(5); // JZB cut at 150
|
57 |
buchmann |
1.1 |
flippedNobs.push_back(1); // JZB cut at 200
|
58 |
|
|
flippedNobs.push_back(0); // JZB cut at 250
|
59 |
buchmann |
1.2 |
flippedNpred.push_back(205); // JZB cut at 50
|
60 |
|
|
flippedNpred.push_back(10); // JZB cut at 100
|
61 |
|
|
flippedNpred.push_back(-24); // JZB cut at 150
|
62 |
|
|
flippedNpred.push_back(-15); // JZB cut at 200
|
63 |
|
|
flippedNpred.push_back(-5); // JZB cut at 250
|
64 |
|
|
flippedNprederr.push_back(287.321); // JZB cut at 50
|
65 |
|
|
flippedNprederr.push_back(76.9838); // JZB cut at 100
|
66 |
|
|
flippedNprederr.push_back(20.3468); // JZB cut at 150
|
67 |
|
|
flippedNprederr.push_back(8.79218); // JZB cut at 200
|
68 |
|
|
flippedNprederr.push_back(4.52327); // JZB cut at 250
|
69 |
buchmann |
1.1 |
|
70 |
|
|
|
71 |
buchmann |
1.2 |
luminosity=4980;
|
72 |
|
|
RestrictToMassPeak=0;//defines the type of analysis we're running
|
73 |
buchmann |
1.1 |
|
74 |
|
|
|
75 |
|
|
cout << "Configuration successfully loaded!" << endl;
|
76 |
|
|
|
77 |
|
|
}
|
78 |
|
|
|