muggHpttHptHpbbH_Hmumu Class Reference

#include <HiggsThObservables.h>

Inheritance diagram for muggHpttHptHpbbH_Hmumu:

Detailed Description

Author

HEPfit Collaboration

Copyright

GNU General Public License

Definition at line 4392 of file HiggsThObservables.h.

Public Member Functions
double	computeThValue ()
	muggHpttHptHpbbH_Hmumu (const StandardModel &SM_i, const double sqrt_s_i)
	Constructor.
Public Member Functions inherited from ThObservable
double	getBinMax ()
	A get method to get the maximum value of the bin.
double	getBinMin ()
	A get method to get the minimum value of the bin.
const StandardModel &	getModel ()
	A get method to get the model.
const std::vector< std::string >	getParametersForObservable ()
	A get method to get the parameters for the specific observable.
void	setBinMax (double max)
	A set method to set the maximum value of the bin.
void	setBinMin (double min)
	A set method to set the minimum value of the bin.
void	setParametersForObservable (std::vector< std::string > parametersForObservable_i)
	A set method to get the parameters for the specific observable.
	ThObservable (const StandardModel &SM_i)
	Constructor.
	ThObservable (const ThObservable &orig)
	The copy constructor.
virtual	~ThObservable ()
	The default destructor.

Private Attributes
const NPbase *	myNPbase
const double	sqrt_s

Additional Inherited Members
Protected Attributes inherited from ThObservable
double	max
	the bin maximum.
double	min
	The bin minimum.
std::vector< std::string >	parametersForObservable
	a vector of parameter namesfor the specific observable
const StandardModel &	SM
	A reference to an object of StandardMode class.

Constructor & Destructor Documentation

muggHpttHptHpbbH_Hmumu()

muggHpttHptHpbbH_Hmumu::muggHpttHptHpbbH_Hmumu	(	const StandardModel &	SM_i,
		const double	sqrt_s_i
	)

Constructor.

Parameters

[in]	SM_i	a reference to a StandardModel object or to any extension of it
[in]	sqrt_s_i	the center-of-mass energy in TeV

Definition at line 2871 of file HiggsThObservables.cpp.

: ThObservable(SM_i), sqrt_s(sqrt_s_i)
{
    if ((myNPbase = dynamic_cast<const NPbase*> (&SM)) == NULL)
        throw std::runtime_error("muggHpttHptHpbbH_Hmumu called with a class whose parent is not NPbase");
}

Member Function Documentation

computeThValue()

double muggHpttHptHpbbH_Hmumu::computeThValue ( )

virtual

Returns

Implements ThObservable.

Definition at line 2878 of file HiggsThObservables.cpp.

{
    //VM:Note that these values are valid for 13 TeV, they are not general
    //We should access the SM function that has all the values (for the 
    //different energies). The values are slightly different, we should 
    //check this. Furthermore, the bbH is not included. In the SM this is
    //very suppressed (and it's probably also the case in the SMEFT) but 
    //in some NP models it may not be the case.Unfortunately, bbH is not
    //obtained in the SMEFT, we could just set delta_mubbH to zero in the
    //SMEFT and add here the general expression.
    //Ref: https://www.hepdata.net/record/ins2104706 Figure2a (symmetrized)
    double xsSM_ggHbbH = 44.745;
    double xsSM_ttH    = 0.4998;
    double xsSM_tH     = 0.084769;
    //double xsSM_ggH = myNPbase->computeSigmaggH(sqrt_s);
    //double xsSM_ttH = myNPbase->computeSigmattH(sqrt_s);
    //double xsSM_tH = myNPbase->computeSigmatH(sqrt_s);
    //double xsSM_bbH = myNPbase->computeSigmabbH(sqrt_s);
 
    
    
    if ((this->getModel()).isModelLinearized() || (this->getModel()).isModelNPquadratic()) {
        //AG: Most general expression including quadratic corrections. 
        //    If FlagQuadraticTerms=false, the quadratic pieces become 0 in NPSMEFTd6General
        double muProd1 = ( xsSM_ggHbbH*(myNPbase->delta_muggH_1(sqrt_s)) + xsSM_ttH*(myNPbase->delta_muttH_1(sqrt_s)) + xsSM_tH*(myNPbase->delta_mutH_1(sqrt_s)) )/(xsSM_ggHbbH+xsSM_ttH+xsSM_tH) ;
        double muProd2 = ( xsSM_ggHbbH*(myNPbase->delta_muggH_2(sqrt_s)) + xsSM_ttH*(myNPbase->delta_muttH_2(sqrt_s)) + xsSM_tH*(myNPbase->delta_mutH_2(sqrt_s)) )/(xsSM_ggHbbH+xsSM_ttH+xsSM_tH) ;
        double dGammaR1 = myNPbase->deltaGammaHmumuRatio1();
        double dGammaR2 = myNPbase->deltaGammaHmumuRatio2();
        
        double dGammaRTot1 = myNPbase->deltaGammaTotalRatio1();
        double dGammaRTot2 = myNPbase->deltaGammaTotalRatio2();
        double Br1 = dGammaR1-dGammaRTot1;
        double Br2 = dGammaR2 -dGammaRTot2 - dGammaR1*dGammaRTot1 + pow(dGammaRTot1,2.0);
 
        double mu;
        
        mu = 1.0 + (muProd1 + Br1) + (muProd2 + Br2 + muProd1*Br1);
        
        if (mu < 0) return std::numeric_limits<double>::quiet_NaN();
        
        return mu; 
        
    } else {
        
        //VM: Just in case someone wants to add directly the production*decay 
        //(which is the observable we fit at the end). Furthermore, the Hbb
        //is not added in the original formula, fine for the SM (probably also
        //for the SMEFT) but not for all NP models.
        double NPmuggHpttHptHpbbH_Hmumu = myNPbase->muggHpttHptHpbbH_Hmumu(sqrt_s);
        if(NPmuggHpttHptHpbbH_Hmumu==1.0){
            return ( (xsSM_ggHbbH*(myNPbase->muggH(sqrt_s))+xsSM_ttH*(myNPbase->muttH(sqrt_s))+xsSM_tH*(myNPbase->mutH(sqrt_s))) / (xsSM_ggHbbH+xsSM_ttH+xsSM_tH)  )*(myNPbase->BrHmumuRatio()) ;
        } else {
            return NPmuggHpttHptHpbbH_Hmumu;
        }  
    }
}

Member Data Documentation

myNPbase

const NPbase* muggHpttHptHpbbH_Hmumu::myNPbase

private

Definition at line 4409 of file HiggsThObservables.h.

sqrt_s

const double muggHpttHptHpbbH_Hmumu::sqrt_s

private

Definition at line 4410 of file HiggsThObservables.h.

---

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

• HiggsThObservables.h
• HiggsThObservables.cpp