MPll Class Reference

A class for the \(M \to P l^+ l^-\) decay. More...

#include <MPll.h>

Collaboration diagram for MPll:

Detailed Description

A class for the \(M \to P l^+ l^-\) decay.

Author

HEPfit Collaboration

Copyright

GNU General Public License

This class is used to compute all the functions needed in order to compute the observables relative to the \(M \to P l^+ l^-\) decay. After the parameters are updated in updateParameters() and the cache is checked in checkCache(), the form factor are build in the transverse basis in the functions f_plus(), f_0() and f_T() [26]. Following [93], they are consequentely translated in the helicity basis through the functions V_L(), V_R(), T_L(), T_R(), S_L() and S_R(). Form factors and parameters are combined together in the functions H_V(), H_A(), H_S() and H_P() in order to build the helicity amplitudes, which are consequentely combined to create the angular coefficients in the function I(). Those coefficients are used to create the CP averaged coefficients in the function Sigma() ad the CP asymmetric coefficients in the function Delta(). Form factors, CP averaged and asymmetric coefficients and hadronic contributions are integrated in the functions integrateSigma() and integrateDelta() in order to be further used to build the observables.

Definition at line 41 of file MPll.h.

Public Member Functions
double	beta (double q2)
	The factor \( \beta \) used in the angular coefficients \(I_i\). More...
void	checkCache ()
	The caching method for MPll. More...
double	Delta (int i, double q2)
	The CP asymmetry \( \Delta_{i} \) . More...
double	F (double q2)
	The factor \( F \) used in the angular coefficients \(I_i\). More...
double	f_0 (double q2)
	The form factor \( f_0 \). More...
double	f_plus (double q2)
	The form factor \( f_+ \). More...
double	f_T (double q2)
	The form factor \( f_T \). More...
double	getDelta0 (double q2)
	The CP asymmetry \( \Delta_{1s} \) . More...
double	getDelta2 (double q2)
	The CP asymmetry \( \Delta_{2s} \) . More...
double	getSigma0 (double q2)
	The CP average \( \Sigma_{1s} \) . More...
double	getSigma2 (double q2)
	The CP average \( \Sigma_{2s} \) . More...
gslpp::complex	H (double q2, double m)
	The \( h(q^2,m) \) function involved into \( C_9^{eff}\). More...
gslpp::complex	H_A (double q2, int bar)
	The helicity amplitude \( H_A^{\lambda} \) . More...
gslpp::complex	H_P (double q2, int bar)
	The helicity amplitude \( H_P^{\lambda} \) . More...
gslpp::complex	H_S (double q2, int bar)
	The helicity amplitude \( H_S^{\lambda} \) . More...
gslpp::complex	H_V (double q2, int bar)
	The helicity amplitude \( H_V^{\lambda} \) . More...
double	I (int i, double q2, int bar)
	The angular coefficient \( I_{i} \) . More...
double	integrateDelta (int i, double q_min, double q_max)
	The integral of \( \Delta_{i} \) from \(q_{min}\) to \(q_{max}\). More...
double	integrateSigma (int i, double q_min, double q_max)
	The integral of \( \Sigma_{i} \) from \(q_{min}\) to \(q_{max}\). More...
double	k2 (double q2)
	The square of the 3-momentum of the recoiling meson in the M rest frame, \( k^2 \) . More...
double	lambda (double q2)
	The factor \( \lambda \) used in the angular coefficients \(I_i\). More...
double	LCSR_fit1 (double q2, double r_1, double r_2, double m_fit2)
	The second fit function from arXiv:hep-ph/0412079v1, \( f_2^{LCSR} \). More...
double	LCSR_fit2 (double q2, double r_2, double m_fit2)
	The third fit function from arXiv:hep-ph/0412079v1, \( f_3^{LCSR} \). More...
	MPll (const StandardModel &SM_i, StandardModel::meson meson_i, StandardModel::meson pseudoscalar_i, StandardModel::lepton lep_i)
	Constructor. More...
gslpp::complex	N ()
	The helicity amplitudes normalization factor \( N \) . More...
double	S_L (double q2)
	The helicity form factor \( S_L \). More...
double	S_R (double q2)
	The helicity form factor \( S_R \). More...
double	Sigma (int i, double q2)
	The CP average \( \Sigma_{i} \) . More...
gslpp::complex	T_L (double q2)
	The helicity form factor \( T_L^0 \). More...
gslpp::complex	T_R (double q2)
	The helicity form factor \( T_R^0 \). More...
void	updateParameters ()
	The update parameter method for MPll. More...
gslpp::complex	V_L (double q2)
	The helicity form factor \( V_L^0 \). More...
gslpp::complex	V_R (double q2)
	The helicity form factor \( V_R^0 \). More...
gslpp::complex	Y (double q2)
	The \( Y(q^2) \) function involved into \( C_9^{eff}\). More...
virtual	~MPll ()
	Destructor. More...

Public Attributes
double	ale
gslpp::vector< gslpp::complex > **	allcoeff
gslpp::vector< gslpp::complex > **	allcoeffprime
double	b
gslpp::complex	C_1
gslpp::complex	C_10
gslpp::complex	C_10p
gslpp::complex	C_2
gslpp::complex	C_3
gslpp::complex	C_4
gslpp::complex	C_5
gslpp::complex	C_6
gslpp::complex	C_7
gslpp::complex	C_7p
gslpp::complex	C_9
gslpp::complex	C_9p
gslpp::complex	C_P
gslpp::complex	C_Pp
gslpp::complex	C_S
gslpp::complex	C_Sp
double	GF
gslpp::complex	h_0
gslpp::complex	h_0_1
gslpp::complex	lambda_t
double	m_fit2_f0
double	m_fit2_fplus
double	m_fit2_fT
double	Mb
double	Mc
double	Mlep
double	MM
double	MP
double	Ms
double	mu_b
double	MW
double	q2
double	r_1_fplus
double	r_1_fT
double	r_2_f0
double	r_2_fplus
double	r_2_fT
double	width

Private Attributes
double	avaDelta0
double	avaDelta2
double	avaSigma0
double	avaSigma2
double	beta_cache
unsigned int	beta_updated
gslpp::complex	C_10_cache
unsigned int	C_10_updated
gslpp::complex	C_10p_cache
unsigned int	C_10p_updated
gslpp::complex	C_1_cache
unsigned int	C_1_updated
gslpp::complex	C_2_cache
unsigned int	C_2_updated
gslpp::complex	C_3_cache
unsigned int	C_3_updated
gslpp::complex	C_4_cache
unsigned int	C_4_updated
gslpp::complex	C_5_cache
unsigned int	C_5_updated
gslpp::complex	C_6_cache
unsigned int	C_6_updated
gslpp::complex	C_7_cache
unsigned int	C_7_updated
gslpp::complex	C_7p_cache
unsigned int	C_7p_updated
gslpp::complex	C_9_cache
unsigned int	C_9_updated
gslpp::complex	C_9p_cache
unsigned int	C_9p_updated
gslpp::complex	C_P_cache
unsigned int	C_P_updated
gslpp::complex	C_Pp_cache
unsigned int	C_Pp_updated
gslpp::complex	C_S_cache
unsigned int	C_S_updated
gslpp::complex	C_Sp_cache
unsigned int	C_Sp_updated
std::map< std::pair< double, double >, double >	cacheDelta0
std::map< std::pair< double, double >, double >	cacheDelta2
std::map< std::pair< double, double >, double >	cacheSigma0
std::map< std::pair< double, double >, double >	cacheSigma2
std::map< std::pair< double, double >, unsigned int >	delta0Cached
std::map< std::pair< double, double >, unsigned int >	delta2Cached
double	errDelta0
double	errDelta2
double	errSigma0
double	errSigma2
double	f0_cache
unsigned int	f0_updated
unsigned int	F_updated
gsl_function	FD0
gsl_function	FD2
gslpp::vector< double >	fplus_cache
unsigned int	fplus_updated
gsl_function	FS0
gsl_function	FS2
gslpp::vector< double >	fT_cache
unsigned int	fT_updated
unsigned int	H_A0updated
gslpp::vector< double >	H_P_cache
unsigned int	H_P_updated
gslpp::vector< double >	H_Scache
unsigned int	H_Supdated
gslpp::vector< double >	H_V0cache
gslpp::complex	H_V0Ccache [2]
unsigned int	H_V0updated
unsigned int	I0_updated
unsigned int	I2_updated
unsigned int	I8_updated
gslpp::vector< double >	k2_cache
unsigned int	k2_updated
double	lambda_cache
unsigned int	lambda_updated
StandardModel::lepton	lep
StandardModel::meson	meson
const StandardModel &	mySM
gslpp::vector< double >	N_cache
unsigned int	N_updated
gslpp::complex	Nc_cache
StandardModel::meson	pseudoscalar
std::map< std::pair< double, double >, unsigned int >	sigma0Cached
std::map< std::pair< double, double >, unsigned int >	sigma2Cached
gslpp::vector< double >	SL_cache
unsigned int	SL_updated
unsigned int	SR_updated
unsigned int	TL_updated
unsigned int	TR_updated
unsigned int	VL_updated
unsigned int	VR_updated
gsl_integration_workspace *	w_delta0
gsl_integration_workspace *	w_delta2
gsl_integration_workspace *	w_sigma0
gsl_integration_workspace *	w_sigma2
gslpp::vector< double >	Ycache
unsigned int	Yupdated

Constructor & Destructor Documentation

MPll::MPll	(	const StandardModel &	SM_i,
		StandardModel::meson	meson_i,
		StandardModel::meson	pseudoscalar_i,
		StandardModel::lepton	lep_i
	)

Constructor.

Parameters

[in]	SM_i	a reference to an object of type StandardModel
[in]	meson_i	initial meson of the decay
[in]	pseudoscalar_i	final pseudoscalar meson of the decay
[in]	lep_i	final leptons of the decay

Definition at line 19 of file MPll.cpp.

:       mySM(SM_i),
        fplus_cache(2, 0.),
        fT_cache(2, 0.),
        k2_cache(2, 0.),
        SL_cache(2, 0.),
        N_cache(3, 0.),
        Ycache(3, 0.),
        H_V0cache(2, 0.),
        H_Scache(2, 0.),
        H_P_cache(4, 0.)
{
    lep = lep_i;
    meson = meson_i;
    pseudoscalar = pseudoscalar_i;
    I0_updated = 0;
    I2_updated = 0;
    I8_updated = 0;
    
    w_sigma0 = gsl_integration_workspace_alloc (50);
    w_sigma2 = gsl_integration_workspace_alloc (50);
    
    w_delta0 = gsl_integration_workspace_alloc (50);
    w_delta2 = gsl_integration_workspace_alloc (50);
}

MPll::~MPll ( )

virtual

Destructor.

Definition at line 46 of file MPll.cpp.

{}

Member Function Documentation

double MPll::beta

(

double

q2

)

The factor \( \beta \) used in the angular coefficients \(I_i\).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \beta \)

Definition at line 528 of file MPll.cpp.

{
    return sqrt(1. - 4.*Mlep*Mlep/q2);
}

void MPll::checkCache

(

)

The caching method for MPll.

Definition at line 119 of file MPll.cpp.

{
    
    if(r_1_fplus == fplus_cache(0) && r_2_fplus == fplus_cache(1)) {
        fplus_updated = 1;
    } else {
        fplus_updated = 0;
        fplus_cache(0) = r_1_fplus;
        fplus_cache(1) = r_2_fplus;
    }
    
    if(r_1_fT == fT_cache(0) && r_2_fT == fT_cache(1)) {
        fT_updated = 1;
    } else {
        fT_updated = 0;
        fT_cache(0) = r_1_fT;
        fT_cache(1) = r_2_fT;
    }
    
    if(r_2_f0 == f0_cache) {
        f0_updated = 1;
    } else {
        f0_updated = 0;
        f0_cache = r_2_f0;
    }
    
    if (MM == k2_cache(0) && MP == k2_cache(1) ) {
        k2_updated = 1;
    } else {
        k2_updated = 0;
        k2_cache(0) = MM;
        k2_cache(1) = MP;
    }
    
    if (Mlep == beta_cache) {
        beta_updated = 1;
    } else {
        beta_updated = 0;
        beta_cache = Mlep;
    }
    
    if (MP == lambda_cache) {
        lambda_updated = k2_updated;
        F_updated = lambda_updated * beta_updated;
    } else {
        lambda_updated = 0;
        F_updated = 0;
        lambda_cache = MP;
    }
    
    VL_updated = k2_updated * fplus_updated;
    
    VR_updated = VL_updated;
    
    TL_updated = k2_updated * fT_updated;
    
    TR_updated = TL_updated;
    
    if (Mb == SL_cache(0) && Ms == SL_cache(1) ){
        SL_updated = k2_updated * f0_updated;
        SR_updated = SL_updated;
    } else {
        SL_updated = 0;
        SR_updated = SL_updated;
        SL_cache(0) = Mb;
        SL_cache(1) = Ms;
    }
    
    if (GF == N_cache(0) && ale == N_cache(1) && MM == N_cache(2) && lambda_t == Nc_cache) {
        N_updated = 1;
    } else {
        N_updated = 0;
        N_cache(0) = GF;
        N_cache(1) = ale;
        N_cache(2) = MM;
        Nc_cache = lambda_t;
    }
    
    if (C_1 == C_1_cache) {
        C_1_updated = 1;
    } else {
        C_1_updated = 0;
        C_1_cache = C_1;
    }
    
    if (C_2 == C_2_cache) {
        C_2_updated = 1;
    } else {
        C_2_updated = 0;
        C_2_cache = C_2;
    }
    
    if (C_3 == C_3_cache) {
        C_3_updated = 1;
    } else {
        C_3_updated = 0;
        C_3_cache = C_3;
    }
    
    if (C_4 == C_4_cache) {
        C_4_updated = 1;
    } else {
        C_4_updated = 0;
        C_4_cache = C_4;
    }
    
    if (C_5 == C_5_cache) {
        C_5_updated = 1;
    } else {
        C_5_updated = 0;
        C_5_cache = C_5;
    }
    
    if (C_6 == C_6_cache) {
        C_6_updated = 1;
    } else {
        C_6_updated = 0;
        C_6_cache = C_6;
    }
    
    if (C_7 == C_7_cache) {
        C_7_updated = 1;
    } else {
        C_7_updated = 0;
        C_7_cache = C_7;
    }
    
    if (C_9 == C_9_cache) {
        C_9_updated = 1;
    } else {
        C_9_updated = 0;
        C_9_cache = C_9;
    }
    
    if (C_10 == C_10_cache) {
        C_10_updated = 1;
    } else {
        C_10_updated = 0;
        C_10_cache = C_10;
    }
    
    if (C_S == C_S_cache) {
        C_S_updated = 1;
    } else {
        C_S_updated = 0;
        C_S_cache = C_S;
    }
    
    if (C_P == C_P_cache) {
        C_P_updated = 1;
    } else {
        C_P_updated = 0;
        C_P_cache = C_P;
    }
    
    if (C_7p == C_7p_cache) {
        C_7p_updated = 1;
    } else {
        C_7p_updated = 0;
        C_7p_cache = C_7p;
    }
    
    if (C_9p == C_9p_cache) {
        C_9p_updated = 1;
    } else {
        C_9p_updated = 0;
        C_9p_cache = C_9p;
    }
    
    if (C_10p == C_10p_cache) {
        C_10p_updated = 1;
    } else {
        C_10p_updated = 0;
        C_10p_cache = C_10p;
    }
    
    if (C_Sp == C_Sp_cache) {
        C_Sp_updated = 1;
    } else {
        C_Sp_updated = 0;
        C_Sp_cache = C_Sp;
    }
    
    if (C_Pp == C_Pp_cache) {
        C_Pp_updated = 1;
    } else {
        C_Pp_updated = 0;
        C_Pp_cache = C_Pp;
    }
    
    if (Mb == Ycache(0) && Mc == Ycache(1) ) {
        Yupdated = C_1_updated * C_2_updated * C_3_updated * C_4_updated * C_5_updated * C_6_updated;
    } else {
        Yupdated = 0;
        Ycache(0) = Mb;
        Ycache(1) = Mc;
    }
    
    if (MM == H_V0cache(0) && Mb == H_V0cache(1) && h_0 == H_V0Ccache[0] && h_0_1 == H_V0Ccache[1]) {
        H_V0updated = N_updated * C_9_updated * Yupdated * VL_updated * C_9p_updated * VR_updated * C_7_updated * TL_updated * C_7p_updated * TR_updated;
    } else {
        H_V0updated = 0;
        H_V0cache(0) = MM;
        H_V0cache(1) = Mb;
        H_V0Ccache[0] = h_0;
        H_V0Ccache[1] = h_0_1;
    }
    
    H_A0updated = N_updated * C_10_updated * VL_updated * C_10p_updated * VR_updated;
    
    if (Mb == H_Scache(0) && MW == H_Scache(1)) {
        H_Supdated = N_updated * C_S_updated * SL_updated * C_Sp_updated * SR_updated;
    } else {
        H_Supdated = 0;
        H_Scache(0) = Mb;
        H_Scache(1) = MW;
    }
    
    if (Mb == H_P_cache(0) && MW == H_P_cache(1) && Mlep == H_P_cache(2) && Ms == H_P_cache(3)) {
        H_P_updated = N_updated * C_P_updated * SL_updated * C_Pp_updated * SR_updated * C_10_updated * C_10p_updated;
    } else {
        H_P_updated = 0;
        H_P_cache(0) = Mb;
        H_P_cache(1) = MW;
        H_P_cache(2) = Mlep;
        H_P_cache(3) = Ms;   
    }
    
    I0_updated = F_updated * H_V0updated * H_A0updated * H_P_updated * beta_updated * H_Supdated;
    I2_updated = F_updated * beta_updated * H_V0updated * H_A0updated;
    I8_updated = F_updated * beta_updated * H_Supdated * H_V0updated;
    
}

double MPll::Delta	(	int	i,
		double	q2
	)

The CP asymmetry \( \Delta_{i} \) .

Parameters

[in]	i	index of the angular coefficient \( I_{i} \)
[in]	q2	\(q^2\) of the decay

Returns

\( \Delta_{i} \)

Definition at line 570 of file MPll.cpp.

{
    return (I(i, q2, 0) - I(i, q2, 1))/2.;
}

double MPll::F

(

double

q2

)

The factor \( F \) used in the angular coefficients \(I_i\).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( F \)

Definition at line 538 of file MPll.cpp.

{
    return sqrt(lambda(q2))*beta(q2)*q2/(96.*M_PI*M_PI*M_PI*MM*MM*MM);
}

double MPll::f_0

(

double

q2

)

The form factor \( f_0 \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( f_0 \)

Definition at line 377 of file MPll.cpp.

{
    return LCSR_fit2(q2, r_2_f0, m_fit2_f0);
}

double MPll::f_plus

(

double

q2

)

The form factor \( f_+ \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( f_+ \)

Definition at line 367 of file MPll.cpp.

{
    return LCSR_fit1(q2, r_1_fplus, r_2_fplus, m_fit2_fplus);
}

double MPll::f_T

(

double

q2

)

The form factor \( f_T \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( f_T \)

Definition at line 372 of file MPll.cpp.

{
    return LCSR_fit1(q2, r_1_fT, r_2_fT, m_fit2_fT);
}

double MPll::getDelta0 ( double  q2 )

inline

The CP asymmetry \( \Delta_{1s} \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Delta_{1s} \)

Definition at line 352 of file MPll.h.

    {
        return Delta(0, q2);
    };

double MPll::getDelta2 ( double  q2 )

inline

The CP asymmetry \( \Delta_{2s} \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Delta_{2s} \)

Definition at line 362 of file MPll.h.

    {
        return Delta(2, q2);
    };

double MPll::getSigma0 ( double  q2 )

inline

The CP average \( \Sigma_{1s} \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Sigma_{1s} \)

Definition at line 332 of file MPll.h.

    {
        return Sigma(0, q2);
    };

double MPll::getSigma2 ( double  q2 )

inline

The CP average \( \Sigma_{2s} \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Sigma_{2s} \)

Definition at line 342 of file MPll.h.

    {
        return Sigma(2, q2);
    };

gslpp::complex MPll::H	(	double	q2,
		double	m
	)

The \( h(q^2,m) \) function involved into \( C_9^{eff}\).

Parameters

[in]	q2	\(q^2\) of the decay
[in]	m	mass

Returns

\( h(q^2,m) \)

Definition at line 422 of file MPll.cpp.

{
    double x = 4.*m*m/q2;
    gslpp::complex par;
    
    if (x>1.) par = sqrt(x - 1.) * atan( 1. / sqrt(x - 1.) );
    else par = sqrt(1. - x) * ( log( ( 1. + sqrt(1. - x) ) / sqrt(x) ) - gslpp::complex::i()*M_PI/2.);
    
    if (x == 0.) return 8. / 27. * (1. + 3. * gslpp::complex::i() * M_PI);
    else return - 4./9. * ( log( m*m / q2 ) - 2./3. - x ) - 4./9. * (-2. + x) * par;
}

gslpp::complex MPll::H_A	(	double	q2,
		int	bar
	)

The helicity amplitude \( H_A^{\lambda} \) .

Parameters

[in]	q2	\(q^2\) of the decay
[in]	bar	index to choose betwen regular coefficient (bar=0) and conjugated coefficient (bar=1)

Returns

\( H_A^{\lambda} \)

Definition at line 460 of file MPll.cpp.

{
    gslpp::complex n;
    switch(bar){
        case 0:
            n = N();
            break;
        case 1:
            n = N().conjugate();
            break;
        default:
            std::stringstream out;
            out << bar;
            throw std::runtime_error("H_A: index " + out.str() + " not allowed for an Angular Coefficient");
    }
     
    return -gslpp::complex::i()*n*( C_10*V_L(q2) + C_10p*V_R(q2) );
}

gslpp::complex MPll::H_P	(	double	q2,
		int	bar
	)

The helicity amplitude \( H_P^{\lambda} \) .

Parameters

[in]	q2	\(q^2\) of the decay
[in]	bar	index to choose betwen regular coefficient (bar=0) and conjugated coefficient (bar=1)

Returns

\( H_P^{\lambda} \)

Definition at line 498 of file MPll.cpp.

{
    gslpp::complex n;
    switch(bar){
        case 0:
            n = N();
            break;
        case 1:
            n = N().conjugate();
            break;
        default:
            std::stringstream out;
            out << bar;
            throw std::runtime_error("H_S: index " + out.str() + " not allowed for an Angular Coefficient");
    }
     
    return gslpp::complex::i()*n*( Mb/MW*( C_P*S_L(q2) + C_Pp*S_R(q2) ) 
            + 2.*Mlep*Mb/q2*( C_10*( S_L(q2) - Ms/Mb*S_R(q2) ) + C_10p*( S_R(q2) - Ms/Mb*S_L(q2) ) ) );
}

gslpp::complex MPll::H_S	(	double	q2,
		int	bar
	)

The helicity amplitude \( H_S^{\lambda} \) .

Parameters

[in]	q2	\(q^2\) of the decay
[in]	bar	index to choose betwen regular coefficient (bar=0) and conjugated coefficient (bar=1)

Returns

\( H_S^{\lambda} \)

Definition at line 479 of file MPll.cpp.

{
    gslpp::complex n;
    switch(bar){
        case 0:
            n = N();
            break;
        case 1:
            n = N().conjugate();
            break;
        default:
            std::stringstream out;
            out << bar;
            throw std::runtime_error("H_S: index " + out.str() + " not allowed for an Angular Coefficient");
    }
     
    return gslpp::complex::i()*n*Mb/MW*( C_S*S_L(q2) + C_Sp*S_R(q2) );
}

gslpp::complex MPll::H_V	(	double	q2,
		int	bar
	)

The helicity amplitude \( H_V^{\lambda} \) .

Parameters

[in]	q2	\(q^2\) of the decay
[in]	bar	index to choose betwen regular coefficient (bar=0) and conjugated coefficient (bar=1)

Returns

\( H_V^{\lambda} \)

Definition at line 440 of file MPll.cpp.

{
    gslpp::complex n;
    switch(bar){
        case 0:
            n = N();
            break;
        case 1:
            n = N().conjugate();
            break;
        default:
            std::stringstream out;
            out << bar;
            throw std::runtime_error("H_V: index " + out.str() + " not allowed for an Angular Coefficient");
    }
                    
    return -gslpp::complex::i()*n*( ( C_9 + Y(q2) )*V_L(q2) 
            + C_9p*V_R(q2) + MM*MM/q2*( 2*Mb/MM*( C_7*T_L(q2) + C_7p*T_R(q2) ) - 16*M_PI*M_PI*(h_0 + h_0_1 * q2)) );
}

double MPll::I	(	int	i,
		double	q2,
		int	bar
	)

The angular coefficient \( I_{i} \) .

Parameters

[in]	i	index of the angular coefficient: 0 for 1c, 2 for 2c, 8 for 6c
[in]	q2	\(q^2\) of the decay
[in]	bar	index to choose betwen regular coefficient (bar=0) and conjugated coefficient (bar=1)

Returns

\( I_{i} \)

Definition at line 543 of file MPll.cpp.

{

    double Mlep2 = Mlep*Mlep;
    double beta2 = beta(q2)*beta(q2);
    

    switch (i){
        case 0: // I1c
            return F(q2)*( ( H_V(q2,bar).abs2() + H_A(q2,bar).abs2() )/2.  +  H_P(q2,bar).abs2()  +  2.*Mlep2/q2*( H_V(q2,bar).abs2() 
                    - H_A(q2,bar).abs2() )  + beta2*H_S(q2,bar).abs2() );
        case 2: // I2c
            return -F(q2)*beta2/2.*( H_V(q2,bar).abs2() + H_A(q2,bar).abs2() );
        case 8: // I6c
            return 2.*F(q2)*beta(q2)*Mlep/sqrt(q2)*( H_S(q2,bar).conjugate()*H_V(q2,bar) ).real();
        default:
            std::stringstream out;
            out << i;
            throw std::runtime_error("I: index " + out.str() + " not implemented");
    }
}

double MPll::integrateDelta	(	int	i,
		double	q_min,
		double	q_max
	)

The integral of \( \Delta_{i} \) from \(q_{min}\) to \(q_{max}\).

Parameters

[in]	i	index of the angular coefficient \( I_{i} \)
[in]	q_min	minimum q^2 of the integral
[in]	q_max	maximum q^2 of the integral

Returns

\( <\Delta_{i}> \)

Definition at line 612 of file MPll.cpp.

{
    
    if (mySM.getMyFlavour()->getUpdateFlag(meson, pseudoscalar, lep)){
        updateParameters();
        mySM.getMyFlavour()->setUpdateFlag(meson, pseudoscalar, lep, false);
    }
        
    std::pair<double, double > qbin = std::make_pair(q_min, q_max);
    switch(i){
        case 0:
            if (delta0Cached[qbin] == 0) {
                FD0 = convertToGslFunction( boost::bind( &MPll::getDelta0, &(*this), _1 ) );
                gsl_integration_qags (&FD0, q_min, q_max, 1.e-5, 1.e-3, 50, w_delta0, &avaDelta0, &errDelta0);
                cacheDelta0[qbin] = avaDelta0;
                delta0Cached[qbin] = 1;
            }
            return cacheDelta0[qbin];
            break;
        case 2:
            if (delta2Cached[qbin] == 0) {
                FD2 = convertToGslFunction( boost::bind( &MPll::getDelta2, &(*this), _1 ) );
                gsl_integration_qags (&FD2, q_min, q_max, 1.e-5, 1.e-3, 50, w_delta2, &avaDelta2, &errDelta2);
                cacheDelta2[qbin] = avaDelta2;
                delta2Cached[qbin] = 1;
            }
            return cacheDelta2[qbin];
            break;
        default:
            std::stringstream out;
            out << i;
            throw std::runtime_error("integrateDelta_e: index " + out.str() + " not implemented"); 
    }
}

double MPll::integrateSigma	(	int	i,
		double	q_min,
		double	q_max
	)

The integral of \( \Sigma_{i} \) from \(q_{min}\) to \(q_{max}\).

Parameters

[in]	i	index of the angular coefficient \( I_{i} \)
[in]	q_min	minimum q^2 of the integral
[in]	q_max	maximum q^2 of the integral

Returns

\( <\Sigma_{i}> \)

Definition at line 575 of file MPll.cpp.

{
    
    if (mySM.getMyFlavour()->getUpdateFlag(meson, pseudoscalar, lep)){
        updateParameters();
        mySM.getMyFlavour()->setUpdateFlag(meson, pseudoscalar, lep, false);
    }
    
    std::pair<double, double > qbin = std::make_pair(q_min, q_max);
    
    switch(i){
        case 0:
            if (sigma0Cached[qbin] == 0) {
                FS0 = convertToGslFunction( boost::bind( &MPll::getSigma0, &(*this), _1 ) );
                gsl_integration_qags (&FS0, q_min, q_max, 1.e-5, 1.e-3, 50, w_sigma0, &avaSigma0, &errSigma0);
                cacheSigma0[qbin] = avaSigma0;
                sigma0Cached[qbin] = 1;
            }
            return cacheSigma0[qbin];
            break;
        case 2:
            if (sigma2Cached[qbin] == 0) {
                FS2 = convertToGslFunction( boost::bind( &MPll::getSigma2, &(*this), _1 ) );
                gsl_integration_qags (&FS2, q_min, q_max, 1.e-5, 1.e-3, 50, w_sigma2, &avaSigma2, &errSigma2);
                cacheSigma2[qbin] = avaSigma2;
                sigma2Cached[qbin] = 1;
            }
            return cacheSigma2[qbin];
            break;
        default:
            std::stringstream out;
            out << i;
            throw std::runtime_error("MVll::integrateSigma_e: index " + out.str() + " not implemented");
    }
    
}

double MPll::k2

(

double

q2

)

The square of the 3-momentum of the recoiling meson in the M rest frame, \( k^2 \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( k^2 \)

Definition at line 523 of file MPll.cpp.

{
    return (pow(MM,4.) + q2*q2 + pow(MP,4.) -2.*MP*MP*q2 -2.*MM*MM*(q2 + MP*MP))/(4.*MM*MM);
}

double MPll::lambda

(

double

q2

)

The factor \( \lambda \) used in the angular coefficients \(I_i\).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \lambda \)

Definition at line 533 of file MPll.cpp.

{
    return 4.*MM*MM*k2(q2);
}

double MPll::LCSR_fit1	(	double	q2,
		double	r_1,
		double	r_2,
		double	m_fit2
	)

The second fit function from arXiv:hep-ph/0412079v1, \( f_2^{LCSR} \).

Parameters

[in]	q2	\(q^2\) of the decay
[in]	r_1	fit parameter
[in]	r_2	fit parameter
[in]	m_fit2	fit parameter

Returns

\( f_2^{LCSR} \)

Definition at line 356 of file MPll.cpp.

{
    return r_1/( 1. - q2/m_fit2 ) + r_2/pow( ( 1. - q2/m_fit2 ) , 2.) ;

}

double MPll::LCSR_fit2	(	double	q2,
		double	r_2,
		double	m_fit2
	)

The third fit function from arXiv:hep-ph/0412079v1, \( f_3^{LCSR} \).

Parameters

[in]	q2	\(q^2\) of the decay
[in]	r_2	fit parameter
[in]	m_fit2	fit parameter

Returns

\( f_3^{LCSR} \)

Definition at line 362 of file MPll.cpp.

{
    return r_2/( 1. - q2/m_fit2 ) ; 
}

gslpp::complex MPll::N

(

)

The helicity amplitudes normalization factor \( N \) .

Returns

\( N \)

Definition at line 417 of file MPll.cpp.

{
    return -(4.*GF*MM*ale*lambda_t)/(sqrt(2.)*4.*M_PI);
}

double MPll::S_L

(

double

q2

)

The helicity form factor \( S_L \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( S_L \)

Definition at line 402 of file MPll.cpp.

{
    return -( MM*MM - MP*MP )/(2*MM*(Mb + Ms)) * ( 1 + Ms/Mb )/( 1 - Ms/Mb ) * f_0(q2);
}

double MPll::S_R

(

double

q2

)

The helicity form factor \( S_R \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( S_R \)

Definition at line 407 of file MPll.cpp.

{
    return -S_L(q2);
}

double MPll::Sigma	(	int	i,
		double	q2
	)

The CP average \( \Sigma_{i} \) .

Parameters

[in]	i	index of the angular coefficient \( I_{i} \)
[in]	q2	\(q^2\) of the decay

Returns

\( \Sigma_{i} \)

Definition at line 565 of file MPll.cpp.

{
    return (I(i, q2, 0) + I(i, q2, 1))/2.;
}

gslpp::complex MPll::T_L

(

double

q2

)

The helicity form factor \( T_L^0 \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( T_L^{\lambda} \)

Definition at line 392 of file MPll.cpp.

{
    return gslpp::complex::i()  * sqrt(lambda(q2)*q2) / (2.*MM*MM*(MM+MP)) * f_T(q2);
}

gslpp::complex MPll::T_R

(

double

q2

)

The helicity form factor \( T_R^0 \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( T_R^{\lambda} \)

Definition at line 397 of file MPll.cpp.

{
    return -T_L(q2);
}

void MPll::updateParameters

(

)

The update parameter method for MPll.

Definition at line 49 of file MPll.cpp.

{
    GF = mySM.getGF();
    ale=mySM.getAle();
    Mlep=mySM.getLeptons(lep).getMass();
    MM=mySM.getMesons(meson).getMass();
    MP=mySM.getMesons(pseudoscalar).getMass();
    Mb=mySM.getQuarks(QCD::BOTTOM).getMass();    // add the PS b mass
    Mc=mySM.getQuarks(QCD::CHARM).getMass();
    Ms=mySM.getQuarks(QCD::STRANGE).getMass();
    MW=mySM.Mw();
    lambda_t=mySM.computelamt_s();
    mu_b = mySM.getMub();
    width = mySM.getMesons(meson).computeWidth();
    
    switch(pseudoscalar){
        case StandardModel::K_P :
            r_1_fplus = mySM.getr_1_fplus();
            r_2_fplus = mySM.getr_2_fplus();
            m_fit2_fplus = mySM.getm_fit2_fplus();
            r_1_fT = mySM.getr_1_fT();
            r_2_fT = mySM.getr_2_fT();
            m_fit2_fT = mySM.getm_fit2_fT();
            r_2_f0 = mySM.getr_2_f0();
            m_fit2_f0 = mySM.getm_fit2_f0();
            break;
        default:
            std::stringstream out;
            out << pseudoscalar;
            throw std::runtime_error("MPll: pseudoscalar " + out.str() + " not implemented");
    }
    
    h_0 = mySM.geth_0_MP();
    h_0_1 = mySM.geth_0_1_MP();
    
    b = 1.;                           //please check
    
    allcoeff = mySM.getMyFlavour()->ComputeCoeffBMll(mu_b);   //check the mass scale, scheme fixed to NDR
    allcoeffprime = mySM.getMyFlavour()->ComputeCoeffprimeBMll(mu_b);   //check the mass scale, scheme fixed to NDR

    C_1 = (*(allcoeff[LO]))(0) + (*(allcoeff[NLO]))(0);
    C_2 =  ((*(allcoeff[LO]))(1) + (*(allcoeff[NLO]))(1));
    C_3 = ((*(allcoeff[LO]))(2) + (*(allcoeff[NLO]))(2));
    C_4 = (*(allcoeff[LO]))(3) + (*(allcoeff[NLO]))(3);
    C_5 =  ((*(allcoeff[LO]))(4) + (*(allcoeff[NLO]))(4));
    C_6 = ((*(allcoeff[LO]))(5) + (*(allcoeff[NLO]))(5));
    C_7 = (*(allcoeff[LO]))(6) + (*(allcoeff[NLO]))(6);
    C_9 = (*(allcoeff[LO]))(8) + (*(allcoeff[NLO]))(8);
    C_10 = (*(allcoeff[LO]))(9) + (*(allcoeff[NLO]))(9);
    C_S = (*(allcoeff[LO]))(10) + (*(allcoeff[NLO]))(10);
    C_P = (*(allcoeff[LO]))(11) + (*(allcoeff[NLO]))(11);
    
    C_7p = (*(allcoeffprime[LO]))(6) + (*(allcoeffprime[NLO]))(6);
    C_9p = (*(allcoeffprime[LO]))(8) + (*(allcoeffprime[NLO]))(8);
    C_10p = (*(allcoeffprime[LO]))(9) + (*(allcoeffprime[NLO]))(9);
    C_Sp = (*(allcoeffprime[LO]))(10) + (*(allcoeffprime[NLO]))(10);
    C_Pp = (*(allcoeffprime[LO]))(11) + (*(allcoeffprime[NLO]))(11);
    
    checkCache();
    
    std::map<std::pair<double, double>, unsigned int >::iterator it;
    
    if (I0_updated == 0) for (it = sigma0Cached.begin(); it != sigma0Cached.end(); ++it) it->second = 0;
    if (I2_updated == 0) for (it = sigma2Cached.begin(); it != sigma2Cached.end(); ++it) it->second = 0;
    
    if (I0_updated == 0) for (it = delta0Cached.begin(); it != delta0Cached.end(); ++it) it->second = 0;
    if (I2_updated == 0) for (it = delta2Cached.begin(); it != delta2Cached.end(); ++it) it->second = 0;
    
}

gslpp::complex MPll::V_L

(

double

q2

)

The helicity form factor \( V_L^0 \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( V_L^{\lambda} \)

Definition at line 382 of file MPll.cpp.

{
    return gslpp::complex::i() * sqrt(lambda(q2)) / (2*MM*sqrt(q2)) * f_plus(q2);
}

gslpp::complex MPll::V_R

(

double

q2

)

The helicity form factor \( V_R^0 \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( V_R^{\lambda} \)

Definition at line 387 of file MPll.cpp.

{
    return -V_L(q2);
}

gslpp::complex MPll::Y

(

double

q2

)

The \( Y(q^2) \) function involved into \( C_9^{eff}\).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( Y(q^2) \)

Definition at line 434 of file MPll.cpp.

{
    return 4./3. * C_3 + 64./9. * C_5 + 64./27. * C_6 - 1./2. * H(q2,0.) * ( C_3 + 4./3.*C_4 + 16. * C_5 + 64./3.*C_6 )
            + H(q2, Mc) * ( 4./3.*C_1 + C_2 + 6.*C_3 + 60.*C_5 ) - 1./2. * H(q2, Mb) * ( 7.*C_3 + 4./3.*C_4 + 76.*C_5 + 64./3.*C_6 );
}

Member Data Documentation

double MPll::ale

alpha electromagnetic

Definition at line 68 of file MPll.h.

gslpp::vector<gslpp::complex>** MPll::allcoeff

vector that contains the Wilson coeffients

Definition at line 95 of file MPll.h.

gslpp::vector<gslpp::complex>** MPll::allcoeffprime

vector that contains the primed Wilson coeffients

Definition at line 96 of file MPll.h.

double MPll::avaDelta0

private

Gsl integral variable

Definition at line 508 of file MPll.h.

double MPll::avaDelta2

private

Gsl integral variable

Definition at line 509 of file MPll.h.

double MPll::avaSigma0

private

Gsl integral variable

Definition at line 502 of file MPll.h.

double MPll::avaSigma2

private

Gsl integral variable

Definition at line 503 of file MPll.h.

double MPll::b

BF of the decay V -> final states

Definition at line 79 of file MPll.h.

double MPll::beta_cache

private

Cache variable

Definition at line 405 of file MPll.h.

unsigned int MPll::beta_updated

private

Cache variable

Definition at line 404 of file MPll.h.

gslpp::complex MPll::C_1

Wilson coeffients \(C_1\)

Definition at line 98 of file MPll.h.

gslpp::complex MPll::C_10

Wilson coeffients \(C_{10}\)

Definition at line 106 of file MPll.h.

gslpp::complex MPll::C_10_cache

private

Cache variable

Definition at line 454 of file MPll.h.

unsigned int MPll::C_10_updated

private

Cache variable

Definition at line 453 of file MPll.h.

gslpp::complex MPll::C_10p

Wilson coeffients \(C_{10}'\)

Definition at line 112 of file MPll.h.

gslpp::complex MPll::C_10p_cache

private

Cache variable

Definition at line 463 of file MPll.h.

unsigned int MPll::C_10p_updated

private

Cache variable

Definition at line 462 of file MPll.h.

gslpp::complex MPll::C_1_cache

private

Cache variable

Definition at line 430 of file MPll.h.

unsigned int MPll::C_1_updated

private

Cache variable

Definition at line 429 of file MPll.h.

gslpp::complex MPll::C_2

Wilson coeffients \(C_2\)

Definition at line 99 of file MPll.h.

gslpp::complex MPll::C_2_cache

private

Cache variable

Definition at line 433 of file MPll.h.

unsigned int MPll::C_2_updated

private

Cache variable

Definition at line 432 of file MPll.h.

gslpp::complex MPll::C_3

Wilson coeffients \(C_3\)

Definition at line 100 of file MPll.h.

gslpp::complex MPll::C_3_cache

private

Cache variable

Definition at line 436 of file MPll.h.

unsigned int MPll::C_3_updated

private

Cache variable

Definition at line 435 of file MPll.h.

gslpp::complex MPll::C_4

Wilson coeffients \(C_4\)

Definition at line 101 of file MPll.h.

gslpp::complex MPll::C_4_cache

private

Cache variable

Definition at line 439 of file MPll.h.

unsigned int MPll::C_4_updated

private

Cache variable

Definition at line 438 of file MPll.h.

gslpp::complex MPll::C_5

Wilson coeffients \(C_5\)

Definition at line 102 of file MPll.h.

gslpp::complex MPll::C_5_cache

private

Cache variable

Definition at line 442 of file MPll.h.

unsigned int MPll::C_5_updated

private

Cache variable

Definition at line 441 of file MPll.h.

gslpp::complex MPll::C_6

Wilson coeffients \(C_6\)

Definition at line 103 of file MPll.h.

gslpp::complex MPll::C_6_cache

private

Cache variable

Definition at line 445 of file MPll.h.

unsigned int MPll::C_6_updated

private

Cache variable

Definition at line 444 of file MPll.h.

gslpp::complex MPll::C_7

Wilson coeffients \(C_7\)

Definition at line 104 of file MPll.h.

gslpp::complex MPll::C_7_cache

private

Cache variable

Definition at line 448 of file MPll.h.

unsigned int MPll::C_7_updated

private

Cache variable

Definition at line 447 of file MPll.h.

gslpp::complex MPll::C_7p

Wilson coeffients \(C_7'\)

Definition at line 110 of file MPll.h.

gslpp::complex MPll::C_7p_cache

private

Cache variable

Definition at line 457 of file MPll.h.

unsigned int MPll::C_7p_updated

private

Cache variable

Definition at line 456 of file MPll.h.

gslpp::complex MPll::C_9

Wilson coeffients \(C_9\)

Definition at line 105 of file MPll.h.

gslpp::complex MPll::C_9_cache

private

Cache variable

Definition at line 451 of file MPll.h.

unsigned int MPll::C_9_updated

private

Cache variable

Definition at line 450 of file MPll.h.

gslpp::complex MPll::C_9p

Wilson coeffients \(C_9'\)

Definition at line 111 of file MPll.h.

gslpp::complex MPll::C_9p_cache

private

Cache variable

Definition at line 460 of file MPll.h.

unsigned int MPll::C_9p_updated

private

Cache variable

Definition at line 459 of file MPll.h.

gslpp::complex MPll::C_P

Wilson coeffients \(C_P\)

Definition at line 108 of file MPll.h.

gslpp::complex MPll::C_P_cache

private

Cache variable

Definition at line 469 of file MPll.h.

unsigned int MPll::C_P_updated

private

Cache variable

Definition at line 468 of file MPll.h.

gslpp::complex MPll::C_Pp

Wilson coeffients \(C_P'\)

Definition at line 114 of file MPll.h.

gslpp::complex MPll::C_Pp_cache

private

Cache variable

Definition at line 475 of file MPll.h.

unsigned int MPll::C_Pp_updated

private

Cache variable

Definition at line 474 of file MPll.h.

gslpp::complex MPll::C_S

Wilson coeffients \(C_S\)

Definition at line 107 of file MPll.h.

gslpp::complex MPll::C_S_cache

private

Cache variable

Definition at line 466 of file MPll.h.

unsigned int MPll::C_S_updated

private

Cache variable

Definition at line 465 of file MPll.h.

gslpp::complex MPll::C_Sp

Wilson coeffients \(C_S'\)

Definition at line 113 of file MPll.h.

gslpp::complex MPll::C_Sp_cache

private

Cache variable

Definition at line 472 of file MPll.h.

unsigned int MPll::C_Sp_updated

private

Cache variable

Definition at line 471 of file MPll.h.

std::map<std::pair<double, double>, double > MPll::cacheDelta0

private

Gsl integral variable

Definition at line 529 of file MPll.h.

std::map<std::pair<double, double>, double > MPll::cacheDelta2

private

Gsl integral variable

Definition at line 530 of file MPll.h.

std::map<std::pair<double, double>, double > MPll::cacheSigma0

private

Gsl integral variable

Definition at line 526 of file MPll.h.

std::map<std::pair<double, double>, double > MPll::cacheSigma2

private

Gsl integral variable

Definition at line 527 of file MPll.h.

std::map<std::pair<double, double>, unsigned int > MPll::delta0Cached

private

Cache variable

Definition at line 499 of file MPll.h.

std::map<std::pair<double, double>, unsigned int > MPll::delta2Cached

private

Cache variable

Definition at line 500 of file MPll.h.

double MPll::errDelta0

private

Gsl integral variable

Definition at line 511 of file MPll.h.

double MPll::errDelta2

private

Gsl integral variable

Definition at line 512 of file MPll.h.

double MPll::errSigma0

private

Gsl integral variable

Definition at line 505 of file MPll.h.

double MPll::errSigma2

private

Gsl integral variable

Definition at line 506 of file MPll.h.

double MPll::f0_cache

private

Cache variable

Definition at line 399 of file MPll.h.

unsigned int MPll::f0_updated

private

Cache variable

Definition at line 398 of file MPll.h.

unsigned int MPll::F_updated

private

Cache variable

Definition at line 410 of file MPll.h.

gsl_function MPll::FD0

private

Gsl integral variable

Definition at line 517 of file MPll.h.

gsl_function MPll::FD2

private

Gsl integral variable

Definition at line 518 of file MPll.h.

gslpp::vector<double> MPll::fplus_cache

private

Cache variable

Definition at line 393 of file MPll.h.

unsigned int MPll::fplus_updated

private

Cache variable

Definition at line 392 of file MPll.h.

gsl_function MPll::FS0

private

Gsl integral variable

Definition at line 514 of file MPll.h.

gsl_function MPll::FS2

private

Gsl integral variable

Definition at line 515 of file MPll.h.

gslpp::vector<double> MPll::fT_cache

private

Cache variable

Definition at line 396 of file MPll.h.

unsigned int MPll::fT_updated

private

Cache variable

Definition at line 395 of file MPll.h.

double MPll::GF

Fermi constant

Definition at line 67 of file MPll.h.

gslpp::complex MPll::h_0

parameter that contains the contribution from the hadronic hamiltonian

Definition at line 80 of file MPll.h.

gslpp::complex MPll::h_0_1

parameter that contains the contribution from the hadronic hamiltonian

Definition at line 81 of file MPll.h.

unsigned int MPll::H_A0updated

private

Cache variable

Definition at line 484 of file MPll.h.

gslpp::vector<double> MPll::H_P_cache

private

Cache variable

Definition at line 490 of file MPll.h.

unsigned int MPll::H_P_updated

private

Cache variable

Definition at line 489 of file MPll.h.

gslpp::vector<double> MPll::H_Scache

private

Cache variable

Definition at line 487 of file MPll.h.

unsigned int MPll::H_Supdated

private

Cache variable

Definition at line 486 of file MPll.h.

gslpp::vector<double> MPll::H_V0cache

private

Cache variable

Definition at line 481 of file MPll.h.

gslpp::complex MPll::H_V0Ccache[2]

private

Cache variable

Definition at line 482 of file MPll.h.

unsigned int MPll::H_V0updated

private

Cache variable

Definition at line 480 of file MPll.h.

unsigned int MPll::I0_updated

private

Cache variable

Definition at line 492 of file MPll.h.

unsigned int MPll::I2_updated

private

Cache variable

Definition at line 493 of file MPll.h.

unsigned int MPll::I8_updated

private

Cache variable

Definition at line 494 of file MPll.h.

gslpp::vector<double> MPll::k2_cache

private

Cache variable

Definition at line 402 of file MPll.h.

unsigned int MPll::k2_updated

private

Cache variable

Definition at line 401 of file MPll.h.

double MPll::lambda_cache

private

Cache variable

Definition at line 408 of file MPll.h.

gslpp::complex MPll::lambda_t

Vckm factor

Definition at line 78 of file MPll.h.

unsigned int MPll::lambda_updated

private

Cache variable

Definition at line 407 of file MPll.h.

StandardModel::lepton MPll::lep

private

Final leptons type

Definition at line 388 of file MPll.h.

double MPll::m_fit2_f0

LCSR fit parameter

Definition at line 92 of file MPll.h.

double MPll::m_fit2_fplus

LCSR fit parameter

Definition at line 87 of file MPll.h.

double MPll::m_fit2_fT

LCSR fit parameter

Definition at line 90 of file MPll.h.

double MPll::Mb

b quark mass

Definition at line 72 of file MPll.h.

double MPll::Mc

c quark mass

Definition at line 74 of file MPll.h.

StandardModel::meson MPll::meson

private

Initial meson type

Definition at line 389 of file MPll.h.

double MPll::Mlep

muon mass

Definition at line 69 of file MPll.h.

double MPll::MM

initial meson mass

Definition at line 70 of file MPll.h.

double MPll::MP

final pseudoscalar meson mass

Definition at line 71 of file MPll.h.

double MPll::Ms

s quark mass

Definition at line 75 of file MPll.h.

double MPll::mu_b

b mass scale

Definition at line 73 of file MPll.h.

double MPll::MW

W boson mass

Definition at line 77 of file MPll.h.

const StandardModel& MPll::mySM

private

Model type

Definition at line 387 of file MPll.h.

gslpp::vector<double> MPll::N_cache

private

Cache variable

Definition at line 426 of file MPll.h.

unsigned int MPll::N_updated

private

Cache variable

Definition at line 425 of file MPll.h.

gslpp::complex MPll::Nc_cache

private

Cache variable

Definition at line 427 of file MPll.h.

StandardModel::meson MPll::pseudoscalar

private

Final pseudoscalar meson type

Definition at line 390 of file MPll.h.

double MPll::q2

\(q^2\) of the decay

Definition at line 82 of file MPll.h.

double MPll::r_1_fplus

LCSR fit parameter

Definition at line 85 of file MPll.h.

double MPll::r_1_fT

LCSR fit parameter

Definition at line 88 of file MPll.h.

double MPll::r_2_f0

LCSR fit parameter

Definition at line 91 of file MPll.h.

double MPll::r_2_fplus

LCSR fit parameter

Definition at line 86 of file MPll.h.

double MPll::r_2_fT

LCSR fit parameter

Definition at line 89 of file MPll.h.

std::map<std::pair<double, double>, unsigned int > MPll::sigma0Cached

private

Cache variable

Definition at line 496 of file MPll.h.

std::map<std::pair<double, double>, unsigned int > MPll::sigma2Cached

private

Cache variable

Definition at line 497 of file MPll.h.

gslpp::vector<double> MPll::SL_cache

private

Cache variable

Definition at line 421 of file MPll.h.

unsigned int MPll::SL_updated

private

Cache variable

Definition at line 420 of file MPll.h.

unsigned int MPll::SR_updated

private

Cache variable

Definition at line 423 of file MPll.h.

unsigned int MPll::TL_updated

private

Cache variable

Definition at line 416 of file MPll.h.

unsigned int MPll::TR_updated

private

Cache variable

Definition at line 418 of file MPll.h.

unsigned int MPll::VL_updated

private

Cache variable

Definition at line 412 of file MPll.h.

unsigned int MPll::VR_updated

private

Cache variable

Definition at line 414 of file MPll.h.

gsl_integration_workspace* MPll::w_delta0

private

Gsl integral variable

Definition at line 523 of file MPll.h.

gsl_integration_workspace* MPll::w_delta2

private

Gsl integral variable

Definition at line 524 of file MPll.h.

gsl_integration_workspace* MPll::w_sigma0

private

Gsl integral variable

Definition at line 520 of file MPll.h.

gsl_integration_workspace* MPll::w_sigma2

private

Gsl integral variable

Definition at line 521 of file MPll.h.

double MPll::width

initial meson width

Definition at line 76 of file MPll.h.

gslpp::vector<double> MPll::Ycache

private

Cache variable

Definition at line 478 of file MPll.h.

unsigned int MPll::Yupdated

private

Cache variable

Definition at line 477 of file MPll.h.

---

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

• MPll.h
• MPll.cpp