3202 bool flag_Quadratic=
false;
3207 double sigma_ttz_bin_0_40_madgraph = 0.064177;
3217 double total = (sigma_ttz_bin_0_40_madgraph -0.010478*
ewgc(
"CG")-0.005804*
ewgc(
"CHq1R",2,2)
3218 +0.005759*
ewgc(
"CHq3R",2,2)+0.002967*
ewgc(
"CHuR",2,2)+0.001126*
ewgc(
"Cqd8R",2,2,0,0)
3219 +0.000122*
ewgc(
"Cqd8R",2,2,1,1)
3220 +0.008742*
ewgc(
"Cqq1R",0,2,2,0)+0.000211*
ewgc(
"Cqq1R",1,2,2,1)
3221 +0.045418*
ewgc(
"Cqq3R",0,2,2,0)+0.004850*
ewgc(
"Cqq3R",1,2,2,1)
3222 +0.005860*
ewgc(
"Cqu8R",0,0,2,2)+0.000532*
ewgc(
"Cqu8R",1,1,2,2)
3223 +0.000976*
ewgc(
"Cqu8R",2,2,0,0)+0.000001*
ewgc(
"Cqu8R",2,2,1,1)
3224 -0.000170*
ewgc(
"CuBR",2,2)+0.000666*
ewgc(
"Cud8R",2,2,0,0)+0.000058*
ewgc(
"Cud8R",2,2,1,1)
3225 -0.018690*
ewgc(
"CuGR",2,2)+0.002251*
ewgc(
"CuuR",0,2,2,0)
3226 +0.000036*
ewgc(
"CuuR",1,2,2,1)-0.000098*
ewgc(
"CuWR",2,2)
3227 +0.000365*
ewgc(
"CHD")+-0.003957*
ewgc(
"CHl3R",0,0)+-0.003957*
ewgc(
"CHl3R",1,1)+-0.000032*
ewgc(
"CHW")
3228 +0.002458*
ewgc(
"CHWB")+0.003876*
ewgc(
"CllR",0,1,1,0)+-0.000058*
ewgc(
"Cqd1R",2,2,0,0)+-0.000030*
ewgc(
"Cqd1R",2,2,1,1)
3229 +-0.001074*
ewgc(
"Cqq1R",0,0,2,2)+-0.000307*
ewgc(
"Cqq1R",1,1,2,2)+0.004072*
ewgc(
"Cqq3R",0,0,2,2)
3230 +0.000317*
ewgc(
"Cqq3R",1,1,2,2)+-0.000065*
ewgc(
"Cqu1R",0,0,2,2)+-0.000069*
ewgc(
"Cqu1R",1,1,2,2)
3231 +0.000058*
ewgc(
"Cqu1R",2,2,0,0)+-0.000030*
ewgc(
"Cqu1R",2,2,1,1)+-0.000072*
ewgc(
"Cud1R",2,2,0,0)
3232 +-0.000031*
ewgc(
"Cud1R",2,2,1,1)+0.000196*
ewgc(
"CuuR",0,0,2,2)+-0.000029*
ewgc(
"CuuR",1,1,2,2))
3233 *(SM_sigma_ttz_bin_0_40/sigma_ttz_bin_0_40_madgraph);
3235 if (total < 0)
return std::numeric_limits<double>::quiet_NaN();
3243 double sigma_ttz_bin_40_70_madgraph = 0.096642;
3253 double total = (sigma_ttz_bin_40_70_madgraph -0.018049*
ewgc(
"CG")-0.009212*
ewgc(
"CHq1R",2,2)
3254 +0.009177*
ewgc(
"CHq3R",2,2)+0.005078*
ewgc(
"CHuR",2,2)+0.001840*
ewgc(
"Cqd8R",2,2,0,0)
3255 +0.000234*
ewgc(
"Cqd8R",2,2,1,1)
3256 +0.013449*
ewgc(
"Cqq1R",0,2,2,0)+0.000372*
ewgc(
"Cqq1R",1,2,2,1)
3257 +0.065355*
ewgc(
"Cqq3R",0,2,2,0)+0.006633*
ewgc(
"Cqq3R",1,2,2,1)
3258 +0.008494*
ewgc(
"Cqu8R",0,0,2,2)+0.000775*
ewgc(
"Cqu8R",1,1,2,2)
3259 +0.001722*
ewgc(
"Cqu8R",2,2,0,0)+0.000055*
ewgc(
"Cqu8R",2,2,1,1)
3260 -0.000190*
ewgc(
"CuBR",2,2)+0.001097*
ewgc(
"Cud8R",2,2,0,0)
3261 +0.000148*
ewgc(
"Cud8R",2,2,1,1)+-0.029042*
ewgc(
"CuGR",2,2)
3262 +0.003988*
ewgc(
"CuuR",0,2,2,0)+0.000126*
ewgc(
"CuuR",1,2,2,1)
3263 -0.000106*
ewgc(
"CuWR",2,2)
3264 +0.000406*
ewgc(
"CHD")+-0.005906*
ewgc(
"CHl3R",0,0)+-0.005906*
ewgc(
"CHl3R",1,1)+0.000019*
ewgc(
"CHW")
3265 +0.003564*
ewgc(
"CHWB")+0.005878*
ewgc(
"CllR",0,1,1,0)+-0.000048*
ewgc(
"Cqd1R",2,2,0,0)+-0.000002*
ewgc(
"Cqd1R",2,2,1,1)
3266 +-0.001543*
ewgc(
"Cqq1R",0,0,2,2)+-0.000396*
ewgc(
"Cqq1R",1,1,2,2)+0.005957*
ewgc(
"Cqq3R",0,0,2,2)
3267 +0.000477*
ewgc(
"Cqq3R",1,1,2,2)+0.000037*
ewgc(
"Cqu1R",0,0,2,2)+-0.000046*
ewgc(
"Cqu1R",1,1,2,2)
3268 +0.000093*
ewgc(
"Cqu1R",2,2,0,0)+0.000007*
ewgc(
"Cqu1R",2,2,1,1)+-0.000040*
ewgc(
"Cud1R",2,2,0,0)
3269 +0.000012*
ewgc(
"Cud1R",2,2,1,1)+0.000312*
ewgc(
"CuuR",0,0,2,2)+0.000027*
ewgc(
"CuuR",1,1,2,2))
3270 *(SM_sigma_ttz_bin_40_70/sigma_ttz_bin_40_70_madgraph);
3272 if (total < 0)
return std::numeric_limits<double>::quiet_NaN();
3281 double sigma_ttz_bin_70_110_madgraph = 0.127419;
3291 double total = (sigma_ttz_bin_70_110_madgraph-0.028383*
ewgc(
"CG")+-0.012875*
ewgc(
"CHq1R",2,2)+0.012884*
ewgc(
"CHq3R",2,2)
3292 +0.007778*
ewgc(
"CHuR",2,2)+0.002468*
ewgc(
"Cqd8R",2,2,0,0)+0.000298*
ewgc(
"Cqd8R",2,2,1,1)
3293 +0.017859*
ewgc(
"Cqq1R",0,2,2,0)+0.000464*
ewgc(
"Cqq1R",1,2,2,1)
3294 +0.080694*
ewgc(
"Cqq3R",0,2,2,0)+0.007774*
ewgc(
"Cqq3R",1,2,2,1)+0.010371*
ewgc(
"Cqu8R",0,0,2,2)
3295 +0.000893*
ewgc(
"Cqu8R",1,1,2,2)+0.002610*
ewgc(
"Cqu8R",2,2,0,0)+0.000063*
ewgc(
"Cqu8R",2,2,1,1)
3296 -0.000135*
ewgc(
"CuBR",2,2)+0.001481*
ewgc(
"Cud8R",2,2,0,0)+0.000183*
ewgc(
"Cud8R",2,2,1,1)
3297 -0.039573*
ewgc(
"CuGR",2,2)+0.006317*
ewgc(
"CuuR",0,2,2,0)+0.000171*
ewgc(
"CuuR",1,2,2,1)
3298 -0.000281*
ewgc(
"CuWR",2,2)
3299 +0.000076*
ewgc(
"CHD")+-0.007774*
ewgc(
"CHl3R",0,0)+-0.007774*
ewgc(
"CHl3R",1,1)+0.000013*
ewgc(
"CHW")
3300 +0.004172*
ewgc(
"CHWB")+0.007759*
ewgc(
"CllR",0,1,1,0)+-0.000123*
ewgc(
"Cqd1R",2,2,0,0)+-0.000006*
ewgc(
"Cqd1R",2,2,1,1)
3301 +-0.001720*
ewgc(
"Cqq1R",0,0,2,2)+-0.000435*
ewgc(
"Cqq1R",1,1,2,2)+0.007293*
ewgc(
"Cqq3R",0,0,2,2)
3302 +0.000589*
ewgc(
"Cqq3R",1,1,2,2)+-0.000047*
ewgc(
"Cqu1R",0,0,2,2)+-0.000039*
ewgc(
"Cqu1R",1,1,2,2)
3303 +0.000144*
ewgc(
"Cqu1R",2,2,0,0)+-0.000001*
ewgc(
"Cqu1R",2,2,1,1)+-0.000117*
ewgc(
"Cud1R",2,2,0,0)
3304 +-0.000012*
ewgc(
"Cud1R",2,2,1,1)+0.000456*
ewgc(
"CuuR",0,0,2,2)+0.000004*
ewgc(
"CuuR",1,1,2,2))
3305 *(SM_sigma_ttz_bin_70_110/sigma_ttz_bin_70_110_madgraph);
3307 if (total < 0)
return std::numeric_limits<double>::quiet_NaN();
3316 double sigma_ttz_bin_110_160_madgraph = 0.118917;
3326 double total = (sigma_ttz_bin_110_160_madgraph -0.032594*
ewgc(
"CG")-0.012662*
ewgc(
"CHq1R",2,2)
3327 +0.012643*
ewgc(
"CHq3R",2,2)+0.008236*
ewgc(
"CHuR",2,2)+0.002540*
ewgc(
"Cqd8R",2,2,0,0)
3328 +0.000253*
ewgc(
"Cqd8R",2,2,1,1)+0.017566*
ewgc(
"Cqq1R",0,2,2,0)+0.000417*
ewgc(
"Cqq1R",1,2,2,1)
3329 +0.071863*
ewgc(
"Cqq3R",0,2,2,0)+0.006279*
ewgc(
"Cqq3R",1,2,2,1)+0.009224*
ewgc(
"Cqu8R",0,0,2,2)
3330 +0.000682*
ewgc(
"Cqu8R",1,1,2,2)+0.002931*
ewgc(
"Cqu8R",2,2,0,0)+0.000049*
ewgc(
"Cqu8R",2,2,1,1)
3331 +-0.000017*
ewgc(
"CuBR",2,2)+0.001474*
ewgc(
"Cud8R",2,2,0,0)+0.000150*
ewgc(
"Cud8R",2,2,1,1)
3332 -0.038866*
ewgc(
"CuGR",2,2)+0.007314*
ewgc(
"CuuR",0,2,2,0)
3333 +0.000171*
ewgc(
"CuuR",1,2,2,1)+-0.000599*
ewgc(
"CuWR",2,2)
3334 +-0.000370*
ewgc(
"CHD")+-0.007261*
ewgc(
"CHl3R",0,0)+-0.007261*
ewgc(
"CHl3R",1,1)+-0.000003*
ewgc(
"CHW")
3335 +0.003479*
ewgc(
"CHWB")+0.007219*
ewgc(
"CllR",0,1,1,0)+-0.000093*
ewgc(
"Cqd1R",2,2,0,0)+-0.000028*
ewgc(
"Cqd1R",2,2,1,1)
3336 +-0.001164*
ewgc(
"Cqq1R",0,0,2,2)+-0.000361*
ewgc(
"Cqq1R",1,1,2,2)+0.006786*
ewgc(
"Cqq3R",0,0,2,2)
3337 +0.000480*
ewgc(
"Cqq3R",1,1,2,2)+0.000040*
ewgc(
"Cqu1R",0,0,2,2)+-0.000053*
ewgc(
"Cqu1R",1,1,2,2)
3338 +0.000170*
ewgc(
"Cqu1R",2,2,0,0)+-0.000016*
ewgc(
"Cqu1R",2,2,1,1)+-0.000090*
ewgc(
"Cud1R",2,2,0,0)
3339 +-0.000038*
ewgc(
"Cud1R",2,2,1,1)+0.000486*
ewgc(
"CuuR",0,0,2,2)+-0.000013*
ewgc(
"CuuR",1,1,2,2))
3340 *(SM_sigma_ttz_bin_110_160/sigma_ttz_bin_110_160_madgraph);
3342 if (total < 0)
return std::numeric_limits<double>::quiet_NaN();
3351 double sigma_ttz_bin_160_220_madgraph = 0.086169;
3361 double total = (sigma_ttz_bin_160_220_madgraph -0.028727*
ewgc(
"CG")+-0.009451*
ewgc(
"CHq1R",2,2)
3362 +0.009451*
ewgc(
"CHq3R",2,2)+0.006505*
ewgc(
"CHuR",2,2)+0.002109*
ewgc(
"Cqd8R",2,2,0,0)
3363 +0.000210*
ewgc(
"Cqd8R",2,2,1,1)+0.014508*
ewgc(
"Cqq1R",0,2,2,0)
3364 +0.000323*
ewgc(
"Cqq1R",1,2,2,1)+0.054025*
ewgc(
"Cqq3R",0,2,2,0)+0.004204*
ewgc(
"Cqq3R",1,2,2,1)
3365 +0.006853*
ewgc(
"Cqu8R",0,0,2,2)+0.000457*
ewgc(
"Cqu8R",1,1,2,2)+0.002647*
ewgc(
"Cqu8R",2,2,0,0)
3366 +0.000061*
ewgc(
"Cqu8R",2,2,1,1)+0.000070*
ewgc(
"CuBR",2,2)+0.001274*
ewgc(
"Cud8R",2,2,0,0)
3367 +0.000130*
ewgc(
"Cud8R",2,2,1,1)+-0.029532*
ewgc(
"CuGR",2,2)
3368 +0.006980*
ewgc(
"CuuR",0,2,2,0)+0.000163*
ewgc(
"CuuR",1,2,2,1)+-0.000537*
ewgc(
"CuWR",2,2)
3369 +-0.000506*
ewgc(
"CHD")+-0.005252*
ewgc(
"CHl3R",0,0)+-0.005252*
ewgc(
"CHl3R",1,1)+0.000016*
ewgc(
"CHW")
3370 +0.002257*
ewgc(
"CHWB")+0.005257*
ewgc(
"CllR",0,1,1,0)+-0.000081*
ewgc(
"Cqd1R",2,2,0,0)+-0.000017*
ewgc(
"Cqd1R",2,2,1,1)
3371 +-0.000799*
ewgc(
"Cqq1R",0,0,2,2)+-0.000244*
ewgc(
"Cqq1R",1,1,2,2)+0.005197*
ewgc(
"Cqq3R",0,0,2,2)
3372 +0.000326*
ewgc(
"Cqq3R",1,1,2,2)+0.000178*
ewgc(
"Cqu1R",0,0,2,2)+-0.000014*
ewgc(
"Cqu1R",1,1,2,2)
3373 +0.000180*
ewgc(
"Cqu1R",2,2,0,0)+0.000002*
ewgc(
"Cqu1R",2,2,1,1)+-0.000112*
ewgc(
"Cud1R",2,2,0,0)
3374 +-0.000009*
ewgc(
"Cud1R",2,2,1,1)+0.000539*
ewgc(
"CuuR",0,0,2,2)+0.000010*
ewgc(
"CuuR",1,1,2,2))
3375 *(SM_sigma_ttz_bin_160_220/sigma_ttz_bin_160_220_madgraph);
3377 if (total < 0)
return std::numeric_limits<double>::quiet_NaN();
3386 double sigma_ttz_bin_220_290_madgraph = 0.051619;
3396 double total = (sigma_ttz_bin_220_290_madgraph -0.020633*
ewgc(
"CG")+-0.005708*
ewgc(
"CHq1R",2,2)
3397 +0.005730*
ewgc(
"CHq3R",2,2)+0.004111*
ewgc(
"CHuR",2,2)+0.001568*
ewgc(
"Cqd8R",2,2,0,0)
3398 +0.000156*
ewgc(
"Cqd8R",2,2,1,1)+0.010800*
ewgc(
"Cqq1R",0,2,2,0)
3399 +0.000238*
ewgc(
"Cqq1R",1,2,2,1)+0.036604*
ewgc(
"Cqq3R",0,2,2,0)+0.002537*
ewgc(
"Cqq3R",1,2,2,1)
3400 +0.004688*
ewgc(
"Cqu8R",0,0,2,2)+0.000287*
ewgc(
"Cqu8R",1,1,2,2)+0.002098*
ewgc(
"Cqu8R",2,2,0,0)
3401 +0.000060*
ewgc(
"Cqu8R",2,2,1,1)+0.000066*
ewgc(
"CuBR",2,2)+0.000986*
ewgc(
"Cud8R",2,2,0,0)
3402 +0.000107*
ewgc(
"Cud8R",2,2,1,1)-0.018706*
ewgc(
"CuGR",2,2)
3403 +0.005776*
ewgc(
"CuuR",0,2,2,0)+0.000142*
ewgc(
"CuuR",1,2,2,1)-0.000356*
ewgc(
"CuWR",2,2)
3404 +-0.000396*
ewgc(
"CHD")+-0.003125*
ewgc(
"CHl3R",0,0)+-0.003125*
ewgc(
"CHl3R",1,1)+0.000030*
ewgc(
"CHW")
3405 +0.001241*
ewgc(
"CHWB")+0.003151*
ewgc(
"CllR",0,1,1,0)+-0.000024*
ewgc(
"Cqd1R",2,2,0,0)+0.000006*
ewgc(
"Cqd1R",2,2,1,1)
3406 +-0.000280*
ewgc(
"Cqq1R",0,0,2,2)+-0.000122*
ewgc(
"Cqq1R",1,1,2,2)+0.003698*
ewgc(
"Cqq3R",0,0,2,2)
3407 +0.000221*
ewgc(
"Cqq3R",1,1,2,2)+0.000151*
ewgc(
"Cqu1R",0,0,2,2)+0.000009*
ewgc(
"Cqu1R",1,1,2,2)
3408 +0.000165*
ewgc(
"Cqu1R",2,2,0,0)+0.000020*
ewgc(
"Cqu1R",2,2,1,1)+-0.000062*
ewgc(
"Cud1R",2,2,0,0)
3409 +0.000012*
ewgc(
"Cud1R",2,2,1,1)+0.000461*
ewgc(
"CuuR",0,0,2,2)+0.000027*
ewgc(
"CuuR",1,1,2,2))
3410 *(SM_sigma_ttz_bin_220_290/sigma_ttz_bin_220_290_madgraph);
3412 if (total < 0)
return std::numeric_limits<double>::quiet_NaN();
3421 double sigma_ttz_bin_290_400_madgraph = 0.032507;
3431 double total = (sigma_ttz_bin_290_400_madgraph -0.015551*
ewgc(
"CG")+-0.003576*
ewgc(
"CHq1R",2,2)
3432 +0.003588*
ewgc(
"CHq3R",2,2)+0.002640*
ewgc(
"CHuR",2,2)+0.001358*
ewgc(
"Cqd8R",2,2,0,0)
3433 +0.000110*
ewgc(
"Cqd8R",2,2,1,1)+0.009422*
ewgc(
"Cqq1R",0,2,2,0)
3434 +0.000177*
ewgc(
"Cqq1R",1,2,2,1)+0.029289*
ewgc(
"Cqq3R",0,2,2,0)+0.001743*
ewgc(
"Cqq3R",1,2,2,1)
3435 +0.003812*
ewgc(
"Cqu8R",0,0,2,2)+0.000192*
ewgc(
"Cqu8R",1,1,2,2)+0.001922*
ewgc(
"Cqu8R",2,2,0,0)
3436 +0.000038*
ewgc(
"Cqu8R",2,2,1,1)+0.000016*
ewgc(
"CuBR",2,2)+0.000882*
ewgc(
"Cud8R",2,2,0,0)
3437 +0.000076*
ewgc(
"Cud8R",2,2,1,1)-0.012764*
ewgc(
"CuGR",2,2)
3438 +0.005527*
ewgc(
"CuuR",0,2,2,0)+0.000108*
ewgc(
"CuuR",1,2,2,1)+-0.000264*
ewgc(
"CuWR",2,2)
3439 +-0.000295*
ewgc(
"CHD")+-0.001978*
ewgc(
"CHl3R",0,0)+-0.001978*
ewgc(
"CHl3R",1,1)+0.000020*
ewgc(
"CHW")
3440 +0.000740*
ewgc(
"CHWB")+0.001985*
ewgc(
"CllR",0,1,1,0)+-0.000037*
ewgc(
"Cqd1R",2,2,0,0)+-0.000001*
ewgc(
"Cqd1R",2,2,1,1)
3441 +-0.000103*
ewgc(
"Cqq1R",0,0,2,2)+-0.000084*
ewgc(
"Cqq1R",1,1,2,2)+0.003015*
ewgc(
"Cqq3R",0,0,2,2)
3442 +0.000159*
ewgc(
"Cqq3R",1,1,2,2)+0.000151*
ewgc(
"Cqu1R",0,0,2,2)+-0.000000*
ewgc(
"Cqu1R",1,1,2,2)
3443 +0.000122*
ewgc(
"Cqu1R",2,2,0,0)+0.000006*
ewgc(
"Cqu1R",2,2,1,1)+-0.000060*
ewgc(
"Cud1R",2,2,0,0)
3444 +-0.000001*
ewgc(
"Cud1R",2,2,1,1)+0.000432*
ewgc(
"CuuR",0,0,2,2)+0.000012*
ewgc(
"CuuR",1,1,2,2))
3445 *(SM_sigma_ttz_bin_290_400/sigma_ttz_bin_290_400_madgraph);
3447 if (total < 0)
return std::numeric_limits<double>::quiet_NaN();
3454 else throw std::runtime_error(
"wrong bin choice in sigma_ttz_diff_LO_ATLAS_210312603");
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