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

Test Observable. More...

#include <TopQuarkObservables.h>

+ Inheritance diagram for charge_asymmetry_tt_diff_mtt_LO:

Detailed Description

Test Observable.

Definition at line 409 of file TopQuarkObservables.h.

Public Member Functions

 charge_asymmetry_tt_diff_mtt_LO (const StandardModel &SM_i)
 charge_asymmetry_tt_diff_mtt_LO constructor. More...
 
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
 
- 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

◆ charge_asymmetry_tt_diff_mtt_LO()

charge_asymmetry_tt_diff_mtt_LO::charge_asymmetry_tt_diff_mtt_LO ( const StandardModel &  SM_i)

charge_asymmetry_tt_diff_mtt_LO constructor.

Definition at line 1569 of file TopQuarkObservables.cpp.

1570: ThObservable(SM_i), mytopobs(static_cast<const NPSMEFTd6General&> (SM))
1571{
1572
1573 setParametersForObservable(make_vector<std::string>() << "SM_charge_asymmetry_bin_mtt_0_500" << "SM_charge_asymmetry_deno_bin_mtt_0_500"
1574 << "SM_charge_asymmetry_bin_mtt_500_750" << "SM_charge_asymmetry_deno_bin_mtt_500_750" << "SM_charge_asymmetry_bin_mtt_750_1000"
1575 << "SM_charge_asymmetry_deno_bin_mtt_750_1000" << "SM_charge_asymmetry_bin_mtt_1000_1500" << "SM_charge_asymmetry_deno_bin_mtt_1000_1500"
1576 << "SM_charge_asymmetry_bin_mtt_1500_3000" << "SM_charge_asymmetry_deno_bin_mtt_1500_3000");
1577 //"SM_sigma_tt_bin_2300_2600" << "SM_sigma_tt_bin_2600_3000" << "SM_sigma_tt_bin_3000_3500" << "SM_sigma_tt_bin_3500_4000");
1578
1579}
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
charge_asymmetry_tt_diff_mtt_LO::mytopobs
const TopQuarkObservables mytopobs
Definition: TopQuarkObservables.h:446

Member Function Documentation

◆ computeThValue()

double charge_asymmetry_tt_diff_mtt_LO::computeThValue ( )
virtual
Returns
The value of charge_asymmetry_tt_diff_mtt_LO

Implements ThObservable.

Definition at line 1582 of file TopQuarkObservables.cpp.

1583{
1584
1585 b_min = getBinMin();
1586 b_max = getBinMax();
1587
1588 bool flag_Quadratic= false; //Needs to be properly defined
1589
1590
1591
1592
1593
1594
1595
1596 if(b_min == 0 && b_max == 500){
1597
1598 double SM_charge_asymmetry_bin_mtt_0_500 = SM.getOptionalParameter("SM_charge_asymmetry_bin_mtt_0_500");
1599 double SM_charge_asymmetry_deno_bin_mtt_0_500 = SM.getOptionalParameter("SM_charge_asymmetry_deno_bin_mtt_0_500");
1600 double SM_charge_asymmetry_num_bin_mtt_0_500 = SM_charge_asymmetry_bin_mtt_0_500*SM_charge_asymmetry_deno_bin_mtt_0_500;
1601
1602 double SM_sigma_pos_bin_mtt_0_500 =0.5*(SM_charge_asymmetry_num_bin_mtt_0_500+SM_charge_asymmetry_deno_bin_mtt_0_500);
1603 double SM_sigma_neg_bin_mtt_0_500 =0.5*(SM_charge_asymmetry_deno_bin_mtt_0_500-SM_charge_asymmetry_num_bin_mtt_0_500);
1604
1605
1606 if(flag_Quadratic){
1607
1608 return 0.;
1609
1610 }
1611 else{
1612
1613
1614
1615 double sigma_pos_bin_mtt_0_500_madgraph = 135.533000;
1616 double sigma_pos_bin_mtt_0_500_NP = -6.478050*ewgc("CG")+0.364812*ewgc("Cqd8R",2,2,0,0)
1617 +0.065870*ewgc("Cqd8R",2,2,1,1)+4.159410*ewgc("Cqq1R",0,2,2,0)+0.184853*ewgc("Cqq1R",1,2,2,1)
1618 +8.461010*ewgc("Cqq3R",0,2,2,0)+0.893405*ewgc("Cqq3R",1,2,2,1)+1.126020*ewgc("Cqu8R",0,0,2,2)
1619 -0.010419*ewgc("Cqu8R",1,1,2,2)+0.588354*ewgc("Cqu8R",2,2,0,0)+-0.002524*ewgc("Cqu8R",2,2,1,1)
1620 +0.455566*ewgc("Cud8R",2,2,0,0)+0.131078*ewgc("Cud8R",2,2,1,1)-39.064300*ewgc("CuGR",2,2)
1621 +4.072180*ewgc("CuuR",0,2,2,0)+0.097783*ewgc("CuuR",1,2,2,1)
1622 +0.022580*ewgc("Cqd1R",2,2,0,0)+0.004491*ewgc("Cqd1R",2,2,1,1)+0.359419*ewgc("Cqq1R",0,0,2,2)
1623 +-0.028395*ewgc("Cqq1R",1,1,2,2)+0.861834*ewgc("Cqq3R",0,0,2,2)+0.132996*ewgc("Cqq3R",1,1,2,2)
1624 +-0.067518*ewgc("Cqu1R",0,0,2,2)+0.020733*ewgc("Cqu1R",1,1,2,2)+-0.041823*ewgc("Cqu1R",2,2,0,0)
1625 +-0.014534*ewgc("Cqu1R",2,2,1,1)+0.134580*ewgc("Cud1R",2,2,0,0)+-0.004280*ewgc("Cud1R",2,2,1,1)
1626 +0.096208*ewgc("CuuR",0,0,2,2)+0.038326*ewgc("CuuR",1,1,2,2);
1627 double sigma_neg_bin_mtt_0_500_madgraph = 135.300000;
1628 double sigma_neg_bin_mtt_0_500_NP = -6.358540*ewgc("CG")+0.411874*ewgc("Cqd8R",2,2,0,0)
1629 +-0.006404*ewgc("Cqd8R",2,2,1,1)+2.642460*ewgc("Cqq1R",0,2,2,0)+0.063511*ewgc("Cqq1R",1,2,2,1)
1630 +5.970760*ewgc("Cqq3R",0,2,2,0)+0.714859*ewgc("Cqq3R",1,2,2,1)+1.615130*ewgc("Cqu8R",0,0,2,2)
1631 -0.079004*ewgc("Cqu8R",1,1,2,2)+0.944519*ewgc("Cqu8R",2,2,0,0)+0.028645*ewgc("Cqu8R",2,2,1,1)
1632 +0.301135*ewgc("Cud8R",2,2,0,0)+-0.078071*ewgc("Cud8R",2,2,1,1)-39.340100*ewgc("CuGR",2,2)
1633 +2.559810*ewgc("CuuR",0,2,2,0)+-0.060753*ewgc("CuuR",1,2,2,1)
1634 +-0.015693*ewgc("Cqd1R",2,2,0,0)+-0.008030*ewgc("Cqd1R",2,2,1,1)+0.123273*ewgc("Cqq1R",0,0,2,2)
1635 +0.052135*ewgc("Cqq1R",1,1,2,2)+0.387403*ewgc("Cqq3R",0,0,2,2)+-0.093296*ewgc("Cqq3R",1,1,2,2)
1636 +-0.012692*ewgc("Cqu1R",0,0,2,2)+-0.019275*ewgc("Cqu1R",1,1,2,2)+-0.015188*ewgc("Cqu1R",2,2,0,0)
1637 +0.021649*ewgc("Cqu1R",2,2,1,1)+0.001446*ewgc("Cud1R",2,2,0,0)+0.011419*ewgc("Cud1R",2,2,1,1)
1638 +-0.088095*ewgc("CuuR",0,0,2,2)+-0.031337*ewgc("CuuR",1,1,2,2);
1639
1640 double sigma_pos_bin_mtt_0_500_NP_Corrected = SM_sigma_pos_bin_mtt_0_500*sigma_pos_bin_mtt_0_500_NP/sigma_pos_bin_mtt_0_500_madgraph;
1641 double sigma_neg_bin_mtt_0_500_NP_Corrected = SM_sigma_neg_bin_mtt_0_500*sigma_neg_bin_mtt_0_500_NP/sigma_neg_bin_mtt_0_500_madgraph;
1642
1643 double NP_charge_asymmetry_num_bin_mtt_0_500 = sigma_pos_bin_mtt_0_500_NP_Corrected - sigma_neg_bin_mtt_0_500_NP_Corrected;
1644 double NP_charge_asymmetry_deno_bin_mtt_0_500 = sigma_pos_bin_mtt_0_500_NP_Corrected + sigma_neg_bin_mtt_0_500_NP_Corrected;;
1645
1646 //return SM_charge_asymmetry_bin_mtt_0_500*(1+(NP_charge_asymmetry_num_bin_mtt_0_500-NP_charge_asymmetry_deno_bin_mtt_0_500)/SM_charge_asymmetry_deno_bin_mtt_0_500);
1647
1648 return SM_charge_asymmetry_bin_mtt_0_500 + (NP_charge_asymmetry_num_bin_mtt_0_500
1649 - SM_charge_asymmetry_bin_mtt_0_500*NP_charge_asymmetry_deno_bin_mtt_0_500
1650 )/SM_charge_asymmetry_deno_bin_mtt_0_500;
1651
1652 }
1653
1654 } else if(b_min == 500 && b_max == 750){
1655
1656 double SM_charge_asymmetry_bin_mtt_500_750 = SM.getOptionalParameter("SM_charge_asymmetry_bin_mtt_500_750");
1657 double SM_charge_asymmetry_deno_bin_mtt_500_750 = SM.getOptionalParameter("SM_charge_asymmetry_deno_bin_mtt_500_750");
1658 double SM_charge_asymmetry_num_bin_mtt_500_750 = SM_charge_asymmetry_bin_mtt_500_750*SM_charge_asymmetry_deno_bin_mtt_500_750;
1659
1660
1661 double SM_sigma_pos_bin_mtt_500_750 =0.5*(SM_charge_asymmetry_num_bin_mtt_500_750+SM_charge_asymmetry_deno_bin_mtt_500_750);
1662 double SM_sigma_neg_bin_mtt_500_750 =0.5*(SM_charge_asymmetry_deno_bin_mtt_500_750-SM_charge_asymmetry_num_bin_mtt_500_750);
1663
1664
1665 if(flag_Quadratic){
1666
1667 return 0.;
1668
1669 }
1670 else{
1671
1672
1673
1674
1675 double sigma_pos_bin_mtt_500_750_madgraph = 80.589400;
1676 double sigma_pos_bin_mtt_500_750_NP = -6.999250*ewgc("CG")+0.386926*ewgc("Cqd8R",2,2,0,0)
1677 +0.078147*ewgc("Cqd8R",2,2,1,1)+4.115120*ewgc("Cqq1R",0,2,2,0)+0.207149*ewgc("Cqq1R",1,2,2,1)
1678 +8.367150*ewgc("Cqq3R",0,2,2,0)+0.726869*ewgc("Cqq3R",1,2,2,1)+0.814576*ewgc("Cqu8R",0,0,2,2)
1679 +0.182552*ewgc("Cqu8R",1,1,2,2)+0.471441*ewgc("Cqu8R",2,2,0,0)+0.020673*ewgc("Cqu8R",2,2,1,1)
1680 +0.603346*ewgc("Cud8R",2,2,0,0)+0.063159*ewgc("Cud8R",2,2,1,1)-20.819700*ewgc("CuGR",2,2)
1681 +3.980290*ewgc("CuuR",0,2,2,0)+0.169850*ewgc("CuuR",1,2,2,1)
1682 +0.047499*ewgc("Cqd1R",2,2,0,0)+0.051676*ewgc("Cqd1R",2,2,1,1)+0.322849*ewgc("Cqq1R",0,0,2,2)
1683 +0.058428*ewgc("Cqq1R",1,1,2,2)+0.909203*ewgc("Cqq3R",0,0,2,2)+0.028468*ewgc("Cqq3R",1,1,2,2)
1684 +0.083136*ewgc("Cqu1R",0,0,2,2)+0.009166*ewgc("Cqu1R",1,1,2,2)+0.095883*ewgc("Cqu1R",2,2,0,0)
1685 +0.046037*ewgc("Cqu1R",2,2,1,1)+-0.029524*ewgc("Cud1R",2,2,0,0)+0.062062*ewgc("Cud1R",2,2,1,1)
1686 +0.373091*ewgc("CuuR",0,0,2,2)+0.067338*ewgc("CuuR",1,1,2,2);
1687 double sigma_neg_bin_mtt_500_750_madgraph = 80.891300;
1688 double sigma_neg_bin_mtt_500_750_NP = -7.057070*ewgc("CG")+0.573031*ewgc("Cqd8R",2,2,0,0)
1689 -0.004786*ewgc("Cqd8R",2,2,1,1)+2.052850*ewgc("Cqq1R",0,2,2,0)+0.125976*ewgc("Cqq1R",1,2,2,1)
1690 +4.836950*ewgc("Cqq3R",0,2,2,0)+0.711224*ewgc("Cqq3R",1,2,2,1)+1.470530*ewgc("Cqu8R",0,0,2,2)
1691 +0.069589*ewgc("Cqu8R",1,1,2,2)+0.938294*ewgc("Cqu8R",2,2,0,0)+-0.045546*ewgc("Cqu8R",2,2,1,1)
1692 +0.347501*ewgc("Cud8R",2,2,0,0)+0.011230*ewgc("Cud8R",2,2,1,1)+-20.768900*ewgc("CuGR",2,2)
1693 +1.942090*ewgc("CuuR",0,2,2,0)+0.161845*ewgc("CuuR",1,2,2,1)
1694 +-0.030151*ewgc("Cqd1R",2,2,0,0)+-0.051728*ewgc("Cqd1R",2,2,1,1)+0.012765*ewgc("Cqq1R",0,0,2,2)
1695 +-0.097104*ewgc("Cqq1R",1,1,2,2)+0.424513*ewgc("Cqq3R",0,0,2,2)+-0.015512*ewgc("Cqq3R",1,1,2,2)
1696 +0.014810*ewgc("Cqu1R",0,0,2,2)+-0.010594*ewgc("Cqu1R",1,1,2,2)+0.012532*ewgc("Cqu1R",2,2,0,0)
1697 +-0.046820*ewgc("Cqu1R",2,2,1,1)+-0.078305*ewgc("Cud1R",2,2,0,0)+-0.072728*ewgc("Cud1R",2,2,1,1)
1698 +0.132974*ewgc("CuuR",0,0,2,2)+-0.067169*ewgc("CuuR",1,1,2,2);
1699
1700 double sigma_pos_bin_mtt_500_750_NP_Corrected = SM_sigma_pos_bin_mtt_500_750*sigma_pos_bin_mtt_500_750_NP/sigma_pos_bin_mtt_500_750_madgraph;
1701 double sigma_neg_bin_mtt_500_750_NP_Corrected = SM_sigma_neg_bin_mtt_500_750*sigma_neg_bin_mtt_500_750_NP/sigma_neg_bin_mtt_500_750_madgraph;
1702
1703 double NP_charge_asymmetry_num_bin_mtt_500_750 = sigma_pos_bin_mtt_500_750_NP_Corrected - sigma_neg_bin_mtt_500_750_NP_Corrected;
1704 double NP_charge_asymmetry_deno_bin_mtt_500_750 = sigma_pos_bin_mtt_500_750_NP_Corrected + sigma_neg_bin_mtt_500_750_NP_Corrected;;
1705
1706 //return SM_charge_asymmetry_bin_mtt_500_750*(1+(NP_charge_asymmetry_num_bin_mtt_500_750-NP_charge_asymmetry_deno_bin_mtt_500_750)/SM_charge_asymmetry_deno_bin_mtt_500_750);
1707
1708 return SM_charge_asymmetry_bin_mtt_500_750 + (NP_charge_asymmetry_num_bin_mtt_500_750
1709 - SM_charge_asymmetry_bin_mtt_500_750*NP_charge_asymmetry_deno_bin_mtt_500_750
1710 )/SM_charge_asymmetry_deno_bin_mtt_500_750;
1711
1712
1713
1714
1715 }
1716
1717 } else if(b_min == 750 && b_max == 1000){
1718
1719
1720 double SM_charge_asymmetry_deno_bin_mtt_750_1000 = SM.getOptionalParameter("SM_charge_asymmetry_deno_bin_mtt_750_1000");
1721 double SM_charge_asymmetry_bin_mtt_750_1000 = SM.getOptionalParameter("SM_charge_asymmetry_bin_mtt_750_1000");
1722 double SM_charge_asymmetry_num_bin_mtt_750_1000 = SM_charge_asymmetry_bin_mtt_750_1000*SM_charge_asymmetry_deno_bin_mtt_750_1000;
1723
1724 double SM_sigma_pos_bin_mtt_750_1000 =0.5*(SM_charge_asymmetry_num_bin_mtt_750_1000+SM_charge_asymmetry_deno_bin_mtt_750_1000);
1725 double SM_sigma_neg_bin_mtt_750_1000 =0.5*(SM_charge_asymmetry_deno_bin_mtt_750_1000-SM_charge_asymmetry_num_bin_mtt_750_1000);
1726
1727
1728 if(flag_Quadratic){
1729
1730 return 0.;
1731
1732 }
1733 else{
1734
1735
1736
1737
1738 double sigma_pos_bin_mtt_750_1000_madgraph = 16.986900;
1739 double sigma_pos_bin_mtt_750_1000_NP = -1.924360*ewgc("CG")+0.144619*ewgc("Cqd8R",2,2,0,0)
1740 +0.015897*ewgc("Cqd8R",2,2,1,1)+1.863020*ewgc("Cqq1R",0,2,2,0)+0.031168*ewgc("Cqq1R",1,2,2,1)
1741 +3.676260*ewgc("Cqq3R",0,2,2,0)+0.245826*ewgc("Cqq3R",1,2,2,1)+0.316581*ewgc("Cqu8R",0,0,2,2)
1742 +0.022092*ewgc("Cqu8R",1,1,2,2)+0.190067*ewgc("Cqu8R",2,2,0,0)+0.021485*ewgc("Cqu8R",2,2,1,1)
1743 +0.258731*ewgc("Cud8R",2,2,0,0)+0.013146*ewgc("Cud8R",2,2,1,1)-4.307000*ewgc("CuGR",2,2)
1744 +1.792720*ewgc("CuuR",0,2,2,0)+0.019910*ewgc("CuuR",1,2,2,1)
1745 +-0.016893*ewgc("Cqd1R",2,2,0,0)+-0.008497*ewgc("Cqd1R",2,2,1,1)+0.141702*ewgc("Cqq1R",0,0,2,2)
1746 +-0.012729*ewgc("Cqq1R",1,1,2,2)+0.438894*ewgc("Cqq3R",0,0,2,2)+0.019810*ewgc("Cqq3R",1,1,2,2)
1747 +0.021260*ewgc("Cqu1R",0,0,2,2)+-0.004833*ewgc("Cqu1R",1,1,2,2)+0.011461*ewgc("Cqu1R",2,2,0,0)
1748 +-0.003083*ewgc("Cqu1R",2,2,1,1)+-0.036348*ewgc("Cud1R",2,2,0,0)+-0.008414*ewgc("Cud1R",2,2,1,1)
1749 +0.152046*ewgc("CuuR",0,0,2,2)+-0.006108*ewgc("CuuR",1,1,2,2);
1750
1751 double sigma_neg_bin_mtt_750_1000_madgraph = 16.965400;
1752 double sigma_neg_bin_mtt_750_1000_NP = -1.909310*ewgc("CG")+0.253028*ewgc("Cqd8R",2,2,0,0)
1753 +0.026868*ewgc("Cqd8R",2,2,1,1)+0.775801*ewgc("Cqq1R",0,2,2,0)+0.021341*ewgc("Cqq1R",1,2,2,1)
1754 +1.850950*ewgc("Cqq3R",0,2,2,0)+0.208119*ewgc("Cqq3R",1,2,2,1)+0.666539*ewgc("Cqu8R",0,0,2,2)
1755 +0.038662*ewgc("Cqu8R",1,1,2,2)+0.437622*ewgc("Cqu8R",2,2,0,0)+0.027679*ewgc("Cqu8R",2,2,1,1)
1756 +0.143565*ewgc("Cud8R",2,2,0,0)+0.028559*ewgc("Cud8R",2,2,1,1)+-4.304670*ewgc("CuGR",2,2)
1757 +0.748939*ewgc("CuuR",0,2,2,0)+0.029156*ewgc("CuuR",1,2,2,1)
1758 +-0.016259*ewgc("Cqd1R",2,2,0,0)+0.007655*ewgc("Cqd1R",2,2,1,1)+0.028036*ewgc("Cqq1R",0,0,2,2)
1759 +-0.000705*ewgc("Cqq1R",1,1,2,2)+0.186361*ewgc("Cqq3R",0,0,2,2)+0.044037*ewgc("Cqq3R",1,1,2,2)
1760 +0.037482*ewgc("Cqu1R",0,0,2,2)+0.005345*ewgc("Cqu1R",1,1,2,2)+0.040132*ewgc("Cqu1R",2,2,0,0)
1761 +0.002865*ewgc("Cqu1R",2,2,1,1)+-0.012822*ewgc("Cud1R",2,2,0,0)+0.007088*ewgc("Cud1R",2,2,1,1)
1762 +0.055812*ewgc("CuuR",0,0,2,2)+0.007647*ewgc("CuuR",1,1,2,2);
1763
1764 double sigma_pos_bin_mtt_750_1000_NP_Corrected = SM_sigma_pos_bin_mtt_750_1000*sigma_pos_bin_mtt_750_1000_NP/sigma_pos_bin_mtt_750_1000_madgraph;
1765 double sigma_neg_bin_mtt_750_1000_NP_Corrected = SM_sigma_neg_bin_mtt_750_1000*sigma_neg_bin_mtt_750_1000_NP/sigma_neg_bin_mtt_750_1000_madgraph;
1766
1767 double NP_charge_asymmetry_num_bin_mtt_750_1000 = sigma_pos_bin_mtt_750_1000_NP_Corrected - sigma_neg_bin_mtt_750_1000_NP_Corrected;
1768 double NP_charge_asymmetry_deno_bin_mtt_750_1000 = sigma_pos_bin_mtt_750_1000_NP_Corrected + sigma_neg_bin_mtt_750_1000_NP_Corrected;;
1769
1770 //return SM_charge_asymmetry_bin_mtt_750_1000*(1+(NP_charge_asymmetry_num_bin_mtt_750_1000-NP_charge_asymmetry_deno_bin_mtt_750_1000)/SM_charge_asymmetry_deno_bin_mtt_750_1000);
1771
1772 return SM_charge_asymmetry_bin_mtt_750_1000 + (NP_charge_asymmetry_num_bin_mtt_750_1000
1773 - SM_charge_asymmetry_bin_mtt_750_1000*NP_charge_asymmetry_deno_bin_mtt_750_1000
1774 )/SM_charge_asymmetry_deno_bin_mtt_750_1000;
1775
1776
1777
1778 }
1779
1780 } else if(b_min == 1000 && b_max == 1500){
1781
1782
1783 double SM_charge_asymmetry_deno_bin_mtt_1000_1500 = SM.getOptionalParameter("SM_charge_asymmetry_deno_bin_mtt_1000_1500");
1784 double SM_charge_asymmetry_bin_mtt_1000_1500 = SM.getOptionalParameter("SM_charge_asymmetry_bin_mtt_1000_1500");
1785 double SM_charge_asymmetry_num_bin_mtt_1000_1500 = SM_charge_asymmetry_bin_mtt_1000_1500*SM_charge_asymmetry_deno_bin_mtt_1000_1500;
1786
1787
1788 double SM_sigma_pos_bin_mtt_1000_1500 =0.5*(SM_charge_asymmetry_num_bin_mtt_1000_1500+SM_charge_asymmetry_deno_bin_mtt_1000_1500);
1789 double SM_sigma_neg_bin_mtt_1000_1500 =0.5*(SM_charge_asymmetry_deno_bin_mtt_1000_1500-SM_charge_asymmetry_num_bin_mtt_1000_1500);
1790
1791
1792 if(flag_Quadratic){
1793
1794 return 0.;
1795
1796 }
1797 else{
1798
1799
1800 double sigma_pos_bin_mtt_1000_1500_madgraph = 5.965330;
1801 double sigma_pos_bin_mtt_1000_1500_NP = -0.732765*ewgc("CG")+0.095774*ewgc("Cqd8R",2,2,0,0)
1802 +0.009861*ewgc("Cqd8R",2,2,1,1)+1.407820*ewgc("Cqq1R",0,2,2,0)+0.023055*ewgc("Cqq1R",1,2,2,1)
1803 +2.789740*ewgc("Cqq3R",0,2,2,0)+0.124751*ewgc("Cqq3R",1,2,2,1)+0.226047*ewgc("Cqu8R",0,0,2,2)
1804 +0.012025*ewgc("Cqu8R",1,1,2,2)+0.132009*ewgc("Cqu8R",2,2,0,0)+0.002076*ewgc("Cqu8R",2,2,1,1)
1805 +0.176967*ewgc("Cud8R",2,2,0,0)+0.014820*ewgc("Cud8R",2,2,1,1)-1.511020*ewgc("CuGR",2,2)
1806 +1.381670*ewgc("CuuR",0,2,2,0)+0.022990*ewgc("CuuR",1,2,2,1)
1807 +0.000480*ewgc("Cqd1R",2,2,0,0)+-0.000506*ewgc("Cqd1R",2,2,1,1)+0.119062*ewgc("Cqq1R",0,0,2,2)
1808 +-0.007047*ewgc("Cqq1R",1,1,2,2)+0.356119*ewgc("Cqq3R",0,0,2,2)+0.012398*ewgc("Cqq3R",1,1,2,2)
1809 +0.012939*ewgc("Cqu1R",0,0,2,2)+0.000035*ewgc("Cqu1R",1,1,2,2)+0.003252*ewgc("Cqu1R",2,2,0,0)
1810 +-0.002001*ewgc("Cqu1R",2,2,1,1)+-0.013101*ewgc("Cud1R",2,2,0,0)+-0.000576*ewgc("Cud1R",2,2,1,1)
1811 +0.117162*ewgc("CuuR",0,0,2,2)+0.001190*ewgc("CuuR",1,1,2,2);
1812 double sigma_neg_bin_mtt_1000_1500_madgraph = 5.963920;
1813 double sigma_neg_bin_mtt_1000_1500_NP = -0.734104*ewgc("CG")+0.177167*ewgc("Cqd8R",2,2,0,0)
1814 +0.016263*ewgc("Cqd8R",2,2,1,1)+0.546178*ewgc("Cqq1R",0,2,2,0)+0.024062*ewgc("Cqq1R",1,2,2,1)
1815 +1.297260*ewgc("Cqq3R",0,2,2,0)+0.100018*ewgc("Cqq3R",1,2,2,1)+0.517319*ewgc("Cqu8R",0,0,2,2)
1816 +0.016952*ewgc("Cqu8R",1,1,2,2)+0.336176*ewgc("Cqu8R",2,2,0,0)+0.003311*ewgc("Cqu8R",2,2,1,1)
1817 +0.095957*ewgc("Cud8R",2,2,0,0)+0.011525*ewgc("Cud8R",2,2,1,1)-1.511200*ewgc("CuGR",2,2)
1818 +0.534920*ewgc("CuuR",0,2,2,0)+0.024264*ewgc("CuuR",1,2,2,1)
1819 +0.000672*ewgc("Cqd1R",2,2,0,0)+0.000477*ewgc("Cqd1R",2,2,1,1)+0.023176*ewgc("Cqq1R",0,0,2,2)
1820 +-0.006814*ewgc("Cqq1R",1,1,2,2)+0.147804*ewgc("Cqq3R",0,0,2,2)+0.009854*ewgc("Cqq3R",1,1,2,2)
1821 +0.029905*ewgc("Cqu1R",0,0,2,2)+-0.000840*ewgc("Cqu1R",1,1,2,2)+0.013292*ewgc("Cqu1R",2,2,0,0)
1822 +0.001980*ewgc("Cqu1R",2,2,1,1)+-0.004801*ewgc("Cud1R",2,2,0,0)+0.001036*ewgc("Cud1R",2,2,1,1)
1823 +0.045432*ewgc("CuuR",0,0,2,2)+0.000336*ewgc("CuuR",1,1,2,2);
1824
1825
1826 double sigma_pos_bin_mtt_1000_1500_NP_Corrected = SM_sigma_pos_bin_mtt_1000_1500*sigma_pos_bin_mtt_1000_1500_NP/sigma_pos_bin_mtt_1000_1500_madgraph;
1827 double sigma_neg_bin_mtt_1000_1500_NP_Corrected = SM_sigma_neg_bin_mtt_1000_1500*sigma_neg_bin_mtt_1000_1500_NP/sigma_neg_bin_mtt_1000_1500_madgraph;
1828
1829 double NP_charge_asymmetry_num_bin_mtt_1000_1500 = sigma_pos_bin_mtt_1000_1500_NP_Corrected - sigma_neg_bin_mtt_1000_1500_NP_Corrected;
1830 double NP_charge_asymmetry_deno_bin_mtt_1000_1500 = sigma_pos_bin_mtt_1000_1500_NP_Corrected + sigma_neg_bin_mtt_1000_1500_NP_Corrected;;
1831
1832 //return SM_charge_asymmetry_bin_mtt_1000_1500*(1+(NP_charge_asymmetry_num_bin_mtt_1000_1500-NP_charge_asymmetry_deno_bin_mtt_1000_1500)/SM_charge_asymmetry_deno_bin_mtt_1000_1500);
1833
1834 return SM_charge_asymmetry_bin_mtt_1000_1500 + (NP_charge_asymmetry_num_bin_mtt_1000_1500
1835 - SM_charge_asymmetry_bin_mtt_1000_1500*NP_charge_asymmetry_deno_bin_mtt_1000_1500
1836 )/SM_charge_asymmetry_deno_bin_mtt_1000_1500;
1837
1838
1839
1840 }
1841
1842 } else if(b_min == 1500 && b_max == 3000){
1843
1844 double SM_charge_asymmetry_bin_mtt_1500_3000 = SM.getOptionalParameter("SM_charge_asymmetry_bin_mtt_1500_3000");
1845 double SM_charge_asymmetry_deno_bin_mtt_1500_3000 = SM.getOptionalParameter("SM_charge_asymmetry_deno_bin_mtt_1500_3000");
1846 double SM_charge_asymmetry_num_bin_mtt_1500_3000 = SM_charge_asymmetry_bin_mtt_1500_3000*SM_charge_asymmetry_deno_bin_mtt_1500_3000;
1847
1848 double SM_sigma_pos_bin_mtt_1500_3000 =0.5*(SM_charge_asymmetry_num_bin_mtt_1500_3000+SM_charge_asymmetry_deno_bin_mtt_1500_3000);
1849 double SM_sigma_neg_bin_mtt_1500_3000 =0.5*(SM_charge_asymmetry_deno_bin_mtt_1500_3000-SM_charge_asymmetry_num_bin_mtt_1500_3000);
1850
1851
1852
1853 if(flag_Quadratic){
1854
1855 return 0.;
1856
1857 }
1858 else{
1859
1860
1861
1862
1863 double sigma_pos_bin_mtt_1500_3000_madgraph = 0.910259;
1864 double sigma_pos_bin_mtt_1500_3000_NP = -0.116441*ewgc("CG")+0.045409*ewgc("Cqd8R",2,2,0,0)
1865 +0.003096*ewgc("Cqd8R",2,2,1,1)+0.711655*ewgc("Cqq1R",0,2,2,0)+0.007175*ewgc("Cqq1R",1,2,2,1)
1866 +1.350330*ewgc("Cqq3R",0,2,2,0)+0.039100*ewgc("Cqq3R",1,2,2,1)+0.110131*ewgc("Cqu8R",0,0,2,2)
1867 +0.005191*ewgc("Cqu8R",1,1,2,2)+0.064112*ewgc("Cqu8R",2,2,0,0)+0.002451*ewgc("Cqu8R",2,2,1,1)
1868 +0.084721*ewgc("Cud8R",2,2,0,0)+0.004134*ewgc("Cud8R",2,2,1,1)-0.231007*ewgc("CuGR",2,2)
1869 +0.691483*ewgc("CuuR",0,2,2,0)+0.006844*ewgc("CuuR",1,2,2,1)
1870 +-0.001536*ewgc("Cqd1R",2,2,0,0)+0.000002*ewgc("Cqd1R",2,2,1,1)+0.063201*ewgc("Cqq1R",0,0,2,2)
1871 +-0.001288*ewgc("Cqq1R",1,1,2,2)+0.178152*ewgc("Cqq3R",0,0,2,2)+0.004233*ewgc("Cqq3R",1,1,2,2)
1872 +0.005775*ewgc("Cqu1R",0,0,2,2)+0.000328*ewgc("Cqu1R",1,1,2,2)+0.003439*ewgc("Cqu1R",2,2,0,0)
1873 +-0.000036*ewgc("Cqu1R",2,2,1,1)+-0.007397*ewgc("Cud1R",2,2,0,0)+-0.000466*ewgc("Cud1R",2,2,1,1)
1874 +0.059131*ewgc("CuuR",0,0,2,2)+0.000437*ewgc("CuuR",1,1,2,2);
1875 double sigma_neg_bin_mtt_1500_3000_madgraph = 0.909569;
1876 double sigma_neg_bin_mtt_1500_3000_NP = -0.116979*ewgc("CG")+0.084970*ewgc("Cqd8R",2,2,0,0)
1877 +0.004067*ewgc("Cqd8R",2,2,1,1)+0.274254*ewgc("Cqq1R",0,2,2,0)+0.007000*ewgc("Cqq1R",1,2,2,1)
1878 +0.618997*ewgc("Cqq3R",0,2,2,0)+0.032143*ewgc("Cqq3R",1,2,2,1)+0.253632*ewgc("Cqu8R",0,0,2,2)
1879 +0.006071*ewgc("Cqu8R",1,1,2,2)+0.167569*ewgc("Cqu8R",2,2,0,0)+0.002373*ewgc("Cqu8R",2,2,1,1)
1880 +0.045316*ewgc("Cud8R",2,2,0,0)+0.003020*ewgc("Cud8R",2,2,1,1)-0.231248*ewgc("CuGR",2,2)
1881 +0.265616*ewgc("CuuR",0,2,2,0)+0.006959*ewgc("CuuR",1,2,2,1)
1882 +-0.002375*ewgc("Cqd1R",2,2,0,0)+-0.000039*ewgc("Cqd1R",2,2,1,1)+0.015305*ewgc("Cqq1R",0,0,2,2)
1883 +-0.001083*ewgc("Cqq1R",1,1,2,2)+0.075644*ewgc("Cqq3R",0,0,2,2)+0.003511*ewgc("Cqq3R",1,1,2,2)
1884 +0.011743*ewgc("Cqu1R",0,0,2,2)+0.000548*ewgc("Cqu1R",1,1,2,2)+0.008199*ewgc("Cqu1R",2,2,0,0)
1885 +0.000069*ewgc("Cqu1R",2,2,1,1)+-0.003571*ewgc("Cud1R",2,2,0,0)+-0.000568*ewgc("Cud1R",2,2,1,1)
1886 +0.022219*ewgc("CuuR",0,0,2,2)+0.000486*ewgc("CuuR",1,1,2,2);
1887
1888
1889
1890 double sigma_pos_bin_mtt_1500_3000_NP_Corrected = SM_sigma_pos_bin_mtt_1500_3000*sigma_pos_bin_mtt_1500_3000_NP/sigma_pos_bin_mtt_1500_3000_madgraph;
1891 double sigma_neg_bin_mtt_1500_3000_NP_Corrected = SM_sigma_neg_bin_mtt_1500_3000*sigma_neg_bin_mtt_1500_3000_NP/sigma_neg_bin_mtt_1500_3000_madgraph;
1892
1893 double NP_charge_asymmetry_num_bin_mtt_1500_3000 = sigma_pos_bin_mtt_1500_3000_NP_Corrected - sigma_neg_bin_mtt_1500_3000_NP_Corrected;
1894 double NP_charge_asymmetry_deno_bin_mtt_1500_3000 = sigma_pos_bin_mtt_1500_3000_NP_Corrected + sigma_neg_bin_mtt_1500_3000_NP_Corrected;;
1895
1896 //return SM_charge_asymmetry_bin_mtt_1500_3000*(1+(NP_charge_asymmetry_num_bin_mtt_1500_3000-NP_charge_asymmetry_deno_bin_mtt_1500_3000)/SM_charge_asymmetry_deno_bin_mtt_1500_3000);
1897
1898
1899 return SM_charge_asymmetry_bin_mtt_1500_3000 + (NP_charge_asymmetry_num_bin_mtt_1500_3000
1900 - SM_charge_asymmetry_bin_mtt_1500_3000*NP_charge_asymmetry_deno_bin_mtt_1500_3000
1901 )/SM_charge_asymmetry_deno_bin_mtt_1500_3000;
1902
1903 }
1904
1905 }
1906
1907
1908 else {
1909 throw std::runtime_error("\nERROR: Please specify a correct binning range for charge_asymmetry_tt_diff_mtt_LO.\n");
1910 }
1911
1912}
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
charge_asymmetry_tt_diff_mtt_LO::b_min
double b_min
Definition: TopQuarkObservables.h:442
charge_asymmetry_tt_diff_mtt_LO::ewgc
double ewgc(const std::string name) const
Definition: TopQuarkObservables.h:425
charge_asymmetry_tt_diff_mtt_LO::b_max
double b_max
Definition: TopQuarkObservables.h:443

◆ ewgc() [1/3]

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

Definition at line 425 of file TopQuarkObservables.h.

426 {
427 return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name)*1000000;
428 }
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 charge_asymmetry_tt_diff_mtt_LO::ewgc ( const std::string  name,
int  i,
int  j 
) const
inline

Definition at line 430 of file TopQuarkObservables.h.

431 {
432 return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j)*1000000;
433 }

◆ ewgc() [3/3]

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

Definition at line 435 of file TopQuarkObservables.h.

436 {
437 return mytopobs.GetmyNPSMEFTd6General().getSMEFTCoeffEW(name, i, j, k, l)*1000000;
438 }

Member Data Documentation

◆ b_max

double charge_asymmetry_tt_diff_mtt_LO::b_max = 0.
private

Definition at line 443 of file TopQuarkObservables.h.

◆ b_min

double charge_asymmetry_tt_diff_mtt_LO::b_min = 0.
private

Definition at line 442 of file TopQuarkObservables.h.

◆ mytopobs

const TopQuarkObservables charge_asymmetry_tt_diff_mtt_LO::mytopobs
private

Definition at line 446 of file TopQuarkObservables.h.

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