302 bool flag_Quadratic=
false;
309 double sigma_tt_bin_300_380_madgraph = 54.003200;
318 double total = SM_sigma_tt_bin_300_380 + (+-1.080720*
ewgc(
"CG")+0.015822*
ewgc(
"Cqd1R",2,2,0,0)
319 +0.027498*
ewgc(
"Cqd1R",2,2,1,1)+0.239796*
ewgc(
"Cqd8R",2,2,0,0)+0.073740*
ewgc(
"Cqd8R",2,2,1,1)
320 +0.089304*
ewgc(
"Cqq1R",0,0,2,2)+1.256970*
ewgc(
"Cqq1R",0,2,2,0)+0.033486*
ewgc(
"Cqq1R",1,1,2,2)
321 +0.108318*
ewgc(
"Cqq1R",1,2,2,1)+0.230184*
ewgc(
"Cqq3R",0,0,2,2)+2.777490*
ewgc(
"Cqq3R",0,2,2,0)
322 +0.067932*
ewgc(
"Cqq3R",1,1,2,2)+0.457110*
ewgc(
"Cqq3R",1,2,2,1)+0.071850*
ewgc(
"Cqu1R",0,0,2,2)
323 +0.029322*
ewgc(
"Cqu1R",1,1,2,2)+0.080706*
ewgc(
"Cqu1R",2,2,0,0)+0.034674*
ewgc(
"Cqu1R",2,2,1,1)
324 +0.535692*
ewgc(
"Cqu8R",0,0,2,2)+0.091128*
ewgc(
"Cqu8R",1,1,2,2)+0.347364*
ewgc(
"Cqu8R",2,2,0,0)
325 +0.057720*
ewgc(
"Cqu8R",2,2,1,1)+0.014940*
ewgc(
"Cud1R",2,2,0,0)+0.030534*
ewgc(
"Cud1R",2,2,1,1)
326 +0.241590*
ewgc(
"Cud8R",2,2,0,0)+0.076716*
ewgc(
"Cud8R",2,2,1,1)+-18.341500*
ewgc(
"CuGR",2,2)
327 +0.081258*
ewgc(
"CuuR",0,0,2,2)+1.238740*
ewgc(
"CuuR",0,2,2,0)+0.036636*
ewgc(
"CuuR",1,1,2,2)
328 +0.108228*
ewgc(
"CuuR",1,2,2,1)
329 )*(SM_sigma_tt_bin_300_380/sigma_tt_bin_300_380_madgraph);
331 if (total < 0)
return std::numeric_limits<double>::quiet_NaN();
341 double sigma_tt_bin_380_470_madgraph = 173.423000;
350 double total = SM_sigma_tt_bin_380_470 + (+-8.359080*
ewgc(
"CG")+-0.042078*
ewgc(
"Cqd1R",2,2,0,0)+0.004590*
ewgc(
"Cqd1R",2,2,1,1)
351 +0.649842*
ewgc(
"Cqd8R",2,2,0,0)+0.144468*
ewgc(
"Cqd8R",2,2,1,1)+0.260772*
ewgc(
"Cqq1R",0,0,2,2)+4.262270*
ewgc(
"Cqq1R",0,2,2,0)
352 +-0.029244*
ewgc(
"Cqq1R",1,1,2,2)+0.230286*
ewgc(
"Cqq1R",1,2,2,1)+0.780282*
ewgc(
"Cqq3R",0,0,2,2)+9.255620*
ewgc(
"Cqq3R",0,2,2,0)
353 +0.084678*
ewgc(
"Cqq3R",1,1,2,2)+1.257790*
ewgc(
"Cqq3R",1,2,2,1)+0.071382*
ewgc(
"Cqu1R",0,0,2,2)+0.007008*
ewgc(
"Cqu1R",1,1,2,2)
354 +0.097896*
ewgc(
"Cqu1R",2,2,0,0)+0.018918*
ewgc(
"Cqu1R",2,2,1,1)+1.653840*
ewgc(
"Cqu8R",0,0,2,2)+0.194502*
ewgc(
"Cqu8R",1,1,2,2)
355 +1.000630*
ewgc(
"Cqu8R",2,2,0,0)+0.057138*
ewgc(
"Cqu8R",2,2,1,1)+-0.038910*
ewgc(
"Cud1R",2,2,0,0)+0.004434*
ewgc(
"Cud1R",2,2,1,1)
356 +0.648138*
ewgc(
"Cud8R",2,2,0,0)+0.144918*
ewgc(
"Cud8R",2,2,1,1)+-51.556900*
ewgc(
"CuGR",2,2)+0.308838*
ewgc(
"CuuR",0,0,2,2)
357 +4.176160*
ewgc(
"CuuR",0,2,2,0)+0.030888*
ewgc(
"CuuR",1,1,2,2)+0.226812*
ewgc(
"CuuR",1,2,2,1)
358 )*(SM_sigma_tt_bin_380_470/sigma_tt_bin_380_470_madgraph);
360 if (total < 0)
return std::numeric_limits<double>::quiet_NaN();
369 double sigma_tt_bin_470_620_madgraph = 152.893000;
378 double total = SM_sigma_tt_bin_470_620 + (+-12.015500*
ewgc(
"CG")+-0.014436*
ewgc(
"Cqd1R",2,2,0,0)+-0.010944*
ewgc(
"Cqd1R",2,2,1,1)
379 +0.738569*
ewgc(
"Cqd8R",2,2,0,0)+0.127001*
ewgc(
"Cqd8R",2,2,1,1)+0.294306*
ewgc(
"Cqq1R",0,0,2,2)+4.783330*
ewgc(
"Cqq1R",0,2,2,0)
380 +-0.041070*
ewgc(
"Cqq1R",1,1,2,2)+0.199307*
ewgc(
"Cqq1R",1,2,2,1)+1.030540*
ewgc(
"Cqq3R",0,0,2,2)+10.384200*
ewgc(
"Cqq3R",0,2,2,0)
381 +0.084312*
ewgc(
"Cqq3R",1,1,2,2)+1.174270*
ewgc(
"Cqq3R",1,2,2,1)+0.134514*
ewgc(
"Cqu1R",0,0,2,2)+-0.007254*
ewgc(
"Cqu1R",1,1,2,2)
382 +0.095286*
ewgc(
"Cqu1R",2,2,0,0)+0.000101*
ewgc(
"Cqu1R",2,2,1,1)+1.890170*
ewgc(
"Cqu8R",0,0,2,2)+0.177748*
ewgc(
"Cqu8R",1,1,2,2)
383 +1.147850*
ewgc(
"Cqu8R",2,2,0,0)+0.048700*
ewgc(
"Cqu8R",2,2,1,1)+-0.042132*
ewgc(
"Cud1R",2,2,0,0)+-0.013392*
ewgc(
"Cud1R",2,2,1,1)
384 +0.739715*
ewgc(
"Cud8R",2,2,0,0)+0.126949*
ewgc(
"Cud8R",2,2,1,1)+-42.458000*
ewgc(
"CuGR",2,2)+0.326724*
ewgc(
"CuuR",0,0,2,2)
385 +4.680410*
ewgc(
"CuuR",0,2,2,0)+0.007338*
ewgc(
"CuuR",1,1,2,2)+0.194285*
ewgc(
"CuuR",1,2,2,1)
386 )*(SM_sigma_tt_bin_470_620/sigma_tt_bin_470_620_madgraph);
388 if (total < 0)
return std::numeric_limits<double>::quiet_NaN();
397 double sigma_tt_bin_620_820_madgraph = 66.954300;
406 double total = SM_sigma_tt_bin_620_820 + (+-6.708550*
ewgc(
"CG")+-0.017060*
ewgc(
"Cqd1R",2,2,0,0)+0.003028*
ewgc(
"Cqd1R",2,2,1,1)
407 +0.511089*
ewgc(
"Cqd8R",2,2,0,0)+0.076736*
ewgc(
"Cqd8R",2,2,1,1)+0.224652*
ewgc(
"Cqq1R",0,0,2,2)+3.446780*
ewgc(
"Cqq1R",0,2,2,0)
408 +-0.023712*
ewgc(
"Cqq1R",1,1,2,2)+0.115517*
ewgc(
"Cqq1R",1,2,2,1)+0.804056*
ewgc(
"Cqq3R",0,0,2,2)+7.363740*
ewgc(
"Cqq3R",0,2,2,0)
409 +0.059872*
ewgc(
"Cqq3R",1,1,2,2)+0.658982*
ewgc(
"Cqq3R",1,2,2,1)+0.068678*
ewgc(
"Cqu1R",0,0,2,2)+0.006047*
ewgc(
"Cqu1R",1,1,2,2)
410 +0.049867*
ewgc(
"Cqu1R",2,2,0,0)+0.005234*
ewgc(
"Cqu1R",2,2,1,1)+1.331380*
ewgc(
"Cqu8R",0,0,2,2)+0.102676*
ewgc(
"Cqu8R",1,1,2,2)
411 +0.823487*
ewgc(
"Cqu8R",2,2,0,0)+0.029953*
ewgc(
"Cqu8R",2,2,1,1)+-0.042581*
ewgc(
"Cud1R",2,2,0,0)+-0.001945*
ewgc(
"Cud1R",2,2,1,1)
412 +0.511913*
ewgc(
"Cud8R",2,2,0,0)+0.077032*
ewgc(
"Cud8R",2,2,1,1)+-18.574200*
ewgc(
"CuGR",2,2)+0.262504*
ewgc(
"CuuR",0,0,2,2)
413 +3.359280*
ewgc(
"CuuR",0,2,2,0)+0.010662*
ewgc(
"CuuR",1,1,2,2)+0.112864*
ewgc(
"CuuR",1,2,2,1)
414 )*(SM_sigma_tt_bin_620_820/sigma_tt_bin_620_820_madgraph);
416 if (total < 0)
return std::numeric_limits<double>::quiet_NaN();
425 double sigma_tt_bin_820_1100_madgraph = 24.163600;
434 double total = SM_sigma_tt_bin_820_1100 + (-2.802570*
ewgc(
"CG")+-0.006901*
ewgc(
"Cqd1R",2,2,0,0)+0.000237*
ewgc(
"Cqd1R",2,2,1,1)
435 +0.333880*
ewgc(
"Cqd8R",2,2,0,0)+0.037791*
ewgc(
"Cqd8R",2,2,1,1)+0.161210*
ewgc(
"Cqq1R",0,0,2,2)+2.286240*
ewgc(
"Cqq1R",0,2,2,0)
436 +-0.012073*
ewgc(
"Cqq1R",1,1,2,2)+0.059007*
ewgc(
"Cqq1R",1,2,2,1)+0.566288*
ewgc(
"Cqq3R",0,0,2,2)+4.802780*
ewgc(
"Cqq3R",0,2,2,0)
437 +0.033808*
ewgc(
"Cqq3R",1,1,2,2)+0.333648*
ewgc(
"Cqq3R",1,2,2,1)+0.042755*
ewgc(
"Cqu1R",0,0,2,2)+0.002326*
ewgc(
"Cqu1R",1,1,2,2)
438 +0.035514*
ewgc(
"Cqu1R",2,2,0,0)+0.001443*
ewgc(
"Cqu1R",2,2,1,1)+0.875141*
ewgc(
"Cqu8R",0,0,2,2)+0.051098*
ewgc(
"Cqu8R",1,1,2,2)
439 +0.542868*
ewgc(
"Cqu8R",2,2,0,0)+0.013234*
ewgc(
"Cqu8R",2,2,1,1)+-0.023700*
ewgc(
"Cud1R",2,2,0,0)+-0.001051*
ewgc(
"Cud1R",2,2,1,1)
440 +0.332700*
ewgc(
"Cud8R",2,2,0,0)+0.037501*
ewgc(
"Cud8R",2,2,1,1)-6.912790*
ewgc(
"CuGR",2,2)+0.178817*
ewgc(
"CuuR",0,0,2,2)
441 +2.227770*
ewgc(
"CuuR",0,2,2,0)+0.006347*
ewgc(
"CuuR",1,1,2,2)+0.057377*
ewgc(
"CuuR",1,2,2,1)
442 )*(SM_sigma_tt_bin_820_1100/sigma_tt_bin_820_1100_madgraph);
444 if (total < 0)
return std::numeric_limits<double>::quiet_NaN();
452 double SM_sigma_tt_bin_1100_1500 =
SM.
getOptionalParameter(
"SM_sigma_tt_bin_mtt_1100_1500_CMS_181106625");
453 double sigma_tt_bin_1100_1500_madgraph = 6.983800;
462 double total = SM_sigma_tt_bin_1100_1500 + (+-0.863200*
ewgc(
"CG")+-0.004510*
ewgc(
"Cqd1R",2,2,0,0)+0.000022*
ewgc(
"Cqd1R",2,2,1,1)
463 +0.191513*
ewgc(
"Cqd8R",2,2,0,0)+0.016414*
ewgc(
"Cqd8R",2,2,1,1)+0.103366*
ewgc(
"Cqq1R",0,0,2,2)+1.351660*
ewgc(
"Cqq1R",0,2,2,0)
464 +-0.005690*
ewgc(
"Cqq1R",1,1,2,2)+0.026285*
ewgc(
"Cqq1R",1,2,2,1)+0.355618*
ewgc(
"Cqq3R",0,0,2,2)+2.795780*
ewgc(
"Cqq3R",0,2,2,0)
465 +0.014960*
ewgc(
"Cqq3R",1,1,2,2)+0.147160*
ewgc(
"Cqq3R",1,2,2,1)+0.024812*
ewgc(
"Cqu1R",0,0,2,2)+0.001067*
ewgc(
"Cqu1R",1,1,2,2)
466 +0.017626*
ewgc(
"Cqu1R",2,2,0,0)+0.000632*
ewgc(
"Cqu1R",2,2,1,1)+0.509667*
ewgc(
"Cqu8R",0,0,2,2)+0.022537*
ewgc(
"Cqu8R",1,1,2,2)
467 +0.318022*
ewgc(
"Cqu8R",2,2,0,0)+0.006561*
ewgc(
"Cqu8R",2,2,1,1)+-0.016528*
ewgc(
"Cud1R",2,2,0,0)+-0.001045*
ewgc(
"Cud1R",2,2,1,1)
468 +0.191209*
ewgc(
"Cud8R",2,2,0,0)+0.016353*
ewgc(
"Cud8R",2,2,1,1)+-2.105990*
ewgc(
"CuGR",2,2)+0.112485*
ewgc(
"CuuR",0,0,2,2)
469 +1.312690*
ewgc(
"CuuR",0,2,2,0)+0.002397*
ewgc(
"CuuR",1,1,2,2)+0.025537*
ewgc(
"CuuR",1,2,2,1)
470 )*(SM_sigma_tt_bin_1100_1500/sigma_tt_bin_1100_1500_madgraph);
472 if (total < 0)
return std::numeric_limits<double>::quiet_NaN();
480 double SM_sigma_tt_bin_1500_2500 =
SM.
getOptionalParameter(
"SM_sigma_tt_bin_mtt_1500_2500_CMS_181106625");
481 double sigma_tt_bin_1500_2500_madgraph = 1.821290;
490 double total = SM_sigma_tt_bin_1500_2500 + (+-0.231926*
ewgc(
"CG")+-0.003089*
ewgc(
"Cqd1R",2,2,0,0)+-0.000047*
ewgc(
"Cqd1R",2,2,1,1)
491 +0.130652*
ewgc(
"Cqd8R",2,2,0,0)+0.007679*
ewgc(
"Cqd8R",2,2,1,1)+0.078637*
ewgc(
"Cqq1R",0,0,2,2)+0.983935*
ewgc(
"Cqq1R",0,2,2,0)
492 +-0.002654*
ewgc(
"Cqq1R",1,1,2,2)+0.014376*
ewgc(
"Cqq1R",1,2,2,1)+0.252179*
ewgc(
"Cqq3R",0,0,2,2)+1.968040*
ewgc(
"Cqq3R",0,2,2,0)
493 +0.007838*
ewgc(
"Cqq3R",1,1,2,2)+0.071898*
ewgc(
"Cqq3R",1,2,2,1)+0.017883*
ewgc(
"Cqu1R",0,0,2,2)+0.000657*
ewgc(
"Cqu1R",1,1,2,2)
494 +0.011678*
ewgc(
"Cqu1R",2,2,0,0)+0.000359*
ewgc(
"Cqu1R",2,2,1,1)+0.363022*
ewgc(
"Cqu8R",0,0,2,2)+0.011088*
ewgc(
"Cqu8R",1,1,2,2)
495 +0.232613*
ewgc(
"Cqu8R",2,2,0,0)+0.003450*
ewgc(
"Cqu8R",2,2,1,1)+-0.011469*
ewgc(
"Cud1R",2,2,0,0)+-0.000549*
ewgc(
"Cud1R",2,2,1,1)
496 +0.130711*
ewgc(
"Cud8R",2,2,0,0)+0.007725*
ewgc(
"Cud8R",2,2,1,1)+-0.607749*
ewgc(
"CuGR",2,2)+0.081595*
ewgc(
"CuuR",0,0,2,2)
497 +0.955959*
ewgc(
"CuuR",0,2,2,0)+0.001335*
ewgc(
"CuuR",1,1,2,2)+0.014025*
ewgc(
"CuuR",1,2,2,1)
498 )*(SM_sigma_tt_bin_1500_2500/sigma_tt_bin_1500_2500_madgraph);
500 if (total < 0)
return std::numeric_limits<double>::quiet_NaN();
509 throw std::runtime_error(
"\nERROR: Please specify a correct binning range for sigma_tt_diff_mtt_LO_CMS_181106625.\n");
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.
double ewgc(const std::string name) const