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.06*pt-3.58612)"; |
23 |
< |
mcjzb="(jzb[1]+0.04*pt-2.32211)"; |
22 |
> |
datajzb="((jzb[1]+0.053*pt)-2.95688)"; |
23 |
> |
mcjzb="((jzb[1]+0.043*pt)-3.05937)"; |
24 |
|
|
25 |
|
|
26 |
|
|
27 |
< |
MCPeakError=0.0990011; |
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 |
34 |
|
jzb_cuts.push_back(250); // JZB cut at 250 |
35 |
|
|
36 |
|
|
37 |
< |
Nobs.push_back(408); // JZB cut at 50 |
38 |
< |
Nobs.push_back(88); // JZB cut at 100 |
39 |
< |
Nobs.push_back(21); // JZB cut at 150 |
40 |
< |
Nobs.push_back(5); // JZB cut at 200 |
41 |
< |
Nobs.push_back(3); // JZB cut at 250 |
42 |
< |
Npred.push_back(407.667); // JZB cut at 50 |
43 |
< |
Npred.push_back(89.3333); // JZB cut at 100 |
44 |
< |
Npred.push_back(22); // JZB cut at 150 |
45 |
< |
Npred.push_back(8); // JZB cut at 200 |
46 |
< |
Npred.push_back(2); // JZB cut at 250 |
47 |
< |
Nprederr.push_back(61.1314); // JZB cut at 50 |
48 |
< |
Nprederr.push_back(13.6275); // JZB cut at 100 |
49 |
< |
Nprederr.push_back(4.25245); // JZB cut at 150 |
50 |
< |
Nprederr.push_back(2.10159); // JZB cut at 200 |
51 |
< |
Nprederr.push_back(0.874007); // JZB cut at 250 |
52 |
< |
|
53 |
< |
|
54 |
< |
flippedNobs.push_back(149); // JZB cut at 50 |
55 |
< |
flippedNobs.push_back(11); // JZB cut at 100 |
56 |
< |
flippedNobs.push_back(3); // JZB cut at 150 |
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(149.333); // JZB cut at 50 |
60 |
< |
flippedNpred.push_back(9.66667); // JZB cut at 100 |
61 |
< |
flippedNpred.push_back(2); // JZB cut at 150 |
62 |
< |
flippedNpred.push_back(-2); // JZB cut at 200 |
63 |
< |
flippedNpred.push_back(1); // JZB cut at 250 |
64 |
< |
flippedNprederr.push_back(114.072); // JZB cut at 50 |
65 |
< |
flippedNprederr.push_back(27.1872); // JZB cut at 100 |
66 |
< |
flippedNprederr.push_back(7.9425); // JZB cut at 150 |
67 |
< |
flippedNprederr.push_back(3.14907); // JZB cut at 200 |
68 |
< |
flippedNprederr.push_back(2.08); // 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=4653.74; |
71 |
> |
luminosity=4980; |
72 |
> |
RestrictToMassPeak=0;//defines the type of analysis we're running |
73 |
|
|
74 |
|
|
75 |
|
cout << "Configuration successfully loaded!" << endl; |