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

#include <TopQuarkObservables.h>

+ Inheritance diagram for sigma_tta_diff_LO_CMS_dilepton_220107301:

Detailed Description

Definition at line 728 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_tta_diff_LO_CMS_dilepton_220107301 (const StandardModel &SM_i)
 FL 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_tta_diff_LO_CMS_dilepton_220107301()

sigma_tta_diff_LO_CMS_dilepton_220107301::sigma_tta_diff_LO_CMS_dilepton_220107301 ( const StandardModel &  SM_i)

FL constructor.

Definition at line 2727 of file TopQuarkObservables.cpp.

2728: ThObservable(SM_i), mytopobs(static_cast<const NPSMEFTd6General&> (SM))
2729{
2730
2731 setParametersForObservable(make_vector<std::string>() << "SM_sigma_tta_bin_20_35_CMS_dilepton" << "SM_sigma_tta_bin_35_50_CMS_dilepton"
2732 << "SM_sigma_tta_bin_50_70_CMS_dilepton" << "SM_sigma_tta_bin_70_130_CMS_dilepton"
2733 << "SM_sigma_tta_bin_130_200_CMS_dilepton" << "SM_sigma_tta_bin_200_300_CMS_dilepton");
2734
2735
2736}
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_tta_diff_LO_CMS_dilepton_220107301::mytopobs
const TopQuarkObservables mytopobs
Definition: TopQuarkObservables.h:765

Member Function Documentation

◆ computeThValue()

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

Implements ThObservable.

Definition at line 2739 of file TopQuarkObservables.cpp.

2740{
2741
2742 b_min = getBinMin();
2743 b_max = getBinMax();
2744
2745 bool flag_Quadratic= false; //Needs to be properly defined
2746
2747
2748 if(b_min == 20 && b_max == 35){
2749
2750 double SM_sigma_tta_bin_20_35 = SM.getOptionalParameter("SM_sigma_tta_bin_20_35_CMS_dilepton");
2751 double sigma_tta_bin_20_35_madgraph = 0.583675; //pb
2752
2753
2754 if(flag_Quadratic){
2755
2756 return 0.;
2757
2758 }
2759 else{
2760
2761 //There is no dependence on ewgc("Cqd8R",2,2,2,2) neither on ewgc("Cud8R",2,2,2,2) with the precision considered
2762
2763 double total = SM_sigma_tta_bin_20_35 +(-0.034106*ewgc("CG")+0.004794*ewgc("Cqd8R",2,2,0,0)
2764 +0.000605*ewgc("Cqd8R",2,2,1,1)+0.082493*ewgc("Cqq1R",0,2,2,0)
2765 +0.002529*ewgc("Cqq1R",1,2,2,1)+0.146282*ewgc("Cqq3R",0,2,2,0)+0.010782*ewgc("Cqq3R",1,2,2,1)
2766 +0.036042*ewgc("Cqu8R",0,0,2,2)+0.001578*ewgc("Cqu8R",1,1,2,2)+0.031190*ewgc("Cqu8R",2,2,0,0)
2767 +0.000906*ewgc("Cqu8R",2,2,1,1)+-0.007204*ewgc("CuBR",2,2)+0.008516*ewgc("Cud8R",2,2,0,0)
2768 +0.001057*ewgc("Cud8R",2,2,1,1)+-0.154829*ewgc("CuGR",2,2)
2769 +0.080748*ewgc("CuuR",0,2,2,0)+0.002498*ewgc("CuuR",1,2,2,1)+-0.005006*ewgc("CuWR",2,2)
2770 +-0.061811*ewgc("CHD")+-0.035571*ewgc("CHl3R",0,0)+-0.035571*ewgc("CHl3R",1,1)+-0.132541*ewgc("CHWB")
2771 +0.035506*ewgc("CllR",0,1,1,0)+-0.000081*ewgc("Cqd1R",2,2,0,0)+-0.000027*ewgc("Cqd1R",2,2,1,1)
2772 +0.006399*ewgc("Cqq1R",0,0,2,2)+-0.000210*ewgc("Cqq1R",1,1,2,2)+0.016394*ewgc("Cqq3R",0,0,2,2)
2773 +0.000986*ewgc("Cqq3R",1,1,2,2)+0.002978*ewgc("Cqu1R",0,0,2,2)+0.000069*ewgc("Cqu1R",1,1,2,2)
2774 +0.002181*ewgc("Cqu1R",2,2,0,0)+0.000072*ewgc("Cqu1R",2,2,1,1)+-0.000689*ewgc("Cud1R",2,2,0,0)
2775 +-0.000057*ewgc("Cud1R",2,2,1,1)+0.005823*ewgc("CuuR",0,0,2,2)+0.000189*ewgc("CuuR",1,1,2,2))
2776 *(SM_sigma_tta_bin_20_35/sigma_tta_bin_20_35_madgraph);
2777
2778 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
2779
2780 return total;
2781
2782 }
2783 } else if(b_min == 35 && b_max == 50){
2784
2785
2786 double SM_sigma_tta_bin_35_50 = SM.getOptionalParameter("SM_sigma_tta_bin_35_50_CMS_dilepton");
2787 double sigma_tta_bin_35_50_madgraph = 0.305601; //pb
2788
2789
2790 if(flag_Quadratic){
2791
2792 return 0.;
2793
2794 }
2795 else{
2796
2797 double total = SM_sigma_tta_bin_35_50 +(-0.018555*ewgc("CG")+0.002547*ewgc("Cqd8R",2,2,0,0)
2798 +0.000295*ewgc("Cqd8R",2,2,1,1)+0.040351*ewgc("Cqq1R",0,2,2,0)
2799 +0.001256*ewgc("Cqq1R",1,2,2,1)+0.075386*ewgc("Cqq3R",0,2,2,0)+0.005401*ewgc("Cqq3R",1,2,2,1)
2800 +0.018788*ewgc("Cqu8R",0,0,2,2)+0.000807*ewgc("Cqu8R",1,1,2,2)+0.016262*ewgc("Cqu8R",2,2,0,0)
2801 +0.000470*ewgc("Cqu8R",2,2,1,1)+-0.006535*ewgc("CuBR",2,2)+0.004622*ewgc("Cud8R",2,2,0,0)
2802 +0.000557*ewgc("Cud8R",2,2,1,1)-0.080445*ewgc("CuGR",2,2)
2803 +0.039496*ewgc("CuuR",0,2,2,0)+0.001217*ewgc("CuuR",1,2,2,1)+-0.004069*ewgc("CuWR",2,2)
2804 +-0.032332*ewgc("CHD")+-0.018575*ewgc("CHl3R",0,0)+-0.018575*ewgc("CHl3R",1,1)+-0.069327*ewgc("CHWB")
2805 +0.018597*ewgc("CllR",0,1,1,0)+-0.000057*ewgc("Cqd1R",2,2,0,0)+-0.000002*ewgc("Cqd1R",2,2,1,1)
2806 +0.003004*ewgc("Cqq1R",0,0,2,2)+-0.000132*ewgc("Cqq1R",1,1,2,2)+0.008326*ewgc("Cqq3R",0,0,2,2)
2807 +0.000523*ewgc("Cqq3R",1,1,2,2)+0.001386*ewgc("Cqu1R",0,0,2,2)+0.000057*ewgc("Cqu1R",1,1,2,2)
2808 +0.001019*ewgc("Cqu1R",2,2,0,0)+0.000055*ewgc("Cqu1R",2,2,1,1)+-0.000289*ewgc("Cud1R",2,2,0,0)
2809 +-0.000021*ewgc("Cud1R",2,2,1,1)+0.002900*ewgc("CuuR",0,0,2,2)+0.000133*ewgc("CuuR",1,1,2,2))
2810 *(SM_sigma_tta_bin_35_50/sigma_tta_bin_35_50_madgraph);
2811
2812 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
2813
2814 return total;
2815
2816
2817 }
2818 } else if(b_min == 50 && b_max == 70){
2819
2820
2821 double SM_sigma_tta_bin_50_70 = SM.getOptionalParameter("SM_sigma_tta_bin_50_70_CMS_dilepton");
2822 double sigma_tta_bin_50_70_madgraph = 0.238500; //pb
2823
2824
2825 if(flag_Quadratic){
2826
2827 return 0.;
2828
2829 }
2830 else{
2831
2832 double total = SM_sigma_tta_bin_50_70 +(-0.014859*ewgc("CG")+0.001934*ewgc("Cqd8R",2,2,0,0)
2833 +0.000192*ewgc("Cqd8R",2,2,1,1)+
2834 +0.029858*ewgc("Cqq1R",0,2,2,0)+0.000825*ewgc("Cqq1R",1,2,2,1)
2835 +0.057705*ewgc("Cqq3R",0,2,2,0)+0.004114*ewgc("Cqq3R",1,2,2,1)
2836 +0.014644*ewgc("Cqu8R",0,0,2,2)+0.000602*ewgc("Cqu8R",1,1,2,2)
2837 +0.012747*ewgc("Cqu8R",2,2,0,0)+-0.007537*ewgc("CuBR",2,2)
2838 +0.003692*ewgc("Cud8R",2,2,0,0)+0.000409*ewgc("Cud8R",2,2,1,1)
2839 -0.061847*ewgc("CuGR",2,2)
2840 +0.029209*ewgc("CuuR",0,2,2,0)+0.000813*ewgc("CuuR",1,2,2,1)
2841 +-0.004460*ewgc("CuWR",2,2)
2842 +-0.025253*ewgc("CHD")+-0.014510*ewgc("CHl3R",0,0)+-0.014510*ewgc("CHl3R",1,1)+-0.054091*ewgc("CHWB")
2843 +0.014521*ewgc("CllR",0,1,1,0)+-0.000031*ewgc("Cqd1R",2,2,0,0)+-0.000008*ewgc("Cqd1R",2,2,1,1)
2844 +0.001996*ewgc("Cqq1R",0,0,2,2)+-0.000150*ewgc("Cqq1R",1,1,2,2)+0.006394*ewgc("Cqq3R",0,0,2,2)
2845 +0.000416*ewgc("Cqq3R",1,1,2,2)+0.001142*ewgc("Cqu1R",0,0,2,2)+0.000027*ewgc("Cqu1R",1,1,2,2)
2846 +0.000807*ewgc("Cqu1R",2,2,0,0)+0.000037*ewgc("Cqu1R",2,2,1,1)+-0.000320*ewgc("Cud1R",2,2,0,0)
2847 +-0.000011*ewgc("Cud1R",2,2,1,1)+0.002163*ewgc("CuuR",0,0,2,2)+0.000065*ewgc("CuuR",1,1,2,2))
2848 *(SM_sigma_tta_bin_50_70/sigma_tta_bin_50_70_madgraph);
2849
2850 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
2851
2852 return total;
2853
2854 }
2855 } else if(b_min == 70 && b_max == 130){
2856
2857
2858 double SM_sigma_tta_bin_70_130 = SM.getOptionalParameter("SM_sigma_tta_bin_70_130_CMS_dilepton");
2859 double sigma_tta_bin_70_100_madgraph = 0.198831; //pb
2860 double sigma_tta_bin_100_130_madgraph =0.110866; //pb
2861 double sigma_tta_bin_70_130_madgraph =sigma_tta_bin_70_100_madgraph+sigma_tta_bin_100_130_madgraph; //pb
2862
2863 if(flag_Quadratic){
2864
2865 return 0.;
2866
2867 }
2868 else{
2869
2870 double sigma_tta_bin_70_100_NP = (-0.013285*ewgc("CG")+0.001630*ewgc("Cqd8R",2,2,0,0)
2871 +0.000250*ewgc("Cqd8R",2,2,1,1)
2872 +0.023767*ewgc("Cqq1R",0,2,2,0)+0.000676*ewgc("Cqq1R",1,2,2,1)
2873 +0.047671*ewgc("Cqq3R",0,2,2,0)+0.003336*ewgc("Cqq3R",1,2,2,1)
2874 +0.012379*ewgc("Cqu8R",0,0,2,2)+0.000535*ewgc("Cqu8R",1,1,2,2)
2875 +0.010765*ewgc("Cqu8R",2,2,0,0)+0.000350*ewgc("Cqu8R",2,2,1,1)
2876 -0.008505*ewgc("CuBR",2,2)+0.003205*ewgc("Cud8R",2,2,0,0)
2877 +0.000431*ewgc("Cud8R",2,2,1,1)
2878 -0.050838*ewgc("CuGR",2,2)+0.023185*ewgc("CuuR",0,2,2,0)
2879 +0.000655*ewgc("CuuR",1,2,2,1)-0.004930*ewgc("CuWR",2,2)
2880 +-0.020982*ewgc("CHD")+-0.012029*ewgc("CHl3R",0,0)+-0.012029*ewgc("CHl3R",1,1)+-0.045039*ewgc("CHWB")
2881 +0.012146*ewgc("CllR",0,1,1,0)+-0.000043*ewgc("Cqd1R",2,2,0,0)+0.000040*ewgc("Cqd1R",2,2,1,1)
2882 +0.001365*ewgc("Cqq1R",0,0,2,2)+-0.000101*ewgc("Cqq1R",1,1,2,2)+0.005271*ewgc("Cqq3R",0,0,2,2)
2883 +0.000313*ewgc("Cqq3R",1,1,2,2)+0.000917*ewgc("Cqu1R",0,0,2,2)+0.000089*ewgc("Cqu1R",1,1,2,2)
2884 +0.000787*ewgc("Cqu1R",2,2,0,0)+0.000041*ewgc("Cqu1R",2,2,1,1)+-0.000264*ewgc("Cud1R",2,2,0,0)
2885 +0.000019*ewgc("Cud1R",2,2,1,1)+0.001650*ewgc("CuuR",0,0,2,2)+0.000103*ewgc("CuuR",1,1,2,2));
2886
2887
2888 double sigma_tta_bin_100_130_NP = (-0.008092*ewgc("CG")+0.000916*ewgc("Cqd8R",2,2,0,0)+0.000088*ewgc("Cqd8R",2,2,1,1)
2889 +0.012744*ewgc("Cqq1R",0,2,2,0)+0.000316*ewgc("Cqq1R",1,2,2,1)
2890 +0.026925*ewgc("Cqq3R",0,2,2,0)+0.001792*ewgc("Cqq3R",1,2,2,1)+0.007060*ewgc("Cqu8R",0,0,2,2)
2891 +0.000247*ewgc("Cqu8R",1,1,2,2)+0.006144*ewgc("Cqu8R",2,2,0,0)+0.000162*ewgc("Cqu8R",2,2,1,1)
2892 +-0.005920*ewgc("CuBR",2,2)+0.001849*ewgc("Cud8R",2,2,0,0)+0.000200*ewgc("Cud8R",2,2,1,1)
2893 +-0.027969*ewgc("CuGR",2,2)+0.012428*ewgc("CuuR",0,2,2,0)
2894 +0.000308*ewgc("CuuR",1,2,2,1)+-0.003364*ewgc("CuWR",2,2)
2895 +-0.011717*ewgc("CHD")+-0.006725*ewgc("CHl3R",0,0)+-0.006725*ewgc("CHl3R",1,1)+-0.025133*ewgc("CHWB")
2896 +0.006730*ewgc("CllR",0,1,1,0)+-0.000029*ewgc("Cqd1R",2,2,0,0)+-0.000010*ewgc("Cqd1R",2,2,1,1)
2897 +0.000698*ewgc("Cqq1R",0,0,2,2)+-0.000082*ewgc("Cqq1R",1,1,2,2)+0.002982*ewgc("Cqq3R",0,0,2,2)
2898 +0.000170*ewgc("Cqq3R",1,1,2,2)+0.000509*ewgc("Cqu1R",0,0,2,2)+0.000010*ewgc("Cqu1R",1,1,2,2)
2899 +0.000380*ewgc("Cqu1R",2,2,0,0)+0.000006*ewgc("Cqu1R",2,2,1,1)+-0.000174*ewgc("Cud1R",2,2,0,0)
2900 +-0.000022*ewgc("Cud1R",2,2,1,1)+0.000917*ewgc("CuuR",0,0,2,2)+0.000026*ewgc("CuuR",1,1,2,2));
2901
2902 double total = SM_sigma_tta_bin_70_130 + (sigma_tta_bin_70_100_NP + sigma_tta_bin_100_130_NP)*(SM_sigma_tta_bin_70_130/sigma_tta_bin_70_130_madgraph);
2903
2904 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
2905
2906 return total;
2907
2908 }
2909 } else if(b_min == 130 && b_max == 200){
2910
2911
2912 double SM_sigma_tta_bin_130_200 = SM.getOptionalParameter("SM_sigma_tta_bin_130_200_CMS_dilepton");
2913 double sigma_tta_bin_130_165_madgraph = 0.075543; //pb
2914 double sigma_tta_bin_165_200_madgraph = 0.044687; //pb
2915 double sigma_tta_bin_130_200_madgraph = sigma_tta_bin_130_165_madgraph+sigma_tta_bin_165_200_madgraph;
2916
2917 if(flag_Quadratic){
2918
2919 return 0.;
2920
2921 }
2922 else{
2923
2924 double sigma_tta_bin_130_165_NP = (-0.006183*ewgc("CG")+0.000634*ewgc("Cqd8R",2,2,0,0)+0.000043*ewgc("Cqd8R",2,2,1,1)
2925 +0.008637*ewgc("Cqq1R",0,2,2,0)+0.000173*ewgc("Cqq1R",1,2,2,1)
2926 +0.018848*ewgc("Cqq3R",0,2,2,0)+0.001189*ewgc("Cqq3R",1,2,2,1)+0.005006*ewgc("Cqu8R",0,0,2,2)
2927 +0.000145*ewgc("Cqu8R",1,1,2,2)+0.004346*ewgc("Cqu8R",2,2,0,0)+0.000081*ewgc("Cqu8R",2,2,1,1)
2928 +-0.004429*ewgc("CuBR",2,2)+0.001330*ewgc("Cud8R",2,2,0,0)+0.000119*ewgc("Cud8R",2,2,1,1)
2929 +-0.018971*ewgc("CuGR",2,2)+0.008402*ewgc("CuuR",0,2,2,0)
2930 +0.000171*ewgc("CuuR",1,2,2,1)+-0.002499*ewgc("CuWR",2,2)
2931 +-0.008004*ewgc("CHD")+-0.004609*ewgc("CHl3R",0,0)+-0.004609*ewgc("CHl3R",1,1)+-0.017140*ewgc("CHWB")
2932 +0.004568*ewgc("CllR",0,1,1,0)+-0.000040*ewgc("Cqd1R",2,2,0,0)+-0.000023*ewgc("Cqd1R",2,2,1,1)
2933 +0.000407*ewgc("Cqq1R",0,0,2,2)+-0.000077*ewgc("Cqq1R",1,1,2,2)+0.002110*ewgc("Cqq3R",0,0,2,2)
2934 +0.000087*ewgc("Cqq3R",1,1,2,2)+0.000341*ewgc("Cqu1R",0,0,2,2)+-0.000013*ewgc("Cqu1R",1,1,2,2)
2935 +0.000255*ewgc("Cqu1R",2,2,0,0)+-0.000015*ewgc("Cqu1R",2,2,1,1)+-0.000134*ewgc("Cud1R",2,2,0,0)
2936 +-0.000034*ewgc("Cud1R",2,2,1,1)+0.000617*ewgc("CuuR",0,0,2,2)+-0.000013*ewgc("CuuR",1,1,2,2));
2937
2938 double sigma_tta_bin_165_200_NP = (-0.004285*ewgc("CG")+0.000401*ewgc("Cqd8R",2,2,0,0)+0.000029*ewgc("Cqd8R",2,2,1,1)
2939 +0.005252*ewgc("Cqq1R",0,2,2,0)+0.000104*ewgc("Cqq1R",1,2,2,1)
2940 +0.011756*ewgc("Cqq3R",0,2,2,0)+0.000720*ewgc("Cqq3R",1,2,2,1)+0.003157*ewgc("Cqu8R",0,0,2,2)
2941 +0.000089*ewgc("Cqu8R",1,1,2,2)+0.002752*ewgc("Cqu8R",2,2,0,0)+0.000051*ewgc("Cqu8R",2,2,1,1)
2942 +-0.002716*ewgc("CuBR",2,2)+0.000855*ewgc("Cud8R",2,2,0,0)+0.000072*ewgc("Cud8R",2,2,1,1)
2943 +-0.011315*ewgc("CuGR",2,2)+0.005106*ewgc("CuuR",0,2,2,0)
2944 +0.000098*ewgc("CuuR",1,2,2,1)+-0.001515*ewgc("CuWR",2,2)
2945 +-0.004733*ewgc("CHD")+-0.002719*ewgc("CHl3R",0,0)+-0.002719*ewgc("CHl3R",1,1)+-0.010136*ewgc("CHWB")
2946 +0.002705*ewgc("CllR",0,1,1,0)+-0.000020*ewgc("Cqd1R",2,2,0,0)+-0.000010*ewgc("Cqd1R",2,2,1,1)
2947 +0.000221*ewgc("Cqq1R",0,0,2,2)+-0.000045*ewgc("Cqq1R",1,1,2,2)+0.001322*ewgc("Cqq3R",0,0,2,2)
2948 +0.000061*ewgc("Cqq3R",1,1,2,2)+0.000204*ewgc("Cqu1R",0,0,2,2)+-0.000003*ewgc("Cqu1R",1,1,2,2)
2949 +0.000161*ewgc("Cqu1R",2,2,0,0)+-0.000007*ewgc("Cqu1R",2,2,1,1)+-0.000081*ewgc("Cud1R",2,2,0,0)
2950 +-0.000014*ewgc("Cud1R",2,2,1,1)+0.000411*ewgc("CuuR",0,0,2,2)+-0.000002*ewgc("CuuR",1,1,2,2));
2951
2952 double total = SM_sigma_tta_bin_130_200 + (sigma_tta_bin_130_165_NP+sigma_tta_bin_165_200_NP)*(SM_sigma_tta_bin_130_200/sigma_tta_bin_130_200_madgraph);
2953
2954 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
2955
2956 return total;
2957
2958 }
2959 } else if(b_min == 200 && b_max == 300){
2960
2961
2962 double SM_sigma_tta_bin_200_300 = SM.getOptionalParameter("SM_sigma_tta_bin_200_300_CMS_dilepton");
2963 double sigma_tta_bin_200_250_madgraph = 0.035911; //pb
2964 double sigma_tta_bin_250_300_madgraph = 0.018900; //pb
2965 double sigma_tta_bin_200_300_madgraph = sigma_tta_bin_200_250_madgraph+sigma_tta_bin_250_300_madgraph; //pb
2966
2967
2968 if(flag_Quadratic){
2969
2970 return 0.;
2971
2972 }
2973 else{
2974
2975 double sigma_tta_bin_200_250_NP = (-0.004232*ewgc("CG")+0.000362*ewgc("Cqd8R",2,2,0,0)+0.000028*ewgc("Cqd8R",2,2,1,1)
2976 +0.004454*ewgc("Cqq1R",0,2,2,0)+0.000087*ewgc("Cqq1R",1,2,2,1)
2977 +0.010395*ewgc("Cqq3R",0,2,2,0)+0.000590*ewgc("Cqq3R",1,2,2,1)+0.002802*ewgc("Cqu8R",0,0,2,2)
2978 +0.000078*ewgc("Cqu8R",1,1,2,2)+0.002436*ewgc("Cqu8R",2,2,0,0)+0.000047*ewgc("Cqu8R",2,2,1,1)
2979 -0.002135*ewgc("CuBR",2,2)+0.000773*ewgc("Cud8R",2,2,0,0)+0.000064*ewgc("Cud8R",2,2,1,1)
2980 -0.009208*ewgc("CuGR",2,2)+0.004354*ewgc("CuuR",0,2,2,0)
2981 +0.000086*ewgc("CuuR",1,2,2,1)+-0.001220*ewgc("CuWR",2,2)
2982 +-0.003802*ewgc("CHD")+-0.002185*ewgc("CHl3R",0,0)+-0.002185*ewgc("CHl3R",1,1)+-0.008144*ewgc("CHWB")
2983 +0.002181*ewgc("CllR",0,1,1,0)+-0.000016*ewgc("Cqd1R",2,2,0,0)+-0.000004*ewgc("Cqd1R",2,2,1,1)
2984 +0.000176*ewgc("Cqq1R",0,0,2,2)+-0.000031*ewgc("Cqq1R",1,1,2,2)+0.001156*ewgc("Cqq3R",0,0,2,2)
2985 +0.000055*ewgc("Cqq3R",1,1,2,2)+0.000179*ewgc("Cqu1R",0,0,2,2)+0.000003*ewgc("Cqu1R",1,1,2,2)
2986 +0.000131*ewgc("Cqu1R",2,2,0,0)+0.000000*ewgc("Cqu1R",2,2,1,1)+-0.000073*ewgc("Cud1R",2,2,0,0)
2987 +-0.000011*ewgc("Cud1R",2,2,1,1)+0.000320*ewgc("CuuR",0,0,2,2)+0.000005*ewgc("CuuR",1,1,2,2));
2988
2989 double sigma_tta_bin_250_300_NP = (-0.002904*ewgc("CG")+0.000216*ewgc("Cqd8R",2,2,0,0)+0.000014*ewgc("Cqd8R",2,2,1,1)
2990 +0.002603*ewgc("Cqq1R",0,2,2,0)+0.000045*ewgc("Cqq1R",1,2,2,1)
2991 +0.006131*ewgc("Cqq3R",0,2,2,0)+0.000331*ewgc("Cqq3R",1,2,2,1)+0.001689*ewgc("Cqu8R",0,0,2,2)
2992 +0.000042*ewgc("Cqu8R",1,1,2,2)+0.001470*ewgc("Cqu8R",2,2,0,0)+0.000026*ewgc("Cqu8R",2,2,1,1)
2993 -0.001064*ewgc("CuBR",2,2)+0.000469*ewgc("Cud8R",2,2,0,0)+0.000035*ewgc("Cud8R",2,2,1,1)
2994 +-0.004994*ewgc("CuGR",2,2)+0.002538*ewgc("CuuR",0,2,2,0)
2995 +0.000044*ewgc("CuuR",1,2,2,1)+-0.000597*ewgc("CuWR",2,2)
2996 +-0.001999*ewgc("CHD")+-0.001149*ewgc("CHl3R",0,0)+-0.001149*ewgc("CHl3R",1,1)+-0.004281*ewgc("CHWB")
2997 +0.001146*ewgc("CllR",0,1,1,0)+-0.000008*ewgc("Cqd1R",2,2,0,0)+-0.000003*ewgc("Cqd1R",2,2,1,1)
2998 +0.000087*ewgc("Cqq1R",0,0,2,2)+-0.000021*ewgc("Cqq1R",1,1,2,2)+0.000688*ewgc("Cqq3R",0,0,2,2)
2999 +0.000029*ewgc("Cqq3R",1,1,2,2)+0.000105*ewgc("Cqu1R",0,0,2,2)+-0.000003*ewgc("Cqu1R",1,1,2,2)
3000 +0.000073*ewgc("Cqu1R",2,2,0,0)+-0.000001*ewgc("Cqu1R",2,2,1,1)+-0.000043*ewgc("Cud1R",2,2,0,0)
3001 +-0.000008*ewgc("Cud1R",2,2,1,1)+0.000197*ewgc("CuuR",0,0,2,2)+0.000001*ewgc("CuuR",1,1,2,2));
3002
3003 double total = SM_sigma_tta_bin_200_300 + (sigma_tta_bin_200_250_NP + sigma_tta_bin_250_300_NP)*(SM_sigma_tta_bin_200_300/sigma_tta_bin_200_300_madgraph);
3004
3005 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
3006
3007 return total;
3008
3009 }
3010 } else{
3011 throw std::runtime_error("\nERROR: Please specify a correct binning range for sigma_tta_diff_LO_CMS_dilepton.\n");
3012 }
3013
3014}
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_tta_diff_LO_CMS_dilepton_220107301::ewgc
double ewgc(const std::string name) const
Definition: TopQuarkObservables.h:744
sigma_tta_diff_LO_CMS_dilepton_220107301::b_min
double b_min
Definition: TopQuarkObservables.h:761
sigma_tta_diff_LO_CMS_dilepton_220107301::b_max
double b_max
Definition: TopQuarkObservables.h:762

◆ ewgc() [1/3]

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

Definition at line 744 of file TopQuarkObservables.h.

745 {
746 return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
747 }
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_tta_diff_LO_CMS_dilepton_220107301::ewgc ( const std::string  name,
int  i,
int  j 
) const
inline

Definition at line 749 of file TopQuarkObservables.h.

750 {
751 return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
752 }

◆ ewgc() [3/3]

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

Definition at line 754 of file TopQuarkObservables.h.

755 {
756 return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
757 }

Member Data Documentation

◆ b_max

double sigma_tta_diff_LO_CMS_dilepton_220107301::b_max = 0.
private

Definition at line 762 of file TopQuarkObservables.h.

◆ b_min

double sigma_tta_diff_LO_CMS_dilepton_220107301::b_min = 0.
private

Definition at line 761 of file TopQuarkObservables.h.

◆ mytopobs

const TopQuarkObservables sigma_tta_diff_LO_CMS_dilepton_220107301::mytopobs
private

Definition at line 765 of file TopQuarkObservables.h.

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