a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
Public Member Functions | Private Attributes | List of all members
charge_asymmetry_tt_diff_mtt_NLO Class Reference

#include <NPSMEFT6dtopquark.h>

+ Inheritance diagram for charge_asymmetry_tt_diff_mtt_NLO:

Detailed Description

Definition at line 1302 of file NPSMEFT6dtopquark.h.

Public Member Functions

 charge_asymmetry_tt_diff_mtt_NLO (const StandardModel &SM_i)
 charge_asymmetry_tt_diff_mtt_LO constructor. More...
 
double computeThValue ()
 
- Public Member Functions inherited from ThObservable
double getBinMax ()
 A get method to get the maximum value of the bin. More...
 
double getBinMin ()
 A get method to get the minimum value of the bin. More...
 
const StandardModel & getModel ()
 A get method to get the model. More...
 
const std::vector< std::string > getParametersForObservable ()
 A get method to get the parameters for the specific observable. More...
 
void setBinMax (double max)
 A set method to set the maximum value of the bin. More...
 
void setBinMin (double min)
 A set method to set the minimum value of the bin. More...
 
void setParametersForObservable (std::vector< std::string > parametersForObservable_i)
 A set method to get the parameters for the specific observable. More...
 
 ThObservable (const StandardModel &SM_i)
 Constructor. More...
 
 ThObservable (const ThObservable &orig)
 The copy constructor. More...
 
virtual ~ThObservable ()
 The default destructor. More...
 

Private Attributes

double b_max = 0.
 
double b_min = 0.
 
const NPSMEFT6dtopquark & myNPSMEFT6dtopquark
 

Additional Inherited Members

- Protected Attributes inherited from ThObservable
double max
 the bin maximum. More...
 
double min
 The bin minimum. More...
 
std::vector< std::string > parametersForObservable
 a vector of parameter namesfor the specific observable More...
 
const StandardModel & SM
 A reference to an object of StandardMode class. More...
 

Constructor & Destructor Documentation

◆ charge_asymmetry_tt_diff_mtt_NLO()

charge_asymmetry_tt_diff_mtt_NLO::charge_asymmetry_tt_diff_mtt_NLO ( const StandardModel &  SM_i)

charge_asymmetry_tt_diff_mtt_LO constructor.

Definition at line 4334 of file NPSMEFT6dtopquark.cpp.

4335: ThObservable(SM_i), myNPSMEFT6dtopquark(static_cast<const NPSMEFT6dtopquark&> (SM_i))
4336{
4337
4338 setParametersForObservable(make_vector<std::string>() << "SM_charge_asymmetry_bin_mtt_0_500" << "SM_charge_asymmetry_deno_bin_mtt_0_500"
4339 << "SM_charge_asymmetry_bin_mtt_500_750" << "SM_charge_asymmetry_deno_bin_mtt_500_750" << "SM_charge_asymmetry_bin_mtt_750_1000"
4340 << "SM_charge_asymmetry_deno_bin_mtt_750_1000" << "SM_charge_asymmetry_bin_mtt_1000_1500" << "SM_charge_asymmetry_deno_bin_mtt_1000_1500"
4341 << "SM_charge_asymmetry_bin_mtt_1500_3000" << "SM_charge_asymmetry_deno_bin_mtt_1500_3000" << "SM_charge_asymmetry_bin_mtt_1500_2000"
4342 << "SM_charge_asymmetry_deno_bin_mtt_1500_2000" << "SM_charge_asymmetry_bin_mtt_2000_2500" << "SM_charge_asymmetry_deno_bin_mtt_2000_2500"
4343 << "SM_charge_asymmetry_bin_mtt_2500_3000" << "SM_charge_asymmetry_deno_bin_mtt_2500_3000");
4344
4345}
ThObservable::setParametersForObservable
void setParametersForObservable(std::vector< std::string > parametersForObservable_i)
A set method to get the parameters for the specific observable.
Definition: ThObservable.h:109
ThObservable::ThObservable
ThObservable(const StandardModel &SM_i)
Constructor.
Definition: ThObservable.h:32
charge_asymmetry_tt_diff_mtt_NLO::myNPSMEFT6dtopquark
const NPSMEFT6dtopquark & myNPSMEFT6dtopquark
Definition: NPSMEFT6dtopquark.h:1323

Member Function Documentation

◆ computeThValue()

double charge_asymmetry_tt_diff_mtt_NLO::computeThValue ( )
virtual
Returns
The value of charge_asymmetry_tt_diff_mtt_LO

Implements ThObservable.

Definition at line 4347 of file NPSMEFT6dtopquark.cpp.

4348{
4349
4350 b_min = getBinMin();
4351 b_max = getBinMax();
4352
4353 double C_tG = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_tG();
4354 double C_Qd8 = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_Qd8();
4355 double C_Qd1 = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_Qd1();
4356 double C_Qu1 = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_Qu1();
4357 double C_Qu8 = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_Qu8();
4358 double C_td1 = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_td1();
4359 double C_td8 = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_td8();
4360 double C_tq1 = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_tq1();
4361 double C_tq8 = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_tq8();
4362 double C_tu1 = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_tu1();
4363 double C_tu8 = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_tu8();
4364 double C_Qq11 = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_Qq11();
4365 double C_Qq31 = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_Qq31();
4366 double C_Qq18 = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_Qq18();
4367 double C_Qq38 = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_Qq38();
4368 bool flag_Quadratic = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_flag_Quadratic();
4369
4370 //double SM_Charge_Asymmetry_bin_tt_0_500 = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_SM_Charge_Asymmetry_bin_tt_0_500();
4371
4372
4373 if(b_min == 0 && b_max == 500){
4374
4375 double SM_charge_asymmetry_bin_mtt_0_500 = SM.getOptionalParameter("SM_charge_asymmetry_bin_mtt_0_500");
4376 double SM_charge_asymmetry_deno_bin_mtt_0_500 = SM.getOptionalParameter("SM_charge_asymmetry_deno_bin_mtt_0_500");
4377 double SM_charge_asymmetry_num_bin_mtt_0_500 = SM_charge_asymmetry_bin_mtt_0_500*SM_charge_asymmetry_deno_bin_mtt_0_500;
4378
4379 double SM_sigma_pos_bin_mtt_0_500 =0.5*(SM_charge_asymmetry_num_bin_mtt_0_500+SM_charge_asymmetry_deno_bin_mtt_0_500);
4380 double SM_sigma_neg_bin_mtt_0_500 =0.5*(SM_charge_asymmetry_deno_bin_mtt_0_500-SM_charge_asymmetry_num_bin_mtt_0_500);
4381
4382
4383 double sigma_pos_bin_mtt_0_500_madgraph_NLO = 206522.71; //fb
4384 double sigma_neg_bin_mtt_0_500_madgraph_NLO = 205736.11; //fb
4385
4386
4387 if(flag_Quadratic){
4388
4389
4390 double sigma_pos_bin_mtt_0_500_NP_quad = 10856.61*C_tG*C_tG + 186.94*C_Qd1*C_Qd1 + 81.37*C_Qd8*C_Qd8 + 655.08*C_Qq11*C_Qq11 +
4391 712.66*C_Qq31*C_Qq31 + 589.45*C_Qq18*C_Qq18 + 121.35*C_Qq38*C_Qq38 + 198.11*C_Qu1*C_Qu1 + 83.7*C_Qu8*C_Qu8 + 244.73*C_td1*C_td1 +
4392 -27.74*C_td8*C_td8 + 499.26*C_tq1*C_tq1 + 35.53*C_tq8*C_tq8 + 380.22*C_tu1*C_tu1 + 128.01*C_tu8*C_tu8;
4393 double sigma_neg_bin_mtt_0_500_NP_quad = 10844.05*C_tG*C_tG + 124.97*C_Qd1*C_Qd1 + -57.71*C_Qd8*C_Qd8 + 415.01*C_Qq11*C_Qq11 +
4394 348.52*C_Qq31*C_Qq31 + 575.55*C_Qq18*C_Qq18 + 148.51*C_Qq38*C_Qq38 + 373.47*C_Qu1*C_Qu1 + 31.75*C_Qu8*C_Qu8 + 258.97*C_td1*C_td1 +
4395 1.84*C_td8*C_td8 + 602.71*C_tq1*C_tq1 + 100.72*C_tq8*C_tq8 + 275.67*C_tu1*C_tu1 + -42.82*C_tu8*C_tu8;
4396
4397
4398 double sigma_pos_bin_mtt_0_500_NP_lin = 69255.55*C_tG + 86.03*C_Qd1 + 631.63*C_Qd8 + -171.25*C_Qq11 + 11.28*C_Qq31 + 977.78*C_Qq18 +
4399 383.52*C_Qq38 + 191.56*C_Qu1 + 1018.38*C_Qu8 + -191.45*C_td1 + 677.13*C_td8 + 273.38*C_tq1 + 1321.58*C_tq8 + -7.5*C_tu1 + 1119.97*C_tu8;
4400 double sigma_neg_bin_mtt_0_500_NP_lin = 68845.04*C_tG + 203.17*C_Qd1 + 727.79*C_Qd8 + -445.67*C_Qq11 + -122.59*C_Qq31 + 423.69*C_Qq18 +
4401 72.21*C_Qq38 + -65.45*C_Qu1 + 760.75*C_Qu8 + -267.01*C_td1 + 508.03*C_td8 + 176.2*C_tq1 + 1589.47*C_tq8 + -336.04*C_tu1 + 610.15*C_tu8;
4402
4403
4404 double sigma_pos_bin_mtt_0_500_NP = sigma_pos_bin_mtt_0_500_NP_lin + sigma_pos_bin_mtt_0_500_NP_quad;
4405 double sigma_neg_bin_mtt_0_500_NP = sigma_neg_bin_mtt_0_500_NP_lin + sigma_neg_bin_mtt_0_500_NP_quad;
4406
4407
4408 double sigma_pos_bin_mtt_0_500_NP_Corrected = SM_sigma_pos_bin_mtt_0_500*sigma_pos_bin_mtt_0_500_NP/sigma_pos_bin_mtt_0_500_madgraph_NLO;
4409 double sigma_neg_bin_mtt_0_500_NP_Corrected = SM_sigma_neg_bin_mtt_0_500*sigma_neg_bin_mtt_0_500_NP/sigma_neg_bin_mtt_0_500_madgraph_NLO;
4410
4411 double NP_charge_asymmetry_num_bin_mtt_0_500 = sigma_pos_bin_mtt_0_500_NP_Corrected - sigma_neg_bin_mtt_0_500_NP_Corrected;
4412 double NP_charge_asymmetry_deno_bin_mtt_0_500 = sigma_pos_bin_mtt_0_500_NP_Corrected + sigma_neg_bin_mtt_0_500_NP_Corrected;
4413
4414 //return SM_charge_asymmetry_bin_mtt_0_500*(1+(NP_charge_asymmetry_num_bin_mtt_0_500-NP_charge_asymmetry_deno_bin_mtt_0_500)/SM_charge_asymmetry_deno_bin_mtt_0_500);
4415
4416
4417 return SM_charge_asymmetry_bin_mtt_0_500 + (NP_charge_asymmetry_num_bin_mtt_0_500
4418 - SM_charge_asymmetry_bin_mtt_0_500*NP_charge_asymmetry_deno_bin_mtt_0_500
4419 )/SM_charge_asymmetry_deno_bin_mtt_0_500;
4420
4421 }
4422 else{
4423
4424
4425 double sigma_pos_bin_mtt_0_500_NP = 69255.55*C_tG + 86.03*C_Qd1 + 631.63*C_Qd8 + -171.25*C_Qq11 + 11.28*C_Qq31 + 977.78*C_Qq18 +
4426 383.52*C_Qq38 + 191.56*C_Qu1 + 1018.38*C_Qu8 + -191.45*C_td1 + 677.13*C_td8 + 273.38*C_tq1 + 1321.58*C_tq8 + -7.5*C_tu1 + 1119.97*C_tu8;
4427 double sigma_neg_bin_mtt_0_500_NP = 68845.04*C_tG + 203.17*C_Qd1 + 727.79*C_Qd8 + -445.67*C_Qq11 + -122.59*C_Qq31 + 423.69*C_Qq18 +
4428 72.21*C_Qq38 + -65.45*C_Qu1 + 760.75*C_Qu8 + -267.01*C_td1 + 508.03*C_td8 + 176.2*C_tq1 + 1589.47*C_tq8 + -336.04*C_tu1 + 610.15*C_tu8;
4429
4430;
4431
4432 double sigma_pos_bin_mtt_0_500_NP_Corrected = SM_sigma_pos_bin_mtt_0_500*sigma_pos_bin_mtt_0_500_NP/sigma_pos_bin_mtt_0_500_madgraph_NLO;
4433 double sigma_neg_bin_mtt_0_500_NP_Corrected = SM_sigma_neg_bin_mtt_0_500*sigma_neg_bin_mtt_0_500_NP/sigma_neg_bin_mtt_0_500_madgraph_NLO;
4434
4435 double NP_charge_asymmetry_num_bin_mtt_0_500 = sigma_pos_bin_mtt_0_500_NP_Corrected - sigma_neg_bin_mtt_0_500_NP_Corrected;
4436 double NP_charge_asymmetry_deno_bin_mtt_0_500 = sigma_pos_bin_mtt_0_500_NP_Corrected + sigma_neg_bin_mtt_0_500_NP_Corrected;
4437
4438 //return SM_charge_asymmetry_bin_mtt_0_500*(1+(NP_charge_asymmetry_num_bin_mtt_0_500-NP_charge_asymmetry_deno_bin_mtt_0_500)/SM_charge_asymmetry_deno_bin_mtt_0_500);
4439
4440
4441 return SM_charge_asymmetry_bin_mtt_0_500 + (NP_charge_asymmetry_num_bin_mtt_0_500
4442 - SM_charge_asymmetry_bin_mtt_0_500*NP_charge_asymmetry_deno_bin_mtt_0_500
4443 )/SM_charge_asymmetry_deno_bin_mtt_0_500;
4444
4445
4446 }
4447
4448 } else if(b_min == 500 && b_max == 750){
4449
4450 double SM_charge_asymmetry_bin_mtt_500_750 = SM.getOptionalParameter("SM_charge_asymmetry_bin_mtt_500_750");
4451 double SM_charge_asymmetry_deno_bin_mtt_500_750 = SM.getOptionalParameter("SM_charge_asymmetry_deno_bin_mtt_500_750");
4452 double SM_charge_asymmetry_num_bin_mtt_500_750 = SM_charge_asymmetry_bin_mtt_500_750*SM_charge_asymmetry_deno_bin_mtt_500_750;
4453
4454
4455 double SM_sigma_pos_bin_mtt_500_750 =0.5*(SM_charge_asymmetry_num_bin_mtt_500_750+SM_charge_asymmetry_deno_bin_mtt_500_750);
4456 double SM_sigma_neg_bin_mtt_500_750 =0.5*(SM_charge_asymmetry_deno_bin_mtt_500_750-SM_charge_asymmetry_num_bin_mtt_500_750);
4457
4458
4459 double sigma_pos_bin_mtt_500_750_madgraph_NLO = 98972.12;
4460 double sigma_neg_bin_mtt_500_750_madgraph_NLO = 98471.35;
4461
4462 if(flag_Quadratic){
4463
4464
4465 double sigma_pos_bin_mtt_500_750_NP_quad = 5493.57*C_tG*C_tG + 315.65*C_Qd1*C_Qd1 + 52.9*C_Qd8*C_Qd8 + 1274.19*C_Qq11*C_Qq11 +
4466 1247.82*C_Qq31*C_Qq31 + 234.29*C_Qq18*C_Qq18 + 220.5*C_Qq38*C_Qq38 + 375.02*C_Qu1*C_Qu1 + 81.75*C_Qu8*C_Qu8 + 449.86*C_td1*C_td1 +
4467 67.51*C_td8*C_td8 + 652.68*C_tq1*C_tq1 + 147.74*C_tq8*C_tq8 + 776.43*C_tu1*C_tu1 + 140.43*C_tu8*C_tu8;
4468 double sigma_neg_bin_mtt_500_750_NP_quad = 5476.01*C_tG*C_tG + 419.52*C_Qd1*C_Qd1 + 77.02*C_Qd8*C_Qd8 + 616.75*C_Qq11*C_Qq11 +
4469 617.3*C_Qq31*C_Qq31 + 141.47*C_Qq18*C_Qq18 + 76.31*C_Qq38*C_Qq38 + 751.03*C_Qu1*C_Qu1 + 118.22*C_Qu8*C_Qu8 + 276.88*C_td1*C_td1 +
4470 24.92*C_td8*C_td8 + 1237.5*C_tq1*C_tq1 + 200.26*C_tq8*C_tq8 + 359.71*C_tu1*C_tu1 + 52.72*C_tu8*C_tu8;
4471
4472
4473
4474 double sigma_pos_bin_mtt_500_750_NP_lin = 28985.27*C_tG + 32.74*C_Qd1 + 382.94*C_Qd8 + -108.77*C_Qq11 + -48.35*C_Qq31 + 1456.01*C_Qq18 +
4475 341.82*C_Qq38 + 115.28*C_Qu1 + 464.61*C_Qu8 + -46.51*C_td1 + 563.4*C_td8 + 103.98*C_tq1 + 886.74*C_tq8 + -18.49*C_tu1 + 904.87*C_tu8;
4476 double sigma_neg_bin_mtt_500_750_NP_lin = 28804.73*C_tG + 48.33*C_Qd1 + 591.25*C_Qd8 + -189.14*C_Qq11 + -21.89*C_Qq31 + 737.03*C_Qq18 +
4477 69.36*C_Qq38 + -81.18*C_Qu1 + 932.13*C_Qu8 + -82.69*C_td1 + 292.91*C_td8 + -53.01*C_tq1 + 1450.11*C_tq8 + -87.3*C_tu1 + 442.24*C_tu8;
4478
4479
4480
4481 double sigma_pos_bin_mtt_500_750_NP = sigma_pos_bin_mtt_500_750_NP_lin + sigma_pos_bin_mtt_500_750_NP_quad;
4482 double sigma_neg_bin_mtt_500_750_NP = sigma_neg_bin_mtt_500_750_NP_lin + sigma_neg_bin_mtt_500_750_NP_quad;
4483
4484
4485 double sigma_pos_bin_mtt_500_750_NP_Corrected = SM_sigma_pos_bin_mtt_500_750*sigma_pos_bin_mtt_500_750_NP/sigma_pos_bin_mtt_500_750_madgraph_NLO;
4486 double sigma_neg_bin_mtt_500_750_NP_Corrected = SM_sigma_neg_bin_mtt_500_750*sigma_neg_bin_mtt_500_750_NP/sigma_neg_bin_mtt_500_750_madgraph_NLO;
4487
4488 double NP_charge_asymmetry_num_bin_mtt_500_750 = sigma_pos_bin_mtt_500_750_NP_Corrected - sigma_neg_bin_mtt_500_750_NP_Corrected;
4489 double NP_charge_asymmetry_deno_bin_mtt_500_750 = sigma_pos_bin_mtt_500_750_NP_Corrected + sigma_neg_bin_mtt_500_750_NP_Corrected;
4490
4491 //return SM_charge_asymmetry_bin_mtt_500_750*(1+(NP_charge_asymmetry_num_bin_mtt_500_750-NP_charge_asymmetry_deno_bin_mtt_500_750)/SM_charge_asymmetry_deno_bin_mtt_500_750);
4492
4493
4494 return SM_charge_asymmetry_bin_mtt_500_750 + (NP_charge_asymmetry_num_bin_mtt_500_750
4495 - SM_charge_asymmetry_bin_mtt_500_750*NP_charge_asymmetry_deno_bin_mtt_500_750
4496 )/SM_charge_asymmetry_deno_bin_mtt_500_750;
4497
4498
4499 }
4500 else{
4501
4502
4503 double sigma_pos_bin_mtt_500_750_NP = 28985.27*C_tG + 32.74*C_Qd1 + 382.94*C_Qd8 + -108.77*C_Qq11 + -48.35*C_Qq31 + 1456.01*C_Qq18 +
4504 341.82*C_Qq38 + 115.28*C_Qu1 + 464.61*C_Qu8 + -46.51*C_td1 + 563.4*C_td8 + 103.98*C_tq1 + 886.74*C_tq8 + -18.49*C_tu1 + 904.87*C_tu8;
4505 double sigma_neg_bin_mtt_500_750_NP = 28804.73*C_tG + 48.33*C_Qd1 + 591.25*C_Qd8 + -189.14*C_Qq11 + -21.89*C_Qq31 + 737.03*C_Qq18 +
4506 69.36*C_Qq38 + -81.18*C_Qu1 + 932.13*C_Qu8 + -82.69*C_td1 + 292.91*C_td8 + -53.01*C_tq1 + 1450.11*C_tq8 + -87.3*C_tu1 + 442.24*C_tu8;
4507
4508
4509 double sigma_pos_bin_mtt_500_750_NP_Corrected = SM_sigma_pos_bin_mtt_500_750*sigma_pos_bin_mtt_500_750_NP/sigma_pos_bin_mtt_500_750_madgraph_NLO;
4510 double sigma_neg_bin_mtt_500_750_NP_Corrected = SM_sigma_neg_bin_mtt_500_750*sigma_neg_bin_mtt_500_750_NP/sigma_neg_bin_mtt_500_750_madgraph_NLO;
4511
4512 double NP_charge_asymmetry_num_bin_mtt_500_750 = sigma_pos_bin_mtt_500_750_NP_Corrected - sigma_neg_bin_mtt_500_750_NP_Corrected;
4513 double NP_charge_asymmetry_deno_bin_mtt_500_750 = sigma_pos_bin_mtt_500_750_NP_Corrected + sigma_neg_bin_mtt_500_750_NP_Corrected;
4514
4515 //return SM_charge_asymmetry_bin_mtt_500_750*(1+(NP_charge_asymmetry_num_bin_mtt_500_750-NP_charge_asymmetry_deno_bin_mtt_500_750)/SM_charge_asymmetry_deno_bin_mtt_500_750);
4516
4517
4518 return SM_charge_asymmetry_bin_mtt_500_750 + (NP_charge_asymmetry_num_bin_mtt_500_750
4519 - SM_charge_asymmetry_bin_mtt_500_750*NP_charge_asymmetry_deno_bin_mtt_500_750
4520 )/SM_charge_asymmetry_deno_bin_mtt_500_750;
4521
4522
4523
4524 }
4525
4526 } else if(b_min == 750 && b_max == 1000){
4527
4528
4529 double SM_charge_asymmetry_deno_bin_mtt_750_1000 = SM.getOptionalParameter("SM_charge_asymmetry_deno_bin_mtt_750_1000");
4530 double SM_charge_asymmetry_bin_mtt_750_1000 = SM.getOptionalParameter("SM_charge_asymmetry_bin_mtt_750_1000");
4531 double SM_charge_asymmetry_num_bin_mtt_750_1000 = SM_charge_asymmetry_bin_mtt_750_1000*SM_charge_asymmetry_deno_bin_mtt_750_1000;
4532
4533 double SM_sigma_pos_bin_mtt_750_1000 =0.5*(SM_charge_asymmetry_num_bin_mtt_750_1000+SM_charge_asymmetry_deno_bin_mtt_750_1000);
4534 double SM_sigma_neg_bin_mtt_750_1000 =0.5*(SM_charge_asymmetry_deno_bin_mtt_750_1000-SM_charge_asymmetry_num_bin_mtt_750_1000);
4535
4536 double sigma_pos_bin_mtt_750_1000_madgraph_NLO = 18609.62;
4537 double sigma_neg_bin_mtt_750_1000_madgraph_NLO = 18459.50;
4538
4539
4540 if(flag_Quadratic){
4541
4542
4543 double sigma_pos_bin_mtt_750_1000_NP_quad = 1339.65*C_tG*C_tG + 210.5*C_Qd1*C_Qd1 + 43.98*C_Qd8*C_Qd8 + 1017.96*C_Qq11*C_Qq11 +
4544 1008.2*C_Qq31*C_Qq31 + 189.09*C_Qq18*C_Qq18 + 201.45*C_Qq38*C_Qq38 + 275.07*C_Qu1*C_Qu1 + 55.95*C_Qu8*C_Qu8 + 378.54*C_td1*C_td1 +
4545 66.48*C_td8*C_td8 + 458.34*C_tq1*C_tq1 + 108.63*C_tq8*C_tq8 + 655.25*C_tu1*C_tu1 + 134.08*C_tu8*C_tu8;
4546 double sigma_neg_bin_mtt_750_1000_NP_quad = 1343.97*C_tG*C_tG + 370.22*C_Qd1*C_Qd1 + 70.13*C_Qd8*C_Qd8 + 481.64*C_Qq11*C_Qq11 +
4547 481.3*C_Qq31*C_Qq31 + 94.58*C_Qq18*C_Qq18 + 96.39*C_Qq38*C_Qq38 + 653.11*C_Qu1*C_Qu1 + 137.27*C_Qu8*C_Qu8 + 206.12*C_td1*C_td1 +
4548 54.0*C_td8*C_td8 + 1028.88*C_tq1*C_tq1 + 202.31*C_tq8*C_tq8 + 279.84*C_tu1*C_tu1 + 53.02*C_tu8*C_tu8;
4549
4550
4551
4552
4553 double sigma_pos_bin_mtt_750_1000_NP_lin = 5266.33*C_tG + 17.3*C_Qd1 + 124.0*C_Qd8 + 20.27*C_Qq11 + 7.09*C_Qq31 + 594.31*C_Qq18 +
4554 162.54*C_Qq38 + 4.65*C_Qu1 + 167.93*C_Qu8 + 3.88*C_td1 + 210.36*C_td8 + 10.89*C_tq1 + 280.34*C_tq8 + 7.1*C_tu1 + 369.73*C_tu8;
4555 double sigma_neg_bin_mtt_750_1000_NP_lin = 5217.21*C_tG + -13.12*C_Qd1 + 209.09*C_Qd8 + -27.46*C_Qq11 + -7.0*C_Qq31 + 255.72*C_Qq18 +
4556 51.39*C_Qq38 + -36.2*C_Qu1 + 356.72*C_Qu8 + -8.9*C_td1 + 120.85*C_td8 + -42.63*C_tq1 + 568.51*C_tq8 + -21.45*C_tu1 + 151.18*C_tu8;
4557
4558
4559 double sigma_pos_bin_mtt_750_1000_NP = sigma_pos_bin_mtt_750_1000_NP_lin + sigma_pos_bin_mtt_750_1000_NP_quad;
4560 double sigma_neg_bin_mtt_750_1000_NP = sigma_neg_bin_mtt_750_1000_NP_lin + sigma_neg_bin_mtt_750_1000_NP_quad;
4561
4562
4563 double sigma_pos_bin_mtt_750_1000_NP_Corrected = SM_sigma_pos_bin_mtt_750_1000*sigma_pos_bin_mtt_750_1000_NP/sigma_pos_bin_mtt_750_1000_madgraph_NLO;
4564 double sigma_neg_bin_mtt_750_1000_NP_Corrected = SM_sigma_neg_bin_mtt_750_1000*sigma_neg_bin_mtt_750_1000_NP/sigma_neg_bin_mtt_750_1000_madgraph_NLO;
4565
4566 double NP_charge_asymmetry_num_bin_mtt_750_1000 = sigma_pos_bin_mtt_750_1000_NP_Corrected - sigma_neg_bin_mtt_750_1000_NP_Corrected;
4567 double NP_charge_asymmetry_deno_bin_mtt_750_1000 = sigma_pos_bin_mtt_750_1000_NP_Corrected + sigma_neg_bin_mtt_750_1000_NP_Corrected;
4568
4569 //return SM_charge_asymmetry_bin_mtt_750_1000*(1+(NP_charge_asymmetry_num_bin_mtt_750_1000-NP_charge_asymmetry_deno_bin_mtt_750_1000)/SM_charge_asymmetry_deno_bin_mtt_750_1000);
4570
4571 return SM_charge_asymmetry_bin_mtt_750_1000 + (NP_charge_asymmetry_num_bin_mtt_750_1000
4572 - SM_charge_asymmetry_bin_mtt_750_1000*NP_charge_asymmetry_deno_bin_mtt_750_1000
4573 )/SM_charge_asymmetry_deno_bin_mtt_750_1000;
4574
4575 }
4576 else{
4577
4578
4579 double sigma_pos_bin_mtt_750_1000_NP = 5266.33*C_tG + 17.3*C_Qd1 + 124.0*C_Qd8 + 20.27*C_Qq11 + 7.09*C_Qq31 + 594.31*C_Qq18 +
4580 162.54*C_Qq38 + 4.65*C_Qu1 + 167.93*C_Qu8 + 3.88*C_td1 + 210.36*C_td8 + 10.89*C_tq1 + 280.34*C_tq8 + 7.1*C_tu1 + 369.73*C_tu8;
4581 double sigma_neg_bin_mtt_750_1000_NP = 5217.21*C_tG + -13.12*C_Qd1 + 209.09*C_Qd8 + -27.46*C_Qq11 + -7.0*C_Qq31 + 255.72*C_Qq18 +
4582 51.39*C_Qq38 + -36.2*C_Qu1 + 356.72*C_Qu8 + -8.9*C_td1 + 120.85*C_td8 + -42.63*C_tq1 + 568.51*C_tq8 + -21.45*C_tu1 + 151.18*C_tu8;
4583
4584
4585 double sigma_pos_bin_mtt_750_1000_NP_Corrected = SM_sigma_pos_bin_mtt_750_1000*sigma_pos_bin_mtt_750_1000_NP/sigma_pos_bin_mtt_750_1000_madgraph_NLO;
4586 double sigma_neg_bin_mtt_750_1000_NP_Corrected = SM_sigma_neg_bin_mtt_750_1000*sigma_neg_bin_mtt_750_1000_NP/sigma_neg_bin_mtt_750_1000_madgraph_NLO;
4587
4588 double NP_charge_asymmetry_num_bin_mtt_750_1000 = sigma_pos_bin_mtt_750_1000_NP_Corrected - sigma_neg_bin_mtt_750_1000_NP_Corrected;
4589 double NP_charge_asymmetry_deno_bin_mtt_750_1000 = sigma_pos_bin_mtt_750_1000_NP_Corrected + sigma_neg_bin_mtt_750_1000_NP_Corrected;
4590
4591 //return SM_charge_asymmetry_bin_mtt_750_1000*(1+(NP_charge_asymmetry_num_bin_mtt_750_1000-NP_charge_asymmetry_deno_bin_mtt_750_1000)/SM_charge_asymmetry_deno_bin_mtt_750_1000);
4592
4593 return SM_charge_asymmetry_bin_mtt_750_1000 + (NP_charge_asymmetry_num_bin_mtt_750_1000
4594 - SM_charge_asymmetry_bin_mtt_750_1000*NP_charge_asymmetry_deno_bin_mtt_750_1000
4595 )/SM_charge_asymmetry_deno_bin_mtt_750_1000;
4596
4597
4598
4599 }
4600
4601 } else if(b_min == 1000 && b_max == 1500){
4602
4603
4604 double SM_charge_asymmetry_deno_bin_mtt_1000_1500 = SM.getOptionalParameter("SM_charge_asymmetry_deno_bin_mtt_1000_1500");
4605 double SM_charge_asymmetry_bin_mtt_1000_1500 = SM.getOptionalParameter("SM_charge_asymmetry_bin_mtt_1000_1500");
4606 double SM_charge_asymmetry_num_bin_mtt_1000_1500 = SM_charge_asymmetry_bin_mtt_1000_1500*SM_charge_asymmetry_deno_bin_mtt_1000_1500;
4607
4608
4609 double SM_sigma_pos_bin_mtt_1000_1500 =0.5*(SM_charge_asymmetry_num_bin_mtt_1000_1500+SM_charge_asymmetry_deno_bin_mtt_1000_1500);
4610 double SM_sigma_neg_bin_mtt_1000_1500 =0.5*(SM_charge_asymmetry_deno_bin_mtt_1000_1500-SM_charge_asymmetry_num_bin_mtt_1000_1500);
4611
4612 double sigma_pos_bin_mtt_1000_1500_madgraph_NLO = 5702.90;
4613 double sigma_neg_bin_mtt_1000_1500_madgraph_NLO = 5641.72;
4614
4615
4616 if(flag_Quadratic){
4617
4618 double sigma_pos_bin_mtt_1000_1500_NP_quad = 563.15*C_tG*C_tG + 247.32*C_Qd1*C_Qd1 + 51.83*C_Qd8*C_Qd8 + 1345.59*C_Qq11*C_Qq11 +
4619 1346.53*C_Qq31*C_Qq31 + 262.79*C_Qq18*C_Qq18 + 273.66*C_Qq38*C_Qq38 + 355.65*C_Qu1*C_Qu1 + 75.51*C_Qu8*C_Qu8 + 471.12*C_td1*C_td1 +
4620 94.38*C_td8*C_td8 + 595.18*C_tq1*C_tq1 + 127.3*C_tq8*C_tq8 + 880.02*C_tu1*C_tu1 + 182.03*C_tu8*C_tu8;
4621 double sigma_neg_bin_mtt_1000_1500_NP_quad = 558.55*C_tG*C_tG + 470.26*C_Qd1*C_Qd1 + 97.31*C_Qd8*C_Qd8 + 597.83*C_Qq11*C_Qq11 +
4622 598.0*C_Qq31*C_Qq31 + 115.55*C_Qq18*C_Qq18 + 121.64*C_Qq38*C_Qq38 + 881.64*C_Qu1*C_Qu1 + 181.08*C_Qu8*C_Qu8 + 246.57*C_td1*C_td1 +
4623 47.61*C_td8*C_td8 + 1347.14*C_tq1*C_tq1 + 275.14*C_tq8*C_tq8 + 353.28*C_tu1*C_tu1 + 72.25*C_tu8*C_tu8;
4624
4625
4626
4627 double sigma_pos_bin_mtt_1000_1500_NP_lin = 1591.22*C_tG + 1.74*C_Qd1 + 67.33*C_Qd8 + 25.81*C_Qq11 + 4.65*C_Qq31 + 411.2*C_Qq18 +
4628 127.63*C_Qq38 + 1.19*C_Qu1 + 97.73*C_Qu8 + 10.72*C_td1 + 142.96*C_td8 + -2.4*C_tq1 + 167.79*C_tq8 + 22.4*C_tu1 + 269.18*C_tu8;
4629 double sigma_neg_bin_mtt_1000_1500_NP_lin = 1572.62*C_tG + -19.54*C_Qd1 + 131.94*C_Qd8 + -9.81*C_Qq11 + 1.09*C_Qq31 + 159.58*C_Qq18 +
4630 22.5*C_Qq38 + -37.72*C_Qu1 + 245.16*C_Qu8 + -5.86*C_td1 + 68.52*C_td8 + -49.98*C_tq1 + 378.21*C_tq8 + -9.48*C_tu1 + 92.32*C_tu8;
4631
4632
4633 double sigma_pos_bin_mtt_1000_1500_NP = sigma_pos_bin_mtt_1000_1500_NP_lin + sigma_pos_bin_mtt_1000_1500_NP_quad;
4634 double sigma_neg_bin_mtt_1000_1500_NP = sigma_neg_bin_mtt_1000_1500_NP_lin + sigma_neg_bin_mtt_1000_1500_NP_quad;
4635
4636
4637
4638 double sigma_pos_bin_mtt_1000_1500_NP_Corrected = SM_sigma_pos_bin_mtt_1000_1500*sigma_pos_bin_mtt_1000_1500_NP/sigma_pos_bin_mtt_1000_1500_madgraph_NLO;
4639 double sigma_neg_bin_mtt_1000_1500_NP_Corrected = SM_sigma_neg_bin_mtt_1000_1500*sigma_neg_bin_mtt_1000_1500_NP/sigma_neg_bin_mtt_1000_1500_madgraph_NLO;
4640
4641 double NP_charge_asymmetry_num_bin_mtt_1000_1500 = sigma_pos_bin_mtt_1000_1500_NP_Corrected - sigma_neg_bin_mtt_1000_1500_NP_Corrected;
4642 double NP_charge_asymmetry_deno_bin_mtt_1000_1500 = sigma_pos_bin_mtt_1000_1500_NP_Corrected + sigma_neg_bin_mtt_1000_1500_NP_Corrected;
4643
4644 //return SM_charge_asymmetry_bin_mtt_1000_1500*(1+(NP_charge_asymmetry_num_bin_mtt_1000_1500-NP_charge_asymmetry_deno_bin_mtt_1000_1500)/SM_charge_asymmetry_deno_bin_mtt_1000_1500);
4645
4646 return SM_charge_asymmetry_bin_mtt_1000_1500 + (NP_charge_asymmetry_num_bin_mtt_1000_1500
4647 - SM_charge_asymmetry_bin_mtt_1000_1500*NP_charge_asymmetry_deno_bin_mtt_1000_1500
4648 )/SM_charge_asymmetry_deno_bin_mtt_1000_1500;
4649
4650
4651 }
4652 else{
4653
4654
4655
4656 double sigma_pos_bin_mtt_1000_1500_NP = 1591.22*C_tG + 1.74*C_Qd1 + 67.33*C_Qd8 + 25.81*C_Qq11 + 4.65*C_Qq31 + 411.2*C_Qq18 +
4657 127.63*C_Qq38 + 1.19*C_Qu1 + 97.73*C_Qu8 + 10.72*C_td1 + 142.96*C_td8 + -2.4*C_tq1 + 167.79*C_tq8 + 22.4*C_tu1 + 269.18*C_tu8;
4658 double sigma_neg_bin_mtt_1000_1500_NP = 1572.62*C_tG + -19.54*C_Qd1 + 131.94*C_Qd8 + -9.81*C_Qq11 + 1.09*C_Qq31 + 159.58*C_Qq18 +
4659 22.5*C_Qq38 + -37.72*C_Qu1 + 245.16*C_Qu8 + -5.86*C_td1 + 68.52*C_td8 + -49.98*C_tq1 + 378.21*C_tq8 + -9.48*C_tu1 + 92.32*C_tu8;
4660
4661
4662 double sigma_pos_bin_mtt_1000_1500_NP_Corrected = SM_sigma_pos_bin_mtt_1000_1500*sigma_pos_bin_mtt_1000_1500_NP/sigma_pos_bin_mtt_1000_1500_madgraph_NLO;
4663 double sigma_neg_bin_mtt_1000_1500_NP_Corrected = SM_sigma_neg_bin_mtt_1000_1500*sigma_neg_bin_mtt_1000_1500_NP/sigma_neg_bin_mtt_1000_1500_madgraph_NLO;
4664
4665 double NP_charge_asymmetry_num_bin_mtt_1000_1500 = sigma_pos_bin_mtt_1000_1500_NP_Corrected - sigma_neg_bin_mtt_1000_1500_NP_Corrected;
4666 double NP_charge_asymmetry_deno_bin_mtt_1000_1500 = sigma_pos_bin_mtt_1000_1500_NP_Corrected + sigma_neg_bin_mtt_1000_1500_NP_Corrected;
4667
4668 //return SM_charge_asymmetry_bin_mtt_1000_1500*(1+(NP_charge_asymmetry_num_bin_mtt_1000_1500-NP_charge_asymmetry_deno_bin_mtt_1000_1500)/SM_charge_asymmetry_deno_bin_mtt_1000_1500);
4669
4670 return SM_charge_asymmetry_bin_mtt_1000_1500 + (NP_charge_asymmetry_num_bin_mtt_1000_1500
4671 - SM_charge_asymmetry_bin_mtt_1000_1500*NP_charge_asymmetry_deno_bin_mtt_1000_1500
4672 )/SM_charge_asymmetry_deno_bin_mtt_1000_1500;
4673
4674
4675
4676 }
4677
4678 } else if(b_min == 1500 && b_max == 3000){
4679
4680 double SM_charge_asymmetry_bin_mtt_1500_3000 = SM.getOptionalParameter("SM_charge_asymmetry_bin_mtt_1500_3000");
4681 double SM_charge_asymmetry_deno_bin_mtt_1500_3000 = SM.getOptionalParameter("SM_charge_asymmetry_deno_bin_mtt_1500_3000");
4682 double SM_charge_asymmetry_num_bin_mtt_1500_3000 = SM_charge_asymmetry_bin_mtt_1500_3000*SM_charge_asymmetry_deno_bin_mtt_1500_3000;
4683
4684 double SM_sigma_pos_bin_mtt_1500_3000 =0.5*(SM_charge_asymmetry_num_bin_mtt_1500_3000+SM_charge_asymmetry_deno_bin_mtt_1500_3000);
4685 double SM_sigma_neg_bin_mtt_1500_3000 =0.5*(SM_charge_asymmetry_deno_bin_mtt_1500_3000-SM_charge_asymmetry_num_bin_mtt_1500_3000);
4686
4687 double sigma_pos_bin_mtt_1500_3000_madgraph_NLO = 746.10;
4688 double sigma_neg_bin_mtt_1500_3000_madgraph_NLO = 735.27;
4689
4690
4691 if(flag_Quadratic){
4692
4693
4694 double sigma_pos_bin_mtt_1500_3000_NP_quad = 120.81*C_tG*C_tG + 202.85*C_Qd1*C_Qd1 + 42.63*C_Qd8*C_Qd8 + 1495.21*C_Qq11*C_Qq11 +
4695 1494.64*C_Qq31*C_Qq31 + 314.09*C_Qq18*C_Qq18 + 321.03*C_Qq38*C_Qq38 + 449.71*C_Qu1*C_Qu1 + 93.31*C_Qu8*C_Qu8 + 471.54*C_td1*C_td1 +
4696 100.31*C_td8*C_td8 + 650.11*C_tq1*C_tq1 + 135.25*C_tq8*C_tq8 + 1025.45*C_tu1*C_tu1 + 220.28*C_tu8*C_tu8;
4697 double sigma_neg_bin_mtt_1500_3000_NP_quad = 119.3*C_tG*C_tG + 471.71*C_Qd1*C_Qd1 + 93.84*C_Qd8*C_Qd8 + 650.6*C_Qq11*C_Qq11 +
4698 650.82*C_Qq31*C_Qq31 + 132.25*C_Qq18*C_Qq18 + 135.33*C_Qq38*C_Qq38 + 1025.52*C_Qu1*C_Qu1 + 203.34*C_Qu8*C_Qu8 + 202.66*C_td1*C_td1 +
4699 42.05*C_td8*C_td8 + 1494.96*C_tq1*C_tq1 + 296.52*C_tq8*C_tq8 + 448.95*C_tu1*C_tu1 + 92.58*C_tu8*C_tu8;
4700
4701
4702 double sigma_pos_bin_mtt_1500_3000_NP_lin = 204.01*C_tG + -2.08*C_Qd1 + 21.2*C_Qd8 + 23.76*C_Qq11 + 8.49*C_Qq31 + 172.2*C_Qq18 +
4703 59.21*C_Qq38 + -3.52*C_Qu1 + 39.69*C_Qu8 + 7.84*C_td1 + 56.64*C_td8 + -5.58*C_tq1 + 61.85*C_tq8 + 16.14*C_tu1 + 115.53*C_tu8;
4704 double sigma_neg_bin_mtt_1500_3000_NP_lin = 200.5*C_tG + -11.33*C_Qd1 + 47.73*C_Qd8 + -0.59*C_Qq11 + -1.24*C_Qq31 + 63.93*C_Qq18 +
4705 18.88*C_Qq38 + -24.09*C_Qu1 + 98.88*C_Qu8 + 0.43*C_td1 + 22.17*C_td8 + -36.32*C_tq1 + 147.64*C_tq8 + -0.96*C_tu1 + 40.92*C_tu8;
4706
4707 double sigma_pos_bin_mtt_1500_3000_NP = sigma_pos_bin_mtt_1500_3000_NP_lin + sigma_pos_bin_mtt_1500_3000_NP_quad;
4708 double sigma_neg_bin_mtt_1500_3000_NP = sigma_neg_bin_mtt_1500_3000_NP_lin + sigma_neg_bin_mtt_1500_3000_NP_quad;
4709
4710
4711 double sigma_pos_bin_mtt_1500_3000_NP_Corrected = SM_sigma_pos_bin_mtt_1500_3000*sigma_pos_bin_mtt_1500_3000_NP/sigma_pos_bin_mtt_1500_3000_madgraph_NLO;
4712 double sigma_neg_bin_mtt_1500_3000_NP_Corrected = SM_sigma_neg_bin_mtt_1500_3000*sigma_neg_bin_mtt_1500_3000_NP/sigma_neg_bin_mtt_1500_3000_madgraph_NLO;
4713
4714 double NP_charge_asymmetry_num_bin_mtt_1500_3000 = sigma_pos_bin_mtt_1500_3000_NP_Corrected - sigma_neg_bin_mtt_1500_3000_NP_Corrected;
4715 double NP_charge_asymmetry_deno_bin_mtt_1500_3000 = sigma_pos_bin_mtt_1500_3000_NP_Corrected + sigma_neg_bin_mtt_1500_3000_NP_Corrected;
4716
4717 //return SM_charge_asymmetry_bin_mtt_1500_3000*(1+(NP_charge_asymmetry_num_bin_mtt_1500_3000-NP_charge_asymmetry_deno_bin_mtt_1500_3000)/SM_charge_asymmetry_deno_bin_mtt_1500_3000);
4718
4719 return SM_charge_asymmetry_bin_mtt_1500_3000 + (NP_charge_asymmetry_num_bin_mtt_1500_3000
4720 - SM_charge_asymmetry_bin_mtt_1500_3000*NP_charge_asymmetry_deno_bin_mtt_1500_3000
4721 )/SM_charge_asymmetry_deno_bin_mtt_1500_3000;
4722
4723
4724 }
4725 else{
4726
4727
4728 double sigma_pos_bin_mtt_1500_3000_NP = 204.01*C_tG + -2.08*C_Qd1 + 21.2*C_Qd8 + 23.76*C_Qq11 + 8.49*C_Qq31 + 172.2*C_Qq18 +
4729 59.21*C_Qq38 + -3.52*C_Qu1 + 39.69*C_Qu8 + 7.84*C_td1 + 56.64*C_td8 + -5.58*C_tq1 + 61.85*C_tq8 + 16.14*C_tu1 + 115.53*C_tu8;
4730 double sigma_neg_bin_mtt_1500_3000_NP = 200.5*C_tG + -11.33*C_Qd1 + 47.73*C_Qd8 + -0.59*C_Qq11 + -1.24*C_Qq31 + 63.93*C_Qq18 +
4731 18.88*C_Qq38 + -24.09*C_Qu1 + 98.88*C_Qu8 + 0.43*C_td1 + 22.17*C_td8 + -36.32*C_tq1 + 147.64*C_tq8 + -0.96*C_tu1 + 40.92*C_tu8;
4732
4733
4734 double sigma_pos_bin_mtt_1500_3000_NP_Corrected = SM_sigma_pos_bin_mtt_1500_3000*sigma_pos_bin_mtt_1500_3000_NP/sigma_pos_bin_mtt_1500_3000_madgraph_NLO;
4735 double sigma_neg_bin_mtt_1500_3000_NP_Corrected = SM_sigma_neg_bin_mtt_1500_3000*sigma_neg_bin_mtt_1500_3000_NP/sigma_neg_bin_mtt_1500_3000_madgraph_NLO;
4736
4737 double NP_charge_asymmetry_num_bin_mtt_1500_3000 = sigma_pos_bin_mtt_1500_3000_NP_Corrected - sigma_neg_bin_mtt_1500_3000_NP_Corrected;
4738 double NP_charge_asymmetry_deno_bin_mtt_1500_3000 = sigma_pos_bin_mtt_1500_3000_NP_Corrected + sigma_neg_bin_mtt_1500_3000_NP_Corrected;
4739
4740 //return SM_charge_asymmetry_bin_mtt_1500_3000*(1+(NP_charge_asymmetry_num_bin_mtt_1500_3000-NP_charge_asymmetry_deno_bin_mtt_1500_3000)/SM_charge_asymmetry_deno_bin_mtt_1500_3000);
4741
4742 return SM_charge_asymmetry_bin_mtt_1500_3000 + (NP_charge_asymmetry_num_bin_mtt_1500_3000
4743 - SM_charge_asymmetry_bin_mtt_1500_3000*NP_charge_asymmetry_deno_bin_mtt_1500_3000
4744 )/SM_charge_asymmetry_deno_bin_mtt_1500_3000;
4745
4746 }
4747
4748 } else if(b_min == 1500 && b_max == 2000){
4749
4750
4751
4752 double SM_charge_asymmetry_bin_mtt_1500_2000 = SM.getOptionalParameter("SM_charge_asymmetry_bin_mtt_1500_2000");
4753 double SM_charge_asymmetry_deno_bin_mtt_1500_2000 = SM.getOptionalParameter("SM_charge_asymmetry_deno_bin_mtt_1500_2000");
4754 double SM_charge_asymmetry_num_bin_mtt_1500_2000 = SM_charge_asymmetry_bin_mtt_1500_2000*SM_charge_asymmetry_deno_bin_mtt_1500_2000;
4755
4756 double SM_sigma_pos_bin_mtt_1500_2000 =0.5*(SM_charge_asymmetry_num_bin_mtt_1500_2000+SM_charge_asymmetry_deno_bin_mtt_1500_2000);
4757 double SM_sigma_neg_bin_mtt_1500_2000 =0.5*(SM_charge_asymmetry_deno_bin_mtt_1500_2000-SM_charge_asymmetry_num_bin_mtt_1500_2000);
4758
4759
4760 //double Delta_Madgraph = -1.400000000000091;
4761 //double Total_Madgraph = 1519.0;
4762 //double pos_Madgraph = 758.8;
4763 //double neg_Madgraph = 760.2;
4764
4765 double sigma_pos_bin_mtt_1500_2000_madgraph = 623.06;
4766 double sigma_neg_bin_mtt_1500_2000_madgraph = 613.37;
4767
4768
4769
4770
4771 if(flag_Quadratic){
4772
4773
4774 //double Delta = Delta_Madgraph +0.6*C_tG +127.3*C_td1*C_td1 +26.33*C_td8 +28.28*C_td8*C_td8 +322.4*C_tu1*C_tu1 +65.61*C_tu8 +71.49*C_tu8*C_tu8 +449.8*C_Qq11*C_Qq11 +449.6*C_Qq31*C_Qq31 +96.74*C_Qq18 +99.97*C_Qq18*C_Qq18 +46.49*C_Qq38 +99.92*C_Qq38*C_Qq38 +2.194*C_tG*C_td8 +6.783*C_tG*C_tu8 +9.33*C_tG*C_Qq18 +3.806*C_tG*C_Qq38 +0.67*C_Qd1*C_td1 +0.163*C_Qu8*C_tu8 +389.5*C_Qq11*C_Qq31 +86.86*C_Qq18*C_Qq38 ;
4775 //double Total = Total_Madgraph +365.6*C_tG +193.64*C_tG*C_tG +81.69*C_Qd8 +86.77*C_Qd8*C_Qd8 +390.8*C_Qd1*C_Qd1 +729.9*C_Qu1*C_Qu1 +161.18*C_Qu8 +162.26*C_Qu8*C_Qu8 +390.7*C_td1*C_td1 +81.07*C_td8 +86.88*C_td8*C_td8 +1118.8*C_tq1*C_tq1 +237.84*C_tq8 +248.78*C_tq8*C_tq8 +729.6*C_tu1*C_tu1 +150.79*C_tu8 +162.11*C_tu8*C_tu8 +1119.2*C_Qq11*C_Qq11 +1119.2*C_Qq31*C_Qq31 +230.46*C_Qq18 +248.63*C_Qq18*C_Qq18 +72.81*C_Qq38 +248.68*C_Qq38*C_Qq38 +12.235*C_tG*C_Qd8 +21.962*C_tG*C_Qu8 +11.538*C_tG*C_td8 +34.19*C_tG*C_tq8 +22.517*C_tG*C_tu8 +34.89*C_tG*C_Qq18 +10.308*C_tG*C_Qq38 +5.379*C_Qd8*C_td8 +23.45*C_Qd1*C_td1 +45.23*C_Qu1*C_tu1 +9.957*C_Qu8*C_tu8 +68.33*C_tq1*C_Qq11 +21.29*C_tq1*C_Qq31 +15.024*C_tq8*C_Qq18 +4.593*C_tq8*C_Qq38 +681.1*C_Qq11*C_Qq31 +151.74*C_Qq18*C_Qq38 ;
4776
4777
4778 double sigma_pos_bin_mtt_1500_2000_NP_quad = 92.57*C_tG*C_tG + 118.16*C_Qd1*C_Qd1 + 24.32*C_Qd8*C_Qd8 + 734.31*C_Qq11*C_Qq11 +
4779 732.63*C_Qq31*C_Qq31 + 148.8*C_Qq18*C_Qq18 + 152.52*C_Qq38*C_Qq38 + 201.03*C_Qu1*C_Qu1 + 41.1*C_Qu8*C_Qu8
4780 + 245.73*C_td1*C_td1 + 50.24*C_td8*C_td8 + 317.54*C_tq1*C_tq1 + 65.51*C_tq8*C_tq8 + 489.18*C_tu1*C_tu1 + 104.44*C_tu8*C_tu8;
4781 double sigma_neg_bin_mtt_1500_2000_NP_quad = 91.95*C_tG*C_tG + 245.59*C_Qd1*C_Qd1 + 48.46*C_Qd8*C_Qd8 + 318.95*C_Qq11*C_Qq11 +
4782 320.29*C_Qq31*C_Qq31 + 63.03*C_Qq18*C_Qq18 + 65.21*C_Qq38*C_Qq38 + 490.42*C_Qu1*C_Qu1 + 96.83*C_Qu8*C_Qu8 +
4783 118.57*C_td1*C_td1 + 24.0*C_td8*C_td8 + 733.93*C_tq1*C_tq1 + 145.07*C_tq8*C_tq8 + 202.0*C_tu1*C_tu1 + 40.84*C_tu8*C_tu8;
4784
4785
4786 double sigma_pos_bin_mtt_1500_2000_NP_lin = 170.81*C_tG + -1.77*C_Qd1 + 15.41*C_Qd8 + 12.18*C_Qq11 + 2.58*C_Qq31 + 114.92*C_Qq18 +
4787 37.21*C_Qq38 + -1.56*C_Qu1 + 25.46*C_Qu8 + 5.03*C_td1 + 39.19*C_td8 + -3.63*C_tq1 + 42.28*C_tq8 + 10.0*C_tu1 + 75.73*C_tu8;
4788 double sigma_neg_bin_mtt_1500_2000_NP_lin = 167.48*C_tG + -6.39*C_Qd1 + 32.97*C_Qd8 + -0.51*C_Qq11 + 0.27*C_Qq31 + 42.69*C_Qq18 +
4789 9.01*C_Qq38 + -15.76*C_Qu1 + 65.88*C_Qu8 + 0.08*C_td1 + 16.22*C_td8 + -23.19*C_tq1 + 100.18*C_tq8 + -2.51*C_tu1 + 25.44*C_tu8;
4790
4791
4792 double sigma_pos_bin_mtt_1500_2000_NP = sigma_pos_bin_mtt_1500_2000_NP_lin + sigma_pos_bin_mtt_1500_2000_NP_quad;
4793 double sigma_neg_bin_mtt_1500_2000_NP = sigma_neg_bin_mtt_1500_2000_NP_lin + sigma_neg_bin_mtt_1500_2000_NP_quad;
4794
4795
4796 double sigma_pos_bin_mtt_1500_2000_NP_Corrected = SM_sigma_pos_bin_mtt_1500_2000*sigma_pos_bin_mtt_1500_2000_NP/sigma_pos_bin_mtt_1500_2000_madgraph;
4797 double sigma_neg_bin_mtt_1500_2000_NP_Corrected = SM_sigma_neg_bin_mtt_1500_2000*sigma_neg_bin_mtt_1500_2000_NP/sigma_neg_bin_mtt_1500_2000_madgraph;
4798
4799 double NP_charge_asymmetry_num_bin_mtt_1500_2000 = sigma_pos_bin_mtt_1500_2000_NP_Corrected - sigma_neg_bin_mtt_1500_2000_NP_Corrected;
4800 double NP_charge_asymmetry_deno_bin_mtt_1500_2000 = sigma_pos_bin_mtt_1500_2000_NP_Corrected + sigma_neg_bin_mtt_1500_2000_NP_Corrected;
4801
4802 //return SM_charge_asymmetry_bin_mtt_1500_2000*(1+(NP_charge_asymmetry_num_bin_mtt_1500_2000-NP_charge_asymmetry_deno_bin_mtt_1500_2000)/SM_charge_asymmetry_deno_bin_mtt_1500_2000);
4803
4804 return SM_charge_asymmetry_bin_mtt_1500_2000 + (NP_charge_asymmetry_num_bin_mtt_1500_2000
4805 - SM_charge_asymmetry_bin_mtt_1500_2000*NP_charge_asymmetry_deno_bin_mtt_1500_2000
4806 )/SM_charge_asymmetry_deno_bin_mtt_1500_2000;
4807
4808 }
4809 else{
4810
4811
4812 double sigma_pos_bin_mtt_1500_2000_NP = 170.81*C_tG + -1.77*C_Qd1 + 15.41*C_Qd8 + 12.18*C_Qq11 + 2.58*C_Qq31 + 114.92*C_Qq18 +
4813 37.21*C_Qq38 + -1.56*C_Qu1 + 25.46*C_Qu8 + 5.03*C_td1 + 39.19*C_td8 + -3.63*C_tq1 + 42.28*C_tq8 + 10.0*C_tu1 + 75.73*C_tu8;
4814 double sigma_neg_bin_mtt_1500_2000_NP = 167.48*C_tG + -6.39*C_Qd1 + 32.97*C_Qd8 + -0.51*C_Qq11 + 0.27*C_Qq31 + 42.69*C_Qq18 +
4815 9.01*C_Qq38 + -15.76*C_Qu1 + 65.88*C_Qu8 + 0.08*C_td1 + 16.22*C_td8 + -23.19*C_tq1 + 100.18*C_tq8 + -2.51*C_tu1 + 25.44*C_tu8;
4816
4817
4818 double sigma_pos_bin_mtt_1500_2000_NP_Corrected = SM_sigma_pos_bin_mtt_1500_2000*sigma_pos_bin_mtt_1500_2000_NP/sigma_pos_bin_mtt_1500_2000_madgraph;
4819 double sigma_neg_bin_mtt_1500_2000_NP_Corrected = SM_sigma_neg_bin_mtt_1500_2000*sigma_neg_bin_mtt_1500_2000_NP/sigma_neg_bin_mtt_1500_2000_madgraph;
4820
4821 double NP_charge_asymmetry_num_bin_mtt_1500_2000 = sigma_pos_bin_mtt_1500_2000_NP_Corrected - sigma_neg_bin_mtt_1500_2000_NP_Corrected;
4822 double NP_charge_asymmetry_deno_bin_mtt_1500_2000 = sigma_pos_bin_mtt_1500_2000_NP_Corrected + sigma_neg_bin_mtt_1500_2000_NP_Corrected;
4823
4824 //return SM_charge_asymmetry_bin_mtt_1500_2000*(1+(NP_charge_asymmetry_num_bin_mtt_1500_2000-NP_charge_asymmetry_deno_bin_mtt_1500_2000)/SM_charge_asymmetry_deno_bin_mtt_1500_2000);
4825
4826 return SM_charge_asymmetry_bin_mtt_1500_2000 + (NP_charge_asymmetry_num_bin_mtt_1500_2000
4827 - SM_charge_asymmetry_bin_mtt_1500_2000*NP_charge_asymmetry_deno_bin_mtt_1500_2000
4828 )/SM_charge_asymmetry_deno_bin_mtt_1500_2000;
4829
4830 }
4831
4832 } else if(b_min == 2000 && b_max == 2500){
4833
4834
4835 double SM_charge_asymmetry_bin_mtt_2000_2500 = SM.getOptionalParameter("SM_charge_asymmetry_bin_mtt_2000_2500");
4836 double SM_charge_asymmetry_deno_bin_mtt_2000_2500 = SM.getOptionalParameter("SM_charge_asymmetry_deno_bin_mtt_2000_2500");
4837 double SM_charge_asymmetry_num_bin_mtt_2000_2500 = SM_charge_asymmetry_bin_mtt_2000_2500*SM_charge_asymmetry_deno_bin_mtt_2000_2500;
4838
4839 double SM_sigma_pos_bin_mtt_2000_2500 =0.5*(SM_charge_asymmetry_num_bin_mtt_2000_2500+SM_charge_asymmetry_deno_bin_mtt_2000_2500);
4840 double SM_sigma_neg_bin_mtt_2000_2500 =0.5*(SM_charge_asymmetry_deno_bin_mtt_2000_2500-SM_charge_asymmetry_num_bin_mtt_2000_2500);
4841
4842
4843 //double Delta_Madgraph = 9.899999999999977;
4844 //double Total_Madgraph = 280.5;
4845 //double pos_Madgraph = 145.2;
4846 //double neg_Madgraph = 135.3;
4847
4848 double sigma_pos_bin_mtt_2000_2500_madgraph = 91.28;
4849 double sigma_neg_bin_mtt_2000_2500_madgraph = 89.45;
4850
4851
4852 if(flag_Quadratic){
4853
4854
4855 double sigma_pos_bin_mtt_2000_2500_NP_quad = 19.16*C_tG*C_tG + 52.13*C_Qd1*C_Qd1 + 11.31*C_Qd8*C_Qd8 + 380.79*C_Qq11*C_Qq11 +
4856 380.23*C_Qq31*C_Qq31 + 80.68*C_Qq18*C_Qq18 + 83.66*C_Qq38*C_Qq38 + 111.34*C_Qu1*C_Qu1 + 24.39*C_Qu8*C_Qu8 +
4857 121.22*C_td1*C_td1 + 26.15*C_td8*C_td8 + 162.15*C_tq1*C_tq1 + 35.82*C_tq8*C_tq8 + 259.37*C_tu1*C_tu1 + 55.42*C_tu8*C_tu8;
4858 double sigma_neg_bin_mtt_2000_2500_NP_quad = 18.96*C_tG*C_tG + 122.29*C_Qd1*C_Qd1 + 25.74*C_Qd8*C_Qd8 + 161.29*C_Qq11*C_Qq11 +
4859 160.46*C_Qq31*C_Qq31 + 33.11*C_Qq18*C_Qq18 + 33.92*C_Qq38*C_Qq38 + 259.13*C_Qu1*C_Qu1 + 54.48*C_Qu8*C_Qu8 +
4860 52.33*C_td1*C_td1 + 11.05*C_td8*C_td8 + 379.61*C_tq1*C_tq1 + 79.62*C_tq8*C_tq8 + 110.6*C_tu1*C_tu1 + 22.87*C_tu8*C_tu8;
4861
4862
4863 double sigma_pos_bin_mtt_2000_2500_NP_lin = 24.76*C_tG + -0.52*C_Qd1 + 3.67*C_Qd8 + 3.98*C_Qq11 + -0.13*C_Qq31 + 36.66*C_Qq18 +
4864 12.61*C_Qq38 + -0.59*C_Qu1 + 8.95*C_Qu8 + 1.25*C_td1 + 11.32*C_td8 + -0.29*C_tq1 + 13.01*C_tq8 + 1.89*C_tu1 + 24.59*C_tu8;
4865 double sigma_neg_bin_mtt_2000_2500_NP_lin = 24.28*C_tG + -2.52*C_Qd1 + 10.16*C_Qd8 + 0.47*C_Qq11 + 0.78*C_Qq31 + 13.11*C_Qq18 +
4866 4.73*C_Qq38 + -4.64*C_Qu1 + 21.13*C_Qu8 + 0.9*C_td1 + 4.33*C_td8 + -7.25*C_tq1 + 31.39*C_tq8 + 1.68*C_tu1 + 8.81*C_tu8;
4867
4868 double sigma_pos_bin_mtt_2000_2500_NP = sigma_pos_bin_mtt_2000_2500_NP_lin + sigma_pos_bin_mtt_2000_2500_NP_quad;
4869 double sigma_neg_bin_mtt_2000_2500_NP = sigma_neg_bin_mtt_2000_2500_NP_lin + sigma_neg_bin_mtt_2000_2500_NP_quad;
4870
4871
4872 double sigma_pos_bin_mtt_2000_2500_NP_Corrected = SM_sigma_pos_bin_mtt_2000_2500*sigma_pos_bin_mtt_2000_2500_NP/sigma_pos_bin_mtt_2000_2500_madgraph;
4873 double sigma_neg_bin_mtt_2000_2500_NP_Corrected = SM_sigma_neg_bin_mtt_2000_2500*sigma_neg_bin_mtt_2000_2500_NP/sigma_neg_bin_mtt_2000_2500_madgraph;
4874
4875 double NP_charge_asymmetry_num_bin_mtt_2000_2500 = sigma_pos_bin_mtt_2000_2500_NP_Corrected - sigma_neg_bin_mtt_2000_2500_NP_Corrected;
4876 double NP_charge_asymmetry_deno_bin_mtt_2000_2500 = sigma_pos_bin_mtt_2000_2500_NP_Corrected + sigma_neg_bin_mtt_2000_2500_NP_Corrected;
4877
4878 //return SM_charge_asymmetry_bin_mtt_2000_2500*(1+(NP_charge_asymmetry_num_bin_mtt_2000_2500-NP_charge_asymmetry_deno_bin_mtt_2000_2500)/SM_charge_asymmetry_deno_bin_mtt_2000_2500);
4879
4880 return SM_charge_asymmetry_bin_mtt_2000_2500 + (NP_charge_asymmetry_num_bin_mtt_2000_2500
4881 - SM_charge_asymmetry_bin_mtt_2000_2500*NP_charge_asymmetry_deno_bin_mtt_2000_2500
4882 )/SM_charge_asymmetry_deno_bin_mtt_2000_2500;
4883
4884 }
4885 else{
4886
4887
4888 double sigma_pos_bin_mtt_2000_2500_NP = 24.76*C_tG + -0.52*C_Qd1 + 3.67*C_Qd8 + 3.98*C_Qq11 + -0.13*C_Qq31 + 36.66*C_Qq18 +
4889 12.61*C_Qq38 + -0.59*C_Qu1 + 8.95*C_Qu8 + 1.25*C_td1 + 11.32*C_td8 + -0.29*C_tq1 + 13.01*C_tq8 + 1.89*C_tu1 + 24.59*C_tu8;
4890 double sigma_neg_bin_mtt_2000_2500_NP = 24.28*C_tG + -2.52*C_Qd1 + 10.16*C_Qd8 + 0.47*C_Qq11 + 0.78*C_Qq31 + 13.11*C_Qq18 +
4891 4.73*C_Qq38 + -4.64*C_Qu1 + 21.13*C_Qu8 + 0.9*C_td1 + 4.33*C_td8 + -7.25*C_tq1 + 31.39*C_tq8 + 1.68*C_tu1 + 8.81*C_tu8;
4892
4893
4894
4895 double sigma_pos_bin_mtt_2000_2500_NP_Corrected = SM_sigma_pos_bin_mtt_2000_2500*sigma_pos_bin_mtt_2000_2500_NP/sigma_pos_bin_mtt_2000_2500_madgraph;
4896 double sigma_neg_bin_mtt_2000_2500_NP_Corrected = SM_sigma_neg_bin_mtt_2000_2500*sigma_neg_bin_mtt_2000_2500_NP/sigma_neg_bin_mtt_2000_2500_madgraph;
4897
4898 double NP_charge_asymmetry_num_bin_mtt_2000_2500 = sigma_pos_bin_mtt_2000_2500_NP_Corrected - sigma_neg_bin_mtt_2000_2500_NP_Corrected;
4899 double NP_charge_asymmetry_deno_bin_mtt_2000_2500 = sigma_pos_bin_mtt_2000_2500_NP_Corrected + sigma_neg_bin_mtt_2000_2500_NP_Corrected;
4900
4901 //return SM_charge_asymmetry_bin_mtt_2000_2500*(1+(NP_charge_asymmetry_num_bin_mtt_2000_2500-NP_charge_asymmetry_deno_bin_mtt_2000_2500)/SM_charge_asymmetry_deno_bin_mtt_2000_2500);
4902
4903 return SM_charge_asymmetry_bin_mtt_2000_2500 + (NP_charge_asymmetry_num_bin_mtt_2000_2500
4904 - SM_charge_asymmetry_bin_mtt_2000_2500*NP_charge_asymmetry_deno_bin_mtt_2000_2500
4905 )/SM_charge_asymmetry_deno_bin_mtt_2000_2500;
4906
4907
4908
4909 }
4910
4911 } else if(b_min == 2500 && b_max == 3000){
4912
4913
4914 double SM_charge_asymmetry_bin_mtt_2500_3000 = SM.getOptionalParameter("SM_charge_asymmetry_bin_mtt_2500_3000");
4915 double SM_charge_asymmetry_deno_bin_mtt_2500_3000 = SM.getOptionalParameter("SM_charge_asymmetry_deno_bin_mtt_2500_3000");
4916 double SM_charge_asymmetry_num_bin_mtt_2500_3000 = SM_charge_asymmetry_bin_mtt_2500_3000*SM_charge_asymmetry_deno_bin_mtt_2500_3000;
4917
4918 double SM_sigma_pos_bin_mtt_2500_3000 =0.5*(SM_charge_asymmetry_num_bin_mtt_2500_3000+SM_charge_asymmetry_deno_bin_mtt_2500_3000);
4919 double SM_sigma_neg_bin_mtt_2500_3000 =0.5*(SM_charge_asymmetry_deno_bin_mtt_2500_3000-SM_charge_asymmetry_num_bin_mtt_2500_3000);
4920
4921
4922
4923
4924 //double Delta_Madgraph = -2.8999999999999986;
4925 //double Total_Madgraph = 61.12;
4926 //double pos_Madgraph = 29.11;
4927 //double neg_Madgraph = 32.01;
4928
4929
4930 double sigma_pos_bin_mtt_2500_3000_madgraph = 17.54;
4931 double sigma_neg_bin_mtt_2500_3000_madgraph = 17.24;
4932
4933
4934 if(flag_Quadratic){
4935
4936 //double Delta = Delta_Madgraph +51.68*C_td1*C_td1 +3.262*C_td8 +11.478*C_td8*C_td8 +120.07*C_tu1*C_tu1 +9.458*C_tu8 +26.67*C_tu8*C_tu8 +171.8*C_Qq11*C_Qq11 +172.2*C_Qq31*C_Qq31 +13.305*C_Qq18 +38.22*C_Qq18*C_Qq18 +5.559*C_Qq38 +38.21*C_Qq38*C_Qq38 +0.564*C_tG*C_td8 +0.901*C_tG*C_tu8 +1.423*C_tG*C_Qq18 +0.386*C_tG*C_Qq38 +0.006*C_Qd8*C_td8 +0.624*C_Qu1*C_tu1 +0.232*C_Qu8*C_tu8 +0.243*C_tq1*C_Qq11 +0.001*C_tq8*C_Qq38 +135.8*C_Qq11*C_Qq31 +30.27*C_Qq18*C_Qq38 ;
4937 //double Total = Total_Madgraph +14.606*C_tG +15.074*C_tG*C_tG +11.4*C_Qd8 +29.416*C_Qd8*C_Qd8 +132.52*C_Qd1*C_Qd1 +302.15*C_Qu1*C_Qu1 +21.356*C_Qu8 +67.13*C_Qu8*C_Qu8 +132.38*C_td1*C_td1 +9.424*C_td8 +29.482*C_td8*C_td8 +434.2*C_tq1*C_tq1 +38.76*C_tq8 +96.46*C_tq8*C_tq8 +302.13*C_tu1*C_tu1 +27.042*C_tu8 +67.15*C_tu8*C_tu8 +434.2*C_Qq11*C_Qq11 +434.4*C_Qq31*C_Qq31 +32.115*C_Qq18 +96.56*C_Qq18*C_Qq18 +12.195*C_Qq38 +96.51*C_Qq38*C_Qq38 +1.869*C_tG*C_Qd8 +3.491*C_tG*C_Qu8 +1.346*C_tG*C_td8 +5.49*C_tG*C_tq8 +3.599*C_tG*C_tu8 +5.263*C_tG*C_Qq18 +2.061*C_tG*C_Qq38 +0.653*C_Qd8*C_td8 +3.018*C_Qd1*C_td1 +7.298*C_Qu1*C_tu1 +1.591*C_Qu8*C_tu8 +10.329*C_tq1*C_Qq11 +3.613*C_tq1*C_Qq31 +2.332*C_tq8*C_Qq18 +0.958*C_tq8*C_Qq38 +339.2*C_Qq11*C_Qq31 +75.31*C_Qq18*C_Qq38 ;
4938
4939
4940
4941 double sigma_pos_bin_mtt_2500_3000_NP_quad = 4.98*C_tG*C_tG + 21.28*C_Qd1*C_Qd1 + 4.74*C_Qd8*C_Qd8 + 192.55*C_Qq11*C_Qq11 +
4942 194.11*C_Qq31*C_Qq31 + 41.97*C_Qq18*C_Qq18 + 43.67*C_Qq38*C_Qq38 + 61.34*C_Qu1*C_Qu1 + 13.52*C_Qu8*C_Qu8 +
4943 58.14*C_td1*C_td1 + 12.9*C_td8*C_td8 + 82.9*C_tq1*C_tq1 + 18.17*C_tq8*C_tq8 + 134.43*C_tu1*C_tu1 + 29.76*C_tu8*C_tu8;
4944 double sigma_neg_bin_mtt_2500_3000_NP_quad = 4.85*C_tG*C_tG + 58.14*C_Qd1*C_Qd1 + 12.44*C_Qd8*C_Qd8 + 83.33*C_Qq11*C_Qq11 +
4945 82.69*C_Qq31*C_Qq31 + 17.09*C_Qq18*C_Qq18 + 17.76*C_Qq38*C_Qq38 + 135.1*C_Qu1*C_Qu1 + 28.74*C_Qu8*C_Qu8 +
4946 21.41*C_td1*C_td1 + 4.69*C_td8*C_td8 + 194.18*C_tq1*C_tq1 + 40.96*C_tq8*C_tq8 + 62.11*C_tu1*C_tu1 + 12.69*C_tu8*C_tu8;
4947
4948
4949 double sigma_pos_bin_mtt_2500_3000_NP_lin = 4.73*C_tG + -0.44*C_Qd1 + 1.06*C_Qd8 + 0.43*C_Qq11 + 0.7*C_Qq31 + 12.57*C_Qq18 +
4950 4.94*C_Qq38 + -0.45*C_Qu1 + 3.36*C_Qu8 + 0.91*C_td1 + 3.77*C_td8 + 0.43*C_tq1 + 4.39*C_tq8 + 1.14*C_tu1 + 8.57*C_tu8;
4951 double sigma_neg_bin_mtt_2500_3000_NP_lin = 4.61*C_tG + -0.79*C_Qd1 + 3.23*C_Qd8 + 0.53*C_Qq11 + -0.51*C_Qq31 + 4.9*C_Qq18 +
4952 2.42*C_Qq38 + -1.35*C_Qu1 + 7.23*C_Qu8 + -0.29*C_td1 + 1.18*C_td8 + -2.82*C_tq1 + 10.79*C_tq8 + 0.59*C_tu1 + 3.24*C_tu8;
4953
4954
4955 double sigma_pos_bin_mtt_2500_3000_NP = sigma_pos_bin_mtt_2500_3000_NP_lin + sigma_pos_bin_mtt_2500_3000_NP_quad;
4956 double sigma_neg_bin_mtt_2500_3000_NP = sigma_neg_bin_mtt_2500_3000_NP_lin + sigma_neg_bin_mtt_2500_3000_NP_quad;
4957
4958
4959 double sigma_pos_bin_mtt_2500_3000_NP_Corrected = SM_sigma_pos_bin_mtt_2500_3000*sigma_pos_bin_mtt_2500_3000_NP/sigma_pos_bin_mtt_2500_3000_madgraph;
4960 double sigma_neg_bin_mtt_2500_3000_NP_Corrected = SM_sigma_neg_bin_mtt_2500_3000*sigma_neg_bin_mtt_2500_3000_NP/sigma_neg_bin_mtt_2500_3000_madgraph;
4961
4962 double NP_charge_asymmetry_num_bin_mtt_2500_3000 = sigma_pos_bin_mtt_2500_3000_NP_Corrected - sigma_neg_bin_mtt_2500_3000_NP_Corrected;
4963 double NP_charge_asymmetry_deno_bin_mtt_2500_3000 = sigma_pos_bin_mtt_2500_3000_NP_Corrected + sigma_neg_bin_mtt_2500_3000_NP_Corrected;
4964
4965 //return SM_charge_asymmetry_bin_mtt_2500_3000*(1+(NP_charge_asymmetry_num_bin_mtt_2500_3000-NP_charge_asymmetry_deno_bin_mtt_2500_3000)/SM_charge_asymmetry_deno_bin_mtt_2500_3000);
4966
4967 return SM_charge_asymmetry_bin_mtt_2500_3000 + (NP_charge_asymmetry_num_bin_mtt_2500_3000
4968 - SM_charge_asymmetry_bin_mtt_2500_3000*NP_charge_asymmetry_deno_bin_mtt_2500_3000
4969 )/SM_charge_asymmetry_deno_bin_mtt_2500_3000;
4970
4971
4972 }
4973 else{
4974
4975
4976 double sigma_pos_bin_mtt_2500_3000_NP = 4.73*C_tG + -0.44*C_Qd1 + 1.06*C_Qd8 + 0.43*C_Qq11 + 0.7*C_Qq31 + 12.57*C_Qq18 +
4977 4.94*C_Qq38 + -0.45*C_Qu1 + 3.36*C_Qu8 + 0.91*C_td1 + 3.77*C_td8 + 0.43*C_tq1 + 4.39*C_tq8 + 1.14*C_tu1 + 8.57*C_tu8;
4978 double sigma_neg_bin_mtt_2500_3000_NP = 4.61*C_tG + -0.79*C_Qd1 + 3.23*C_Qd8 + 0.53*C_Qq11 + -0.51*C_Qq31 + 4.9*C_Qq18 +
4979 2.42*C_Qq38 + -1.35*C_Qu1 + 7.23*C_Qu8 + -0.29*C_td1 + 1.18*C_td8 + -2.82*C_tq1 + 10.79*C_tq8 + 0.59*C_tu1 + 3.24*C_tu8;
4980
4981
4982
4983 double sigma_pos_bin_mtt_2500_3000_NP_Corrected = SM_sigma_pos_bin_mtt_2500_3000*sigma_pos_bin_mtt_2500_3000_NP/sigma_pos_bin_mtt_2500_3000_madgraph;
4984 double sigma_neg_bin_mtt_2500_3000_NP_Corrected = SM_sigma_neg_bin_mtt_2500_3000*sigma_neg_bin_mtt_2500_3000_NP/sigma_neg_bin_mtt_2500_3000_madgraph;
4985
4986 double NP_charge_asymmetry_num_bin_mtt_2500_3000 = sigma_pos_bin_mtt_2500_3000_NP_Corrected - sigma_neg_bin_mtt_2500_3000_NP_Corrected;
4987 double NP_charge_asymmetry_deno_bin_mtt_2500_3000 = sigma_pos_bin_mtt_2500_3000_NP_Corrected + sigma_neg_bin_mtt_2500_3000_NP_Corrected;
4988
4989 //return SM_charge_asymmetry_bin_mtt_2500_3000*(1+(NP_charge_asymmetry_num_bin_mtt_2500_3000-NP_charge_asymmetry_deno_bin_mtt_2500_3000)/SM_charge_asymmetry_deno_bin_mtt_2500_3000);
4990
4991
4992 return SM_charge_asymmetry_bin_mtt_2500_3000 + (NP_charge_asymmetry_num_bin_mtt_2500_3000
4993 - SM_charge_asymmetry_bin_mtt_2500_3000*NP_charge_asymmetry_deno_bin_mtt_2500_3000
4994 )/SM_charge_asymmetry_deno_bin_mtt_2500_3000;
4995
4996
4997 }
4998 }
4999
5000
5001 else {
5002 throw std::runtime_error("\nERROR: Please specify a correct binning range for charge_asymmetry_tt_diff_mtt_NLO.\n");
5003 }
5004
5005}
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_flag_Quadratic
double getNPSMEFT6dtopquark_flag_Quadratic() const
Definition: NPSMEFT6dtopquark.h:220
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_C_Qd1
double getNPSMEFT6dtopquark_C_Qd1() const
Definition: NPSMEFT6dtopquark.h:164
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_C_Qq38
double getNPSMEFT6dtopquark_C_Qq38() const
Definition: NPSMEFT6dtopquark.h:147
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_C_Qq31
double getNPSMEFT6dtopquark_C_Qq31() const
Definition: NPSMEFT6dtopquark.h:182
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_C_Qq18
double getNPSMEFT6dtopquark_C_Qq18() const
Definition: NPSMEFT6dtopquark.h:141
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_C_tq1
double getNPSMEFT6dtopquark_C_tq1() const
Definition: NPSMEFT6dtopquark.h:176
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_C_tG
double getNPSMEFT6dtopquark_C_tG() const
Definition: NPSMEFT6dtopquark.h:101
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_C_Qu1
double getNPSMEFT6dtopquark_C_Qu1() const
Definition: NPSMEFT6dtopquark.h:167
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_C_tu8
double getNPSMEFT6dtopquark_C_tu8() const
Definition: NPSMEFT6dtopquark.h:134
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_C_Qu8
double getNPSMEFT6dtopquark_C_Qu8() const
Definition: NPSMEFT6dtopquark.h:150
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_C_td8
double getNPSMEFT6dtopquark_C_td8() const
Definition: NPSMEFT6dtopquark.h:138
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_C_Qd8
double getNPSMEFT6dtopquark_C_Qd8() const
Definition: NPSMEFT6dtopquark.h:153
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_C_tq8
double getNPSMEFT6dtopquark_C_tq8() const
Definition: NPSMEFT6dtopquark.h:144
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_C_Qq11
double getNPSMEFT6dtopquark_C_Qq11() const
Definition: NPSMEFT6dtopquark.h:179
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_C_td1
double getNPSMEFT6dtopquark_C_td1() const
Definition: NPSMEFT6dtopquark.h:170
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_C_tu1
double getNPSMEFT6dtopquark_C_tu1() const
Definition: NPSMEFT6dtopquark.h:173
QCD::getOptionalParameter
const double getOptionalParameter(std::string name) const
A method to get parameters that are specific to only one set of observables.
Definition: QCD.h:450
ThObservable::getBinMin
double getBinMin()
A get method to get the minimum value of the bin.
Definition: ThObservable.h:82
ThObservable::getBinMax
double getBinMax()
A get method to get the maximum value of the bin.
Definition: ThObservable.h:91
ThObservable::SM
const StandardModel & SM
A reference to an object of StandardMode class.
Definition: ThObservable.h:121
charge_asymmetry_tt_diff_mtt_NLO::b_min
double b_min
Definition: NPSMEFT6dtopquark.h:1319
charge_asymmetry_tt_diff_mtt_NLO::b_max
double b_max
Definition: NPSMEFT6dtopquark.h:1320

Member Data Documentation

◆ b_max

double charge_asymmetry_tt_diff_mtt_NLO::b_max = 0.
private

Definition at line 1320 of file NPSMEFT6dtopquark.h.

◆ b_min

double charge_asymmetry_tt_diff_mtt_NLO::b_min = 0.
private

Definition at line 1319 of file NPSMEFT6dtopquark.h.

◆ myNPSMEFT6dtopquark

const NPSMEFT6dtopquark& charge_asymmetry_tt_diff_mtt_NLO::myNPSMEFT6dtopquark
private

Definition at line 1323 of file NPSMEFT6dtopquark.h.

The documentation for this class was generated from the following files: