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

Differential cross section of pp->ttbar w.r.t. the invariant mass of the ttbar pair with the binning of 2108.02803. More...

#include <TopQuarkObservables.h>

+ Inheritance diagram for sigma_tt_diff_mtt_CMS_LO:

Detailed Description

Differential cross section of pp->ttbar w.r.t. the invariant mass of the ttbar pair with the binning of 2108.02803.

Definition at line 172 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_tt_diff_mtt_CMS_LO (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_tt_diff_mtt_CMS_LO()

sigma_tt_diff_mtt_CMS_LO::sigma_tt_diff_mtt_CMS_LO ( const StandardModel &  SM_i)

FL constructor.

Definition at line 520 of file TopQuarkObservables.cpp.

521: ThObservable(SM_i), mytopobs(static_cast<const NPSMEFTd6General&> (SM))
522{
523
524 setParametersForObservable(make_vector<std::string>() << "SM_sigma_tt_CMS_bin_mtt_250_400" << "SM_sigma_tt_CMS_bin_mtt_400_480" <<
525 "SM_sigma_tt_CMS_bin_mtt_480_560" << "SM_sigma_tt_CMS_bin_mtt_560_640" << "SM_sigma_tt_CMS_bin_mtt_640_720" <<
526 "SM_sigma_tt_CMS_bin_mtt_720_800" << "SM_sigma_tt_CMS_bin_mtt_800_900" << "SM_sigma_tt_CMS_bin_mtt_900_1000" <<
527 "SM_sigma_tt_CMS_bin_mtt_1000_1150" << "SM_sigma_tt_CMS_bin_mtt_1150_1300" << "SM_sigma_tt_CMS_bin_mtt_1300_1500" <<
528 "SM_sigma_tt_CMS_bin_mtt_1500_1700" << "SM_sigma_tt_CMS_bin_mtt_1700_2000" << "SM_sigma_tt_CMS_bin_mtt_2000_2300" <<
529 "SM_sigma_tt_CMS_bin_mtt_2300_3500");
530 //"SM_sigma_tt_bin_2300_2600" << "SM_sigma_tt_bin_2600_3000" << "SM_sigma_tt_bin_3000_3500" << "SM_sigma_tt_bin_3500_4000");
531
532}
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_tt_diff_mtt_CMS_LO::mytopobs
const TopQuarkObservables mytopobs
Definition: TopQuarkObservables.h:210

Member Function Documentation

◆ computeThValue()

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

Implements ThObservable.

Definition at line 534 of file TopQuarkObservables.cpp.

535{
536
537 b_min = getBinMin();
538 b_max = getBinMax();
539
540 bool flag_Quadratic= false; //Needs to be properly defined
541
542
543 //double C_tG = ewgc("CuGR",2,2);
544
545
546
547
548 if(b_min == 250 && b_max == 400){
549
550 double SM_sigma_tt_bin_250_400 = SM.getOptionalParameter("SM_sigma_tt_CMS_bin_mtt_250_400");
551 //double sigma_tt_bin_250_400_madgraph = 105600.0;//fb maybe over the bin width? Check!
552 double sigma_tt_bin_250_400_madgraph = 95.612100;//pb
553
554 if(flag_Quadratic){
555
556 return 0.;
557
558 }
559 else{
560
561
562
563 double total = SM_sigma_tt_bin_250_400 + (-2.422920*ewgc("CG")+0.401694*ewgc("Cqd8R",2,2,0,0)+0.123768*ewgc("Cqd8R",2,2,1,1)
564 +2.248510*ewgc("Cqq1R",0,2,2,0)+0.184272*ewgc("Cqq1R",1,2,2,1)+4.951130*ewgc("Cqq3R",0,2,2,0)
565 +0.780684*ewgc("Cqq3R",1,2,2,1)+0.926298*ewgc("Cqu8R",0,0,2,2)+0.152874*ewgc("Cqu8R",1,1,2,2)
566 +0.575544*ewgc("Cqu8R",2,2,0,0)+0.079620*ewgc("Cqu8R",2,2,1,1)+0.404088*ewgc("Cud8R",2,2,0,0)
567 +0.127740*ewgc("Cud8R",2,2,1,1)-31.370400*ewgc("CuGR",2,2)+2.216290*ewgc("CuuR",0,2,2,0)
568 +0.180834*ewgc("CuuR",1,2,2,1)+
569 0.007578*ewgc("Cqd1R",2,2,0,0)+0.042222*ewgc("Cqd1R",2,2,1,1)+0.159660*ewgc("Cqq1R",0,0,2,2)
570 +0.024942*ewgc("Cqq1R",1,1,2,2)+0.397038*ewgc("Cqq3R",0,0,2,2)+0.088098*ewgc("Cqq3R",1,1,2,2)
571 +0.100374*ewgc("Cqu1R",0,0,2,2)+0.039066*ewgc("Cqu1R",1,1,2,2)+0.110028*ewgc("Cqu1R",2,2,0,0)
572 +0.051414*ewgc("Cqu1R",2,2,1,1)+0.023658*ewgc("Cud1R",2,2,0,0)+0.044238*ewgc("Cud1R",2,2,1,1)
573 +0.146250*ewgc("CuuR",0,0,2,2)+0.059376*ewgc("CuuR",1,1,2,2)
574 )*(SM_sigma_tt_bin_250_400/sigma_tt_bin_250_400_madgraph);
575
576 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
577
578 return total;
579
580
581 }
582
583 } else if(b_min == 400 && b_max == 480){
584
585 double SM_sigma_tt_bin_400_480 = SM.getOptionalParameter("SM_sigma_tt_CMS_bin_mtt_400_480");
586 double sigma_tt_bin_400_480_madgraph = 147.218000;//pb
587
588 if(flag_Quadratic){
589
590 return 0.;
591
592 }
593 else{
594
595 double total = SM_sigma_tt_bin_400_480 + (-8.063800*ewgc("CG")+0.537774*ewgc("Cqd8R",2,2,0,0)+0.102036*ewgc("Cqd8R",2,2,1,1)
596 +3.659260*ewgc("Cqq1R",0,2,2,0)+0.168744*ewgc("Cqq1R",1,2,2,1)+7.940650*ewgc("Cqq3R",0,2,2,0)
597 +1.036730*ewgc("Cqq3R",1,2,2,1)+1.420910*ewgc("Cqu8R",0,0,2,2)+0.147192*ewgc("Cqu8R",1,1,2,2)
598 +0.859650*ewgc("Cqu8R",2,2,0,0)+0.036930*ewgc("Cqu8R",2,2,1,1)+0.537786*ewgc("Cud8R",2,2,0,0)
599 +0.101004*ewgc("Cud8R",2,2,1,1)-42.879900*ewgc("CuGR",2,2)+3.585550*ewgc("CuuR",0,2,2,0)
600 +0.166566*ewgc("CuuR",1,2,2,1)
601 -0.030360*ewgc("Cqd1R",2,2,0,0)+-0.014292*ewgc("Cqd1R",2,2,1,1)+0.202368*ewgc("Cqq1R",0,0,2,2)
602 +-0.028410*ewgc("Cqq1R",1,1,2,2)+0.682230*ewgc("Cqq3R",0,0,2,2)+0.060234*ewgc("Cqq3R",1,1,2,2)
603 +0.056178*ewgc("Cqu1R",0,0,2,2)+-0.007908*ewgc("Cqu1R",1,1,2,2)+0.073824*ewgc("Cqu1R",2,2,0,0)
604 +-0.001914*ewgc("Cqu1R",2,2,1,1)+-0.048564*ewgc("Cud1R",2,2,0,0)+-0.013968*ewgc("Cud1R",2,2,1,1)
605 +0.267738*ewgc("CuuR",0,0,2,2)+0.001314*ewgc("CuuR",1,1,2,2)
606 )*(SM_sigma_tt_bin_400_480/sigma_tt_bin_400_480_madgraph);
607
608 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
609
610 return total;
611
612 }
613
614 } else if(b_min == 480 && b_max == 560){
615
616 double SM_sigma_tt_bin_480_560 = SM.getOptionalParameter("SM_sigma_tt_CMS_bin_mtt_480_560");
617 double sigma_tt_bin_480_560_madgraph = 93.593900;//pb
618
619 if(flag_Quadratic){
620
621 return 0.;
622
623 }
624 else{
625
626 double total = SM_sigma_tt_bin_480_560+(-7.025760*ewgc("CG")+0.449310*ewgc("Cqd8R",2,2,0,0)+0.079662*ewgc("Cqd8R",2,2,1,1)
627 +2.769450*ewgc("Cqq1R",0,2,2,0)+0.129120*ewgc("Cqq1R",1,2,2,1)+6.058480*ewgc("Cqq3R",0,2,2,0)
628 +0.711270*ewgc("Cqq3R",1,2,2,1)+1.107320*ewgc("Cqu8R",0,0,2,2)+0.113052*ewgc("Cqu8R",1,1,2,2)
629 +0.676488*ewgc("Cqu8R",2,2,0,0)+0.029088*ewgc("Cqu8R",2,2,1,1)+0.448860*ewgc("Cud8R",2,2,0,0)
630 +0.080076*ewgc("Cud8R",2,2,1,1)-26.003200*ewgc("CuGR",2,2)+2.714080*ewgc("CuuR",0,2,2,0)
631 +0.128328*ewgc("CuuR",1,2,2,1)
632 -0.010302*ewgc("Cqd1R",2,2,0,0)+-0.001752*ewgc("Cqd1R",2,2,1,1)+0.150552*ewgc("Cqq1R",0,0,2,2)
633 +-0.018192*ewgc("Cqq1R",1,1,2,2)+0.575952*ewgc("Cqq3R",0,0,2,2)+0.060726*ewgc("Cqq3R",1,1,2,2)
634 +0.086460*ewgc("Cqu1R",0,0,2,2)+0.002484*ewgc("Cqu1R",1,1,2,2)+0.062640*ewgc("Cqu1R",2,2,0,0)
635 +0.001020*ewgc("Cqu1R",2,2,1,1)+-0.025974*ewgc("Cud1R",2,2,0,0)+-0.005118*ewgc("Cud1R",2,2,1,1)
636 +0.175404*ewgc("CuuR",0,0,2,2)+0.010116*ewgc("CuuR",1,1,2,2)
637 )*(SM_sigma_tt_bin_480_560/sigma_tt_bin_480_560_madgraph);
638
639 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
640
641 return total;
642
643 }
644
645 } else if(b_min == 560 && b_max == 640){
646
647 double SM_sigma_tt_bin_560_640 = SM.getOptionalParameter("SM_sigma_tt_CMS_bin_mtt_560_640");
648 double sigma_tt_bin_560_640_madgraph = 55.041700;//pb
649
650 if(flag_Quadratic){
651
652 return 0.;
653
654 }
655 else{
656
657 double total = SM_sigma_tt_bin_560_640+(-4.951430*ewgc("CG")+0.305342*ewgc("Cqd8R",2,2,0,0)+0.048212*ewgc("Cqd8R",2,2,1,1)
658 +2.089870*ewgc("Cqq1R",0,2,2,0)+0.073246*ewgc("Cqq1R",1,2,2,1)+4.469060*ewgc("Cqq3R",0,2,2,0)
659 +0.460417*ewgc("Cqq3R",1,2,2,1)+0.799277*ewgc("Cqu8R",0,0,2,2)+0.062208*ewgc("Cqu8R",1,1,2,2)
660 +0.489781*ewgc("Cqu8R",2,2,0,0)+0.018157*ewgc("Cqu8R",2,2,1,1)+0.304480*ewgc("Cud8R",2,2,0,0)
661 +0.048640*ewgc("Cud8R",2,2,1,1)-15.174600*ewgc("CuGR",2,2)+2.032110*ewgc("CuuR",0,2,2,0)
662 +0.070646*ewgc("CuuR",1,2,2,1)
663 -0.015037*ewgc("Cqd1R",2,2,0,0)+-0.004797*ewgc("Cqd1R",2,2,1,1)+0.156666*ewgc("Cqq1R",0,0,2,2)
664 +-0.019126*ewgc("Cqq1R",1,1,2,2)+0.483078*ewgc("Cqq3R",0,0,2,2)+0.035353*ewgc("Cqq3R",1,1,2,2)
665 +0.040200*ewgc("Cqu1R",0,0,2,2)+-0.004762*ewgc("Cqu1R",1,1,2,2)+0.026975*ewgc("Cqu1R",2,2,0,0)
666 +0.000032*ewgc("Cqu1R",2,2,1,1)+-0.024794*ewgc("Cud1R",2,2,0,0)+-0.003859*ewgc("Cud1R",2,2,1,1)
667 +0.150933*ewgc("CuuR",0,0,2,2)+0.002163*ewgc("CuuR",1,1,2,2)
668 )*(SM_sigma_tt_bin_560_640/sigma_tt_bin_560_640_madgraph);
669
670 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
671
672 return total;
673
674 }
675
676 } else if(b_min == 640 && b_max == 720){
677
678 double SM_sigma_tt_bin_640_720 = SM.getOptionalParameter("SM_sigma_tt_CMS_bin_mtt_640_720");
679 double sigma_tt_bin_640_720_madgraph = 32.530700;//pb
680
681 if(flag_Quadratic){
682
683 return 0.;
684
685 }
686 else{
687
688 double total = SM_sigma_tt_bin_640_720+(-3.219570*ewgc("CG")+0.232042*ewgc("Cqd8R",2,2,0,0)+0.037100*ewgc("Cqd8R",2,2,1,1)
689 +1.550940*ewgc("Cqq1R",0,2,2,0)+0.053066*ewgc("Cqq1R",1,2,2,1)+3.301500*ewgc("Cqq3R",0,2,2,0)
690 +0.307033*ewgc("Cqq3R",1,2,2,1)+0.604172*ewgc("Cqu8R",0,0,2,2)+0.049607*ewgc("Cqu8R",1,1,2,2)
691 +0.375725*ewgc("Cqu8R",2,2,0,0)+0.014703*ewgc("Cqu8R",2,2,1,1)+0.232665*ewgc("Cud8R",2,2,0,0)
692 +0.036439*ewgc("Cud8R",2,2,1,1)-8.998020*ewgc("CuGR",2,2)+1.513220*ewgc("CuuR",0,2,2,0)
693 +0.052060*ewgc("CuuR",1,2,2,1)
694 -0.002628*ewgc("Cqd1R",2,2,0,0)+0.001228*ewgc("Cqd1R",2,2,1,1)+0.109033*ewgc("Cqq1R",0,0,2,2)
695 +-0.011971*ewgc("Cqq1R",1,1,2,2)+0.357790*ewgc("Cqq3R",0,0,2,2)+0.027758*ewgc("Cqq3R",1,1,2,2)
696 +0.031403*ewgc("Cqu1R",0,0,2,2)+0.002604*ewgc("Cqu1R",1,1,2,2)+0.030589*ewgc("Cqu1R",2,2,0,0)
697 +0.004275*ewgc("Cqu1R",2,2,1,1)+-0.015530*ewgc("Cud1R",2,2,0,0)+-0.001006*ewgc("Cud1R",2,2,1,1)
698 +0.127340*ewgc("CuuR",0,0,2,2)+0.005749*ewgc("CuuR",1,1,2,2)
699 )*(SM_sigma_tt_bin_640_720/sigma_tt_bin_640_720_madgraph);
700
701 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
702
703 return total;
704
705
706 }
707
708 } else if(b_min == 720 && b_max == 800){
709
710 double SM_sigma_tt_bin_720_800 = SM.getOptionalParameter("SM_sigma_tt_CMS_bin_mtt_720_800");
711 double sigma_tt_bin_720_800_madgraph = 19.672800;//pb
712
713 //PROBABLY THE OPERATORS PURELY FROM THE THIRD FAMILY NEED TO BE NEGLECTED SINCE THEIR CONTRIBUTION IS REALLY SMALL
714
715 if(flag_Quadratic){
716
717 return 0.;
718
719 }
720 else{
721
722 double total = SM_sigma_tt_bin_720_800+(-2.083310*ewgc("CG")+0.175894*ewgc("Cqd8R",2,2,0,0)+0.025732*ewgc("Cqd8R",2,2,1,1)
723 +1.182780*ewgc("Cqq1R",0,2,2,0)+0.038678*ewgc("Cqq1R",1,2,2,1)+2.528020*ewgc("Cqq3R",0,2,2,0)
724 +0.209536*ewgc("Cqq3R",1,2,2,1)+0.457521*ewgc("Cqu8R",0,0,2,2)+0.034597*ewgc("Cqu8R",1,1,2,2)
725 +0.282292*ewgc("Cqu8R",2,2,0,0)+0.011773*ewgc("Cqu8R",2,2,1,1)+0.176461*ewgc("Cud8R",2,2,0,0)
726 +0.026183*ewgc("Cud8R",2,2,1,1)-5.490750*ewgc("CuGR",2,2)+1.152010*ewgc("CuuR",0,2,2,0)
727 +0.037207*ewgc("CuuR",1,2,2,1)
728 -0.005403*ewgc("Cqd1R",2,2,0,0)+0.001667*ewgc("Cqd1R",2,2,1,1)+0.074296*ewgc("Cqq1R",0,0,2,2)
729 +-0.006382*ewgc("Cqq1R",1,1,2,2)+0.282862*ewgc("Cqq3R",0,0,2,2)+0.020603*ewgc("Cqq3R",1,1,2,2)
730 +0.024782*ewgc("Cqu1R",0,0,2,2)+0.003390*ewgc("Cqu1R",1,1,2,2)+0.019139*ewgc("Cqu1R",2,2,0,0)
731 +0.003552*ewgc("Cqu1R",2,2,1,1)+-0.011411*ewgc("Cud1R",2,2,0,0)+0.000128*ewgc("Cud1R",2,2,1,1)
732 +0.087248*ewgc("CuuR",0,0,2,2)+0.006031*ewgc("CuuR",1,1,2,2)
733 )*(SM_sigma_tt_bin_720_800/sigma_tt_bin_720_800_madgraph);
734
735 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
736
737 return total;
738
739 }
740
741 } else if(b_min == 800 && b_max == 900){
742
743 double SM_sigma_tt_bin_800_900 = SM.getOptionalParameter("SM_sigma_tt_CMS_bin_mtt_800_900");
744 double sigma_tt_bin_800_900_madgraph = 14.503700;//pb
745
746 if(flag_Quadratic){
747
748 return 0.;
749
750 }
751 else{
752
753 double total = SM_sigma_tt_bin_800_900+(-1.630550*ewgc("CG")+0.163468*ewgc("Cqd8R",2,2,0,0)+0.019984*ewgc("Cqd8R",2,2,1,1)
754 +1.108320*ewgc("Cqq1R",0,2,2,0)+0.030141*ewgc("Cqq1R",1,2,2,1)+2.335190*ewgc("Cqq3R",0,2,2,0)
755 +0.177932*ewgc("Cqq3R",1,2,2,1)+0.423182*ewgc("Cqu8R",0,0,2,2)+0.026737*ewgc("Cqu8R",1,1,2,2)
756 +0.263074*ewgc("Cqu8R",2,2,0,0)+0.007495*ewgc("Cqu8R",2,2,1,1)+0.162380*ewgc("Cud8R",2,2,0,0)
757 +0.019683*ewgc("Cud8R",2,2,1,1)-4.096130*ewgc("CuGR",2,2)+1.080140*ewgc("CuuR",0,2,2,0)
758 +0.029855*ewgc("CuuR",1,2,2,1)
759 -0.003016*ewgc("Cqd1R",2,2,0,0)+-0.000014*ewgc("Cqd1R",2,2,1,1)+0.082508*ewgc("Cqq1R",0,0,2,2)
760 +-0.006467*ewgc("Cqq1R",1,1,2,2)+0.284167*ewgc("Cqq3R",0,0,2,2)+0.016937*ewgc("Cqq3R",1,1,2,2)
761 +0.024673*ewgc("Cqu1R",0,0,2,2)+0.001100*ewgc("Cqu1R",1,1,2,2)+0.017210*ewgc("Cqu1R",2,2,0,0)
762 +0.000256*ewgc("Cqu1R",2,2,1,1)+-0.012320*ewgc("Cud1R",2,2,0,0)+-0.001016*ewgc("Cud1R",2,2,1,1)
763 +0.092839*ewgc("CuuR",0,0,2,2)+0.002999*ewgc("CuuR",1,1,2,2)
764 )*(SM_sigma_tt_bin_800_900/sigma_tt_bin_800_900_madgraph);
765
766 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
767
768 return total;
769
770 }
771
772 } else if(b_min == 900 && b_max == 1000){
773
774 double SM_sigma_tt_bin_900_1000 = SM.getOptionalParameter("SM_sigma_tt_CMS_bin_mtt_900_1000");
775 double sigma_tt_bin_900_1000_madgraph = 8.319110;
776
777 if(flag_Quadratic){
778
779 return 0.;
780
781 }
782 else{
783
784 double total = SM_sigma_tt_bin_900_1000+(-0.968369*ewgc("CG")+0.119704*ewgc("Cqd8R",2,2,0,0)+0.014122*ewgc("Cqd8R",2,2,1,1)
785 +0.816250*ewgc("Cqq1R",0,2,2,0)+0.021400*ewgc("Cqq1R",1,2,2,1)+1.722240*ewgc("Cqq3R",0,2,2,0)
786 +0.118511*ewgc("Cqq3R",1,2,2,1)+0.313127*ewgc("Cqu8R",0,0,2,2)+0.019055*ewgc("Cqu8R",1,1,2,2)
787 +0.193743*ewgc("Cqu8R",2,2,0,0)+0.005861*ewgc("Cqu8R",2,2,1,1)+0.118769*ewgc("Cud8R",2,2,0,0)
788 +0.014256*ewgc("Cud8R",2,2,1,1)-2.388600*ewgc("CuGR",2,2)+0.792650*ewgc("CuuR",0,2,2,0)
789 +0.020774*ewgc("CuuR",1,2,2,1)
790 -0.003293*ewgc("Cqd1R",2,2,0,0)+0.000536*ewgc("Cqd1R",2,2,1,1)+0.050591*ewgc("Cqq1R",0,0,2,2)
791 +-0.004069*ewgc("Cqq1R",1,1,2,2)+0.198784*ewgc("Cqq3R",0,0,2,2)+0.012664*ewgc("Cqq3R",1,1,2,2)
792 +0.014075*ewgc("Cqu1R",0,0,2,2)+0.000829*ewgc("Cqu1R",1,1,2,2)+0.012454*ewgc("Cqu1R",2,2,0,0)
793 +0.001399*ewgc("Cqu1R",2,2,1,1)+-0.009767*ewgc("Cud1R",2,2,0,0)+-0.000207*ewgc("Cud1R",2,2,1,1)
794 +0.062268*ewgc("CuuR",0,0,2,2)+0.002893*ewgc("CuuR",1,1,2,2)
795 )*(SM_sigma_tt_bin_900_1000/sigma_tt_bin_900_1000_madgraph);
796
797
798 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
799
800 return total;
801
802 }
803
804 } else if(b_min == 1000 && b_max == 1150){
805
806 double SM_sigma_tt_bin_1000_1150 = SM.getOptionalParameter("SM_sigma_tt_CMS_bin_mtt_1000_1150");
807 double sigma_tt_bin_1000_1150_madgraph = 6.639310;
808
809 if(flag_Quadratic){
810
811 return 0.;
812
813 }
814 else{
815
816 double total = SM_sigma_tt_bin_1000_1150+(-0.804304*ewgc("CG")+0.123780*ewgc("Cqd8R",2,2,0,0)+0.012627*ewgc("Cqd8R",2,2,1,1)
817 +0.861600*ewgc("Cqq1R",0,2,2,0)+0.019534*ewgc("Cqq1R",1,2,2,1)+1.803180*ewgc("Cqq3R",0,2,2,0)
818 +0.110947*ewgc("Cqq3R",1,2,2,1)+0.330188*ewgc("Cqu8R",0,0,2,2)+0.016923*ewgc("Cqu8R",1,1,2,2)
819 +0.204389*ewgc("Cqu8R",2,2,0,0)+0.004738*ewgc("Cqu8R",2,2,1,1)+0.124652*ewgc("Cud8R",2,2,0,0)
820 +0.012582*ewgc("Cud8R",2,2,1,1)-1.942700*ewgc("CuGR",2,2)+0.840043*ewgc("CuuR",0,2,2,0)
821 +0.018953*ewgc("CuuR",1,2,2,1)
822 -0.003310*ewgc("Cqd1R",2,2,0,0)+-0.000267*ewgc("Cqd1R",2,2,1,1)+0.067012*ewgc("Cqq1R",0,0,2,2)
823 +-0.004102*ewgc("Cqq1R",1,1,2,2)+0.217415*ewgc("Cqq3R",0,0,2,2)+0.011236*ewgc("Cqq3R",1,1,2,2)
824 +0.013693*ewgc("Cqu1R",0,0,2,2)+0.001001*ewgc("Cqu1R",1,1,2,2)+0.009635*ewgc("Cqu1R",2,2,0,0)
825 +0.000624*ewgc("Cqu1R",2,2,1,1)+-0.011255*ewgc("Cud1R",2,2,0,0)+-0.000678*ewgc("Cud1R",2,2,1,1)
826 +0.068914*ewgc("CuuR",0,0,2,2)+0.001922*ewgc("CuuR",1,1,2,2)
827 )*(SM_sigma_tt_bin_1000_1150/sigma_tt_bin_1000_1150_madgraph);
828
829 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
830
831 return total;
832
833
834 }
835
836 } else if(b_min == 1150 && b_max == 1300){
837
838 double SM_sigma_tt_bin_1150_1300 = SM.getOptionalParameter("SM_sigma_tt_CMS_bin_mtt_1150_1300");
839 double sigma_tt_bin_1150_1300_madgraph = 3.246500;
840
841 if(flag_Quadratic){
842
843 return 0.;
844
845 }
846 else{
847
848 double total = SM_sigma_tt_bin_1150_1300+(-0.401966*ewgc("CG")+0.083879*ewgc("Cqd8R",2,2,0,0)+0.007322*ewgc("Cqd8R",2,2,1,1)
849 +0.586193*ewgc("Cqq1R",0,2,2,0)+0.011479*ewgc("Cqq1R",1,2,2,1)+1.212270*ewgc("Cqq3R",0,2,2,0)
850 +0.066123*ewgc("Cqq3R",1,2,2,1)+0.220114*ewgc("Cqu8R",0,0,2,2)+0.009786*ewgc("Cqu8R",1,1,2,2)
851 +0.138518*ewgc("Cqu8R",2,2,0,0)+0.002758*ewgc("Cqu8R",2,2,1,1)+0.083424*ewgc("Cud8R",2,2,0,0)
852 +0.007362*ewgc("Cud8R",2,2,1,1)-0.974239*ewgc("CuGR",2,2)+0.569482*ewgc("CuuR",0,2,2,0)
853 +0.011161*ewgc("CuuR",1,2,2,1)
854 -0.001702*ewgc("Cqd1R",2,2,0,0)+-0.000120*ewgc("Cqd1R",2,2,1,1)+0.040512*ewgc("Cqq1R",0,0,2,2)
855 +-0.002375*ewgc("Cqq1R",1,1,2,2)+0.151994*ewgc("Cqq3R",0,0,2,2)+0.007076*ewgc("Cqq3R",1,1,2,2)
856 +0.013447*ewgc("Cqu1R",0,0,2,2)+0.000342*ewgc("Cqu1R",1,1,2,2)+0.008418*ewgc("Cqu1R",2,2,0,0)
857 +0.000261*ewgc("Cqu1R",2,2,1,1)+-0.006759*ewgc("Cud1R",2,2,0,0)+-0.000548*ewgc("Cud1R",2,2,1,1)
858 +0.048311*ewgc("CuuR",0,0,2,2)+0.001036*ewgc("CuuR",1,1,2,2)
859 )*(SM_sigma_tt_bin_1150_1300/sigma_tt_bin_1150_1300_madgraph);
860
861 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
862
863 return total;
864
865 }
866
867 } else if(b_min == 1300 && b_max == 1500){
868
869 double SM_sigma_tt_bin_1300_1500 = SM.getOptionalParameter("SM_sigma_tt_CMS_bin_mtt_1300_1500");
870 double sigma_tt_bin_1300_1500_madgraph = 2.043090;
871
872 if(flag_Quadratic){
873
874 return 0.;
875
876 }
877 else{
878
879 double total = SM_sigma_tt_bin_1300_1500+(-0.257644*ewgc("CG")+0.071800*ewgc("Cqd8R",2,2,0,0)+0.005583*ewgc("Cqd8R",2,2,1,1)
880 +0.514089*ewgc("Cqq1R",0,2,2,0)+0.009040*ewgc("Cqq1R",1,2,2,1)+1.057920*ewgc("Cqq3R",0,2,2,0)
881 +0.049771*ewgc("Cqq3R",1,2,2,1)+0.193778*ewgc("Cqu8R",0,0,2,2)+0.007659*ewgc("Cqu8R",1,1,2,2)
882 +0.121016*ewgc("Cqu8R",2,2,0,0)+0.002184*ewgc("Cqu8R",2,2,1,1)+0.071700*ewgc("Cud8R",2,2,0,0)
883 +0.005503*ewgc("Cud8R",2,2,1,1)-0.630898*ewgc("CuGR",2,2)+0.499461*ewgc("CuuR",0,2,2,0)
884 +0.008718*ewgc("CuuR",1,2,2,1)
885 -0.001690*ewgc("Cqd1R",2,2,0,0)+0.000018*ewgc("Cqd1R",2,2,1,1)+0.040412*ewgc("Cqq1R",0,0,2,2)
886 +-0.002151*ewgc("Cqq1R",1,1,2,2)+0.135654*ewgc("Cqq3R",0,0,2,2)+0.004878*ewgc("Cqq3R",1,1,2,2)
887 +0.008707*ewgc("Cqu1R",0,0,2,2)+0.000370*ewgc("Cqu1R",1,1,2,2)+0.005981*ewgc("Cqu1R",2,2,0,0)
888 +0.000084*ewgc("Cqu1R",2,2,1,1)+-0.006194*ewgc("Cud1R",2,2,0,0)+-0.000432*ewgc("Cud1R",2,2,1,1)
889 +0.043381*ewgc("CuuR",0,0,2,2)+0.000671*ewgc("CuuR",1,1,2,2)
890 )*(SM_sigma_tt_bin_1300_1500/sigma_tt_bin_1300_1500_madgraph);
891
892 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
893
894 return total;
895
896 }
897
898 } else if(b_min == 1500 && b_max == 1700){
899
900 double SM_sigma_tt_bin_1500_1700 = SM.getOptionalParameter("SM_sigma_tt_CMS_bin_mtt_1500_1700");
901 double sigma_tt_bin_1500_1700_madgraph = 0.914445;
902
903 if(flag_Quadratic){
904
905 return 0.;
906
907 }
908 else{
909
910 double total = SM_sigma_tt_bin_1500_1700+(-0.117785*ewgc("CG")+0.045019*ewgc("Cqd8R",2,2,0,0)+0.003048*ewgc("Cqd8R",2,2,1,1)
911 +0.324342*ewgc("Cqq1R",0,2,2,0)+0.005084*ewgc("Cqq1R",1,2,2,1)+0.662237*ewgc("Cqq3R",0,2,2,0)
912 +0.028217*ewgc("Cqq3R",1,2,2,1)+0.121878*ewgc("Cqu8R",0,0,2,2)+0.004270*ewgc("Cqu8R",1,1,2,2)
913 +0.077000*ewgc("Cqu8R",2,2,0,0)+0.001206*ewgc("Cqu8R",2,2,1,1)+0.045112*ewgc("Cud8R",2,2,0,0)
914 +0.003088*ewgc("Cud8R",2,2,1,1)-0.292794*ewgc("CuGR",2,2)+0.315186*ewgc("CuuR",0,2,2,0)
915 +0.004977*ewgc("CuuR",1,2,2,1)
916 -0.001201*ewgc("Cqd1R",2,2,0,0)+-0.000002*ewgc("Cqd1R",2,2,1,1)+0.021836*ewgc("Cqq1R",0,0,2,2)
917 +-0.001153*ewgc("Cqq1R",1,1,2,2)+0.081725*ewgc("Cqq3R",0,0,2,2)+0.003047*ewgc("Cqq3R",1,1,2,2)
918 +0.006181*ewgc("Cqu1R",0,0,2,2)+0.000329*ewgc("Cqu1R",1,1,2,2)+0.004044*ewgc("Cqu1R",2,2,0,0)
919 +0.000197*ewgc("Cqu1R",2,2,1,1)+-0.003978*ewgc("Cud1R",2,2,0,0)+-0.000220*ewgc("Cud1R",2,2,1,1)
920 +0.026489*ewgc("CuuR",0,0,2,2)+0.000505*ewgc("CuuR",1,1,2,2)
921 )*(SM_sigma_tt_bin_1500_1700/sigma_tt_bin_1500_1700_madgraph);
922
923 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
924
925 return total;
926
927 }
928
929 } else if(b_min == 1700 && b_max == 2000){
930
931 double SM_sigma_tt_bin_1700_2000 = SM.getOptionalParameter("SM_sigma_tt_CMS_bin_mtt_1700_2000");
932 double sigma_tt_bin_1700_2000_madgraph = 0.562282;
933
934 if(flag_Quadratic){
935
936 return 0.;
937
938 }
939 else{
940
941 double total = SM_sigma_tt_bin_1700_2000+(-0.071302*ewgc("CG")+0.039170*ewgc("Cqd8R",2,2,0,0)+0.002352*ewgc("Cqd8R",2,2,1,1)
942 +0.289186*ewgc("Cqq1R",0,2,2,0)+0.004238*ewgc("Cqq1R",1,2,2,1)+0.584013*ewgc("Cqq3R",0,2,2,0)
943 +0.021609*ewgc("Cqq3R",1,2,2,1)+0.107387*ewgc("Cqu8R",0,0,2,2)+0.003355*ewgc("Cqu8R",1,1,2,2)
944 +0.068348*ewgc("Cqu8R",2,2,0,0)+0.001051*ewgc("Cqu8R",2,2,1,1)+0.039148*ewgc("Cud8R",2,2,0,0)
945 +0.002369*ewgc("Cud8R",2,2,1,1)-0.187937*ewgc("CuGR",2,2)+0.281115*ewgc("CuuR",0,2,2,0)
946 +0.004135*ewgc("CuuR",1,2,2,1)
947 -0.001093*ewgc("Cqd1R",2,2,0,0)+0.000001*ewgc("Cqd1R",2,2,1,1)+0.027035*ewgc("Cqq1R",0,0,2,2)
948 +-0.000869*ewgc("Cqq1R",1,1,2,2)+0.077631*ewgc("Cqq3R",0,0,2,2)+0.002319*ewgc("Cqq3R",1,1,2,2)
949 +0.005833*ewgc("Cqu1R",0,0,2,2)+0.000155*ewgc("Cqu1R",1,1,2,2)+0.003598*ewgc("Cqu1R",2,2,0,0)
950 +0.000098*ewgc("Cqu1R",2,2,1,1)+-0.003566*ewgc("Cud1R",2,2,0,0)+-0.000132*ewgc("Cud1R",2,2,1,1)
951 +0.023764*ewgc("CuuR",0,0,2,2)+0.000360*ewgc("CuuR",1,1,2,2)
952 )*(SM_sigma_tt_bin_1700_2000/sigma_tt_bin_1700_2000_madgraph);
953
954 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
955
956 return total;
957
958
959 }
960
961 } else if(b_min == 2000 && b_max == 2300){
962
963 double SM_sigma_tt_bin_2000_2300 = SM.getOptionalParameter("SM_sigma_tt_CMS_bin_mtt_2000_2300");
964 double sigma_tt_bin_2000_2300_madgraph = 0.204084;
965
966 if(flag_Quadratic){
967
968 return 0.;
969
970 }
971 else{
972
973 double total = SM_sigma_tt_bin_2000_2300+(-0.025926*ewgc("CG")+0.020907*ewgc("Cqd8R",2,2,0,0)+0.001111*ewgc("Cqd8R",2,2,1,1)
974 +0.158780*ewgc("Cqq1R",0,2,2,0)+0.002134*ewgc("Cqq1R",1,2,2,1)+0.315561*ewgc("Cqq3R",0,2,2,0)
975 +0.010496*ewgc("Cqq3R",1,2,2,1)+0.058108*ewgc("Cqu8R",0,0,2,2)+0.001605*ewgc("Cqu8R",1,1,2,2)
976 +0.037313*ewgc("Cqu8R",2,2,0,0)+0.000506*ewgc("Cqu8R",2,2,1,1)+0.020857*ewgc("Cud8R",2,2,0,0)
977 +0.001103*ewgc("Cud8R",2,2,1,1)-0.072063*ewgc("CuGR",2,2)+0.154056*ewgc("CuuR",0,2,2,0)
978 +0.002073*ewgc("CuuR",1,2,2,1)
979 -0.000189*ewgc("Cqd1R",2,2,0,0)+-0.000024*ewgc("Cqd1R",2,2,1,1)+0.011560*ewgc("Cqq1R",0,0,2,2)
980 +-0.000319*ewgc("Cqq1R",1,1,2,2)+0.039980*ewgc("Cqq3R",0,0,2,2)+0.001166*ewgc("Cqq3R",1,1,2,2)
981 +0.002430*ewgc("Cqu1R",0,0,2,2)+0.000093*ewgc("Cqu1R",1,1,2,2)+0.001730*ewgc("Cqu1R",2,2,0,0)
982 +0.000036*ewgc("Cqu1R",2,2,1,1)+-0.001642*ewgc("Cud1R",2,2,0,0)+-0.000091*ewgc("Cud1R",2,2,1,1)
983 +0.013720*ewgc("CuuR",0,0,2,2)+0.000218*ewgc("CuuR",1,1,2,2)
984 )*(SM_sigma_tt_bin_2000_2300/sigma_tt_bin_2000_2300_madgraph);
985
986 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
987
988 return total;
989
990
991 }
992
993 } else if(b_min == 2300 && b_max == 3500){
994
995 double SM_sigma_tt_bin_2300_3500 = SM.getOptionalParameter("SM_sigma_tt_CMS_bin_mtt_2300_3500");
996 double sigma_tt_bin_2300_3500_madgraph = 0.140478;
997
998 if(flag_Quadratic){
999
1000 return 0.;
1001
1002 }
1003 else{
1004
1005 double total = SM_sigma_tt_bin_2300_3500+(-0.016913*ewgc("CG")+0.025555*ewgc("Cqd8R",2,2,0,0)+0.001168*ewgc("Cqd8R",2,2,1,1)
1006 +0.211627*ewgc("Cqq1R",0,2,2,0)+0.002920*ewgc("Cqq1R",1,2,2,1)+0.406226*ewgc("Cqq3R",0,2,2,0)
1007 +0.011576*ewgc("Cqq3R",1,2,2,1)+0.075650*ewgc("Cqu8R",0,0,2,2)+0.001858*ewgc("Cqu8R",1,1,2,2)
1008 +0.049952*ewgc("Cqu8R",2,2,0,0)+0.000687*ewgc("Cqu8R",2,2,1,1)+0.025594*ewgc("Cud8R",2,2,0,0)
1009 +0.001165*ewgc("Cud8R",2,2,1,1)-0.054955*ewgc("CuGR",2,2)+0.205602*ewgc("CuuR",0,2,2,0)
1010 +0.002840*ewgc("CuuR",1,2,2,1)
1011 -0.000606*ewgc("Cqd1R",2,2,0,0)+-0.000021*ewgc("Cqd1R",2,2,1,1)+0.018206*ewgc("Cqq1R",0,0,2,2)
1012 +-0.000313*ewgc("Cqq1R",1,1,2,2)+0.052844*ewgc("Cqq3R",0,0,2,2)+0.001307*ewgc("Cqq3R",1,1,2,2)
1013 +0.003438*ewgc("Cqu1R",0,0,2,2)+0.000081*ewgc("Cqu1R",1,1,2,2)+0.002306*ewgc("Cqu1R",2,2,0,0)
1014 +0.000029*ewgc("Cqu1R",2,2,1,1)+-0.002283*ewgc("Cud1R",2,2,0,0)+-0.000106*ewgc("Cud1R",2,2,1,1)
1015 +0.017622*ewgc("CuuR",0,0,2,2)+0.000252*ewgc("CuuR",1,1,2,2)
1016 )*(SM_sigma_tt_bin_2300_3500/sigma_tt_bin_2300_3500_madgraph);
1017
1018 if (total < 0) return std::numeric_limits<double>::quiet_NaN();
1019
1020 return total;
1021
1022
1023
1024 }
1025
1026 } else {
1027 throw std::runtime_error("\nERROR: Please specify a correct binning range for sigma_tt_diff_LO.\n");
1028 }
1029
1030}
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_tt_diff_mtt_CMS_LO::b_max
double b_max
Definition: TopQuarkObservables.h:207
sigma_tt_diff_mtt_CMS_LO::ewgc
double ewgc(const std::string name) const
Definition: TopQuarkObservables.h:189
sigma_tt_diff_mtt_CMS_LO::b_min
double b_min
Definition: TopQuarkObservables.h:206

◆ ewgc() [1/3]

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

Definition at line 189 of file TopQuarkObservables.h.

190 {
191 return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
192 }
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_tt_diff_mtt_CMS_LO::ewgc ( const std::string  name,
int  i,
int  j 
) const
inline

Definition at line 194 of file TopQuarkObservables.h.

195 {
196 return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
197 }

◆ ewgc() [3/3]

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

Definition at line 199 of file TopQuarkObservables.h.

200 {
201 return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
202 }

Member Data Documentation

◆ b_max

double sigma_tt_diff_mtt_CMS_LO::b_max = 0.
private

Definition at line 207 of file TopQuarkObservables.h.

◆ b_min

double sigma_tt_diff_mtt_CMS_LO::b_min = 0.
private

Definition at line 206 of file TopQuarkObservables.h.

◆ mytopobs

const TopQuarkObservables sigma_tt_diff_mtt_CMS_LO::mytopobs
private

Definition at line 210 of file TopQuarkObservables.h.

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