a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
sigma_tHj_diff_LO_mth Class Reference

#include <NPSMEFT6dtopquark.h>

+ Inheritance diagram for sigma_tHj_diff_LO_mth:

Detailed Description

Definition at line 3818 of file NPSMEFT6dtopquark.h.

Public Member Functions

double computeThValue ()
 
 sigma_tHj_diff_LO_mth (const StandardModel &SM_i)
 sigma_tHj_diff_LO_mth constructor. More...
 
- 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 StandardModelgetModel ()
 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 NPSMEFT6dtopquarkmyNPSMEFT6dtopquark
 

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 StandardModelSM
 A reference to an object of StandardMode class. More...
 

Constructor & Destructor Documentation

◆ sigma_tHj_diff_LO_mth()

sigma_tHj_diff_LO_mth::sigma_tHj_diff_LO_mth ( const StandardModel SM_i)

sigma_tHj_diff_LO_mth constructor.

Definition at line 18470 of file NPSMEFT6dtopquark.cpp.

18471: ThObservable(SM_i),myNPSMEFT6dtopquark(static_cast<const NPSMEFT6dtopquark&> (SM_i))
18472{
18473
18474 setParametersForObservable(make_vector<std::string>() << "SM_sigma_tHj_diff_LO_mth_bin_200_340"
18475 << "SM_sigma_tHj_diff_LO_mth_bin_340_424" << "SM_sigma_tHj_diff_LO_mth_bin_424_620"
18476 << "SM_sigma_tHj_diff_LO_mth_bin_620_1600" );
18477
18478}
void setParametersForObservable(std::vector< std::string > parametersForObservable_i)
A set method to get the parameters for the specific observable.
Definition: ThObservable.h:109
ThObservable(const StandardModel &SM_i)
Constructor.
Definition: ThObservable.h:32
const NPSMEFT6dtopquark & myNPSMEFT6dtopquark

Member Function Documentation

◆ computeThValue()

double sigma_tHj_diff_LO_mth::computeThValue ( )
virtual
Returns
The value of sigma_tHj_diff_LO_mth

Implements ThObservable.

Definition at line 18480 of file NPSMEFT6dtopquark.cpp.

18481{
18482
18483 b_min = getBinMin();
18484 b_max = getBinMax();
18485
18486
18489
18490
18491 //double cHG = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_phiG();
18492 //double cHGtil = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_phiGtil();
18497 //double ctGRe = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_tG();
18498 //double ctGIm = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_tGIm();
18501
18502
18503
18504 //Expressions are written in the basis of SMEFTsim
18505 //let's change to the basis of dim6top (same as
18506 //SMEFT@NLO except for the gs on CtG). These
18507 //models define the covariance derivative
18508 //with a different convention
18509
18510 if(flag_LHC_WG_Basis){
18511// double ctZRe = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_tZ();
18512// double ctZIm = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_tZIm();
18513// sw2 = 0.22305;
18514// cw2 = 1 - sw2;
18515// tw2 = sw2/cw2;
18516// sw = sqrt(sw2);
18517// tw = sqrt(tw2);
18518// ctBRe = ctZRe/sw - ctWRe/tw;
18519// ctBIm = ctZIm/sw - ctWIm/tw;
18520 ctWRe = -ctWRe;
18521 ctWIm = -ctWIm;
18522// ctGRe = -ctGRe;
18523// ctGIm = -ctGIm;
18524 }
18525// else{
18526// double ctBRe = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_tB();
18527// double ctBIm = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_tBIm();
18528// }
18529
18530
18531
18532 if(b_min == 200 && b_max == 340){
18533
18534 double SM_sigma_tHj_diff_LO_mth_bin_200_340 = SM.getOptionalParameter("SM_sigma_tHj_diff_LO_mth_bin_200_340");
18535 double sigma_tHj_diff_LO_mth_bin_200_340_madgraph = 4.8312389333236;
18536 double total;
18537
18538
18539 if(flag_Quadratic){
18540
18541 total = SM_sigma_tHj_diff_LO_mth_bin_200_340 + (2.044625422571*cHW + -0.0019697321068240514*cHWtil + -0.0002770594241113811*ctHIm +
18542 -0.9319723301799*ctHRe + -0.00476895114764897*ctWIm + -4.0956850184244*ctWRe +
18543 0.2928532476119*cHW*cHW + 0.09546992012006428*cHWtil*cHWtil + 0.01454612189076554*ctHIm*ctHIm +
18544 0.1209362969154*ctHRe*ctHRe + 2.7343139034392*ctWIm*ctWIm + 2.6995685969779997*ctWRe*ctWRe +
18545 -0.005934498577654642*cHW*cHWtil + -0.0018803268831102064*cHW*ctHIm + -0.2194966254055*cHW*ctHRe +
18546 0.003760700951024876*cHW*ctWIm + -1.2677935653004*cHW*ctWRe + -0.02264745989419263*cHWtil*ctHIm +
18547 -0.000901401797751721*cHWtil*ctHRe + -0.5172919603443*cHWtil*ctWIm + 0.0005135710164916318*cHWtil*ctWRe +
18548 0.001622458841596397*ctHIm*ctHRe + 0.1287365371493*ctHIm*ctWIm + -0.005794226903136551*ctHIm*ctWRe +
18549 -0.0036847586393026423*ctHRe*ctWIm + 0.37264722012039997*ctHRe*ctWRe + 0.00219928660555313*ctWIm*ctWRe
18550 )*(SM_sigma_tHj_diff_LO_mth_bin_200_340/sigma_tHj_diff_LO_mth_bin_200_340_madgraph);
18551 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
18552 return total;
18553
18554 }
18555 else{
18556
18557 total = SM_sigma_tHj_diff_LO_mth_bin_200_340 + (2.044625422571*cHW +
18558 -0.9319723301799*ctHRe + -4.0956850184244*ctWRe
18559 )*(SM_sigma_tHj_diff_LO_mth_bin_200_340/sigma_tHj_diff_LO_mth_bin_200_340_madgraph);
18560 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
18561 return total;
18562
18563 }
18564
18565 } else if(b_min == 340 && b_max == 424){
18566
18567 double SM_sigma_tHj_diff_LO_mth_bin_340_424 = SM.getOptionalParameter("SM_sigma_tHj_diff_LO_mth_bin_340_424");
18568 double sigma_tHj_diff_LO_mth_bin_340_424_madgraph = 15.7487145559461;
18569 double total;
18570
18571
18572 if(flag_Quadratic){
18573
18574 total = SM_sigma_tHj_diff_LO_mth_bin_340_424 + (4.9981313432903995*cHW + -0.00433369039049003*cHWtil + 0.00039224905580559366*ctHIm +
18575 -1.0434688823926*ctHRe + 0.003958581262702987*ctWIm + -11.327199630097*ctWRe +
18576 1.3719672914796999*cHW*cHW + 0.6464375218769001*cHWtil*cHWtil + 0.06497929284932091*ctHIm*ctHIm +
18577 0.2766726496071*ctHRe*ctHRe + 10.297289637753002*ctWIm*ctWIm + 10.2481785950402*ctWRe*ctWRe +
18578 -0.015556904787081122*cHW*cHWtil + 0.00012779846093879946*cHW*ctHIm + -0.5974518855705999*cHW*ctHRe +
18579 0.01289784306163222*cHW*ctWIm + -4.6693468708844*cHW*ctWRe + -0.07474505546859601*cHWtil*ctHIm +
18580 -0.0017800266970389987*cHWtil*ctHRe + -2.2400298269615*cHWtil*ctWIm + -0.018183520573615563*cHWtil*ctWRe +
18581 -0.0003375385250080826*ctHIm*ctHRe + 0.5066845255559*ctHIm*ctWIm + 0.02001141484961333*ctHIm*ctWRe +
18582 0.02423769252751784*ctHRe*ctWIm + 1.0459736936309*ctHRe*ctWRe + 0.0640694541759329*ctWIm*ctWRe
18583 )*(SM_sigma_tHj_diff_LO_mth_bin_340_424/sigma_tHj_diff_LO_mth_bin_340_424_madgraph);
18584 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
18585 return total;
18586
18587 }
18588 else{
18589
18590 total = SM_sigma_tHj_diff_LO_mth_bin_340_424 + (4.9981313432903995*cHW +
18591 -1.0434688823926*ctHRe + -11.327199630097*ctWRe
18592 )*(SM_sigma_tHj_diff_LO_mth_bin_340_424/sigma_tHj_diff_LO_mth_bin_340_424_madgraph);
18593 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
18594 return total;
18595
18596
18597 }
18598
18599
18600 } else if(b_min == 424 && b_max == 620){
18601
18602 double SM_sigma_tHj_diff_LO_mth_bin_424_620 = SM.getOptionalParameter("SM_sigma_tHj_diff_LO_mth_bin_424_620");
18603 double sigma_tHj_diff_LO_mth_bin_424_620_madgraph = 10.9357549691572;
18604 double total;
18605
18606
18607 if(flag_Quadratic){
18608
18609 total = SM_sigma_tHj_diff_LO_mth_bin_424_620 + (1.5666660066113*cHW + 0.009873964296414073*cHWtil + -0.003017092917417808*ctHIm +
18610 0.08256721872489103*ctHRe + -0.008516060692950744*ctWIm + -4.136280835257099*ctWRe +
18611 1.0890377895256*cHW*cHW + 0.7007065924062*cHWtil*cHWtil + 0.0390975212028466*ctHIm*ctHIm +
18612 0.09289750818528009*ctHRe*ctHRe + 7.2496001463796*ctWIm*ctWIm + 7.2836843498583*ctWRe*ctWRe +
18613 0.009983943811839311*cHW*cHWtil + -0.0012669630490372552*cHW*ctHIm + -0.1910815830423*cHW*ctHRe +
18614 -0.0110446887577412*cHW*ctWIm + -2.7163356332474002*cHW*ctWRe + -0.02914841005248382*cHWtil*ctHIm +
18615 -0.0011939154296200202*cHWtil*ctHRe + -1.6033076535583*cHWtil*ctWIm + -0.003923502799203852*cHWtil*ctWRe +
18616 0.0007781238255122336*ctHIm*ctHRe + 0.2845742436244*ctHIm*ctWIm + -0.006589865014594622*ctHIm*ctWRe +
18617 -0.014562482561483663*ctHRe*ctWIm + 0.4522172735994*ctHRe*ctWRe + -0.03856803764934983*ctWIm*ctWRe
18618 )*(SM_sigma_tHj_diff_LO_mth_bin_424_620/sigma_tHj_diff_LO_mth_bin_424_620_madgraph);
18619 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
18620 return total;
18621
18622 }
18623 else{
18624
18625 total = SM_sigma_tHj_diff_LO_mth_bin_424_620 + (1.5666660066113*cHW +
18626 0.08256721872489103*ctHRe -4.136280835257099*ctWRe
18627 )*(SM_sigma_tHj_diff_LO_mth_bin_424_620/sigma_tHj_diff_LO_mth_bin_424_620_madgraph);
18628 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
18629 return total;
18630
18631 }
18632
18633
18634 } else if(b_min == 620 && b_max == 1600){
18635
18636 double SM_sigma_tHj_diff_LO_mth_bin_620_1600 = SM.getOptionalParameter("SM_sigma_tHj_diff_LO_mth_bin_620_1600");
18637 double sigma_tHj_diff_LO_mth_bin_620_1600_madgraph = 2.0062389761137;
18638 double total;
18639
18640
18641 if(flag_Quadratic){
18642
18643 total = SM_sigma_tHj_diff_LO_mth_bin_620_1600 + (0.1111424096413*cHW + -0.0014860638211933177*cHWtil + 0.0003135970091938667*ctHIm +
18644 0.02440742106196714*ctHRe + 0.0029584200534060745*ctWIm + -0.3390724228498*ctWRe +
18645 0.29202717282470003*cHW*cHW + 0.2430289940472*cHWtil*cHWtil + 0.005608877718665806*ctHIm*ctHIm +
18646 0.008411593347051882*ctHRe*ctHRe + 2.6193177426385*ctWIm*ctWIm + 2.6275190748015*ctWRe*ctWRe +
18647 0.0003166402726267892*cHW*cHWtil + 0.0003708859349638161*cHW*ctHIm + -0.010811319724541921*cHW*ctHRe +
18648 -0.004126259538358132*cHW*ctWIm + -0.4598013457796*cHW*ctWRe + -0.0019298379837962567*cHWtil*ctHIm +
18649 0.0005878696857919419*cHWtil*ctHRe + -0.3198592839751*cHWtil*ctWIm + 0.0026820414812528397*cHWtil*ctWRe +
18650 -8.366197674275488e-05*ctHIm*ctHRe + 0.045894456459680266*ctHIm*ctWIm + -0.0014537499925816463*ctHIm*ctWRe +
18651 -0.0007806602952328085*ctHRe*ctWIm + 0.05445691902987437*ctHRe*ctWRe + -0.0006867750649745359*ctWIm*ctWRe
18652 )*(SM_sigma_tHj_diff_LO_mth_bin_620_1600/sigma_tHj_diff_LO_mth_bin_620_1600_madgraph);
18653 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
18654 return total;
18655
18656 }
18657 else{
18658
18659 total = SM_sigma_tHj_diff_LO_mth_bin_620_1600 + (0.1111424096413*cHW +
18660 0.02440742106196714*ctHRe + -0.3390724228498*ctWRe
18661 )*(SM_sigma_tHj_diff_LO_mth_bin_620_1600/sigma_tHj_diff_LO_mth_bin_620_1600_madgraph);
18662 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
18663 return total;
18664
18665 }
18666
18667
18668 } else{
18669 throw std::runtime_error("\nERROR: Please specify a correct binning for sigma_tHj_diff_LO_mth. \n");
18670 }
18671
18672}
double getNPSMEFT6dtopquark_flag_Quadratic() const
double getNPSMEFT6dtopquark_C_phiWtil() const
double getNPSMEFT6dtopquark_C_phiW() const
double getNPSMEFT6dtopquark_C_tphi() const
double getNPSMEFT6dtopquark_C_tWIm() const
double getNPSMEFT6dtopquark_C_tW() const
double getNPSMEFT6dtopquark_flag_LHC_WG_Basis() const
double getNPSMEFT6dtopquark_C_tphiIm() const
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
double getBinMin()
A get method to get the minimum value of the bin.
Definition: ThObservable.h:82
double getBinMax()
A get method to get the maximum value of the bin.
Definition: ThObservable.h:91
const StandardModel & SM
A reference to an object of StandardMode class.
Definition: ThObservable.h:121

Member Data Documentation

◆ b_max

double sigma_tHj_diff_LO_mth::b_max = 0.
private

Definition at line 3834 of file NPSMEFT6dtopquark.h.

◆ b_min

double sigma_tHj_diff_LO_mth::b_min = 0.
private

Definition at line 3833 of file NPSMEFT6dtopquark.h.

◆ myNPSMEFT6dtopquark

const NPSMEFT6dtopquark& sigma_tHj_diff_LO_mth::myNPSMEFT6dtopquark
private

Definition at line 3836 of file NPSMEFT6dtopquark.h.


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