36 bool flag_Quadratic=
false;
48 double SM_FB_asymmetry_deno_bin_mtt_350_450 =
SM.
getOptionalParameter(
"SM_FB_asymmetry_deno_bin_mtt_350_450");
49 double SM_FB_asymmetry_num_bin_mtt_350_450 = SM_FB_asymmetry_bin_mtt_350_450*SM_FB_asymmetry_deno_bin_mtt_350_450;
51 double SM_sigma_pos_bin_mtt_350_450 =0.5*(SM_FB_asymmetry_num_bin_mtt_350_450+SM_FB_asymmetry_deno_bin_mtt_350_450);
52 double SM_sigma_neg_bin_mtt_350_450 =0.5*(SM_FB_asymmetry_deno_bin_mtt_350_450-SM_FB_asymmetry_num_bin_mtt_350_450);
64 double sigma_pos_bin_mtt_350_450_madgraph = 1.675650;
65 double sigma_pos_bin_mtt_350_450_NP = -0.006110*
ewgc(
"CG")+0.015540*
ewgc(
"Cqd8R",2,2,0,0)
66 +0.304934*
ewgc(
"Cqq1R",0,2,2,0)+0.419488*
ewgc(
"Cqq3R",0,2,2,0)+0.004257*
ewgc(
"Cqq3R",1,2,2,1)
67 +0.087584*
ewgc(
"Cqu8R",0,0,2,2)+0.073317*
ewgc(
"Cqu8R",2,2,0,0)+0.014866*
ewgc(
"Cud8R",2,2,0,0)
68 -0.385947*
ewgc(
"CuGR",2,2)+0.298571*
ewgc(
"CuuR",0,2,2,0)
69 +-0.000866*
ewgc(
"Cqd1R",2,2,0,0)+-0.000177*
ewgc(
"Cqd1R",2,2,1,1)+0.030033*
ewgc(
"Cqq1R",0,0,2,2)
70 +-0.000533*
ewgc(
"Cqq1R",1,1,2,2)+0.040348*
ewgc(
"Cqq3R",0,0,2,2)+-0.000078*
ewgc(
"Cqq3R",1,1,2,2)
71 +0.011868*
ewgc(
"Cqu1R",0,0,2,2)+0.000311*
ewgc(
"Cqu1R",1,1,2,2)+0.007404*
ewgc(
"Cqu1R",2,2,0,0)
72 +0.000365*
ewgc(
"Cqu1R",2,2,1,1)+-0.000412*
ewgc(
"Cud1R",2,2,0,0)+0.000135*
ewgc(
"Cud1R",2,2,1,1)
73 +0.014596*
ewgc(
"CuuR",0,0,2,2)+-0.000309*
ewgc(
"CuuR",1,1,2,2);
74 double sigma_neg_bin_mtt_350_450_madgraph = 1.674510;
75 double sigma_neg_bin_mtt_350_450_NP = -0.007566*
ewgc(
"CG")+0.014887*
ewgc(
"Cqd8R",2,2,0,0)
76 +0.305105*
ewgc(
"Cqq1R",0,2,2,0)+0.418602*
ewgc(
"Cqq3R",0,2,2,0)+0.002955*
ewgc(
"Cqq3R",1,2,2,1)
77 +0.085906*
ewgc(
"Cqu8R",0,0,2,2)+0.073521*
ewgc(
"Cqu8R",2,2,0,0)+0.015697*
ewgc(
"Cud8R",2,2,0,0)
78 -0.383667*
ewgc(
"CuGR",2,2)+0.300446*
ewgc(
"CuuR",0,2,2,0)
79 -0.001188*
ewgc(
"Cqd1R",2,2,0,0)+0.000242*
ewgc(
"Cqd1R",2,2,1,1)+0.028715*
ewgc(
"Cqq1R",0,0,2,2)
80 +0.000588*
ewgc(
"Cqq1R",1,1,2,2)+0.041531*
ewgc(
"Cqq3R",0,0,2,2)+0.000146*
ewgc(
"Cqq3R",1,1,2,2)
81 +0.013026*
ewgc(
"Cqu1R",0,0,2,2)+-0.000294*
ewgc(
"Cqu1R",1,1,2,2)+0.007109*
ewgc(
"Cqu1R",2,2,0,0)
82 +-0.000347*
ewgc(
"Cqu1R",2,2,1,1)+-0.001721*
ewgc(
"Cud1R",2,2,0,0)+-0.000067*
ewgc(
"Cud1R",2,2,1,1)
83 +0.016846*
ewgc(
"CuuR",0,0,2,2)+0.000338*
ewgc(
"CuuR",1,1,2,2);
85 double sigma_pos_bin_mtt_350_450_NP_Corrected = SM_sigma_pos_bin_mtt_350_450*sigma_pos_bin_mtt_350_450_NP/sigma_pos_bin_mtt_350_450_madgraph;
86 double sigma_neg_bin_mtt_350_450_NP_Corrected = SM_sigma_neg_bin_mtt_350_450*sigma_neg_bin_mtt_350_450_NP/sigma_neg_bin_mtt_350_450_madgraph;
88 double NP_FB_asymmetry_num_bin_mtt_350_450 = sigma_pos_bin_mtt_350_450_NP_Corrected - sigma_neg_bin_mtt_350_450_NP_Corrected;
89 double NP_FB_asymmetry_deno_bin_mtt_350_450 = sigma_pos_bin_mtt_350_450_NP_Corrected + sigma_neg_bin_mtt_350_450_NP_Corrected;;
93 return SM_FB_asymmetry_bin_mtt_350_450 + (NP_FB_asymmetry_num_bin_mtt_350_450
94 - SM_FB_asymmetry_bin_mtt_350_450*NP_FB_asymmetry_deno_bin_mtt_350_450
95 )/SM_FB_asymmetry_deno_bin_mtt_350_450;
104 double SM_FB_asymmetry_deno_bin_mtt_450_550 =
SM.
getOptionalParameter(
"SM_FB_asymmetry_deno_bin_mtt_450_550");
105 double SM_FB_asymmetry_num_bin_mtt_450_550 = SM_FB_asymmetry_bin_mtt_450_550*SM_FB_asymmetry_deno_bin_mtt_450_550;
107 double SM_sigma_pos_bin_mtt_450_550 =0.5*(SM_FB_asymmetry_num_bin_mtt_450_550+SM_FB_asymmetry_deno_bin_mtt_450_550);
108 double SM_sigma_neg_bin_mtt_450_550 =0.5*(SM_FB_asymmetry_deno_bin_mtt_450_550-SM_FB_asymmetry_num_bin_mtt_450_550);
120 double sigma_pos_bin_mtt_450_550_madgraph = 0.683803;
121 double sigma_pos_bin_mtt_450_550_NP = -0.003753*
ewgc(
"CG")+0.007844*
ewgc(
"Cqd8R",2,2,0,0)
122 +0.201208*
ewgc(
"Cqq1R",0,2,2,0)+0.257805*
ewgc(
"Cqq3R",0,2,2,0)+0.000709*
ewgc(
"Cqq3R",1,2,2,1)
123 +0.056092*
ewgc(
"Cqu8R",0,0,2,2)+0.047574*
ewgc(
"Cqu8R",2,2,0,0)+0.008085*
ewgc(
"Cud8R",2,2,0,0)
124 -0.168317*
ewgc(
"CuGR",2,2)+0.198755*
ewgc(
"CuuR",0,2,2,0)
125 +-0.000527*
ewgc(
"Cqd1R",2,2,0,0)+0.000144*
ewgc(
"Cqd1R",2,2,1,1)+0.022808*
ewgc(
"Cqq1R",0,0,2,2)
126 +0.000181*
ewgc(
"Cqq1R",1,1,2,2)+0.030815*
ewgc(
"Cqq3R",0,0,2,2)+0.000065*
ewgc(
"Cqq3R",1,1,2,2)
127 +0.005817*
ewgc(
"Cqu1R",0,0,2,2)+0.000027*
ewgc(
"Cqu1R",1,1,2,2)+0.003722*
ewgc(
"Cqu1R",2,2,0,0)
128 +0.000096*
ewgc(
"Cqu1R",2,2,1,1)+0.000104*
ewgc(
"Cud1R",2,2,0,0)+0.000087*
ewgc(
"Cud1R",2,2,1,1)
129 +0.013193*
ewgc(
"CuuR",0,0,2,2)+0.000111*
ewgc(
"CuuR",1,1,2,2);
130 double sigma_neg_bin_mtt_450_550_madgraph = 0.685429;
131 double sigma_neg_bin_mtt_450_550_NP = -0.004203*
ewgc(
"CG")+0.008080*
ewgc(
"Cqd8R",2,2,0,0)
132 +0.201688*
ewgc(
"Cqq1R",0,2,2,0)+0.257507*
ewgc(
"Cqq3R",0,2,2,0)+0.000816*
ewgc(
"Cqq3R",1,2,2,1)
133 +0.056997*
ewgc(
"Cqu8R",0,0,2,2)+0.047589*
ewgc(
"Cqu8R",2,2,0,0)+0.007820*
ewgc(
"Cud8R",2,2,0,0)
134 -0.169079*
ewgc(
"CuGR",2,2)+0.198609*
ewgc(
"CuuR",0,2,2,0)
135 +-0.000152*
ewgc(
"Cqd1R",2,2,0,0)+-0.000219*
ewgc(
"Cqd1R",2,2,1,1)+0.022948*
ewgc(
"Cqq1R",0,0,2,2)
136 +-0.000255*
ewgc(
"Cqq1R",1,1,2,2)+0.030379*
ewgc(
"Cqq3R",0,0,2,2)+-0.000024*
ewgc(
"Cqq3R",1,1,2,2)
137 +0.006338*
ewgc(
"Cqu1R",0,0,2,2)+-0.000058*
ewgc(
"Cqu1R",1,1,2,2)+0.004314*
ewgc(
"Cqu1R",2,2,0,0)
138 +-0.000126*
ewgc(
"Cqu1R",2,2,1,1)+0.000286*
ewgc(
"Cud1R",2,2,0,0)+-0.000164*
ewgc(
"Cud1R",2,2,1,1)
139 +0.013542*
ewgc(
"CuuR",0,0,2,2)+-0.000141*
ewgc(
"CuuR",1,1,2,2);
141 double sigma_pos_bin_mtt_450_550_NP_Corrected = SM_sigma_pos_bin_mtt_450_550*sigma_pos_bin_mtt_450_550_NP/sigma_pos_bin_mtt_450_550_madgraph;
142 double sigma_neg_bin_mtt_450_550_NP_Corrected = SM_sigma_neg_bin_mtt_450_550*sigma_neg_bin_mtt_450_550_NP/sigma_neg_bin_mtt_450_550_madgraph;
144 double NP_FB_asymmetry_num_bin_mtt_450_550 = sigma_pos_bin_mtt_450_550_NP_Corrected - sigma_neg_bin_mtt_450_550_NP_Corrected;
145 double NP_FB_asymmetry_deno_bin_mtt_450_550 = sigma_pos_bin_mtt_450_550_NP_Corrected + sigma_neg_bin_mtt_450_550_NP_Corrected;;
150 return SM_FB_asymmetry_bin_mtt_450_550 + (NP_FB_asymmetry_num_bin_mtt_450_550
151 - SM_FB_asymmetry_bin_mtt_450_550*NP_FB_asymmetry_deno_bin_mtt_450_550
152 )/SM_FB_asymmetry_deno_bin_mtt_450_550;
161 double SM_FB_asymmetry_deno_bin_mtt_550_650 =
SM.
getOptionalParameter(
"SM_FB_asymmetry_deno_bin_mtt_550_650");
162 double SM_FB_asymmetry_num_bin_mtt_550_650 = SM_FB_asymmetry_bin_mtt_550_650*SM_FB_asymmetry_deno_bin_mtt_550_650;
164 double SM_sigma_pos_bin_mtt_550_650 =0.5*(SM_FB_asymmetry_num_bin_mtt_550_650+SM_FB_asymmetry_deno_bin_mtt_550_650);
165 double SM_sigma_neg_bin_mtt_550_650 =0.5*(SM_FB_asymmetry_deno_bin_mtt_550_650-SM_FB_asymmetry_num_bin_mtt_550_650);
177 double sigma_pos_bin_mtt_550_650_madgraph = 0.216436;
178 double sigma_pos_bin_mtt_550_650_NP = -0.001345*
ewgc(
"CG")+0.002993*
ewgc(
"Cqd8R",2,2,0,0)
179 +0.098574*
ewgc(
"Cqq1R",0,2,2,0)+0.119856*
ewgc(
"Cqq3R",0,2,2,0)+0.000356*
ewgc(
"Cqq3R",1,2,2,1)
180 +0.026186*
ewgc(
"Cqu8R",0,0,2,2)+0.023259*
ewgc(
"Cqu8R",2,2,0,0)+0.002973*
ewgc(
"Cud8R",2,2,0,0)
181 -0.057346*
ewgc(
"CuGR",2,2)+0.095602*
ewgc(
"CuuR",0,2,2,0)
182 +0.000089*
ewgc(
"Cqd1R",2,2,0,0)+0.000160*
ewgc(
"Cqd1R",2,2,1,1)+0.012427*
ewgc(
"Cqq1R",0,0,2,2)
183 +0.000182*
ewgc(
"Cqq1R",1,1,2,2)+0.015602*
ewgc(
"Cqq3R",0,0,2,2)+0.000078*
ewgc(
"Cqq3R",1,1,2,2)
184 +0.002898*
ewgc(
"Cqu1R",0,0,2,2)+0.000184*
ewgc(
"Cqu1R",1,1,2,2)+0.001348*
ewgc(
"Cqu1R",2,2,0,0)
185 +0.000237*
ewgc(
"Cqu1R",2,2,1,1)+-0.000559*
ewgc(
"Cud1R",2,2,0,0)+0.000145*
ewgc(
"Cud1R",2,2,1,1)
186 +0.007087*
ewgc(
"CuuR",0,0,2,2)+0.000149*
ewgc(
"CuuR",1,1,2,2);
187 double sigma_neg_bin_mtt_550_650_madgraph = 0.216949;
188 double sigma_neg_bin_mtt_550_650_NP = -0.001414*
ewgc(
"CG")+0.002846*
ewgc(
"Cqd8R",2,2,0,0)
189 +0.098612*
ewgc(
"Cqq1R",0,2,2,0)+0.120103*
ewgc(
"Cqq3R",0,2,2,0)+0.000216*
ewgc(
"Cqq3R",1,2,2,1)
190 +0.025800*
ewgc(
"Cqu8R",0,0,2,2)+0.023527*
ewgc(
"Cqu8R",2,2,0,0)+0.002835*
ewgc(
"Cud8R",2,2,0,0)
191 -0.057275*
ewgc(
"CuGR",2,2)+0.095612*
ewgc(
"CuuR",0,2,2,0)
192 +-0.000225*
ewgc(
"Cqd1R",2,2,0,0)+-0.000155*
ewgc(
"Cqd1R",2,2,1,1)+0.012485*
ewgc(
"Cqq1R",0,0,2,2)
193 +-0.000175*
ewgc(
"Cqq1R",1,1,2,2)+0.015723*
ewgc(
"Cqq3R",0,0,2,2)+-0.000069*
ewgc(
"Cqq3R",1,1,2,2)
194 +0.002745*
ewgc(
"Cqu1R",0,0,2,2)+-0.000178*
ewgc(
"Cqu1R",1,1,2,2)+0.001266*
ewgc(
"Cqu1R",2,2,0,0)
195 +-0.000230*
ewgc(
"Cqu1R",2,2,1,1)+-0.000703*
ewgc(
"Cud1R",2,2,0,0)+-0.000139*
ewgc(
"Cud1R",2,2,1,1)
196 +0.006972*
ewgc(
"CuuR",0,0,2,2)+-0.000143*
ewgc(
"CuuR",1,1,2,2);
198 double sigma_pos_bin_mtt_550_650_NP_Corrected = SM_sigma_pos_bin_mtt_550_650*sigma_pos_bin_mtt_550_650_NP/sigma_pos_bin_mtt_550_650_madgraph;
199 double sigma_neg_bin_mtt_550_650_NP_Corrected = SM_sigma_neg_bin_mtt_550_650*sigma_neg_bin_mtt_550_650_NP/sigma_neg_bin_mtt_550_650_madgraph;
201 double NP_FB_asymmetry_num_bin_mtt_550_650 = sigma_pos_bin_mtt_550_650_NP_Corrected - sigma_neg_bin_mtt_550_650_NP_Corrected;
202 double NP_FB_asymmetry_deno_bin_mtt_550_650 = sigma_pos_bin_mtt_550_650_NP_Corrected + sigma_neg_bin_mtt_550_650_NP_Corrected;;
206 return SM_FB_asymmetry_bin_mtt_550_650 + (NP_FB_asymmetry_num_bin_mtt_550_650
207 - SM_FB_asymmetry_bin_mtt_550_650*NP_FB_asymmetry_deno_bin_mtt_550_650
208 )/SM_FB_asymmetry_deno_bin_mtt_550_650;
216 double SM_FB_asymmetry_deno_bin_mtt_650_750 =
SM.
getOptionalParameter(
"SM_FB_asymmetry_deno_bin_mtt_650_750");
217 double SM_FB_asymmetry_num_bin_mtt_650_750 = SM_FB_asymmetry_bin_mtt_650_750*SM_FB_asymmetry_deno_bin_mtt_650_750;
219 double SM_sigma_pos_bin_mtt_650_750 =0.5*(SM_FB_asymmetry_num_bin_mtt_650_750+SM_FB_asymmetry_deno_bin_mtt_650_750);
220 double SM_sigma_neg_bin_mtt_650_750 =0.5*(SM_FB_asymmetry_deno_bin_mtt_650_750-SM_FB_asymmetry_num_bin_mtt_650_750);
232 double sigma_pos_bin_mtt_650_750_madgraph = 0.097727;
233 double sigma_pos_bin_mtt_650_750_NP = -0.000271*
ewgc(
"CG")+0.001360*
ewgc(
"Cqd8R",2,2,0,0)
234 +0.073375*
ewgc(
"Cqq1R",0,2,2,0)+0.084010*
ewgc(
"Cqq3R",0,2,2,0)+0.000196*
ewgc(
"Cqq3R",1,2,2,1)
235 +0.018728*
ewgc(
"Cqu8R",0,0,2,2)+0.017419*
ewgc(
"Cqu8R",2,2,0,0)+0.001357*
ewgc(
"Cud8R",2,2,0,0)
236 -0.026810*
ewgc(
"CuGR",2,2)+0.071319*
ewgc(
"CuuR",0,2,2,0)
237 +-0.000015*
ewgc(
"Cqd1R",2,2,0,0)+0.000006*
ewgc(
"Cqd1R",2,2,1,1)+0.010420*
ewgc(
"Cqq1R",0,0,2,2)
238 +0.000021*
ewgc(
"Cqq1R",1,1,2,2)+0.012135*
ewgc(
"Cqq3R",0,0,2,2)+0.000029*
ewgc(
"Cqq3R",1,1,2,2)
239 +0.001508*
ewgc(
"Cqu1R",0,0,2,2)+0.000036*
ewgc(
"Cqu1R",1,1,2,2)+0.001196*
ewgc(
"Cqu1R",2,2,0,0)
240 +0.000014*
ewgc(
"Cqu1R",2,2,1,1)+-0.000041*
ewgc(
"Cud1R",2,2,0,0)+0.000016*
ewgc(
"Cud1R",2,2,1,1)
241 +0.005619*
ewgc(
"CuuR",0,0,2,2)+0.000034*
ewgc(
"CuuR",1,1,2,2);
242 double sigma_neg_bin_mtt_650_750_madgraph = 0.097766;
243 double sigma_neg_bin_mtt_650_750_NP = -0.000309*
ewgc(
"CG")+0.001327*
ewgc(
"Cqd8R",2,2,0,0)
244 +0.073361*
ewgc(
"Cqq1R",0,2,2,0)+0.083778*
ewgc(
"Cqq3R",0,2,2,0)+0.000188*
ewgc(
"Cqq3R",1,2,2,1)
245 +0.018756*
ewgc(
"Cqu8R",0,0,2,2)+0.017360*
ewgc(
"Cqu8R",2,2,0,0)+0.001327*
ewgc(
"Cud8R",2,2,0,0)
246 -0.026852*
ewgc(
"CuGR",2,2)+0.071283*
ewgc(
"CuuR",0,2,2,0)
247 +-0.000003*
ewgc(
"Cqd1R",2,2,0,0)+0.000000*
ewgc(
"Cqd1R",2,2,1,1)+0.010473*
ewgc(
"Cqq1R",0,0,2,2)
248 +-0.000013*
ewgc(
"Cqq1R",1,1,2,2)+0.012141*
ewgc(
"Cqq3R",0,0,2,2)+-0.000012*
ewgc(
"Cqq3R",1,1,2,2)
249 +0.001523*
ewgc(
"Cqu1R",0,0,2,2)+-0.000029*
ewgc(
"Cqu1R",1,1,2,2)+0.001164*
ewgc(
"Cqu1R",2,2,0,0)
250 +-0.000007*
ewgc(
"Cqu1R",2,2,1,1)+-0.000073*
ewgc(
"Cud1R",2,2,0,0)+-0.000009*
ewgc(
"Cud1R",2,2,1,1)
251 +0.005643*
ewgc(
"CuuR",0,0,2,2)+-0.000027*
ewgc(
"CuuR",1,1,2,2);
253 double sigma_pos_bin_mtt_650_750_NP_Corrected = SM_sigma_pos_bin_mtt_650_750*sigma_pos_bin_mtt_650_750_NP/sigma_pos_bin_mtt_650_750_madgraph;
254 double sigma_neg_bin_mtt_650_750_NP_Corrected = SM_sigma_neg_bin_mtt_650_750*sigma_neg_bin_mtt_650_750_NP/sigma_neg_bin_mtt_650_750_madgraph;
256 double NP_FB_asymmetry_num_bin_mtt_650_750 = sigma_pos_bin_mtt_650_750_NP_Corrected - sigma_neg_bin_mtt_650_750_NP_Corrected;
257 double NP_FB_asymmetry_deno_bin_mtt_650_750 = sigma_pos_bin_mtt_650_750_NP_Corrected + sigma_neg_bin_mtt_650_750_NP_Corrected;;
261 return SM_FB_asymmetry_bin_mtt_650_750 + (NP_FB_asymmetry_num_bin_mtt_650_750
262 - SM_FB_asymmetry_bin_mtt_650_750*NP_FB_asymmetry_deno_bin_mtt_650_750
263 )/SM_FB_asymmetry_deno_bin_mtt_650_750;
269 throw std::runtime_error(
"\nERROR: Please specify a correct binning range for FB_asymmetry_Tevatron_tt_diff_mtt_LO.\n");
double ewgc(const std::string name) const
const double getOptionalParameter(std::string name) const
A method to get parameters that are specific to only one set of observables.
double getBinMin()
A get method to get the minimum value of the bin.
double getBinMax()
A get method to get the maximum value of the bin.