| 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 |
float DataPeakError=-999;
|
| 19 |
}
|
| 20 |
|
| 21 |
void read_config() {
|
| 22 |
datajzb="((jzb[1]+0.053*pt)-2.95688)";
|
| 23 |
mcjzb="((jzb[1]+0.043*pt)-3.05937)";
|
| 24 |
|
| 25 |
|
| 26 |
|
| 27 |
MCPeakError=0.159709;
|
| 28 |
DataPeakError=0.105771;
|
| 29 |
|
| 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 |
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 |
flippedNobs.push_back(1); // JZB cut at 200
|
| 58 |
flippedNobs.push_back(0); // JZB cut at 250
|
| 59 |
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 |
|
| 70 |
|
| 71 |
luminosity=4980;
|
| 72 |
RestrictToMassPeak=0;//defines the type of analysis we're running
|
| 73 |
|
| 74 |
|
| 75 |
cout << "Configuration successfully loaded!" << endl;
|
| 76 |
|
| 77 |
}
|
| 78 |
|
