| 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 |
}
|
| 19 |
|
| 20 |
void read_config() {
|
| 21 |
datajzb="(jzb[1]+0.06*pt-3.53058)";
|
| 22 |
mcjzb="(jzb[1]+0.04*pt-2.4392)";
|
| 23 |
|
| 24 |
|
| 25 |
|
| 26 |
MCPeakError=0.097537;
|
| 27 |
|
| 28 |
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 |
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 |
|
| 51 |
|
| 52 |
flippedNobs.push_back(145); // JZB cut at 50
|
| 53 |
flippedNobs.push_back(10); // JZB cut at 100
|
| 54 |
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 |
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 |
|
| 68 |
|
| 69 |
luminosity=4653.74;
|
| 70 |
|
| 71 |
|
| 72 |
cout << "Configuration successfully loaded!" << endl;
|
| 73 |
|
| 74 |
}
|
| 75 |
|
