1 |
buchmann |
1.1 |
#include <iostream>
|
2 |
|
|
#include <vector>
|
3 |
buchmann |
1.11 |
#ifndef SampleClassLoaded
|
4 |
|
|
#include "SampleClass.C"
|
5 |
|
|
#endif
|
6 |
|
|
#define SetupLoaded
|
7 |
|
|
#ifndef ResultLibraryClassLoaded
|
8 |
|
|
#include "ResultLibraryClass.C"
|
9 |
|
|
#endif
|
10 |
buchmann |
1.1 |
|
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 |
|
|
}
|
19 |
|
|
|
20 |
|
|
void read_config() {
|
21 |
buchmann |
1.11 |
datajzb="(jzb[1]+0.06*pt-3.53058)";
|
22 |
buchmann |
1.10 |
mcjzb="(jzb[1]+0.04*pt-2.4392)";
|
23 |
buchmann |
1.1 |
|
24 |
buchmann |
1.2 |
|
25 |
|
|
|
26 |
buchmann |
1.11 |
MCPeakError=0.097537;
|
27 |
buchmann |
1.1 |
|
28 |
buchmann |
1.10 |
jzb_cuts.push_back(50); // JZB cut at 50
|
29 |
|
|
jzb_cuts.push_back(100); // JZB cut at 100
|
30 |
|
|
jzb_cuts.push_back(150); // JZB cut at 150
|
31 |
|
|
jzb_cuts.push_back(200); // JZB cut at 200
|
32 |
|
|
jzb_cuts.push_back(250); // JZB cut at 250
|
33 |
|
|
|
34 |
|
|
|
35 |
buchmann |
1.11 |
Nobs.push_back(400); // JZB cut at 50
|
36 |
|
|
Nobs.push_back(83); // JZB cut at 100
|
37 |
|
|
Nobs.push_back(20); // JZB cut at 150
|
38 |
|
|
Nobs.push_back(5); // JZB cut at 200
|
39 |
|
|
Nobs.push_back(3); // JZB cut at 250
|
40 |
|
|
Npred.push_back(377.667); // JZB cut at 50
|
41 |
|
|
Npred.push_back(79.6667); // JZB cut at 100
|
42 |
|
|
Npred.push_back(19.6667); // JZB cut at 150
|
43 |
|
|
Npred.push_back(6); // JZB cut at 200
|
44 |
|
|
Npred.push_back(1.66667); // JZB cut at 250
|
45 |
|
|
Nprederr.push_back(90.4941); // JZB cut at 50
|
46 |
|
|
Nprederr.push_back(22.1861); // JZB cut at 100
|
47 |
|
|
Nprederr.push_back(5.96285); // JZB cut at 150
|
48 |
|
|
Nprederr.push_back(2.24846); // JZB cut at 200
|
49 |
|
|
Nprederr.push_back(0.897527); // JZB cut at 250
|
50 |
buchmann |
1.10 |
|
51 |
|
|
|
52 |
buchmann |
1.11 |
flippedNobs.push_back(145); // JZB cut at 50
|
53 |
|
|
flippedNobs.push_back(10); // JZB cut at 100
|
54 |
buchmann |
1.10 |
flippedNobs.push_back(2); // JZB cut at 150
|
55 |
|
|
flippedNobs.push_back(0); // JZB cut at 200
|
56 |
|
|
flippedNobs.push_back(0); // JZB cut at 250
|
57 |
buchmann |
1.11 |
flippedNpred.push_back(167.333); // JZB cut at 50
|
58 |
|
|
flippedNpred.push_back(13.3333); // JZB cut at 100
|
59 |
|
|
flippedNpred.push_back(2.33333); // JZB cut at 150
|
60 |
|
|
flippedNpred.push_back(-1); // JZB cut at 200
|
61 |
|
|
flippedNpred.push_back(1.33333); // JZB cut at 250
|
62 |
|
|
flippedNprederr.push_back(130.882); // JZB cut at 50
|
63 |
|
|
flippedNprederr.push_back(31.4569); // JZB cut at 100
|
64 |
|
|
flippedNprederr.push_back(8.84904); // JZB cut at 150
|
65 |
|
|
flippedNprederr.push_back(3.40853); // JZB cut at 200
|
66 |
|
|
flippedNprederr.push_back(2.08999); // JZB cut at 250
|
67 |
buchmann |
1.2 |
|
68 |
|
|
|
69 |
buchmann |
1.11 |
luminosity=4653.74;
|
70 |
buchmann |
1.4 |
|
71 |
buchmann |
1.1 |
|
72 |
|
|
cout << "Configuration successfully loaded!" << endl;
|
73 |
|
|
|
74 |
|
|
}
|
75 |
|
|
|