4395{
4396
4397 bool flag_Quadratic= false;
4398
4399 double sigma_ttwp_LO_madgraph = 0.349605;
4400 double sigma_ttwm_LO_madgraph = 0.166990;
4403
4404
4405 double sigma_ttwp_SM = sigma_ttw_SM/(1+1/R_ttw_SM);
4406 double sigma_ttwm_SM = sigma_ttw_SM/(1+R_ttw_SM);
4407
4408
4409
4410
4411 if(flag_Quadratic){
4412 return 0.;
4413 }
4414 else{
4415
4416 double sigma_ttwp_NP = -0.000229*
ewgc(
"CHq1R",2,2)+0.000409*
ewgc(
"CHq3R",2,2)+-0.000105*
ewgc(
"CHuR",2,2)
4417 +0.131627*
ewgc(
"Cqq1R",0,2,2,0)+0.004173*
ewgc(
"Cqq1R",1,2,2,1)+0.622300*
ewgc(
"Cqq3R",0,2,2,0)
4418 +0.029216*
ewgc(
"Cqq3R",1,2,2,1)+0.095061*
ewgc(
"Cqu8R",0,0,2,2)+0.004239*
ewgc(
"Cqu8R",1,1,2,2)
4419 +0.000066*
ewgc(
"CuBR",2,2)-0.085559*
ewgc(
"CuGR",2,2)-0.004044*
ewgc(
"CuWR",2,2)
4420 +-0.000229*
ewgc(
"CHq1R",2,2)+0.000409*
ewgc(
"CHq3R",2,2)+-0.000105*
ewgc(
"CHuR",2,2)+0.131627*
ewgc(
"Cqq1R",0,2,2,0)
4421 +0.004173*
ewgc(
"Cqq1R",1,2,2,1)+0.622300*
ewgc(
"Cqq3R",0,2,2,0)+0.029216*
ewgc(
"Cqq3R",1,2,2,1)
4422 +0.095061*
ewgc(
"Cqu8R",0,0,2,2)+0.004239*
ewgc(
"Cqu8R",1,1,2,2)+0.000066*
ewgc(
"CuBR",2,2)
4423 +-0.085559*
ewgc(
"CuGR",2,2)+-0.004044*
ewgc(
"CuWR",2,2)
4424 +-0.021513*
ewgc(
"CHl3R",0,0)+-0.021513*
ewgc(
"CHl3R",1,1)+0.000171*
ewgc(
"CHW")+-0.000141*
ewgc(
"CHWB")
4425 +0.021524*
ewgc(
"CllR",0,1,1,0)+-0.007924*
ewgc(
"Cqq1R",0,0,2,2)+-0.000761*
ewgc(
"Cqq1R",1,1,2,2)
4426 +0.068930*
ewgc(
"Cqq3R",0,0,2,2)+0.002428*
ewgc(
"Cqq3R",1,1,2,2)+0.003035*
ewgc(
"Cqu1R",0,0,2,2)
4427 +-0.000034*
ewgc(
"Cqu1R",1,1,2,2);
4428 double sigma_ttwm_NP = -0.000104*
ewgc(
"CHq1R",2,2)+0.000188*
ewgc(
"CHq3R",2,2)+-0.000050*
ewgc(
"CHuR",2,2)
4429 +0.051853*
ewgc(
"Cqq1R",0,2,2,0)+0.006487*
ewgc(
"Cqq1R",1,2,2,1)+0.251097*
ewgc(
"Cqq3R",0,2,2,0)
4430 +0.025657*
ewgc(
"Cqq3R",1,2,2,1)+0.038346*
ewgc(
"Cqu8R",0,0,2,2)+0.004035*
ewgc(
"Cqu8R",1,1,2,2)
4431 +0.000030*
ewgc(
"CuBR",2,2)-0.040504*
ewgc(
"CuGR",2,2)-0.001812*
ewgc(
"CuWR",2,2)
4432 +-0.000104*
ewgc(
"CHq1R",2,2)+0.000188*
ewgc(
"CHq3R",2,2)+-0.000050*
ewgc(
"CHuR",2,2)+0.051853*
ewgc(
"Cqq1R",0,2,2,0)
4433 +0.006487*
ewgc(
"Cqq1R",1,2,2,1)+0.251097*
ewgc(
"Cqq3R",0,2,2,0)+0.025657*
ewgc(
"Cqq3R",1,2,2,1)
4434 +0.038346*
ewgc(
"Cqu8R",0,0,2,2)+0.004035*
ewgc(
"Cqu8R",1,1,2,2)+0.000030*
ewgc(
"CuBR",2,2)
4435 +-0.040504*
ewgc(
"CuGR",2,2)+-0.001812*
ewgc(
"CuWR",2,2)
4436 +-0.010266*
ewgc(
"CHl3R",0,0)+-0.010266*
ewgc(
"CHl3R",1,1)+0.000075*
ewgc(
"CHW")+-0.000066*
ewgc(
"CHWB")
4437 +0.010267*
ewgc(
"CllR",0,1,1,0)+-0.003590*
ewgc(
"Cqq1R",0,0,2,2)+-0.000254*
ewgc(
"Cqq1R",1,1,2,2)
4438 +0.025859*
ewgc(
"Cqq3R",0,0,2,2)+0.002285*
ewgc(
"Cqq3R",1,1,2,2)+0.000816*
ewgc(
"Cqu1R",0,0,2,2)
4439 +0.000098*
ewgc(
"Cqu1R",1,1,2,2);
4440
4441 double sigma_ttwp_NP_Corrected = sigma_ttwp_NP*(sigma_ttwp_SM/sigma_ttwp_LO_madgraph);
4442 double sigma_ttwm_NP_Corrected = sigma_ttwm_NP*(sigma_ttwm_SM/sigma_ttwm_LO_madgraph);
4443
4444 double total = R_ttw_SM + (sigma_ttwp_NP_Corrected - sigma_ttwm_NP_Corrected*sigma_ttwp_SM/sigma_ttwm_SM)/sigma_ttwm_SM ;
4445
4446 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
4447
4448 return total;
4449
4450 }
4451
4452}
const double getOptionalParameter(std::string name) const
A method to get parameters that are specific to only one set of observables.
double ewgc(const std::string name) const