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

#include <NPSMEFT6dtopquark.h>

+ Inheritance diagram for b4_ttH_eta_cut_3_LO:

Detailed Description

Definition at line 3891 of file NPSMEFT6dtopquark.h.

Public Member Functions

 b4_ttH_eta_cut_3_LO (const StandardModel &SM_i)
 b4_ttH_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 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

◆ b4_ttH_eta_cut_3_LO()

b4_ttH_eta_cut_3_LO::b4_ttH_eta_cut_3_LO ( const StandardModel SM_i)

b4_ttH_LO constructor.

Definition at line 14192 of file NPSMEFT6dtopquark.cpp.

14193: ThObservable(SM_i),myNPSMEFT6dtopquark(static_cast<const NPSMEFT6dtopquark&> (SM_i))
14194{
14195
14196 setParametersForObservable(make_vector<std::string>() << "SM_b4_ttH_eta_cut_3_LO_bin_m1_m0p5" << "SM_b4_ttH_eta_cut_3_LO_bin_m0p5_0"
14197 << "SM_b4_ttH_eta_cut_3_LO_bin_0_0p5" << "SM_b4_ttH_eta_cut_3_LO_bin_0p5_1");
14198
14199}
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 b4_ttH_eta_cut_3_LO::computeThValue ( )
virtual
Returns
The value of b4_ttH_LO

Implements ThObservable.

Definition at line 14201 of file NPSMEFT6dtopquark.cpp.

14202{
14203
14204 b_min = getBinMin();
14205 b_max = getBinMax();
14206
14207
14210
14211
14214 //double cHW = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_phiW();
14215 //double cHWtil = myNPSMEFT6dtopquark.getNPSMEFT6dtopquark_C_phiWtil();
14222 double ctBRe = 0.;
14223 double ctBIm = 0.;
14224
14225
14226
14227
14228 //Expressions are written in the basis of SMEFTsim
14229 //let's change to the basis of dim6top (same as
14230 //SMEFT@NLO except for the gs on CtG). These
14231 //models define the covariance derivative
14232 //with a different convention
14233
14234 if(flag_LHC_WG_Basis){
14237 double sw2 = 0.22305;
14238 double cw2 = 1 - sw2;
14239 double tw2 = sw2/cw2;
14240 double sw = sqrt(sw2);
14241 double tw = sqrt(tw2);
14242 ctBRe = ctZRe/sw - ctWRe/tw;
14243 ctBIm = ctZIm/sw - ctWIm/tw;
14244 ctWRe = -ctWRe;
14245 ctWIm = -ctWIm;
14246 ctGRe = -ctGRe;
14247 ctGIm = -ctGIm;
14248 }
14249 else{
14252 }
14253
14254
14255
14256 if(b_min == -1 && b_max == -0.5){
14257
14258 double SM_b4_tth_bin_m1_m0p5 = SM.getOptionalParameter("SM_b4_ttH_eta_cut_3_LO_bin_m1_m0p5");
14259 double b4_tth_bin_m1_m0p5_madgraph_LO = 0.9684430541146044 ;//pb
14260 double total;
14261
14262
14263
14264 if(flag_Quadratic){
14265 total = SM_b4_tth_bin_m1_m0p5 + (
14266 +0.8734320303016335 * cHG + -0.0001508759372294 * cHGtil
14267 + -0.0025829369337817 * ctGIm + -1.055582442402413 * ctGRe
14268 + 0.0004214757799254 * ctHIm + -0.1194279109508295 * ctHRe
14269 + -0.0013695034263697 * ctWIm + -0.00093654123036 * ctWRe
14270 + 0.8018553830032952 * cHG * cHG + 0.7689420372074605 * cHGtil * cHGtil
14271 + 0.9073938708776672 * ctGIm * ctGIm + 0.9271480966621598 * ctGRe * ctGRe
14272 + 0.0048101219104801 * ctHIm * ctHIm + 0.0042590787439545 * ctHRe * ctHRe
14273 + 0.1412202101822702 * ctWIm * ctWIm + 0.020697839481579 * ctWRe * ctWRe
14274 + -0.0024331673698921 * cHG * cHGtil + -0.0073214591015668 * cHG * ctGIm
14275 + -1.0644997356352084 * cHG * ctGRe + -0.0004197491449199 * cHG * ctHIm
14276 + -0.0764189877477869 * cHG * ctHRe + 0.0018086744290113 * cHG * ctWIm
14277 + 0.0004350964819894 * cHG * ctWRe + -0.7872729954294924 * cHGtil * ctGIm
14278 + -0.0073411906122446 * cHGtil * ctGRe + 0.0802782937008257 * cHGtil * ctHIm
14279 + -0.0003548796572252 * cHGtil * ctHRe + -0.0046932204616949 * cHGtil * ctWIm
14280 + -0.0018702048906836 * cHGtil * ctWRe + 0.0034097982029186 * ctGIm * ctGRe
14281 + -0.0336097850578753 * ctGIm * ctHIm + 0.0006320373106642 * ctGIm * ctHRe
14282 + -0.0017862474813037 * ctGIm * ctWIm + -0.002642626076949 * ctGIm * ctWRe
14283 + -0.0002110204108775 * ctGRe * ctHIm + 0.072489088712194 * ctGRe * ctHRe
14284 + 0.0376030033326701 * ctGRe * ctWIm + 0.003324753924339 * ctGRe * ctWRe
14285 + -0.0001726324778062 * ctHIm * ctHRe + -0.0002431322793271 * ctHIm * ctWIm
14286 + -0.0004548330423232 * ctHIm * ctWRe + 8.341952703337652e-05 * ctHRe * ctWIm
14287 + -0.0009474934437639 * ctHRe * ctWRe + -0.0022153955316867 * ctWIm * ctWRe
14288 + -0.0008641840283073743 * ctBRe + 0.00202906905524447 * ctBRe * ctBRe+
14289 -0.00014275867242054263 * ctBIm + 0.0005523818643201062 * ctBIm * ctBIm
14290 )*(SM_b4_tth_bin_m1_m0p5/b4_tth_bin_m1_m0p5_madgraph_LO);
14291 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
14292 return total;
14293
14294 }
14295 else{
14296 total = SM_b4_tth_bin_m1_m0p5 + (
14297 +0.8734320303016335 * cHG + -0.0001508759372294 * cHGtil
14298 + -0.0025829369337817 * ctGIm + -1.055582442402413 * ctGRe
14299 + 0.0004214757799254 * ctHIm + -0.1194279109508295 * ctHRe
14300 + -0.0013695034263697 * ctWIm + -0.00093654123036 * ctWRe
14301 -0.0008641840283073743 * ctBRe+ -0.00014275867242054263 * ctBIm
14302 )*(SM_b4_tth_bin_m1_m0p5/b4_tth_bin_m1_m0p5_madgraph_LO);
14303 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
14304 return total;
14305 }
14306
14307 } else if(b_min == -0.5 && b_max == 0.){
14308
14309 double SM_b4_tth_bin_m0p5_0 = SM.getOptionalParameter("SM_b4_ttH_eta_cut_3_LO_bin_m0p5_0");
14310 double b4_tth_bin_m0p5_0_madgraph_LO = 1.86995179928823 ;
14311 double total;
14312
14313
14314 if(flag_Quadratic){
14315 total = SM_b4_tth_bin_m0p5_0 + (
14316 +2.5082006692574814 * cHG + -0.0072777625948975 * cHGtil
14317 + 0.0049523319413338 * ctGIm + -2.7100713913002883 * ctGRe
14318 + 4.601097731871118e-05 * ctHIm + -0.226945305824754 * ctHRe
14319 + -0.0054671458335096 * ctWIm + -0.0012224407671053 * ctWRe
14320 + 1.978937263465199 * cHG * cHG + 1.9077219925760047 * cHGtil * cHGtil
14321 + 2.4711095051488763 * ctGIm * ctGIm + 2.4826970521659635 * ctGRe * ctGRe
14322 + 0.0036758272860543 * ctHIm * ctHIm + 0.0075525314558487 * ctHRe * ctHRe
14323 + 0.4526238518110063 * ctWIm * ctWIm + 0.0647774338725194 * ctWRe * ctWRe
14324 + -0.0096779615662583 * cHG * cHGtil + 0.0179862640902231 * cHG * ctGIm
14325 + -2.2076036907173187 * cHG * ctGRe + 0.0002380301652577 * cHG * ctHIm
14326 + -0.2019167669157224 * cHG * ctHRe + -0.0043505464946114 * cHG * ctWIm
14327 + 0.000680424002248 * cHG * ctWRe + -1.4044727589128492 * cHGtil * ctGIm
14328 + -0.0254389459522563 * cHGtil * ctGRe + 0.1013782106873143 * cHGtil * ctHIm
14329 + -0.0021354305980373 * cHGtil * ctHRe + 0.0024387884638665 * cHGtil * ctWIm
14330 + -0.0043688706516078 * cHGtil * ctWRe + 0.0007127628008351 * ctGIm * ctGRe
14331 + 0.0045711906404675 * ctGIm * ctHIm + -0.0009092135533148 * ctGIm * ctHRe
14332 + -0.0300846877664431 * ctGIm * ctWIm + 0.0057492860309402 * ctGIm * ctWRe
14333 + 1.334308996889244e-06 * ctGRe * ctHIm + 0.1728347517089635 * ctGRe * ctHRe
14334 + 0.0040621684242522 * ctGRe * ctWIm + -0.0040069857939023 * ctGRe * ctWRe
14335 + 3.748064076936064e-06 * ctHIm * ctHRe + -0.0011467866836193 * ctHIm * ctWIm
14336 + -0.0001101854248126 * ctHIm * ctWRe + 0.0002385479337846 * ctHRe * ctWIm
14337 + -0.0012914493313962 * ctHRe * ctWRe + 0.0092816063216226 * ctWIm * ctWRe
14338 -0.003156266259222204 * ctBRe + 0.007933633377850108 * ctBRe * ctBRe
14339 + 0.00013663796676958995 * ctBIm + 0.0036409627044012294 * ctBIm * ctBIm
14340 )*(SM_b4_tth_bin_m0p5_0/b4_tth_bin_m0p5_0_madgraph_LO);
14341 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
14342 return total;
14343
14344 }
14345 else{
14346 total = SM_b4_tth_bin_m0p5_0 + (
14347 +2.5082006692574814 * cHG + -0.0072777625948975 * cHGtil
14348 + 0.0049523319413338 * ctGIm + -2.7100713913002883 * ctGRe
14349 + 4.601097731871118e-05 * ctHIm + -0.226945305824754 * ctHRe
14350 + -0.0054671458335096 * ctWIm + -0.0012224407671053 * ctWRe
14351 -0.003156266259222204 * ctBRe+ 0.00013663796676958995 * ctBIm
14352 )*(SM_b4_tth_bin_m0p5_0/b4_tth_bin_m0p5_0_madgraph_LO);
14353 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
14354 return total;
14355 }
14356
14357 } else if(b_min == 0. && b_max == 0.5){
14358
14359 double SM_b4_tth_bin_0_0p5 = SM.getOptionalParameter("SM_b4_ttH_eta_cut_3_LO_bin_0_0p5");
14360 double b4_tth_bin_0_0p5_madgraph_LO = 2.397551646237114 ;
14361 double total;
14362
14363
14364 if(flag_Quadratic){
14365 total = SM_b4_tth_bin_0_0p5 + (
14366 +3.1202407629736086 * cHG + 0.0025837393308619 * cHGtil
14367 + -0.017115553372471 * ctGIm + -3.4070340346935453 * ctGRe
14368 + -5.958601608919978e-05 * ctHIm + -0.2926733879690542 * ctHRe
14369 + 0.0089327792324831 * ctWIm + -0.0056071336952001 * ctWRe
14370 + 2.432408683805834 * cHG * cHG + 2.3631373800720485 * cHGtil * cHGtil
14371 + 3.102465632717866 * ctGIm * ctGIm + 3.079177599282639 * ctGRe * ctGRe
14372 + 0.0032842988136334 * ctHIm * ctHIm + 0.0096392105518944 * ctHRe * ctHRe
14373 + 0.4518085242347641 * ctWIm * ctWIm + 0.074408951180635 * ctWRe * ctWRe
14374 + 0.0082715886036484 * cHG * cHGtil + -0.0214310465869187 * cHG * ctGIm
14375 + -2.6233186760313467 * cHG * ctGRe + -0.0021758311388906 * cHG * ctHIm
14376 + -0.2435661927057007 * cHG * ctHRe + 0.004265298306129 * cHG * ctWIm
14377 + -0.0009206781113095 * cHG * ctWRe + -1.646242954978188 * cHGtil * ctGIm
14378 + 0.0137368063220981 * cHGtil * ctGRe + 0.1198126060642813 * cHGtil * ctHIm
14379 + 0.0011944815868255 * cHGtil * ctHRe + 0.0045209225044467 * cHGtil * ctWIm
14380 + 0.0008273617570505 * cHGtil * ctWRe + -0.0160790288621262 * ctGIm * ctGRe
14381 + 0.0174699207340045 * ctGIm * ctHIm + -0.0002657591051274 * ctGIm * ctHRe
14382 + 0.0274859955582831 * ctGIm * ctWIm + -0.0001721273656364 * ctGIm * ctWRe
14383 + -0.0002578493077056 * ctGRe * ctHIm + 0.2173921499818756 * ctGRe * ctHRe
14384 + -0.0162476469337616 * ctGRe * ctWIm + -0.0010760912144503 * ctGRe * ctWRe
14385 + -2.382554683892142e-05 * ctHIm * ctHRe + -0.001746144437731 * ctHIm * ctWIm
14386 + 0.0005194431818834 * ctHIm * ctWRe + -0.0003018732773732 * ctHRe * ctWIm
14387 + -0.0026562552988104 * ctHRe * ctWRe + -0.0015871656098833 * ctWIm * ctWRe
14388 -0.005086946128797884 * ctBRe + 0.0117366493701683 * ctBRe * ctBRe+
14389 -0.0001602421534973586 * ctBIm + 0.005570109834694012 * ctBIm * ctBIm
14390 )*(SM_b4_tth_bin_0_0p5/b4_tth_bin_0_0p5_madgraph_LO);
14391 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
14392 return total;
14393
14394 }
14395 else{
14396 total = SM_b4_tth_bin_0_0p5 + (
14397 +3.1202407629736086 * cHG + 0.0025837393308619 * cHGtil
14398 + -0.017115553372471 * ctGIm + -3.4070340346935453 * ctGRe
14399 + -5.958601608919978e-05 * ctHIm + -0.2926733879690542 * ctHRe
14400 + 0.0089327792324831 * ctWIm + -0.0056071336952001 * ctWRe
14401 -0.005086946128797884 * ctBRe+ -0.0001602421534973586 * ctBIm
14402 )*(SM_b4_tth_bin_0_0p5/b4_tth_bin_0_0p5_madgraph_LO);
14403 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
14404 return total;
14405 }
14406
14407 } else if(b_min == 0.5 && b_max == 1.){
14408
14409 double SM_b4_tth_bin_0p5_1 = SM.getOptionalParameter("SM_b4_ttH_eta_cut_3_LO_bin_0p5_1");
14410 double b4_tth_bin_0p5_1_madgraph_LO = 3.870751347545547;//pb
14411 double total;
14412
14413
14414 if(flag_Quadratic){
14415 total = SM_b4_tth_bin_0p5_1 + (
14416 +3.3485789752862267 * cHG + 0.0029440698670099 * cHGtil
14417 + 0.0364904479907755 * ctGIm + -4.263214480363309 * ctGRe
14418 + -0.0002786851476225 * ctHIm + -0.466940753763355 * ctHRe
14419 + 0.0013794115989795 * ctWIm + -0.0215433485065755 * ctWRe
14420 + 2.063155473027622 * cHG * cHG + 2.0219556230466966 * cHGtil * cHGtil
14421 + 2.8206281712618573 * ctGIm * ctGIm + 2.9967795865270896 * ctGRe * ctGRe
14422 + 0.003721402922715 * ctHIm * ctHIm + 0.0140811884327776 * ctHRe * ctHRe
14423 + 0.149365254520435 * ctWIm * ctWIm + 0.0518640741018987 * ctWRe * ctWRe
14424 + 0.0066796844867344 * cHG * cHGtil + 0.0268556472722223 * cHG * ctGIm
14425 + -2.728096152321642 * cHG * ctGRe + 0.0005190885940814 * cHG * ctHIm
14426 + -0.2087011098917203 * cHG * ctHRe + 0.0019505183062117 * cHG * ctWIm
14427 + 0.0015702917342409 * cHG * ctWRe + -1.636616624499878 * cHGtil * ctGIm
14428 + 0.0147962038200352 * cHGtil * ctGRe + 0.1606758436320679 * cHGtil * ctHIm
14429 + -0.0001985789407428 * cHGtil * ctHRe + 0.0010503720540124 * cHGtil * ctWIm
14430 + 0.0068225479068752 * cHGtil * ctWRe + -0.0037893213650285 * ctGIm * ctGRe
14431 + 0.0098655094771906 * ctGIm * ctHIm + 0.0010406217628624 * ctGIm * ctHRe
14432 + 0.0002576492735923 * ctGIm * ctWIm + 0.0021459067606248 * ctGIm * ctWRe
14433 + -0.0001812107388418 * ctGRe * ctHIm + 0.3161364830453751 * ctGRe * ctHRe
14434 + -0.031306745077561 * ctGRe * ctWIm + 0.0050307413068803 * ctGRe * ctWRe
14435 + 7.57357426761085e-05 * ctHIm * ctHRe + -0.0007680728871719 * ctHIm * ctWIm
14436 + 0.0001934901906697 * ctHIm * ctWRe + -4.692781471730093e-05 * ctHRe * ctWIm
14437 + -0.0002216373963714 * ctHRe * ctWRe + -0.002050876442346 * ctWIm * ctWRe
14438 -0.006496243526542184 * ctBRe + 0.014293810877612556 * ctBRe * ctBRe
14439 + 0.00024059095088696525 * ctBIm + 0.004504803676770738 * ctBIm * ctBIm
14440 )*(SM_b4_tth_bin_0p5_1/b4_tth_bin_0p5_1_madgraph_LO);
14441 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
14442 return total;
14443
14444 }
14445 else{
14446 total = SM_b4_tth_bin_0p5_1 + (
14447 +3.3485789752862267 * cHG + 0.0029440698670099 * cHGtil
14448 + 0.0364904479907755 * ctGIm + -4.263214480363309 * ctGRe
14449 + -0.0002786851476225 * ctHIm + -0.466940753763355 * ctHRe
14450 + 0.0013794115989795 * ctWIm + -0.0215433485065755 * ctWRe
14451 -0.006496243526542184 * ctBRe+ 0.00024059095088696525 * ctBIm
14452 )*(SM_b4_tth_bin_0p5_1/b4_tth_bin_0p5_1_madgraph_LO);
14453 //if (total < 0) return std::numeric_limits<double>::quiet_NaN();
14454 return total;
14455 }
14456
14457 } else{
14458 throw std::runtime_error("\nERROR: Please specify a correct binning for b4_tth_eta_cut_3_LO. \n");
14459 }
14460
14461}
double getNPSMEFT6dtopquark_flag_Quadratic() const
double getNPSMEFT6dtopquark_C_tBIm() const
double getNPSMEFT6dtopquark_C_tG() const
double getNPSMEFT6dtopquark_C_phiG() const
double getNPSMEFT6dtopquark_C_phiGtil() const
double getNPSMEFT6dtopquark_C_tphi() const
double getNPSMEFT6dtopquark_C_tZIm() const
double getNPSMEFT6dtopquark_C_tB() const
double getNPSMEFT6dtopquark_C_tZ() const
double getNPSMEFT6dtopquark_C_tWIm() const
double getNPSMEFT6dtopquark_C_tW() const
double getNPSMEFT6dtopquark_C_tGIm() 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 b4_ttH_eta_cut_3_LO::b_max = 0.
private

Definition at line 3907 of file NPSMEFT6dtopquark.h.

◆ b_min

double b4_ttH_eta_cut_3_LO::b_min = 0.
private

Definition at line 3906 of file NPSMEFT6dtopquark.h.

◆ myNPSMEFT6dtopquark

const NPSMEFT6dtopquark& b4_ttH_eta_cut_3_LO::myNPSMEFT6dtopquark
private

Definition at line 3910 of file NPSMEFT6dtopquark.h.


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