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

#include <TopQuarkObservables.h>

+ Inheritance diagram for sigma_ttz_diff_LO_CMS_190711270:

Detailed Description

Definition at line 461 of file TopQuarkObservables.h.

Public Member Functions

double computeThValue ()
 
double ewgc (const std::string name) const
 
double ewgc (const std::string name, int i, int j) const
 
double ewgc (const std::string name, int i, int j, int k, int l) const
 
 sigma_ttz_diff_LO_CMS_190711270 (const StandardModel &SM_i)
 sigma_ttz_diff_LO 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 TopQuarkObservables mytopobs
 

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_ttz_diff_LO_CMS_190711270()

sigma_ttz_diff_LO_CMS_190711270::sigma_ttz_diff_LO_CMS_190711270 ( const StandardModel &  SM_i)

sigma_ttz_diff_LO constructor.

Definition at line 3020 of file TopQuarkObservables.cpp.

3021: ThObservable(SM_i), mytopobs(static_cast<const NPSMEFTd6General&> (SM))
3022{
3023
3024 setParametersForObservable(make_vector<std::string>() << "SM_sigma_ttz_bin_0_75_CMS_190711270"
3025 << "SM_sigma_ttz_bin_75_150_CMS_190711270" << "SM_sigma_ttz_bin_150_250_CMS_190711270"
3026 << "SM_sigma_ttz_bin_250_500_CMS_190711270");
3027
3028
3029}
NPSMEFTd6General
A model class for new physics in the form of the dimension-six effective Lagrangian.
Definition: NPSMEFTd6General.h:593
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
ThObservable::SM
const StandardModel & SM
A reference to an object of StandardMode class.
Definition: ThObservable.h:121
sigma_ttz_diff_LO_CMS_190711270::mytopobs
const TopQuarkObservables mytopobs
Definition: TopQuarkObservables.h:498

Member Function Documentation

◆ computeThValue()

double sigma_ttz_diff_LO_CMS_190711270::computeThValue ( )
virtual
Returns
The value of Test_direct

Implements ThObservable.

Definition at line 3032 of file TopQuarkObservables.cpp.

3033{
3034
3035 b_min = getBinMin();
3036 b_max = getBinMax();
3037
3038 bool flag_Quadratic= false; //Needs to be properly defined
3039
3040 if(b_min == 0 && b_max == 75){
3041
3042 double SM_sigma_ttz_bin_0_75 = SM.getOptionalParameter("SM_sigma_ttz_bin_0_75_CMS_190711270");
3043 double sigma_ttz_bin_0_75_madgraph = 0.177589; //pb
3044
3045
3046 if(flag_Quadratic){
3047
3048 return 0.;
3049
3050 }
3051 else{
3052
3053 double total = SM_sigma_ttz_bin_0_75 + (-0.031991*ewgc("CG")+-0.016710*ewgc("CHq1R",2,2)+0.016573*ewgc("CHq3R",2,2)
3054 +0.009051*ewgc("CHuR",2,2)+0.003299*ewgc("Cqd8R",2,2,0,0)+0.000362*ewgc("Cqd8R",2,2,1,1)
3055 +0.024362*ewgc("Cqq1R",0,2,2,0)+0.000624*ewgc("Cqq1R",1,2,2,1)+0.119583*ewgc("Cqq3R",0,2,2,0)
3056 +0.012655*ewgc("Cqq3R",1,2,2,1)+0.015687*ewgc("Cqu8R",0,0,2,2)+0.001422*ewgc("Cqu8R",1,1,2,2)
3057 +0.002967*ewgc("Cqu8R",2,2,0,0)+0.000052*ewgc("Cqu8R",2,2,1,1)+-0.000318*ewgc("CuBR",2,2)
3058 +0.001835*ewgc("Cud8R",2,2,0,0)+0.000209*ewgc("Cud8R",2,2,1,1)-0.052919*ewgc("CuGR",2,2)
3059 +0.006972*ewgc("CuuR",0,2,2,0)+0.000145*ewgc("CuuR",1,2,2,1)+-0.000317*ewgc("CuWR",2,2)
3060 +0.000764*ewgc("CHD")+-0.010859*ewgc("CHl3R",0,0)+-0.010859*ewgc("CHl3R",1,1)+-0.000041*ewgc("CHW")
3061 +0.006519*ewgc("CHWB")+0.010767*ewgc("CllR",0,1,1,0)+-0.000148*ewgc("Cqd1R",2,2,0,0)+-0.000064*ewgc("Cqd1R",2,2,1,1)
3062 +-0.002423*ewgc("Cqq1R",0,0,2,2)+-0.000764*ewgc("Cqq1R",1,1,2,2)+0.011282*ewgc("Cqq3R",0,0,2,2)
3063 +0.000883*ewgc("Cqq3R",1,1,2,2)+-0.000077*ewgc("Cqu1R",0,0,2,2)+-0.000116*ewgc("Cqu1R",1,1,2,2)
3064 +0.000149*ewgc("Cqu1R",2,2,0,0)+-0.000054*ewgc("Cqu1R",2,2,1,1)+-0.000248*ewgc("Cud1R",2,2,0,0)
3065 +-0.000072*ewgc("Cud1R",2,2,1,1)+0.000417*ewgc("CuuR",0,0,2,2)+-0.000049*ewgc("CuuR",1,1,2,2))
3066 *(SM_sigma_ttz_bin_0_75/sigma_ttz_bin_0_75_madgraph);
3067
3068 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
3069
3070 return total;
3071
3072 }
3073
3074 } else if(b_min == 75 && b_max == 150){
3075
3076 double SM_sigma_ttz_bin_75_150 = SM.getOptionalParameter("SM_sigma_ttz_bin_75_150_CMS_190711270");
3077 double sigma_ttz_bin_75_150_madgraph = 0.209588; //pb
3078
3079
3080 if(flag_Quadratic){
3081
3082 return 0.;
3083
3084 }
3085 else{
3086
3087 double total = SM_sigma_ttz_bin_75_150 + (-0.051622*ewgc("CG")+-0.021757*ewgc("CHq1R",2,2)+0.021813*ewgc("CHq3R",2,2)
3088 +0.013632*ewgc("CHuR",2,2)+0.004262*ewgc("Cqd8R",2,2,0,0)+0.000521*ewgc("Cqd8R",2,2,1,1)
3089 +0.030335*ewgc("Cqq1R",0,2,2,0)+0.000796*ewgc("Cqq1R",1,2,2,1)+0.129627*ewgc("Cqq3R",0,2,2,0)
3090 +0.012107*ewgc("Cqq3R",1,2,2,1)+0.016607*ewgc("Cqu8R",0,0,2,2)+0.001371*ewgc("Cqu8R",1,1,2,2)
3091 +0.004665*ewgc("Cqu8R",2,2,0,0)+0.000163*ewgc("Cqu8R",2,2,1,1)+-0.000026*ewgc("CuBR",2,2)
3092 +0.002538*ewgc("Cud8R",2,2,0,0)+0.000335*ewgc("Cud8R",2,2,1,1)-0.067116*ewgc("CuGR",2,2)
3093 +0.011624*ewgc("CuuR",0,2,2,0)+0.000337*ewgc("CuuR",1,2,2,1)+-0.000735*ewgc("CuWR",2,2)
3094 +-0.000153*ewgc("CHD")+-0.012754*ewgc("CHl3R",0,0)+-0.012754*ewgc("CHl3R",1,1)+0.000067*ewgc("CHW")
3095 +0.006566*ewgc("CHWB")+0.012822*ewgc("CllR",0,1,1,0)+-0.000189*ewgc("Cqd1R",2,2,0,0)+0.000022*ewgc("Cqd1R",2,2,1,1)
3096 +-0.002576*ewgc("Cqq1R",0,0,2,2)+-0.000592*ewgc("Cqq1R",1,1,2,2)+0.012033*ewgc("Cqq3R",0,0,2,2)
3097 +0.000970*ewgc("Cqq3R",1,1,2,2)+0.000170*ewgc("Cqu1R",0,0,2,2)+-0.000031*ewgc("Cqu1R",1,1,2,2)
3098 +0.000301*ewgc("Cqu1R",2,2,0,0)+0.000043*ewgc("Cqu1R",2,2,1,1)+-0.000125*ewgc("Cud1R",2,2,0,0)
3099 +0.000002*ewgc("Cud1R",2,2,1,1)+0.000863*ewgc("CuuR",0,0,2,2)+0.000061*ewgc("CuuR",1,1,2,2))
3100 *(SM_sigma_ttz_bin_75_150/sigma_ttz_bin_75_150_madgraph);
3101
3102 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
3103
3104 return total;
3105
3106 }
3107
3108 } else if(b_min == 150 && b_max == 250){
3109
3110 double SM_sigma_ttz_bin_150_250 = SM.getOptionalParameter("SM_sigma_ttz_bin_150_250_CMS_190711270");
3111 double sigma_ttz_bin_150_250_madgraph = 0.132915; //pb
3112
3113
3114 if(flag_Quadratic){
3115
3116 return 0.;
3117
3118 }
3119 else{
3120
3121 double total = SM_sigma_ttz_bin_150_250 + (-0.044763*ewgc("CG")-0.014558*ewgc("CHq1R",2,2)+0.014567*ewgc("CHq3R",2,2)
3122 +0.010078*ewgc("CHuR",2,2)+0.003298*ewgc("Cqd8R",2,2,0,0)+0.000314*ewgc("Cqd8R",2,2,1,1)
3123 +0.022847*ewgc("Cqq1R",0,2,2,0)+0.000511*ewgc("Cqq1R",1,2,2,1)+0.084393*ewgc("Cqq3R",0,2,2,0)
3124 +0.006512*ewgc("Cqq3R",1,2,2,1)+0.010724*ewgc("Cqu8R",0,0,2,2)+0.000696*ewgc("Cqu8R",1,1,2,2)
3125 +0.004171*ewgc("Cqu8R",2,2,0,0)+0.000098*ewgc("Cqu8R",2,2,1,1)+0.000118*ewgc("CuBR",2,2)
3126 +0.002011*ewgc("Cud8R",2,2,0,0)+0.000196*ewgc("Cud8R",2,2,1,1)-0.045686*ewgc("CuGR",2,2)
3127 +0.011095*ewgc("CuuR",0,2,2,0)+0.000255*ewgc("CuuR",1,2,2,1)+-0.000811*ewgc("CuWR",2,2)
3128 +-0.000798*ewgc("CHD")+-0.008089*ewgc("CHl3R",0,0)+-0.008089*ewgc("CHl3R",1,1)+0.000024*ewgc("CHW")
3129 +0.003455*ewgc("CHWB")+0.008087*ewgc("CllR",0,1,1,0)+-0.000114*ewgc("Cqd1R",2,2,0,0)+-0.000029*ewgc("Cqd1R",2,2,1,1)
3130 +-0.001109*ewgc("Cqq1R",0,0,2,2)+-0.000367*ewgc("Cqq1R",1,1,2,2)+0.008379*ewgc("Cqq3R",0,0,2,2)
3131 +0.000510*ewgc("Cqq3R",1,1,2,2)+0.000203*ewgc("Cqu1R",0,0,2,2)+-0.000032*ewgc("Cqu1R",1,1,2,2)
3132 +0.000246*ewgc("Cqu1R",2,2,0,0)+-0.000005*ewgc("Cqu1R",2,2,1,1)+-0.000148*ewgc("Cud1R",2,2,0,0)
3133 +-0.000031*ewgc("Cud1R",2,2,1,1)+0.000849*ewgc("CuuR",0,0,2,2)+0.000014*ewgc("CuuR",1,1,2,2))
3134 *(SM_sigma_ttz_bin_150_250/sigma_ttz_bin_150_250_madgraph);
3135
3136 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
3137
3138 return total;
3139
3140 }
3141
3142 } else if(b_min == 250 && b_max == 500){
3143
3144 double SM_sigma_ttz_bin_250_500 = SM.getOptionalParameter("SM_sigma_ttz_bin_250_500_CMS_190711270");
3145 double sigma_ttz_bin_250_500_madgraph = 0.067464; //pb
3146
3147
3148 if(flag_Quadratic){
3149
3150 return 0.;
3151
3152 }
3153 else{
3154
3155 double total = SM_sigma_ttz_bin_250_500 + (-0.031862*ewgc("CG")+-0.007411*ewgc("CHq1R",2,2)+0.007439*ewgc("CHq3R",2,2)
3156 +0.005451*ewgc("CHuR",2,2)+0.002779*ewgc("Cqd8R",2,2,0,0)+0.000235*ewgc("Cqd8R",2,2,1,1)
3157 +0.019482*ewgc("Cqq1R",0,2,2,0)+0.000381*ewgc("Cqq1R",1,2,2,1)+0.060750*ewgc("Cqq3R",0,2,2,0)
3158 +0.003621*ewgc("Cqq3R",1,2,2,1)+0.007854*ewgc("Cqu8R",0,0,2,2)+0.000404*ewgc("Cqu8R",1,1,2,2)
3159 +0.003938*ewgc("Cqu8R",2,2,0,0)+0.000089*ewgc("Cqu8R",2,2,1,1)+0.000042*ewgc("CuBR",2,2)
3160 +0.001798*ewgc("Cud8R",2,2,0,0)+0.000162*ewgc("Cud8R",2,2,1,1)-0.026610*ewgc("CuGR",2,2)
3161 +0.011353*ewgc("CuuR",0,2,2,0)+0.000235*ewgc("CuuR",1,2,2,1)+-0.000517*ewgc("CuWR",2,2)
3162 +-0.000594*ewgc("CHD")+-0.004099*ewgc("CHl3R",0,0)+-0.004099*ewgc("CHl3R",1,1)+0.000046*ewgc("CHW")
3163 +0.001552*ewgc("CHWB")+0.004130*ewgc("CllR",0,1,1,0)+-0.000078*ewgc("Cqd1R",2,2,0,0)+0.000009*ewgc("Cqd1R",2,2,1,1)
3164 +-0.000036*ewgc("Cqq1R",0,0,2,2)+-0.000177*ewgc("Cqq1R",1,1,2,2)+0.006452*ewgc("Cqq3R",0,0,2,2)
3165 +0.000324*ewgc("Cqq3R",1,1,2,2)+0.000333*ewgc("Cqu1R",0,0,2,2)+0.000016*ewgc("Cqu1R",1,1,2,2)
3166 +0.000275*ewgc("Cqu1R",2,2,0,0)+0.000020*ewgc("Cqu1R",2,2,1,1)+-0.000138*ewgc("Cud1R",2,2,0,0)
3167 +0.000004*ewgc("Cud1R",2,2,1,1)+0.000886*ewgc("CuuR",0,0,2,2)+0.000029*ewgc("CuuR",1,1,2,2))
3168 *(SM_sigma_ttz_bin_250_500/sigma_ttz_bin_250_500_madgraph);
3169
3170 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
3171
3172 return total;
3173
3174 }
3175
3176 } else throw std::runtime_error("wrong bin choice in sigma_ttz_diff_LO_CMS_190711270");
3177
3178}
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
sigma_ttz_diff_LO_CMS_190711270::b_min
double b_min
Definition: TopQuarkObservables.h:494
sigma_ttz_diff_LO_CMS_190711270::ewgc
double ewgc(const std::string name) const
Definition: TopQuarkObservables.h:477
sigma_ttz_diff_LO_CMS_190711270::b_max
double b_max
Definition: TopQuarkObservables.h:495

◆ ewgc() [1/3]

double sigma_ttz_diff_LO_CMS_190711270::ewgc ( const std::string  name) const
inline

Definition at line 477 of file TopQuarkObservables.h.

478 {
479 return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
480 }
NPSMEFTd6General::getSMEFTCoeffEW
double getSMEFTCoeffEW(const std::string name) const
Definition: NPSMEFTd6General.h:6790
TopQuarkObservables::GetmyNPSMEFTd6General
const NPSMEFTd6General & GetmyNPSMEFTd6General() const
Definition: TopQuarkObservables.h:27

◆ ewgc() [2/3]

double sigma_ttz_diff_LO_CMS_190711270::ewgc ( const std::string  name,
int  i,
int  j 
) const
inline

Definition at line 482 of file TopQuarkObservables.h.

483 {
484 return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
485 }

◆ ewgc() [3/3]

double sigma_ttz_diff_LO_CMS_190711270::ewgc ( const std::string  name,
int  i,
int  j,
int  k,
int  l 
) const
inline

Definition at line 487 of file TopQuarkObservables.h.

488 {
489 return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
490 }

Member Data Documentation

◆ b_max

double sigma_ttz_diff_LO_CMS_190711270::b_max = 0.
private

Definition at line 495 of file TopQuarkObservables.h.

◆ b_min

double sigma_ttz_diff_LO_CMS_190711270::b_min = 0.
private

Definition at line 494 of file TopQuarkObservables.h.

◆ mytopobs

const TopQuarkObservables sigma_ttz_diff_LO_CMS_190711270::mytopobs
private

Definition at line 498 of file TopQuarkObservables.h.

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