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
sigma_tHj_eta_cut_3_diff_LO_mth Class Reference

#include <NPSMEFT6dtopquark.h>

+ Inheritance diagram for sigma_tHj_eta_cut_3_diff_LO_mth:

Detailed Description

Definition at line 4240 of file NPSMEFT6dtopquark.h.

Public Member Functions

double computeThValue ()
 
 sigma_tHj_eta_cut_3_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 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

◆ sigma_tHj_eta_cut_3_diff_LO_mth()

sigma_tHj_eta_cut_3_diff_LO_mth::sigma_tHj_eta_cut_3_diff_LO_mth ( const StandardModel &  SM_i)

sigma_tHj_diff_LO_mth constructor.

Definition at line 16165 of file NPSMEFT6dtopquark.cpp.

16166: ThObservable(SM_i),myNPSMEFT6dtopquark(static_cast<const NPSMEFT6dtopquark&> (SM_i))
16167{
16168
16169 setParametersForObservable(make_vector<std::string>() << "SM_sigma_tHj_eta_cut_3_diff_LO_mth_bin_200_340"
16170 << "SM_sigma_tHj_eta_cut_3_diff_LO_mth_bin_340_424" << "SM_sigma_tHj_eta_cut_3_diff_LO_mth_bin_424_620"
16171 << "SM_sigma_tHj_eta_cut_3_diff_LO_mth_bin_620_1600" );
16172
16173}
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
sigma_tHj_eta_cut_3_diff_LO_mth::myNPSMEFT6dtopquark
const NPSMEFT6dtopquark & myNPSMEFT6dtopquark
Definition: NPSMEFT6dtopquark.h:4258

Member Function Documentation

◆ computeThValue()

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

Implements ThObservable.

Definition at line 16175 of file NPSMEFT6dtopquark.cpp.

16176{
16177
16178 b_min = getBinMin();
16179 b_max = getBinMax();
16180
16181
16182 bool flag_LHC_WG_Basis=myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_flag_LHC_WG_Basis();
16183 bool flag_Quadratic=myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_flag_Quadratic();
16184
16185
16186 //double cHG = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_phiG();
16187 //double cHGtil = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_phiGtil();
16188 double cHW = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_phiW();
16189 double cHWtil = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_phiWtil();
16190 double ctHRe = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_tphi();
16191 double ctHIm = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_tphiIm();
16192 //double ctGRe = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_tG();
16193 //double ctGIm = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_tGIm();
16194 double ctWRe = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_tW();
16195 double ctWIm = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_tWIm();
16196
16197
16198
16199 //Expressions are written in the basis of SMEFTsim
16200 //let's change to the basis of dim6top (same as
16201 //SMEFT@NLO except for the gs on CtG). These
16202 //models define the covariance derivative
16203 //with a different convention
16204
16205 if(flag_LHC_WG_Basis){
16206// double ctZRe = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_tZ();
16207// double ctZIm = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_tZIm();
16208// sw2 = 0.22305;
16209// cw2 = 1 - sw2;
16210// tw2 = sw2/cw2;
16211// sw = sqrt(sw2);
16212// tw = sqrt(tw2);
16213// ctBRe = ctZRe/sw - ctWRe/tw;
16214// ctBIm = ctZIm/sw - ctWIm/tw;
16215 ctWRe = -ctWRe;
16216 ctWIm = -ctWIm;
16217// ctGRe = -ctGRe;
16218// ctGIm = -ctGIm;
16219 }
16220// else{
16221// double ctBRe = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_tB();
16222// double ctBIm = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_tBIm();
16223// }
16224
16225
16226
16227 if(b_min == 200 && b_max == 340){
16228
16229 double SM_sigma_tHj_diff_LO_mth_bin_200_340 = SM.getOptionalParameter("SM_sigma_tHj_eta_cut_3_diff_LO_mth_bin_200_340");
16230 double sigma_tHj_diff_LO_mth_bin_200_340_madgraph = 3.619170460509;
16231 double total;
16232
16233
16234 if(flag_Quadratic){
16235
16236 total = SM_sigma_tHj_diff_LO_mth_bin_200_340 + (
16237 +1.6688407195767 * cHW + -3.5812013807479115e-06 * cHWtil
16238 + -3.4892943340247846e-05 * ctHIm + -0.0006791132796056 * ctHRe
16239 + -0.0010936424908458695 * ctWIm + -0.003224678245972 * ctWRe
16240 + 0.2418823809179 * cHW * cHW + 7.162653900651128e-05 * cHWtil * cHWtil
16241 + 0.008469162500260308 * ctHIm * ctHIm + 6.923755378236109e-05 * ctHRe * ctHRe
16242 + 2.2933642674933 * ctWIm * ctWIm + 0.0022532892306283 * ctWRe * ctWRe
16243 + -0.00261869580730388 * cHW * cHWtil + -3.405171232406982e-07 * cHW * ctHIm
16244 + -0.1670538936151 * cHW * ctHRe + -4.6452008964002846e-07 * cHW * ctWIm
16245 + -1.02348582078 * cHW * ctWRe + -1.682847484238523e-05 * cHWtil * ctHIm
16246 + 0.001594654388675787 * cHWtil * ctHRe + -0.0003923925156083 * cHWtil * ctWIm
16247 + 0.0004999535798966281 * cHWtil * ctWRe + 1.2325714571354007e-06 * ctHIm * ctHRe
16248 + 0.093722219741188 * ctHIm * ctWIm + -2.116668056491424e-06 * ctHIm * ctWRe
16249 + 0.001935814944203917 * ctHRe * ctWIm + 0.0002769065914965 * ctHRe * ctWRe
16250 + 0.00290519535819933 * ctWIm * ctWRe
16251 )*(SM_sigma_tHj_diff_LO_mth_bin_200_340/sigma_tHj_diff_LO_mth_bin_200_340_madgraph);
16252 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
16253 return total;
16254
16255 }
16256 else{
16257
16258 total = SM_sigma_tHj_diff_LO_mth_bin_200_340 + (
16259 +1.6688407195767 * cHW + -3.5812013807479115e-06 * cHWtil
16260 + -3.4892943340247846e-05 * ctHIm + -0.0006791132796056 * ctHRe
16261 + -0.0010936424908458695 * ctWIm + -0.003224678245972 * ctWRe
16262 )*(SM_sigma_tHj_diff_LO_mth_bin_200_340/sigma_tHj_diff_LO_mth_bin_200_340_madgraph);
16263 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
16264 return total;
16265
16266 }
16267
16268 } else if(b_min == 340 && b_max == 424){
16269
16270 double SM_sigma_tHj_diff_LO_mth_bin_340_424 = SM.getOptionalParameter("SM_sigma_tHj_eta_cut_3_diff_LO_mth_bin_340_424");
16271 double sigma_tHj_diff_LO_mth_bin_340_424_madgraph = 10.1005706067801 ;
16272 double total;
16273
16274
16275 if(flag_Quadratic){
16276
16277 total = SM_sigma_tHj_diff_LO_mth_bin_340_424 + (
16278 +4.2197545277171 * cHW + 4.65626312914913e-06 * cHWtil
16279 + 0.00110257955547588 * ctHIm + -0.0009232206777953 * ctHRe
16280 + 0.016350925755151402 * ctWIm + -0.0091507596702359 * ctWRe
16281 + 1.1372943240898 * cHW * cHW + 0.0004991102290958 * cHWtil * cHWtil
16282 + 0.03732598136366172 * ctHIm * ctHIm + 0.0001579255440007 * ctHRe * ctHRe
16283 + 8.641173989847701 * ctWIm * ctWIm + 0.0085860888364268 * ctWRe * ctWRe
16284 + -0.0096821268035852 * cHW * cHWtil + -2.2611564512684646e-07 * cHW * ctHIm
16285 + -0.4676610149002 * cHW * ctHRe + -1.1181980792429467e-05 * cHW * ctWIm
16286 + -3.8410103403863 * cHW * ctWRe + -5.810895365138154e-05 * cHWtil * ctHIm
16287 + -0.002368717938218697 * cHWtil * ctHRe + -0.0017742847988793 * cHWtil * ctWIm
16288 + -0.0027609453798538826 * cHWtil * ctWRe + -3.475760730948959e-07 * ctHIm * ctHRe
16289 + 0.3639693446778 * ctHIm * ctWIm + 1.2234552083004169e-05 * ctHIm * ctWRe
16290 + 0.023309500109203185 * ctHRe * ctWIm + 0.0007873256126956 * ctHRe * ctWRe
16291 + 0.012731860115647309 * ctWIm * ctWRe
16292 )*(SM_sigma_tHj_diff_LO_mth_bin_340_424/sigma_tHj_diff_LO_mth_bin_340_424_madgraph);
16293 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
16294 return total;
16295
16296 }
16297 else{
16298
16299 total = SM_sigma_tHj_diff_LO_mth_bin_340_424 + (
16300 +4.2197545277171 * cHW + 4.65626312914913e-06 * cHWtil
16301 + 0.00110257955547588 * ctHIm + -0.0009232206777953 * ctHRe
16302 + 0.016350925755151402 * ctWIm + -0.0091507596702359 * ctWRe
16303 )*(SM_sigma_tHj_diff_LO_mth_bin_340_424/sigma_tHj_diff_LO_mth_bin_340_424_madgraph);
16304 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
16305 return total;
16306
16307
16308 }
16309
16310
16311 } else if(b_min == 424 && b_max == 620){
16312
16313 double SM_sigma_tHj_diff_LO_mth_bin_424_620 = SM.getOptionalParameter("SM_sigma_tHj_eta_cut_3_diff_LO_mth_bin_424_620");
16314 double sigma_tHj_diff_LO_mth_bin_424_620_madgraph = 5.7714361029407;
16315 double total;
16316
16317
16318 if(flag_Quadratic){
16319
16320 total = SM_sigma_tHj_diff_LO_mth_bin_424_620 + (
16321 +1.389635531438 * cHW + 2.0075927352758427e-06 * cHWtil
16322 + -0.001606695693268101 * ctHIm + -4.878234711728907e-05 * ctHRe
16323 + -0.017851429430962296 * ctWIm + -0.0034886170635417 * ctWRe
16324 + 0.9121486583392 * cHW * cHW + 0.0005627776186415 * cHWtil * cHWtil
16325 + 0.020487480249700707 * ctHIm * ctHIm + 4.743040827425722e-05 * ctHRe * ctHRe
16326 + 6.176908178650599 * ctWIm * ctWIm + 0.0061673172456817 * ctWRe * ctWRe
16327 + 0.009237818024768196 * cHW * cHWtil + -1.8763342924085208e-06 * cHW * ctHIm
16328 + -0.1543347988108 * cHW * ctHRe + -8.639129171580245e-06 * cHW * ctWIm
16329 + -2.2912016280382996 * cHW * ctWRe + -2.565744440081829e-05 * cHWtil * ctHIm
16330 + -0.001054510328439647 * cHWtil * ctHRe + -0.0012993864923725 * cHWtil * ctWIm
16331 + -0.0054575783759468555 * cHWtil * ctWRe + 2.992717955396795e-08 * ctHIm * ctHRe
16332 + 0.1988338348315 * ctHIm * ctWIm + -1.029540757266518e-05 * ctHIm * ctWRe
16333 + -0.00336567140662442 * ctHRe * ctWIm + 0.0003272434816741 * ctHRe * ctWRe
16334 + -0.03863568663995408 * ctWIm * ctWRe
16335 )*(SM_sigma_tHj_diff_LO_mth_bin_424_620/sigma_tHj_diff_LO_mth_bin_424_620_madgraph);
16336 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
16337 return total;
16338
16339 }
16340 else{
16341
16342 total = SM_sigma_tHj_diff_LO_mth_bin_424_620 + (
16343 +1.389635531438 * cHW + 2.0075927352758427e-06 * cHWtil
16344 + -0.001606695693268101 * ctHIm + -4.878234711728907e-05 * ctHRe
16345 + -0.017851429430962296 * ctWIm + -0.0034886170635417 * ctWRe
16346 )*(SM_sigma_tHj_diff_LO_mth_bin_424_620/sigma_tHj_diff_LO_mth_bin_424_620_madgraph);
16347 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
16348 return total;
16349
16350 }
16351
16352
16353 } else if(b_min == 620 && b_max == 1600){
16354
16355 double SM_sigma_tHj_diff_LO_mth_bin_620_1600 = SM.getOptionalParameter("SM_sigma_tHj_eta_cut_3_diff_LO_mth_bin_620_1600");
16356 double sigma_tHj_diff_LO_mth_bin_620_1600_madgraph = 0.8348276093313001;
16357 double total;
16358
16359
16360 if(flag_Quadratic){
16361
16362 total = SM_sigma_tHj_diff_LO_mth_bin_620_1600 + (
16363 +0.10821364459869999 * cHW + -2.8874340988610525e-06 * cHWtil
16364 + 0.00022487637149626322 * ctHIm + 7.364134720744264e-06 * ctHRe
16365 + 0.0002415664477062318 * ctWIm + -0.0003139176772788 * ctWRe
16366 + 0.2451823332766 * cHW * cHW + 0.0001982244893762 * cHWtil * cHWtil
16367 + 0.002821945454561936 * ctHIm * ctHIm + 4.463189174436222e-06 * ctHRe * ctHRe
16368 + 2.2668638066297 * ctWIm * ctWIm + 0.0022708676532284 * ctWRe * ctWRe
16369 + -0.0010133779178814002 * cHW * cHWtil + -1.3685136584326957e-07 * cHW * ctHIm
16370 + -0.009371844326277735 * cHW * ctHRe + -1.1690871514530163e-06 * cHW * ctWIm
16371 + -0.3997185744133 * cHW * ctWRe + -1.5282911600881814e-06 * cHWtil * ctHIm
16372 + -4.7102894794844364e-05 * cHWtil * ctHRe + -0.0002736081199381 * cHWtil * ctWIm
16373 + 0.003457044334853609 * cHWtil * ctWRe + 1.265766582797534e-07 * ctHIm * ctHRe
16374 + 0.03202182013571642 * ctHIm * ctWIm + -3.289807089625314e-07 * ctHIm * ctWRe
16375 + -0.0005056660404694216 * ctHRe * ctWIm + 3.93069423164083e-05 * ctHRe * ctWRe
16376 + -0.005593267455097791 * ctWIm * ctWRe
16377 )*(SM_sigma_tHj_diff_LO_mth_bin_620_1600/sigma_tHj_diff_LO_mth_bin_620_1600_madgraph);
16378 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
16379 return total;
16380
16381 }
16382 else{
16383
16384 total = SM_sigma_tHj_diff_LO_mth_bin_620_1600 + (
16385 +0.10821364459869999 * cHW + -2.8874340988610525e-06 * cHWtil
16386 + 0.00022487637149626322 * ctHIm + 7.364134720744264e-06 * ctHRe
16387 + 0.0002415664477062318 * ctWIm + -0.0003139176772788 * ctWRe
16388 )*(SM_sigma_tHj_diff_LO_mth_bin_620_1600/sigma_tHj_diff_LO_mth_bin_620_1600_madgraph);
16389 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
16390 return total;
16391
16392 }
16393
16394
16395 } else{
16396 throw std::runtime_error("\nERROR: Please specify a correct binning for sigma_tHj_eta_cut_3_diff_LO_mth. \n");
16397 }
16398
16399}
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_flag_Quadratic
double getNPSMEFT6dtopquark_flag_Quadratic() const
Definition: NPSMEFT6dtopquark.h:220
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_C_phiWtil
double getNPSMEFT6dtopquark_C_phiWtil() const
Definition: NPSMEFT6dtopquark.h:197
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_C_phiW
double getNPSMEFT6dtopquark_C_phiW() const
Definition: NPSMEFT6dtopquark.h:194
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_C_tphi
double getNPSMEFT6dtopquark_C_tphi() const
Definition: NPSMEFT6dtopquark.h:77
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_C_tWIm
double getNPSMEFT6dtopquark_C_tWIm() const
Definition: NPSMEFT6dtopquark.h:206
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_C_tW
double getNPSMEFT6dtopquark_C_tW() const
Definition: NPSMEFT6dtopquark.h:65
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_flag_LHC_WG_Basis
double getNPSMEFT6dtopquark_flag_LHC_WG_Basis() const
Definition: NPSMEFT6dtopquark.h:224
NPSMEFT6dtopquark::getNPSMEFT6dtopquark_C_tphiIm
double getNPSMEFT6dtopquark_C_tphiIm() const
Definition: NPSMEFT6dtopquark.h:200
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
sigma_tHj_eta_cut_3_diff_LO_mth::b_min
double b_min
Definition: NPSMEFT6dtopquark.h:4255
sigma_tHj_eta_cut_3_diff_LO_mth::b_max
double b_max
Definition: NPSMEFT6dtopquark.h:4256

Member Data Documentation

◆ b_max

double sigma_tHj_eta_cut_3_diff_LO_mth::b_max = 0.
private

Definition at line 4256 of file NPSMEFT6dtopquark.h.

◆ b_min

double sigma_tHj_eta_cut_3_diff_LO_mth::b_min = 0.
private

Definition at line 4255 of file NPSMEFT6dtopquark.h.

◆ myNPSMEFT6dtopquark

const NPSMEFT6dtopquark& sigma_tHj_eta_cut_3_diff_LO_mth::myNPSMEFT6dtopquark
private

Definition at line 4258 of file NPSMEFT6dtopquark.h.

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