MVll Class Reference

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

#include <MVll.h>

Collaboration diagram for MVll:

Detailed Description

A class for the \(M \to V 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 V 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 V(), A_0(), A_1(), A_2(), T_1() and T_2() using the fit function FF_fit() from [42]. The form factor are consequentely translated in the helicity basis through the functions V_0t(), V_p(), V_m(), T_0t(), T_p(), T_m() and S_L() from [93]. The QCDF corrections to Wilson coefficient \(C_9\) are computed according to [40] . The basic elements are build in the functions Tperpplus(), Tparplus(), Tparminus(), Cperp() and Cpar(); these corrections have to be integrated to be computed, so the final correction is either obtaind through direct integration in the functions DeltaC9_p(), DeltaC9_m() and DeltaC9_0(), or obtained through fitting in the functions fDeltaC9_p(), fDeltaC9_m() and fDeltaC9_0(). Form factors, Wilson coefficients and parameters are combined together in the functions H_V_0(), H_V_p(), H_V_m(), H_A_0(), H_A_p(), H_A_m(), H_S() and H_P() [93] in order to build the helicity aplitudes, which are consequentely combined to create the angular coefficients in the function I_1c(), I_1s(), I_2c(), I_2s(), I_3(), I_4(), I_5(), I_6c(), I_6s(), I_7(), I_8(), I_9(). Those coefficients are used to create the CP averaged coefficients in the functions getSigma1c(), getSigma1s(), getSigma2c(), getSigma2s(), getSigma3(), getSigma4(), getSigma5(), getSigma6c(), getSigma6s(), getSigma7(), getSigma8(), getSigma9(), and the CP asymmetric coefficients in the function Delta(). The CP averaged and asymmetric coefficients are integrated over the \(q^2\) bin in the functions integrateSigma() and integrateDelta(), in order to be further used to build the observables.

Definition at line 87 of file MVll.h.

Public Member Functions
double	beta (double q2)
	The factor \( \beta \) used in the angular coefficients \(I_i\). More...
double	getgtilde_1_im (double q2)
	The immaginary part of \( \tilde{g}^1 \). More...
double	getgtilde_1_re (double q2)
	The real part of \( \tilde{g}^1 \). More...
double	getgtilde_2_im (double q2)
	The immaginary part of \( \tilde{g}^2 \). More...
double	getgtilde_2_re (double q2)
	The real part of \( \tilde{g}^2 \). More...
double	getgtilde_3_im (double q2)
	The imaginary part of \( \tilde{g}^3 \). More...
double	getgtilde_3_re (double q2)
	The real part of \( \tilde{g}^3 \). More...
double	geth_0_im (double q2)
	The imaginary part of \( h_0 \). More...
double	geth_0_re (double q2)
	The real part of \( h_0 \). More...
double	geth_m_im (double q2)
	The imaginary part of \( h_- \). More...
double	geth_m_re (double q2)
	The real part of \( h_- \). More...
double	geth_p_im (double q2)
	The imaginary part of \( h_+ \). More...
double	geth_p_re (double q2)
	The real part of \( h_+ \). More...
double	gethm0_h00_abs ()
	The absolute value of the ratio \( h_-^{(0)}/h_0^{(0)} \). More...
double	gethp0_hm0_abs ()
	The absolute value of the ratio \( h_+^{(0)}/h_0^{(0)} \). More...
double	getMlep ()
	The mass of the lepton l. More...
double	getS (double q2)
	The form factor \( S \) . More...
double	getSigma (int i, double q_2)
	The value of \( \Sigma_{i} \) from \(q_{min}\) to \(q_{max}\). More...
double	getT0 (double q2)
	The form factor \( T_0 \) . More...
double	getTm (double q2)
	The form factor \( T_- \) . More...
double	getTp (double q2)
	The form factor \( T_+ \) . More...
double	getV0 (double q2)
	The form factor \( V_0 \) . More...
double	getVm (double q2)
	The form factor \( V_- \) . More...
double	getVp (double q2)
	The form factor \( V_+ \) . More...
double	getwidth ()
	The width of the meson M. More...
gslpp::complex	H_A_0 (double q2)
	The helicity amplitude \( H_A^0 \) . More...
gslpp::complex	H_A_m (double q2)
	The helicity amplitude \( H_A^- \) . More...
gslpp::complex	H_A_p (double q2)
	The helicity amplitude \( H_A^+ \) . More...
gslpp::complex	H_P (double q2)
	The helicity amplitude \( H_P \) . More...
gslpp::complex	H_S (double q2)
	The helicity amplitude \( H_S \) . More...
gslpp::complex	H_V_0 (double q2)
	The helicity amplitude \( H_V^0 \) . More...
gslpp::complex	H_V_m (double q2)
	The helicity amplitude \( H_V^- \) . More...
gslpp::complex	H_V_p (double q2)
	The helicity amplitude \( H_V^+ \) . 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...
	MVll (const StandardModel &SM_i, StandardModel::meson meson_i, StandardModel::meson vector_i, StandardModel::lepton lep_i)
	Constructor. More...
virtual	~MVll ()
	Destructor. More...

Private Member Functions
double	A_0 (double q2)
	The transverse form factor \( A_0 \). More...
double	A_1 (double q2)
	The transverse form factor \( A_1 \). More...
double	A_2 (double q2)
	The transverse form factor \( A_2 \). More...
double	beta2 (double q2)
	The factor \( \beta^2 \) used in the angular coefficients \(I_i\). More...
void	checkCache ()
	The caching method for MVll. More...
gslpp::complex	Cpar (double q2)
	The correction \( C_{\parallel} \) from [40] . More...
gslpp::complex	Cperp (double q2)
	The correction \( C_{\perp} \) from [40] . More...
double	Delta (int i, double q2)
	The CP asymmetry \( \Delta_{i} \) . More...
gslpp::complex	DeltaC9_0 (double q2)
	The total QCDF correction \( \Delta C_9^0 \) computed integrating over \( u \). More...
gslpp::complex	DeltaC9_m (double q2)
	The total QCDF correction \( \Delta C_9^- \) computed integrating over \( u \). More...
gslpp::complex	DeltaC9_p (double q2)
	The total QCDF correction \( \Delta C_9^+ \) computed integrating over \( u \). More...
gslpp::complex	deltaTpar (double q2)
	The total correction \( \Delta \mathcal{T}^{\parallel} \) from [40] . More...
gslpp::complex	deltaTperp (double q2)
	The total correction \( \Delta \mathcal{T}^{\perp} \) from [40] . More...
double	F (double q2, double b_i)
	The factor \( F \) used in the angular coefficients I_i. More...
gslpp::complex	F19 (double q2)
	The correction \( F_{19} \) from [15]. More...
gslpp::complex	F27 (double q2)
	The correction \( F_{27} \) from [15]. More...
gslpp::complex	F29 (double q2)
	The correction \( F_{29} \) from [15]. More...
gslpp::complex	F87 (double q2)
	The correction \( F_{87} \) from [15]. More...
double	F89 (double q2)
	The correction \( F_{89} \) from [15]. More...
gslpp::complex	fDeltaC9_0 (double q2)
	The total QCDF correction \( \Delta C_9^0 \) computed fitting over \( u \). More...
gslpp::complex	fDeltaC9_m (double q2)
	The total QCDF correction \( \Delta C_9^+ \) computed fitting over \( u \). More...
gslpp::complex	fDeltaC9_p (double q2)
	The total QCDF correction \( \Delta C_9^0 \) computed fitting over \( u \). More...
double	FF_fit (double q2, double a_0, double a_1, double a_2, double MR2)
	The fit function from [42], \( FF^{\rm fit} \). More...
void	fit_DeltaC9_0_ee ()
	The fitting routine for the QCDF correction \( \Delta C_9^0 \) in the electron channel. More...
void	fit_DeltaC9_0_mumu ()
	The fitting routine for the QCDF correction \( \Delta C_9^0 \) in the muon channel. More...
void	fit_DeltaC9_m_ee ()
	The fitting routine for the QCDF correction \( \Delta C_9^- \) in the electron channel. More...
void	fit_DeltaC9_m_mumu ()
	The fitting routine for the QCDF correction \( \Delta C_9^- \) in the muon channel. More...
void	fit_DeltaC9_p_ee ()
	The fitting routine for the QCDF correction \( \Delta C_9^+ \) in the electron channel. More...
void	fit_DeltaC9_p_mumu ()
	The fitting routine for the QCDF correction \( \Delta C_9^+ \) in the muon channel. More...
double	getDelta0 (double q2)
	The CP asymmetry \( \Delta_{1s} \) . More...
double	getDelta1 (double q2)
	The CP asymmetry \( \Delta_{1c} \) . More...
double	getDelta11 (double q2)
	The CP asymmetry \( \Delta_{9} \) . More...
double	getDelta2 (double q2)
	The CP asymmetry \( \Delta_{2s} \) . More...
double	getDelta3 (double q2)
	The CP asymmetry \( \Delta_{2c} \) . More...
double	getDelta7 (double q2)
	The CP asymmetry \( \Delta_{6s} \) . More...
double	getSigma1c (double q2)
	The CP average \( \Sigma_{1s} \) . More...
double	getSigma1s (double q2)
	The CP average \( \Sigma_{1c} \) . More...
double	getSigma2c (double q2)
	The CP average \( \Sigma_{2s} \) . More...
double	getSigma2s (double q2)
	The CP average \( \Sigma_{2c} \) . More...
double	getSigma3 (double q2)
	The CP average \( \Sigma_{3} \) . More...
double	getSigma4 (double q2)
	The CP average \( \Sigma_{4} \) . More...
double	getSigma5 (double q2)
	The CP average \( \Sigma_{5} \) . More...
double	getSigma6c (double q2)
	The CP average \( \Sigma_{7} \) . More...
double	getSigma6s (double q2)
	The CP average \( \Sigma_{6s} \) . More...
double	getSigma7 (double q2)
	The CP average \( \Sigma_{8} \) . More...
double	getSigma8 (double q2)
	The CP average \( \Sigma_{9} \) . More...
double	getSigma9 (double q2)
	The CP average \( \Sigma_{9} \) . More...
gslpp::complex	H_0 (double q2)
	The \( h(q^2,0) \) function involved into \( C_9^{eff}\). More...
gslpp::complex	H_b (double q2)
	The \( h(q^2,m_b) \) function involved into \( C_9^{eff}\). More...
gslpp::complex	H_c (double q2, double mu)
	The \( h(q^2,m_c) \) function involved into \( C_9^{eff}\). More...
gslpp::complex	I1 (double u, double q2)
	The \( I_1 \) function from [40] . More...
double	I_1c (double q2)
	The angular coefficient \( I_{1c} \) . More...
double	I_1s (double q2)
	The angular coefficient \( I_{1s} \) . More...
double	I_2c (double q2)
	The angular coefficient \( I_{2c} \) . More...
double	I_2s (double q2)
	The angular coefficient \( I_{2s} \) . More...
double	I_3 (double q2)
	The angular coefficient \( I_3 \) . More...
double	I_4 (double q2)
	The angular coefficient \( I_4 \) . More...
double	I_5 (double q2)
	The angular coefficient \( I_5 \) . More...
double	I_6c (double q2)
	The angular coefficient \( I_{6c} \) . More...
double	I_6s (double q2)
	The angular coefficient \( I_{6s} \) . More...
double	I_7 (double q2)
	The angular coefficient \( I_7 \) . More...
double	I_8 (double q2)
	The angular coefficient \( I_8 \) . More...
double	I_9 (double q2)
	The angular coefficient \( I_9 \) . More...
double	imDC9fit (double *x, double *p)
	The fit function for the imaginary part of the QCDF correction \( \Delta C_9^{\lambda} \). More...
double	Integrand_ImTpar_pm (double up)
	The sum of Integrand_ImTparplus() and Integrand_ImTparminus(). More...
double	Integrand_ImTparminus (double up)
	The imaginary part of the integral involving \( T^{\parallel}_- \) at fixed \( q^2 \), according to [40] . More...
double	Integrand_ImTparplus (double up)
	The imaginary part of the integral involving \( T^{\parallel}_+ \) at fixed \( q^2 \), according to [40] . More...
double	Integrand_ImTperpplus (double up)
	The imaginary part of the integral involving \( T^{\perp}_+ \) at fixed \( q^2 \), according to [40] . More...
double	Integrand_ReTpar_pm (double up)
	The sum of Integrand_ReTparplus() and Integrand_ReTparminus(). More...
double	Integrand_ReTparminus (double up)
	The real part of the integral involving \( T^{\parallel}_- \) at fixed \( q^2 \), according to [40] . More...
double	Integrand_ReTparplus (double up)
	The real part of the integral involving \( T^{\parallel}_+ \) at fixed \( q^2 \), according to [40] . More...
double	Integrand_ReTperpplus (double up)
	The real part of the integral involving \( T^{\perp}_+ \) at fixed \( q^2 \), according to [40] . 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	reDC9fit (double *x, double *p)
	The fit function for the real part of the QCDF correction \( \Delta C_9^{\lambda} \). More...
double	S_L (double q2)
	The helicity form factor \( S_L \). More...
double	T_0t (double q2)
	The helicity form factor \( \tilde{T}_0 \). More...
double	T_1 (double q2)
	The transverse form factor \( T_1 \). More...
double	T_2 (double q2)
	The transverse form factor \( T_2 \). More...
double	T_m (double q2)
	The helicity form factor \( T_- \). More...
double	T_p (double q2)
	The helicity form factor \( T_+ \). More...
gslpp::complex	Tparminus (double u, double q2)
	The \( T^{\parallel}_- \) function from [40] . More...
gslpp::complex	Tparplus (double u, double q2)
	The \( T^{\parallel}_+ \) function from [40] . More...
gslpp::complex	Tperpplus (double u, double q2)
	The \( T^{\perp}_+ \) function from [40] . More...
void	updateParameters ()
	The update parameter method for MVll. More...
double	V (double q2)
	The transverse form factor \( V \). More...
double	V_0t (double q2)
	The helicity form factor \( \tilde{V}_0 \). More...
double	V_m (double q2)
	The helicity form factor \( V_- \). More...
double	V_p (double q2)
	The helicity form factor \( V_+ \). More...
gslpp::complex	Y (double q2)
	The \( Y(q^2) \) function involved into \( C_9^{eff}\). More...
double	z (double q2)
	The lattice parameter \( z \) from arXiv:1310.3722v3. More...

Private Attributes
gslpp::vector< double >	A0_cache
unsigned int	A0_updated
gslpp::vector< double >	A1_cache
unsigned int	A1_updated
double	a_0A0
double	a_0A1
double	a_0A12
double	a_0A12_LCSR
double	a_0T1
double	a_0T2
double	a_0T23
double	a_0T2_LCSR
double	a_0V
double	a_1A0
double	a_1A1
double	a_1A12
double	a_1T1
double	a_1T2
double	a_1T23
double	a_1V
double	a_2A0
double	a_2A1
double	a_2A12
double	a_2T1
double	a_2T2
double	a_2T23
double	a_2V
double	ale
gslpp::vector< gslpp::complex > **	allcoeff
gslpp::vector< gslpp::complex > **	allcoeffh
gslpp::vector< gslpp::complex > **	allcoeffprime
gslpp::complex	arg1
double	avaDelta
double	avaDTPPR
double	avaSigma
double	b
gslpp::complex	B00
gslpp::complex	B01
double	beta_cache
unsigned int	beta_updated
double	C2_inv
gslpp::complex	C_1
gslpp::complex	C_10
gslpp::complex	C_10_cache
unsigned int	C_10_updated
gslpp::complex	C_10p
gslpp::complex	C_10p_cache
unsigned int	C_10p_updated
gslpp::complex	C_1_cache
unsigned int	C_1_updated
gslpp::complex	C_1L_bar
gslpp::complex	C_1Lh_bar
gslpp::complex	C_2
gslpp::complex	C_2_cache
unsigned int	C_2_updated
gslpp::complex	C_2L_bar
gslpp::complex	C_2Lh_bar
gslpp::complex	C_2Lh_cache
unsigned int	C_2Lh_updated
gslpp::complex	C_3
gslpp::complex	C_3_cache
unsigned int	C_3_updated
gslpp::complex	C_4
gslpp::complex	C_4_cache
unsigned int	C_4_updated
gslpp::complex	C_5
gslpp::complex	C_5_cache
unsigned int	C_5_updated
gslpp::complex	C_6
gslpp::complex	C_6_cache
unsigned int	C_6_updated
gslpp::complex	C_7
gslpp::complex	C_7_cache
unsigned int	C_7_updated
gslpp::complex	C_7p
gslpp::complex	C_7p_cache
unsigned int	C_7p_updated
gslpp::complex	C_8L
gslpp::complex	C_8Lh
gslpp::complex	C_8Lh_cache
unsigned int	C_8Lh_updated
gslpp::complex	C_9
gslpp::complex	C_9_cache
unsigned int	C_9_updated
gslpp::complex	C_9p
gslpp::complex	C_9p_cache
unsigned int	C_9p_updated
gslpp::complex	C_P
gslpp::complex	C_P_cache
unsigned int	C_P_updated
gslpp::complex	C_Pp
gslpp::complex	C_Pp_cache
unsigned int	C_Pp_updated
gslpp::complex	C_S
gslpp::complex	C_S_cache
unsigned int	C_S_updated
gslpp::complex	C_Sp
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 >	cacheDelta1
std::map< std::pair< double, double >, double >	cacheDelta11
std::map< std::pair< double, double >, double >	cacheDelta2
std::map< std::pair< double, double >, double >	cacheDelta3
std::map< std::pair< double, double >, double >	cacheDelta7
std::map< double, gslpp::complex >	cacheDeltaTparm
std::map< double, gslpp::complex >	cacheDeltaTparp
std::map< double, gslpp::complex >	cacheDeltaTperp
std::map< std::pair< double, double >, gslpp::complex >	cacheI1
std::map< std::pair< double, double >, double >	cacheSigma0
std::map< std::pair< double, double >, double >	cacheSigma1
std::map< std::pair< double, double >, double >	cacheSigma10
std::map< std::pair< double, double >, double >	cacheSigma11
std::map< std::pair< double, double >, double >	cacheSigma2
std::map< std::pair< double, double >, double >	cacheSigma3
std::map< std::pair< double, double >, double >	cacheSigma4
std::map< std::pair< double, double >, double >	cacheSigma5
std::map< std::pair< double, double >, double >	cacheSigma6
std::map< std::pair< double, double >, double >	cacheSigma7
std::map< std::pair< double, double >, double >	cacheSigma9
double	CF
std::map< std::pair< double, double >, unsigned int >	delta0Cached
std::map< std::pair< double, double >, unsigned int >	delta11Cached
std::map< std::pair< double, double >, unsigned int >	delta1Cached
std::map< std::pair< double, double >, unsigned int >	delta2Cached
std::map< std::pair< double, double >, unsigned int >	delta3Cached
std::map< std::pair< double, double >, unsigned int >	delta7Cached
double	deltaT_0
double	deltaT_1par
double	deltaT_1perp
std::map< double, unsigned int >	deltaTparmCached
unsigned int	deltaTparmupdated
std::map< double, unsigned int >	deltaTparpCached
unsigned int	deltaTparpupdated
std::map< double, unsigned int >	deltaTperpCached
unsigned int	deltaTperpupdated
gsl_function	DTPPR
double	Ee
double	errDelta
double	errDTPPR
double	errSigma
gslpp::complex	F27_0
gslpp::complex	F27_1
gslpp::complex	F27_2
gslpp::complex	F27_3
gslpp::complex	F27_L1_1
gslpp::complex	F27_L1_2
gslpp::complex	F27_L1_3
gslpp::complex	F29_0
gslpp::complex	F29_1
gslpp::complex	F29_2
gslpp::complex	F29_3
gslpp::complex	F29_L1
gslpp::complex	F29_L1_1
gslpp::complex	F29_L1_2
gslpp::complex	F29_L1_3
gslpp::complex	F87_0
gslpp::complex	F87_1
gslpp::complex	F87_2
gslpp::complex	F87_3
double	F89_0
double	F89_1
double	F89_2
double	F89_3
unsigned int	F_updated
gsl_function	FD
double	fourMb2
double	fourMc2
double	fourMM2
double	fourMV
double	fournineth
gsl_function	FS
double	GF
TGraph	gr1
TGraph	gr2
double	gtilde_1_pre
double	gtilde_2_pre
double	gtilde_3_pre
unsigned int	h0_updated
gslpp::complex	h0Ccache [3]
unsigned int	h1_updated
gslpp::complex	h1Ccache [3]
unsigned int	h2_updated
gslpp::complex	h2Ccache [3]
gslpp::complex	h_0 [3]
gslpp::complex	H_0_pre
gslpp::complex	H_0_WC
gslpp::complex	h_1 [3]
gslpp::complex	h_2 [3]
unsigned int	H_A0updated
unsigned int	H_A1updated
unsigned int	H_A2updated
gslpp::complex	H_b_WC
gslpp::complex	H_c_WC
gslpp::vector< double >	H_Pcache
unsigned int	H_Pupdated
gslpp::vector< double >	H_Scache
unsigned int	H_Supdated
gslpp::vector< double >	H_V0cache
unsigned int	H_V0updated
gslpp::vector< double >	H_V1cache
unsigned int	H_V1updated
gslpp::vector< double >	H_V2cache
unsigned int	H_V2updated
double	half
unsigned int	I0_updated
unsigned int	I10_updated
unsigned int	I11_updated
unsigned int	I1_updated
std::map< std::pair< double, double >, unsigned int >	I1Cached
unsigned int	I2_updated
unsigned int	I3_updated
unsigned int	I4_updated
unsigned int	I5_updated
unsigned int	I6_updated
unsigned int	I7_updated
unsigned int	I8_updated
unsigned int	I9_updated
gslpp::complex	ihalfMPI
std::vector< double >	ImDeltaC9_0_ee
std::vector< double >	ImDeltaC9_0_mumu
std::vector< double >	ImDeltaC9_m_ee
std::vector< double >	ImDeltaC9_m_mumu
std::vector< double >	ImDeltaC9_p_ee
std::vector< double >	ImDeltaC9_p_mumu
TF1	imffit
TFitResultPtr	imfres_0_ee
TFitResultPtr	imfres_0_mumu
TFitResultPtr	imfres_m_ee
TFitResultPtr	imfres_m_mumu
TFitResultPtr	imfres_p_ee
TFitResultPtr	imfres_p_mumu
gslpp::vector< double >	k2_cache
unsigned int	k2_updated
gslpp::complex	L1xm
gslpp::complex	L1xp
gslpp::complex	L1ym
gslpp::complex	L1yp
gslpp::complex	lambda_t
unsigned int	lambda_updated
StandardModel::lepton	lep
double	logMb
double	logMc
double	m4downcharge
double	M_PI2osix
double	Mb
double	Mb2
unsigned int	Mb_Ms_updated
double	MboMW
double	Mc
double	Mc2
StandardModel::meson	meson
double	Mlep
double	Mlep2
double	MM
double	MM2
double	MM2mMV2
double	MM4
double	MM_MMpMV
double	MMmMV
double	MMmMV2
double	MMMV
double	MMpMV
double	MMpMV2
double	MRA0_2
double	MRA12_2
double	MRA1_2
double	MRT1_2
double	MRT23_2
double	MRT2_2
double	MRV_2
double	Ms
double	MsoMb
double	mu_b
double	mu_b2
double	mu_h
double	MV
double	MV2
double	MV4
double	MW
std::vector< double >	myq2
const StandardModel &	mySM
gslpp::vector< double >	N_cache
unsigned int	N_updated
gslpp::complex	Nc_cache
double	ninetysixM_PI3MM3
double	NN
gsl_error_handler_t *	old_handler
double	onepMMoMV
std::vector< double >	ReDeltaC9_0_ee
std::vector< double >	ReDeltaC9_0_mumu
std::vector< double >	ReDeltaC9_m_ee
std::vector< double >	ReDeltaC9_m_mumu
std::vector< double >	ReDeltaC9_p_ee
std::vector< double >	ReDeltaC9_p_mumu
TF1	reffit
TFitResultPtr	refres_0_ee
TFitResultPtr	refres_0_mumu
TFitResultPtr	refres_m_ee
TFitResultPtr	refres_m_mumu
TFitResultPtr	refres_p_ee
TFitResultPtr	refres_p_mumu
double	S_L_pre
std::map< std::pair< double, double >, unsigned int >	sigma0Cached
std::map< std::pair< double, double >, unsigned int >	sigma10Cached
std::map< std::pair< double, double >, unsigned int >	sigma11Cached
std::map< std::pair< double, double >, unsigned int >	sigma1Cached
std::map< std::pair< double, double >, unsigned int >	sigma2Cached
std::map< std::pair< double, double >, unsigned int >	sigma3Cached
std::map< std::pair< double, double >, unsigned int >	sigma4Cached
std::map< std::pair< double, double >, unsigned int >	sigma5Cached
std::map< std::pair< double, double >, unsigned int >	sigma6Cached
std::map< std::pair< double, double >, unsigned int >	sigma7Cached
std::map< std::pair< double, double >, unsigned int >	sigma9Cached
double	sixMMoMb
double	sixteenM_PI2
double	sixteenM_PI2MM2
gslpp::vector< double >	SL_cache
unsigned int	SL_updated
unsigned int	SR_updated
gslpp::vector< double >	T1_cache
unsigned int	T1_updated
gslpp::vector< double >	T2_cache
unsigned int	T2_updated
double	t_0
gslpp::vector< double >	T_cache
double	t_m
double	t_p
unsigned int	T_updated
double	threeGegen0
double	threeGegen1otwo
gslpp::vector< double >	TL0_cache
unsigned int	TL0_updated
unsigned int	TL1_updated
unsigned int	TL2_updated
double	tmpq2
unsigned int	TR0_updated
unsigned int	TR1_updated
unsigned int	TR2_updated
double	twoMboMM
double	twoMc2
double	twoMlepMb
double	twoMM
double	twoMM2
double	twoMM3
double	twoMM_mbpms
double	twoMV2
double	twothird
gslpp::complex	ubar
gslpp::vector< double >	V_cache
unsigned int	V_updated
StandardModel::meson	vectorM
gslpp::vector< double >	VL0_cache
unsigned int	VL0_updated
unsigned int	VL1_updated
unsigned int	VL2_updated
unsigned int	VR0_updated
unsigned int	VR1_updated
unsigned int	VR2_updated
gsl_integration_cquad_workspace *	w_delta
gsl_integration_cquad_workspace *	w_DTPPR
gsl_integration_cquad_workspace *	w_sigma
double	width
gslpp::complex	xm
gslpp::complex	xp
gslpp::vector< double >	Ycache
gslpp::complex	ym
gslpp::complex	yp
unsigned int	Yupdated
double	z_0
unsigned int	z_updated

Constructor & Destructor Documentation

MVll::MVll	(	const StandardModel &	SM_i,
		StandardModel::meson	meson_i,
		StandardModel::meson	vector_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]	vector_i	final vector meson of the decay
[in]	lep_i	final leptons of the decay

Definition at line 16 of file MVll.cpp.

: mySM(SM_i),
N_cache(3, 0.),
V_cache(3, 0.),
A0_cache(3, 0.),
A1_cache(3, 0.),
T1_cache(3, 0.),
T2_cache(3, 0.),
k2_cache(2, 0.),
VL0_cache(3, 0.),
TL0_cache(3, 0.),
SL_cache(2, 0.),
Ycache(2, 0.),
H_V0cache(2, 0.),
H_V1cache(2, 0.),
H_V2cache(2, 0.),
H_Scache(2, 0.),
H_Pcache(4, 0.),
T_cache(5, 0.) 
{

    
    lep = lep_i;
    meson = meson_i;
    vectorM = vector_i;
    I0_updated = 0;
    I1_updated = 0;
    I2_updated = 0;
    I3_updated = 0;
    I4_updated = 0;
    I5_updated = 0;
    I6_updated = 0;
    I7_updated = 0;
    I8_updated = 0;
    I9_updated = 0;
    I10_updated = 0;
    I11_updated = 0;

    VL1_updated = 0;
    VL2_updated = 0;
    TL1_updated = 0;
    TL2_updated = 0;
    VR1_updated = 0;
    VR2_updated = 0;
    TR1_updated = 0;
    TR2_updated = 0;
    VL0_updated = 0;
    TL0_updated = 0;
    VR0_updated = 0;
    TR0_updated = 0;
    SL_updated = 0;
    SR_updated = 0;
    
    deltaTparpupdated = 0;
    deltaTparmupdated = 0;
    deltaTperpupdated = 0;

    w_sigma = gsl_integration_cquad_workspace_alloc (100);   
    w_DTPPR = gsl_integration_cquad_workspace_alloc (100);
    w_delta = gsl_integration_cquad_workspace_alloc (100);

}

MVll::~MVll ( )

virtual

Destructor.

Definition at line 79 of file MVll.cpp.

{
}

Member Function Documentation

double MVll::A_0 ( double  q2 )

private

The transverse form factor \( A_0 \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( A_0 \)

Definition at line 764 of file MVll.cpp.

{
    return FF_fit(q2, a_0A0, a_1A0, a_2A0, MRA0_2);
}

double MVll::A_1 ( double  q2 )

private

The transverse form factor \( A_1 \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( A_1 \)

Definition at line 769 of file MVll.cpp.

{
    return FF_fit(q2, a_0A1, a_1A1, a_2A1, MRA1_2);
}

double MVll::A_2 ( double  q2 )

private

The transverse form factor \( A_2 \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( A_2 \)

Definition at line 774 of file MVll.cpp.

{
    return (MMpMV2 * (MM2mMV2 - q2) * A_1(q2) - 16. * MM * MV2 * MMpMV * FF_fit(q2, a_0A12_LCSR, a_1A12, a_2A12, MRA12_2)) / lambda(q2);
}

double MVll::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 1421 of file MVll.cpp.

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

double MVll::beta2 ( double  q2 )

private

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

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \beta^2 \)

Definition at line 1426 of file MVll.cpp.

{
    return 1. - 4. * Mlep2 / q2;
}

void MVll::checkCache ( )

private

The caching method for MVll.

Definition at line 381 of file MVll.cpp.

{

    if (MM == k2_cache(0) && MV == k2_cache(1)) {
        k2_updated = 1;
        z_updated = 1;
    } else {
        k2_updated = 0;
        z_updated = 0;
        k2_cache(0) = MM;
        k2_cache(1) = MV;
    }

    if (Mlep == beta_cache) {
        beta_updated = 1;
    } else {
        beta_updated = 0;
        beta_cache = Mlep;
    }

    lambda_updated = k2_updated;
    F_updated = lambda_updated * beta_updated;

    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 (a_0V == V_cache(0) && a_1V == V_cache(1) && a_2V == V_cache(2)) {
        V_updated = V_updated * z_updated;
    } else {
        V_updated = 0;
        V_cache(0) = a_0V;
        V_cache(1) = a_1V;
        V_cache(2) = a_2V;
    }

    if (a_0A0 == A0_cache(0) && a_1A0 == A0_cache(1) && a_2A0 == A0_cache(2)) {
        A0_updated = A0_updated * z_updated;
    } else {
        A0_updated = 0;
        A0_cache(0) = a_0A0;
        A0_cache(1) = a_1A0;
        A0_cache(2) = a_2A0;
    }

    if (a_0A1 == A1_cache(0) && a_1A1 == A1_cache(1) && a_2A1 == A1_cache(2)) {
        A1_updated = A1_updated * z_updated;
    } else {
        A1_updated = 0;
        A1_cache(0) = a_0A1;
        A1_cache(1) = a_1A1;
        A1_cache(2) = a_2A1;
    }

    if (a_0T1 == T1_cache(0) && a_1T1 == T1_cache(1) && a_2T1 == T1_cache(2)) {
        T1_updated = T1_updated * z_updated;
    } else {
        T1_updated = 0;
        T1_cache(0) = a_0T1;
        T1_cache(1) = a_1T1;
        T1_cache(2) = a_2T1;
    }

    if (a_0T2 == T2_cache(0) && a_1T2 == T2_cache(1) && a_2T2 == T2_cache(2)) {
        T2_updated = T2_updated * z_updated;
    } else {
        T2_updated = 0;
        T2_cache(0) = a_0T2;
        T2_cache(1) = a_1T2;
        T2_cache(2) = a_2T2;
    }

    VL1_updated = k2_updated * lambda_updated * A1_updated * V_updated;
    VL2_updated = VL1_updated;

    TL1_updated = k2_updated * lambda_updated * T1_updated * T2_updated;
    TL2_updated = TL1_updated;

    VR1_updated = VL2_updated;
    VR2_updated = VL1_updated;

    TR1_updated = TL2_updated;
    TR2_updated = TL1_updated;

    if (Mb == SL_cache(0) && Ms == SL_cache(1)) {
        Mb_Ms_updated = 1;
        SL_updated = lambda_updated * A0_updated;
        SR_updated = SL_updated;
    } else {
        Mb_Ms_updated = 0;
        SL_updated = 0;
        SR_updated = SL_updated;
        SL_cache(0) = Mb;
        SL_cache(1) = Ms;
    }

    if (a_0A12 == VL0_cache(0) && a_1A12 == VL0_cache(1) && a_2A12 == VL0_cache(2)) {
        VL0_updated = VL0_updated * z_updated;
        VR0_updated = VL0_updated;
    } else {
        VL0_updated = 0;
        VR0_updated = VL0_updated;
        VL0_cache(0) = a_0A12;
        VL0_cache(1) = a_1A12;
        VL0_cache(2) = a_2A12;
    }

    if (a_0T23 == TL0_cache(0) && a_1T23 == TL0_cache(1) && a_2T23 == TL0_cache(2)) {
        TL0_updated = TL0_updated * z_updated;
        TR0_updated = TL0_updated;
    } else {
        TL0_updated = 0;
        TR0_updated = TL0_updated;
        TL0_cache(0) = a_0T23;
        TL0_cache(1) = a_1T23;
        TL0_cache(2) = a_2T23;
    }


    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 (C_2Lh_bar == C_2Lh_cache) {
        C_2Lh_updated = 1;
    } else {
        C_2Lh_updated = 0;
        C_2Lh_cache = C_2Lh_bar;
    }

    if (C_8Lh == C_8Lh_cache) {
        C_8Lh_updated = 1;
    } else {
        C_8Lh_updated = 0;
        C_8Lh_cache = C_8Lh;
    }

    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 (h_0[0] == h0Ccache[0] && h_1[0] == h0Ccache[1] && h_2[0] == h0Ccache[2]) {
        h0_updated = 1;
    } else {
        h0_updated = 0;
        h0Ccache[0] = h_0[0];
        h0Ccache[1] = h_1[0];
        h0Ccache[2] = h_2[0];
    }

    if (h_0[1] == h1Ccache[0] && h_1[1] == h1Ccache[1] && h_2[1] == h1Ccache[2]) {
        h1_updated = 1;
    } else {
        h1_updated = 0;
        h1Ccache[0] = h_0[1];
        h1Ccache[1] = h_1[1];
        h1Ccache[2] = h_2[1];
    }

    if (h_0[2] == h2Ccache[0] && h_1[2] == h2Ccache[1] && h_2[2] == h2Ccache[2]) {
        h2_updated = 1;
    } else {
        h2_updated = 0;
        h2Ccache[0] = h_0[2];
        h2Ccache[1] = h_1[2];
        h2Ccache[2] = h_2[2];
    }

    if (MM == H_V0cache(0) && Mb == H_V0cache(1)) {
        H_V0updated = N_updated * C_9_updated * Yupdated * VL0_updated * C_9p_updated * VR0_updated * C_7_updated * TL0_updated * C_7p_updated * TR0_updated * h0_updated;
    } else {
        H_V0updated = 0;
        H_V0cache(0) = MM;
        H_V0cache(1) = Mb;
    }

    if (MM == H_V1cache(0) && Mb == H_V1cache(1)) {
        H_V1updated = N_updated * C_9_updated * Yupdated * VL1_updated * C_9p_updated * VR1_updated * C_7_updated * TL1_updated * C_7p_updated * TR1_updated * h1_updated;
    } else {
        H_V1updated = 0;
        H_V1cache(0) = MM;
        H_V1cache(1) = Mb;
    }

    if (MM == H_V2cache(0) && Mb == H_V2cache(1)) {
        H_V2updated = N_updated * C_9_updated * Yupdated * VL2_updated * C_9p_updated * VR2_updated * C_7_updated * TL2_updated * C_7p_updated * TR2_updated * h2_updated;
    } else {
        H_V2updated = 0;
        H_V2cache(0) = MM;
        H_V2cache(1) = Mb;
    }

    H_A0updated = N_updated * C_10_updated * VL0_updated * C_10p_updated * VR0_updated;
    H_A1updated = N_updated * C_10_updated * VL1_updated * C_10p_updated * VR1_updated;
    H_A2updated = N_updated * C_10_updated * VL2_updated * C_10p_updated * VR2_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_Pcache(0) && MW == H_Pcache(1) && Mlep == H_Pcache(2) && Ms == H_Pcache(3)) {
        H_Pupdated = N_updated * C_P_updated * SL_updated * C_Pp_updated * SR_updated * C_10_updated * C_10p_updated;
    } else {
        H_Pupdated = 0;
        H_Pcache(0) = Mb;
        H_Pcache(1) = MW;
        H_Pcache(2) = Mlep;
        H_Pcache(3) = Ms;

    }
    
        if (MM == T_cache(0) && Mb == T_cache(1) && Mc == T_cache(2) && 
            mySM.getMesons(vectorM).getGegenalpha(0) == T_cache(3) && mySM.getMesons(vectorM).getGegenalpha(1) == T_cache(4) ) {
        T_updated = 1;
    } else {
        T_updated = 0;
        T_cache(0) = MM;
        T_cache(1) = Mb;
        T_cache(2) = Mc;
        T_cache(3) = mySM.getMesons(vectorM).getGegenalpha(0);
        T_cache(4) = mySM.getMesons(vectorM).getGegenalpha(1);
    }

    deltaTparpupdated = C_2Lh_updated * T_updated;
    deltaTparmupdated = C_2Lh_updated * C_8Lh_updated * T_updated;
    deltaTperpupdated = deltaTparpupdated;

    I0_updated = F_updated * H_V0updated * H_A0updated * H_Pupdated * beta_updated * H_Supdated * deltaTparmupdated;
    I1_updated = F_updated * beta_updated * H_V1updated * H_V2updated * H_A1updated * H_A2updated * deltaTparmupdated;
    I2_updated = F_updated * beta_updated * H_V0updated * H_A0updated  * deltaTparmupdated;
    I3_updated = F_updated * H_V1updated * H_V2updated * H_A1updated * H_A2updated * beta_updated * deltaTparmupdated;
    I4_updated = F_updated * H_V1updated * H_V2updated * H_A1updated * H_A2updated * deltaTparmupdated;
    I5_updated = F_updated * H_V0updated * H_V1updated * H_V2updated * H_A0updated * H_A1updated * H_A2updated * beta_updated * deltaTparmupdated;
    I6_updated = F_updated * H_V1updated * H_V2updated * H_A0updated * H_A1updated * H_A2updated * H_V0updated * beta_updated * H_Supdated * deltaTparmupdated;
    I7_updated = I4_updated * beta_updated ;
    I8_updated = F_updated * beta_updated * H_Supdated * H_V0updated * deltaTparmupdated;
    I9_updated = I6_updated;
    I10_updated = I5_updated;
    I11_updated = I7_updated;

}

gslpp::complex MVll::Cpar ( double  q2 )

private

The correction \( C_{\parallel} \) from [40] .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( C_{\parallel} \)

Definition at line 998 of file MVll.cpp.

{
    return 1. / CF * (C_2L_bar * F27(q2) + C_8L * F87(q2) + MM / 2. / Mb * 
            (C_2L_bar * F29(q2) + 2.*C_1L_bar*(F19(q2) + F29(q2)/6.) + C_8L * F89(q2)));
}

gslpp::complex MVll::Cperp ( double  q2 )

private

The correction \( C_{\perp} \) from [40] .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( C_{\perp} \)

Definition at line 992 of file MVll.cpp.

{
    return 1. / CF * (-C_2L_bar * F27(q2) - C_8L * F87(q2) - q2 / 2. / Mb / MM * 
            (C_2L_bar * F29(q2) + 2.*C_1L_bar*(F19(q2) + F29(q2)/6.) + C_8L * F89(q2)));
}

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

private

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 1505 of file MVll.cpp.

{
    return 0; /* FIX CPV */
    //return (I(i, q2,0) - I(i, q2,1))/2;
}

gslpp::complex MVll::DeltaC9_0 ( double  q2 )

private

The total QCDF correction \( \Delta C_9^0 \) computed integrating over \( u \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Delta C_9^0 \)

Definition at line 1064 of file MVll.cpp.

{
    return 1. / 2. / MV / MM / sqrt(q2) * ((MM2mMV2 * (MM2mMV2 - q2) - lambda(q2))* (MM2 - q2) * 
            Mb/MM2/q2 * deltaTperp(q2) - lambda(q2) * (deltaTpar(q2) + deltaTperp(q2))* Mb/MM2mMV2);
}

gslpp::complex MVll::DeltaC9_m ( double  q2 )

private

The total QCDF correction \( \Delta C_9^- \) computed integrating over \( u \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Delta C_9^- \)

Definition at line 1057 of file MVll.cpp.

{
    return 1./q2 * Mb/MM * (MM2mMV2 * (MM2 - q2)/MM2 +
            sqrt(lambda(q2))) * deltaTperp(q2);
}

gslpp::complex MVll::DeltaC9_p ( double  q2 )

private

The total QCDF correction \( \Delta C_9^+ \) computed integrating over \( u \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Delta C_9^+ \)

Definition at line 1051 of file MVll.cpp.

{
    return 1./q2 * Mb/MM * (MM2mMV2 * (MM2 - q2)/MM2 -
            sqrt(lambda(q2))) * deltaTperp(q2);
}

gslpp::complex MVll::deltaTpar ( double  q2 )

private

The total correction \( \Delta \mathcal{T}^{\parallel} \) from [40] .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Delta \mathcal{T}^{\parallel} \)

Definition at line 1024 of file MVll.cpp.

{
    double Lambdaplus = mySM.getMesons(meson).getLambdaM();
    gslpp::complex Lambdamin = exp(-q2 / MM / Lambdaplus) / Lambdaplus * (-gsl_sf_expint_Ei(q2 / MM / Lambdaplus) + gslpp::complex::i() * M_PI);
    double T3q2 = MM2mMV2/lambda(q2) * ( (MM2 + 3.*MV2 - q2) * T_2(q2) - 8.*MM*MV2/MMpMV * FF_fit(q2, a_0T23, a_1T23, a_2T23, MRT23_2) );
    tmpq2 = q2;
    Ee = (MM2 - tmpq2) / twoMM;
    
    if (deltaTparpCached[q2] == 0) {
        
        DTPPR = convertToGslFunction(boost::bind(&MVll::Integrand_ReTpar_pm, &(*this), _1));
        if (gsl_integration_cquad(&DTPPR, 0., 1., 1.e-2, 1.e-1, w_DTPPR, &avaDTPPR, &errDTPPR, NULL) != 0) return std::numeric_limits<double>::quiet_NaN();
        double ReTppint = avaDTPPR;

        DTPPR = convertToGslFunction(boost::bind(&MVll::Integrand_ImTpar_pm, &(*this), _1));
        if (gsl_integration_cquad(&DTPPR, 0., 1., 1.e-2, 1.e-1, w_DTPPR, &avaDTPPR, &errDTPPR, NULL) != 0) return std::numeric_limits<double>::quiet_NaN();
        double ImTppint = avaDTPPR;        

        cacheDeltaTparp[q2] = (ReTppint + gslpp::complex::i() * ImTppint);
        deltaTparpCached[q2] = 1;
    }

    return deltaT_0 * Cpar(q2) + deltaT_1par * MV/Ee / (T_1(q2) - T3q2) * (cacheDeltaTparp[q2]);
}

gslpp::complex MVll::deltaTperp ( double  q2 )

private

The total correction \( \Delta \mathcal{T}^{\perp} \) from [40] .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Delta \mathcal{T}^{\perp} \)

Definition at line 1004 of file MVll.cpp.

{
    tmpq2 = q2;
    if (deltaTperpCached[q2] == 0) {
        
        DTPPR = convertToGslFunction(boost::bind(&MVll::Integrand_ReTperpplus, &(*this), _1));
        if (gsl_integration_cquad(&DTPPR, 0., 1., 1.e-2, 1.e-1, w_DTPPR, &avaDTPPR, &errDTPPR, NULL) != 0) return std::numeric_limits<double>::quiet_NaN();
        double ReTppint = avaDTPPR;
        
        DTPPR = convertToGslFunction(boost::bind(&MVll::Integrand_ImTperpplus, &(*this), _1));
        if (gsl_integration_cquad(&DTPPR, 0., 1., 1.e-2, 1.e-1, w_DTPPR, &avaDTPPR, &errDTPPR, NULL) != 0) return std::numeric_limits<double>::quiet_NaN();
        double ImTppint = avaDTPPR;

        cacheDeltaTperp[q2] = ReTppint + gslpp::complex::i() * ImTppint;
        deltaTperpCached[q2] = 1;
    }

    return deltaT_0 * Cperp(q2) + deltaT_1perp / T_1(q2) / mySM.getMesons(meson).getLambdaM() * cacheDeltaTperp[q2];
}

double MVll::F ( double  q2, double  b_i  )

private

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

Parameters

[in]	q2	\(q^2\) of the decay
[in]	b_i	BF of the decay \( V \to M_1 M_2\)

Returns

\( F \)

Definition at line 1436 of file MVll.cpp.

{
    return sqrt(lambda(q2)) * beta(q2) * q2 * b_i / (ninetysixM_PI3MM3);
}

gslpp::complex MVll::F19 ( double  q2 )

private

The correction \( F_{19} \) from [15].

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( F_{19} \)

Definition at line 946 of file MVll.cpp.

{
    double s = q2 / Mb2;
    double s2 = s*s;
    double Ls = log(s);
    double Lc = log(Mc/Mb);
    double Lm = log(mu_b/Mb);
    gslpp::complex i = gslpp::complex::i();
    return (-1424./729. + 16./243.*i*M_PI + 64./27.*Lc)*Lm - 16./243.*Lm*Ls + (16./1215. - 32./135./Mc2*Mb2)*Lm*s
            +(4./2835. - 8./315./Mc2*Mb2/Mc2*Mb2)*Lm*s2 + (16./76545. - 32./8505/Mc2*Mb2/Mc2*Mb2/Mc2*Mb2)*Lm*s*s2 - 256./243.*Lm*Lm
            +(-11.65 + 0.18223*i + (-24.709 - 0.13474*i) * s + (-43.588 - 0.4738*i) *s2 + (-86.22 - 1.3542 * i) *s*s2 
            + (-0.080959 - 0.051864*i + (-0.036585 + 0.024753*i) * s + (-0.021692 + 0.036925*i) *s2 + (0.013282 + 0.052023 * i) *s*s2)*Ls );
}

gslpp::complex MVll::F27 ( double  q2 )

private

The correction \( F_{27} \) from [15].

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( F_{27} \)

Definition at line 960 of file MVll.cpp.

{
    double s = q2 / Mb2;
    double s2 = s*s;
    double Ls = log(s);
    gslpp::complex i = gslpp::complex::i();
    return F27_0 + F27_1 * s + F27_2 * s2 + F27_3 * s * s2 + F27_L1_1 * Ls * s + F27_L1_2 * Ls * s2 + F27_L1_3 * Ls * s * s2;
}

gslpp::complex MVll::F29 ( double  q2 )

private

The correction \( F_{29} \) from [15].

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( F_{29} \)

Definition at line 969 of file MVll.cpp.

{
    double s = q2 / Mb2;
    double s2 = s*s;
    double Ls = log(s);
    gslpp::complex i = gslpp::complex::i();
    return F29_0 + F29_L1 * Ls + F29_1 * s +F29_2 * s2 +F29_3 * s * s2 + F29_L1_1 * Ls * s + F29_L1_2 * Ls * s2 + F29_L1_3 * Ls * s2 *s;
}

gslpp::complex MVll::F87 ( double  q2 )

private

The correction \( F_{87} \) from [15].

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( F_{87} \)

Definition at line 978 of file MVll.cpp.

{
    double s = q2 / Mb2;
    double s2 = s*s;
    return F87_0 + F87_1 * s + F87_2 * s2 + F87_3 * s * s2 - 0.888889 * log(s)*(1. + s + s2 + s * s2);
}

double MVll::F89 ( double  q2 )

private

The correction \( F_{89} \) from [15].

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( F_{89} \)

Definition at line 985 of file MVll.cpp.

{
    double s = q2 / Mb2;
    double s2 = s*s;
    return F89_0 + F89_1 * s + F89_2 * s2 + F89_3 * s * s2 + 1.77778 * log(s)*(1. + s + s2 + s * s2);
}

gslpp::complex MVll::fDeltaC9_0 ( double  q2 )

private

The total QCDF correction \( \Delta C_9^0 \) computed fitting over \( u \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Delta C_9^+ \)

Definition at line 1320 of file MVll.cpp.

{
    switch (lep) {
        
        case StandardModel::MU:
            if (q2 < SWITCH) return (reDC9fit(&q2, const_cast<double *>(refres_0_mumu->GetParams())) + gslpp::complex::i()*imDC9fit(&q2, const_cast<double *>(imfres_0_mumu->GetParams())));
            else return (reDC9fit(&q2, const_cast<double *>(refres_0_mumu->GetParams())) + gslpp::complex::i()*imDC9fit(&q2, const_cast<double *>(imfres_0_mumu->GetParams())))/q2;
            break;
        case StandardModel::ELECTRON:
            return (reDC9fit(&q2, const_cast<double *>(refres_0_ee->GetParams())) + gslpp::complex::i()*imDC9fit(&q2, const_cast<double *>(imfres_0_ee->GetParams())));
            break;
        default:
            std::stringstream out;
            out << lep;
            throw std::runtime_error("MVll::fDeltaC9_0 : " + out.str() + " not implemented");
    }
}

gslpp::complex MVll::fDeltaC9_m ( double  q2 )

private

The total QCDF correction \( \Delta C_9^+ \) computed fitting over \( u \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Delta C_9^+ \)

Definition at line 1301 of file MVll.cpp.

{
    switch (lep) {
        
        case StandardModel::MU:
            if (q2 < SWITCH) return (reDC9fit(&q2, const_cast<double *>(refres_m_mumu->GetParams())) + gslpp::complex::i()*imDC9fit(&q2, const_cast<double *>(imfres_m_mumu->GetParams())));
            else return (reDC9fit(&q2, const_cast<double *>(refres_m_mumu->GetParams())) + gslpp::complex::i()*imDC9fit(&q2, const_cast<double *>(imfres_m_mumu->GetParams())))/q2;
            break;
        case StandardModel::ELECTRON:
            return (reDC9fit(&q2, const_cast<double *>(refres_m_ee->GetParams())) + gslpp::complex::i()*imDC9fit(&q2, const_cast<double *>(imfres_m_ee->GetParams())));
            break;
        default:
            std::stringstream out;
            out << lep;
            throw std::runtime_error("MVll::fDeltaC9_m : " + out.str() + " not implemented");
    }
}

gslpp::complex MVll::fDeltaC9_p ( double  q2 )

private

The total QCDF correction \( \Delta C_9^0 \) computed fitting over \( u \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Delta C_9^0 \)

Definition at line 1283 of file MVll.cpp.

{
    switch (lep) {
        
        case StandardModel::MU:
            if (q2 < SWITCH) return (reDC9fit(&q2, const_cast<double *>(refres_p_mumu->GetParams())) + gslpp::complex::i()*imDC9fit(&q2, const_cast<double *>(imfres_p_mumu->GetParams())));
            else return (reDC9fit(&q2, const_cast<double *>(refres_p_mumu->GetParams())) + gslpp::complex::i()*imDC9fit(&q2, const_cast<double *>(imfres_p_mumu->GetParams())))/q2;
            break;
        case StandardModel::ELECTRON:
            return (reDC9fit(&q2, const_cast<double *>(refres_p_ee->GetParams())) + gslpp::complex::i()*imDC9fit(&q2, const_cast<double *>(imfres_p_ee->GetParams())));;
            break;
        default:
            std::stringstream out;
            out << lep;
            throw std::runtime_error("MVll::fDeltaC9_p : " + out.str() + " not implemented");
    }
}

double MVll::FF_fit ( double  q2, double  a_0, double  a_1, double  a_2, double  MR2  )

private

The fit function from [42], \( FF^{\rm fit} \).

Parameters

[in]	q2	\(q^2\) of the decay
[in]	a_0	fit parameter
[in]	a_1	fit parameter
[in]	a_2	fit parameter
[in]	MR2	square of the nearest resonance mass

Returns

\( FF^{\rm fit} \)

Definition at line 749 of file MVll.cpp.

{
    return 1. / (1. - q2 / MR_2) * (a_0 + a_1 * (z(q2) - z_0) + a_2 * (z(q2) - z_0) * (z(q2) - z_0));
}

void MVll::fit_DeltaC9_0_ee ( )

private

The fitting routine for the QCDF correction \( \Delta C_9^0 \) in the electron channel.

Definition at line 1247 of file MVll.cpp.

{
    int dim = 0;
    for (double i=0.002; i<0.05; i+=0.005) {
        double q2tmp = i;
        myq2.push_back(q2tmp);
        ReDeltaC9_0_ee.push_back((1. / 2. / MV / MM / sqrt(q2tmp) * ((MM2mMV2 * (MM2mMV2 - q2tmp) - lambda(q2tmp))* (MM2 - q2tmp) *
                                                             Mb/MM2/q2tmp * deltaTperp(q2tmp) - lambda(q2tmp) * (deltaTpar(q2tmp) + deltaTperp(q2tmp))* Mb/MM2mMV2)).real());
        ImDeltaC9_0_ee.push_back((1. / 2. / MV / MM / sqrt(q2tmp) * ((MM2mMV2 * (MM2mMV2 - q2tmp) - lambda(q2tmp))* (MM2 - q2tmp) *
                                                             Mb/MM2/q2tmp * deltaTperp(q2tmp) - lambda(q2tmp) * (deltaTpar(q2tmp) + deltaTperp(q2tmp))* Mb/MM2mMV2)).imag());
        dim++;
    }
    for (double i=0.05; i<1.12; i+=0.05) {
        double q2tmp = i;
        myq2.push_back(q2tmp);
        ReDeltaC9_0_ee.push_back((1. / 2. / MV / MM / sqrt(q2tmp) * ((MM2mMV2 * (MM2mMV2 - q2tmp) - lambda(q2tmp))* (MM2 - q2tmp) *
                                                             Mb/MM2/q2tmp * deltaTperp(q2tmp) - lambda(q2tmp) * (deltaTpar(q2tmp) + deltaTperp(q2tmp))* Mb/MM2mMV2)).real());
        ImDeltaC9_0_ee.push_back((1. / 2. / MV / MM / sqrt(q2tmp) * ((MM2mMV2 * (MM2mMV2 - q2tmp) - lambda(q2tmp))* (MM2 - q2tmp) *
                                                             Mb/MM2/q2tmp * deltaTperp(q2tmp) - lambda(q2tmp) * (deltaTpar(q2tmp) + deltaTperp(q2tmp))* Mb/MM2mMV2)).imag());
        dim++;
    }
    
    gr1 = TGraph(dim, myq2.data(), ReDeltaC9_0_ee.data());
    gr2 = TGraph(dim, myq2.data(), ImDeltaC9_0_ee.data());
    
    reffit = TF1("reffit",this,&MVll::reDC9fit,0,8.1,7,"MVll","reDC9fit");
    imffit = TF1("imffit",this,&MVll::imDC9fit,0,8.1,8,"MVll","imDC9fit");
    
    refres_0_ee = gr1.Fit(&reffit, "SQN0+rob=0.99");
    imfres_0_ee = gr2.Fit(&imffit, "SQN0+rob=0.99");
    
    ReDeltaC9_0_ee.clear();
    ImDeltaC9_0_ee.clear();
    myq2.clear();
}

void MVll::fit_DeltaC9_0_mumu ( )

private

The fitting routine for the QCDF correction \( \Delta C_9^0 \) in the muon channel.

Definition at line 1211 of file MVll.cpp.

{
    int dim = 0;
    for (double i=0.1; i<SWITCH; i+=0.4) {
        double q2tmp = i;
        myq2.push_back(q2tmp);
        ReDeltaC9_0_mumu.push_back((1. / 2. / MV / MM / sqrt(q2tmp) * ((MM2mMV2 * (MM2mMV2 - q2tmp) - lambda(q2tmp))* (MM2 - q2tmp) *
                                                             Mb/MM2/q2tmp * deltaTperp(q2tmp) - lambda(q2tmp) * (deltaTpar(q2tmp) + deltaTperp(q2tmp))* Mb/MM2mMV2)).real());
        ImDeltaC9_0_mumu.push_back((1. / 2. / MV / MM / sqrt(q2tmp) * ((MM2mMV2 * (MM2mMV2 - q2tmp) - lambda(q2tmp))* (MM2 - q2tmp) *
                                                             Mb/MM2/q2tmp * deltaTperp(q2tmp) - lambda(q2tmp) * (deltaTpar(q2tmp) + deltaTperp(q2tmp))* Mb/MM2mMV2)).imag());
        dim++;
    }
    for (double i=SWITCH; i<8.2; i+=0.4) {
        double q2tmp = i;
        myq2.push_back(q2tmp);
        ReDeltaC9_0_mumu.push_back(q2tmp * (1. / 2. / MV / MM / sqrt(q2tmp) * ((MM2mMV2 * (MM2mMV2 - q2tmp) - lambda(q2tmp))* (MM2 - q2tmp) *
                                                             Mb/MM2/q2tmp * deltaTperp(q2tmp) - lambda(q2tmp) * (deltaTpar(q2tmp) + deltaTperp(q2tmp))* Mb/MM2mMV2)).real());
        ImDeltaC9_0_mumu.push_back(q2tmp * (1. / 2. / MV / MM / sqrt(q2tmp) * ((MM2mMV2 * (MM2mMV2 - q2tmp) - lambda(q2tmp))* (MM2 - q2tmp) *
                                                             Mb/MM2/q2tmp * deltaTperp(q2tmp) - lambda(q2tmp) * (deltaTpar(q2tmp) + deltaTperp(q2tmp))* Mb/MM2mMV2)).imag());
        dim++;
    }
    
    gr1 = TGraph(dim, myq2.data(), ReDeltaC9_0_mumu.data());
    gr2 = TGraph(dim, myq2.data(), ImDeltaC9_0_mumu.data());
    
    reffit = TF1("reffit",this,&MVll::reDC9fit,0,8.1,7,"MVll","reDC9fit");
    imffit = TF1("imffit",this,&MVll::imDC9fit,0,8.1,8,"MVll","imDC9fit");
    
    refres_0_mumu = gr1.Fit(&reffit, "SQN0+rob=0.99");
    imfres_0_mumu = gr2.Fit(&imffit, "SQN0+rob=0.99");
    
    ReDeltaC9_0_mumu.clear();
    ImDeltaC9_0_mumu.clear();
    myq2.clear();
}

void MVll::fit_DeltaC9_m_ee ( )

private

The fitting routine for the QCDF correction \( \Delta C_9^- \) in the electron channel.

Definition at line 1184 of file MVll.cpp.

{
    int dim = 0;
    for (double i=0.002; i<1.12; i+=0.05) {
        double q2tmp = i;
        myq2.push_back(q2tmp);
        ReDeltaC9_m_ee.push_back((1./q2tmp * Mb/MM * (MM2mMV2 * (MM2 - q2tmp)/MM2 +
                                              sqrt(lambda(q2tmp))) * deltaTperp(q2tmp)).real());
        ImDeltaC9_m_ee.push_back((1./q2tmp * Mb/MM * (MM2mMV2 * (MM2 - q2tmp)/MM2 +
                                              sqrt(lambda(q2tmp))) * deltaTperp(q2tmp)).imag());
        dim++;
    }
    
    gr1 = TGraph(dim, myq2.data(), ReDeltaC9_m_ee.data());
    gr2 = TGraph(dim, myq2.data(), ImDeltaC9_m_ee.data());
    
    reffit = TF1("reffit",this,&MVll::reDC9fit,0,8.1,7,"MVll","reDC9fit");
    imffit = TF1("imffit",this,&MVll::imDC9fit,0,8.1,8,"MVll","imDC9fit");
    
    refres_m_ee = gr1.Fit(&reffit, "SQN0+rob=0.99");
    imfres_m_ee = gr2.Fit(&imffit, "SQN0+rob=0.99");
    
    ReDeltaC9_m_ee.clear();
    ImDeltaC9_m_ee.clear();
    myq2.clear();
}

void MVll::fit_DeltaC9_m_mumu ( )

private

The fitting routine for the QCDF correction \( \Delta C_9^- \) in the muon channel.

Definition at line 1148 of file MVll.cpp.

{
    int dim = 0;
    for (double i=0.1; i<SWITCH; i+=0.4) {
        double q2tmp = i;
        myq2.push_back(q2tmp);
        ReDeltaC9_m_mumu.push_back((1./q2tmp * Mb/MM * (MM2mMV2 * (MM2 - q2tmp)/MM2 +
                                              sqrt(lambda(q2tmp))) * deltaTperp(q2tmp)).real());
        ImDeltaC9_m_mumu.push_back((1./q2tmp * Mb/MM * (MM2mMV2 * (MM2 - q2tmp)/MM2 +
                                              sqrt(lambda(q2tmp))) * deltaTperp(q2tmp)).imag());
        dim++;
    }
    for (double i=SWITCH; i<8.2; i+=0.4) {
        double q2tmp = i;
        myq2.push_back(q2tmp);
        ReDeltaC9_m_mumu.push_back(q2tmp * (1./q2tmp * Mb/MM * (MM2mMV2 * (MM2 - q2tmp)/MM2 +
                                              sqrt(lambda(q2tmp))) * deltaTperp(q2tmp)).real());
        ImDeltaC9_m_mumu.push_back(q2tmp * (1./q2tmp * Mb/MM * (MM2mMV2 * (MM2 - q2tmp)/MM2 +
                                              sqrt(lambda(q2tmp))) * deltaTperp(q2tmp)).imag());
        dim++;
    }
    
    gr1 = TGraph(dim, myq2.data(), ReDeltaC9_m_mumu.data());
    gr2 = TGraph(dim, myq2.data(), ImDeltaC9_m_mumu.data());
    
    reffit = TF1("reffit",this,&MVll::reDC9fit,0,8.1,7,"MVll","reDC9fit");
    imffit = TF1("imffit",this,&MVll::imDC9fit,0,8.1,8,"MVll","imDC9fit");
    
    refres_m_mumu = gr1.Fit(&reffit, "SQN0+rob=0.99");
    imfres_m_mumu = gr2.Fit(&imffit, "SQN0+rob=0.99");
    
    ReDeltaC9_m_mumu.clear();
    ImDeltaC9_m_mumu.clear();
    myq2.clear();
}

void MVll::fit_DeltaC9_p_ee ( )

private

The fitting routine for the QCDF correction \( \Delta C_9^+ \) in the electron channel.

Definition at line 1121 of file MVll.cpp.

{
    int dim = 0;
    for (double i=0.002; i<1.12; i+=0.05) {
        double q2tmp = i;        
        myq2.push_back(q2tmp);
        ReDeltaC9_p_ee.push_back((1./q2tmp * Mb/MM * (MM2mMV2 * (MM2 - q2tmp)/MM2 -
                                              sqrt(lambda(q2tmp))) * deltaTperp(q2tmp)).real());
        ImDeltaC9_p_ee.push_back((1./q2tmp * Mb/MM * (MM2mMV2 * (MM2 - q2tmp)/MM2 -
                                              sqrt(lambda(q2tmp))) * deltaTperp(q2tmp)).imag());
        dim++;
    }
    
    gr1 =TGraph(dim, myq2.data(), ReDeltaC9_p_ee.data());
    gr2 =TGraph(dim, myq2.data(), ImDeltaC9_p_ee.data());
    
    reffit = TF1("reffit",this,&MVll::reDC9fit,0,8.1,7,"MVll","reDC9fit");
    imffit = TF1("imffit",this,&MVll::imDC9fit,0,8.1,8,"MVll","imDC9fit");
    
    refres_p_ee = gr1.Fit(&reffit, "SQN0+rob=0.99");
    imfres_p_ee = gr2.Fit(&imffit, "SQN0+rob=0.99");
    
    ReDeltaC9_p_ee.clear();
    ImDeltaC9_p_ee.clear();
    myq2.clear();
}

void MVll::fit_DeltaC9_p_mumu ( )

private

The fitting routine for the QCDF correction \( \Delta C_9^+ \) in the muon channel.

Definition at line 1086 of file MVll.cpp.

{
    int dim = 0;
    for (double i=0.1; i<SWITCH; i+=0.4) {
        double q2tmp = i;        
        myq2.push_back(q2tmp);
        ReDeltaC9_p_mumu.push_back((1./q2tmp * Mb/MM * (MM2mMV2 * (MM2 - q2tmp)/MM2 -
                                              sqrt(lambda(q2tmp))) * deltaTperp(q2tmp)).real());
        ImDeltaC9_p_mumu.push_back((1./q2tmp * Mb/MM * (MM2mMV2 * (MM2 - q2tmp)/MM2 -
                                              sqrt(lambda(q2tmp))) * deltaTperp(q2tmp)).imag());
        dim++;
    }
    for (double i=SWITCH; i<8.2; i+=0.4) {
        double q2tmp = i;        
        myq2.push_back(q2tmp);
        ReDeltaC9_p_mumu.push_back(q2tmp * (1./q2tmp * Mb/MM * (MM2mMV2 * (MM2 - q2tmp)/MM2 -
                                              sqrt(lambda(q2tmp))) * deltaTperp(q2tmp)).real());
        ImDeltaC9_p_mumu.push_back(q2tmp * (1./q2tmp * Mb/MM * (MM2mMV2 * (MM2 - q2tmp)/MM2 -
                                              sqrt(lambda(q2tmp))) * deltaTperp(q2tmp)).imag());
        dim++;
    }
    gr1 =TGraph(dim, myq2.data(), ReDeltaC9_p_mumu.data());
    gr2 =TGraph(dim, myq2.data(), ImDeltaC9_p_mumu.data());
    
    reffit = TF1("reffit",this,&MVll::reDC9fit,0.1,8.1,7,"MVll","reDC9fit");
    imffit = TF1("imffit",this,&MVll::imDC9fit,0.1,8.1,8,"MVll","imDC9fit");
    
    refres_p_mumu = gr1.Fit(&reffit, "SQN0+rob=0.99");
    imfres_p_mumu = gr2.Fit(&imffit, "SQN0+rob=0.99");
    
    ReDeltaC9_p_mumu.clear();
    ImDeltaC9_p_mumu.clear();
    myq2.clear();
}

double MVll::getDelta0 ( double  q2 )

inlineprivate

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

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Delta_{1s} \)

Definition at line 1275 of file MVll.h.

    {
        return Delta(0, q2);
    };

double MVll::getDelta1 ( double  q2 )

inlineprivate

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

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Delta_{1c} \)

Definition at line 1285 of file MVll.h.

    {
        return Delta(1, q2);
    };

double MVll::getDelta11 ( double  q2 )

inlineprivate

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

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Delta_{9} \)

Definition at line 1325 of file MVll.h.

    {
        return Delta(11, q2);
    };

double MVll::getDelta2 ( double  q2 )

inlineprivate

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

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Delta_{2s} \)

Definition at line 1295 of file MVll.h.

    {
        return Delta(2, q2);
    };

double MVll::getDelta3 ( double  q2 )

inlineprivate

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

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Delta_{2c} \)

Definition at line 1305 of file MVll.h.

    {
        return Delta(3, q2);
    };

double MVll::getDelta7 ( double  q2 )

inlineprivate

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

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Delta_{6s} \)

Definition at line 1315 of file MVll.h.

    {
        return Delta(7, q2);
    };

double MVll::getgtilde_1_im ( double  q2 )

inline

The immaginary part of \( \tilde{g}^1 \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \mathrm{IM}(\tilde{g}^1) \)

Definition at line 303 of file MVll.h.

    {
        updateParameters();
        return C2_inv * (gtilde_1_pre/(sqrt(lambda(q2)) * V(q2)) * (h_0[2]/q2 + h_1[2] + h_2[2] * q2 - (h_0[1]/q2 + h_1[1] + h_2[1] * q2))).imag();
    }

double MVll::getgtilde_1_re ( double  q2 )

inline

The real part of \( \tilde{g}^1 \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \mathrm{RE}(\tilde{g}^1) \)

Definition at line 292 of file MVll.h.

    {
        updateParameters();
        return C2_inv * (gtilde_1_pre/(sqrt(lambda(q2)) * V(q2)) * (h_0[2]/q2 + h_1[2] + h_2[2] * q2 - (h_0[1]/q2 + h_1[1] + h_2[1] * q2))).real();
    }

double MVll::getgtilde_2_im ( double  q2 )

inline

The immaginary part of \( \tilde{g}^2 \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \mathrm{IM}(\tilde{g}^2) \)

Definition at line 325 of file MVll.h.

    {
        updateParameters();
        return C2_inv * (gtilde_2_pre/A_1(q2) * (h_0[2]/q2 + h_1[2] + h_2[2] * q2 + h_0[1]/q2 + h_1[1] + h_2[1] * q2)).imag();
    }

double MVll::getgtilde_2_re ( double  q2 )

inline

The real part of \( \tilde{g}^2 \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \mathrm{RE}(\tilde{g}^2) \)

Definition at line 314 of file MVll.h.

    {
        updateParameters();
        return C2_inv * (gtilde_2_pre/A_1(q2) * (h_0[2]/q2 + h_1[2] + h_2[2] * q2 + h_0[1]/q2 + h_1[1] + h_2[1] * q2)).real();
    }

double MVll::getgtilde_3_im ( double  q2 )

inline

The imaginary part of \( \tilde{g}^3 \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \mathrm{IM}(\tilde{g}^3) \)

Definition at line 349 of file MVll.h.

    {
        updateParameters();
        return C2_inv * (gtilde_3_pre/(lambda(q2) * A_2(q2)) * (sqrt(q2)*(h_0[0]/q2 + 
                h_1[0] + h_2[0] * q2)-(MM2mMV2 - q2)/(4.*MV) * (h_0[2]/q2 + h_1[2] + 
                h_2[2] * q2 + h_0[1]/q2 + h_1[1] + h_2[1] * q2))).imag();
    }

double MVll::getgtilde_3_re ( double  q2 )

inline

The real part of \( \tilde{g}^3 \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \mathrm{RE}(\tilde{g}^3) \)

Definition at line 336 of file MVll.h.

    {
        updateParameters();
        return C2_inv * (gtilde_3_pre/(lambda(q2) * A_2(q2)) * (sqrt(q2)*(h_0[0]/q2 + 
                h_1[0] + h_2[0] * q2)-(MM2mMV2 - q2)/(4.*MV) * (h_0[2]/q2 + h_1[2] + 
                h_2[2] * q2 + h_0[1]/q2 + h_1[1] + h_2[1] * q2))).real();
    }

double MVll::geth_0_im ( double  q2 )

inline

The imaginary part of \( h_0 \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \mathrm{IM}(h_0) \)

Definition at line 372 of file MVll.h.

    {
        return (sixteenM_PI2MM2 * (h_0[0]/q2 + h_1[0] + h_2[0] * q2)).imag();
    }

double MVll::geth_0_re ( double  q2 )

inline

The real part of \( h_0 \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \mathrm{RE}(h_0) \)

Definition at line 362 of file MVll.h.

    {
        return (sixteenM_PI2MM2 * (h_0[0]/q2 + h_1[0] + h_2[0] * q2)).real();
    }

double MVll::geth_m_im ( double  q2 )

inline

The imaginary part of \( h_- \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \mathrm{IM}(h_-) \)

Definition at line 412 of file MVll.h.

    {
        return (sixteenM_PI2MM2 * (h_0[2]/q2 + h_1[2] + h_2[2] * q2)).imag();
    }

double MVll::geth_m_re ( double  q2 )

inline

The real part of \( h_- \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \mathrm{RE}(h_-) \)

Definition at line 402 of file MVll.h.

    {
        return (sixteenM_PI2MM2 * (h_0[2]/q2 + h_1[2] + h_2[2] * q2)).real();
    }

double MVll::geth_p_im ( double  q2 )

inline

The imaginary part of \( h_+ \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \mathrm{IM}(h_+) \)

Definition at line 392 of file MVll.h.

    {
        return (sixteenM_PI2MM2 * (h_0[1]/q2 + h_1[1] + h_2[1] * q2)).imag();
    }

double MVll::geth_p_re ( double  q2 )

inline

The real part of \( h_+ \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \mathrm{RE}(h_+) \)

Definition at line 382 of file MVll.h.

    {
        return (sixteenM_PI2MM2 * (h_0[1]/q2 + h_1[1] + h_2[1] * q2)).real();
    }

double MVll::gethm0_h00_abs ( )

inline

The absolute value of the ratio \( h_-^{(0)}/h_0^{(0)} \).

Returns

\( h_-^{(0)}/h_0^{(0)} \)

Definition at line 431 of file MVll.h.

    {
        updateParameters();
        return (h_0[2]/h_0[0]).abs();
    }

double MVll::gethp0_hm0_abs ( )

inline

The absolute value of the ratio \( h_+^{(0)}/h_0^{(0)} \).

Returns

\( h_+^{(0)}/h_0^{(0)} \)

Definition at line 421 of file MVll.h.

    {
        updateParameters();
        return (h_0[1]/h_0[2]).abs();
    }

double MVll::getMlep ( )

inline

The mass of the lepton l.

Returns

\( m_l \)

Definition at line 143 of file MVll.h.

                    {
        updateParameters();
        return Mlep;
    }

double MVll::getS ( double  q2 )

inline

The form factor \( S \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( S \)

Definition at line 225 of file MVll.h.

    {
        updateParameters();
        return S_L_pre/sqrt(lambda(q2)) * S_L(q2);
    };

double MVll::getSigma	(	int	i,
		double	q_2
	)

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

Parameters

[in]	i	index of the angular coefficient \( I_{i} \)
[in]

Definition at line 1630 of file MVll.cpp.

{
    updateParameters();

    switch (i) {
        case 0:
            return getSigma1c(q_2);
            break;
        case 1:
            return getSigma1s(q_2);
            break;
        case 2:
            return getSigma2c(q_2);
            break;
        case 3:
            return getSigma2s(q_2);
            break;
        case 4:
            return getSigma3(q_2);
            break;
        case 5:
            return getSigma4(q_2);
            break;
        case 6:
            return getSigma5(q_2);
            break;
        case 7:
            return getSigma6s(q_2);
            break;
        case 9:
            return getSigma7(q_2);
            break;
        case 10:
            return getSigma8(q_2);
            break;
        case 11:
            return getSigma9(q_2);
            break;
        default:
            std::stringstream out;
            out << i;
            throw std::runtime_error("MVll::integrateSigma: index " + out.str() + " not implemented");
    }
}

double MVll::getSigma1c ( double  q2 )

inlineprivate

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

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Sigma_{1s} \)

Definition at line 1155 of file MVll.h.

    {
        return I_1c(q2);
    };

double MVll::getSigma1s ( double  q2 )

inlineprivate

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

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Sigma_{1c} \)

Definition at line 1165 of file MVll.h.

    {
        return I_1s(q2);
    };

double MVll::getSigma2c ( double  q2 )

inlineprivate

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

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Sigma_{2s} \)

Definition at line 1175 of file MVll.h.

    {
        return I_2c(q2);
    };

double MVll::getSigma2s ( double  q2 )

inlineprivate

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

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Sigma_{2c} \)

Definition at line 1185 of file MVll.h.

    {
        return I_2s(q2);
    };

double MVll::getSigma3 ( double  q2 )

inlineprivate

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

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Sigma_{3} \)

Definition at line 1195 of file MVll.h.

    {
        return I_3(q2);
    };

double MVll::getSigma4 ( double  q2 )

inlineprivate

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

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Sigma_{4} \)

Definition at line 1205 of file MVll.h.

    {
        return I_4(q2);
    };

double MVll::getSigma5 ( double  q2 )

inlineprivate

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

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Sigma_{5} \)

Definition at line 1215 of file MVll.h.

    {
        return I_5(q2);
    };

double MVll::getSigma6c ( double  q2 )

inlineprivate

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

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Sigma_{7} \)

Definition at line 1235 of file MVll.h.

    {
        return I_6c(q2);
    };

double MVll::getSigma6s ( double  q2 )

inlineprivate

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

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Sigma_{6s} \)

Definition at line 1225 of file MVll.h.

    {
        return I_6s(q2);
    };

double MVll::getSigma7 ( double  q2 )

inlineprivate

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

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Sigma_{8} \)

Definition at line 1245 of file MVll.h.

    {
        return I_7(q2);
    };

double MVll::getSigma8 ( double  q2 )

inlineprivate

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

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Sigma_{9} \)

Definition at line 1255 of file MVll.h.

    {
        return I_8(q2);
    };

double MVll::getSigma9 ( double  q2 )

inlineprivate

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

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \Sigma_{9} \)

Definition at line 1265 of file MVll.h.

    {
        return I_9(q2);
    };

double MVll::getT0 ( double  q2 )

inline

The form factor \( T_0 \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( T_0 \)

Definition at line 192 of file MVll.h.

    {
        updateParameters();
        return twoMM3/sqrt(q2 * lambda(q2)) * T_0t(q2);
    };

double MVll::getTm ( double  q2 )

inline

The form factor \( T_- \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( T_- \)

Definition at line 214 of file MVll.h.

    {
        updateParameters();
        return T_m(q2);
    };

double MVll::getTp ( double  q2 )

inline

The form factor \( T_+ \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( T_+ \)

Definition at line 203 of file MVll.h.

    {
        updateParameters();
        return T_p(q2);
    };

double MVll::getV0 ( double  q2 )

inline

The form factor \( V_0 \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( V_0 \)

Definition at line 160 of file MVll.h.

    {
        updateParameters();
        return (2. * MM * sqrt(q2))/sqrt(lambda(q2)) * V_0t(q2);
    };

double MVll::getVm ( double  q2 )

inline

The form factor \( V_- \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( V_- \)

Definition at line 182 of file MVll.h.

    {
        return V_m(q2);
    };

double MVll::getVp ( double  q2 )

inline

The form factor \( V_+ \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( V_+ \)

Definition at line 171 of file MVll.h.

    {
        updateParameters();
        return V_p(q2);
    };

double MVll::getwidth ( )

inline

The width of the meson M.

Returns

\( \Gamma_M \)

Definition at line 134 of file MVll.h.

                     {
        updateParameters();
        return width;
    }

gslpp::complex MVll::H_0 ( double  q2 )

private

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

Parameters

[in]

q2

\(q^2\) of the decay

Returns

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

Definition at line 1363 of file MVll.cpp.

{
    return (H_0_pre - fournineth * log(q2 / mu_b2));
}

gslpp::complex MVll::H_A_0

(

double

q2

)

The helicity amplitude \( H_A^0 \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( H_A^0 \)

Definition at line 1388 of file MVll.cpp.

{
    return ( -C_10 + C_10p) * V_0t(q2);
}

gslpp::complex MVll::H_A_m

(

double

q2

)

The helicity amplitude \( H_A^- \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( H_A^- \)

Definition at line 1398 of file MVll.cpp.

{
    return ( -C_10 * V_m(q2) + C_10p * V_p(q2));
}

gslpp::complex MVll::H_A_p

(

double

q2

)

The helicity amplitude \( H_A^+ \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( H_A^+ \)

Definition at line 1393 of file MVll.cpp.

{
    return ( -C_10 * V_p(q2) + C_10p * V_m(q2));
}

gslpp::complex MVll::H_b ( double  q2 )

private

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

Parameters

[in]

q2

\(q^2\) of the decay

Returns

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

Definition at line 1352 of file MVll.cpp.

{
    double x = fourMb2 / 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)) - ihalfMPI);

    return -fournineth * (logMb - twothird - x) - fournineth * (2. + x) * par;
}

gslpp::complex MVll::H_c ( double  q2, double  mu  )

private

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

Parameters

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

Returns

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

Definition at line 1341 of file MVll.cpp.

{
    double x = fourMc2 / 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)) - ihalfMPI);

    return -fournineth * (log(Mc2 / mu2) - twothird - x) - fournineth * (2. + x) * par;
}

gslpp::complex MVll::H_P

(

double

q2

)

The helicity amplitude \( H_P \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( H_P \)

Definition at line 1408 of file MVll.cpp.

{
    return ( MboMW * (C_P - C_Pp) + twoMlepMb / q2 * (C_10 - C_10p) * (1. + MsoMb)) * S_L(q2);
}

gslpp::complex MVll::H_S

(

double

q2

)

The helicity amplitude \( H_S \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( H_S \)

Definition at line 1403 of file MVll.cpp.

{
    return MboMW * (C_S - C_Sp) * S_L(q2);
}

gslpp::complex MVll::H_V_0

(

double

q2

)

The helicity amplitude \( H_V^0 \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( H_V^0 \)

Definition at line 1373 of file MVll.cpp.

{
    return -(((C_9 + fDeltaC9_0(q2) + Y(q2)) - C_9p) * V_0t(q2) + MM2 / q2 * (twoMboMM * (C_7 - C_7p) * T_0t(q2) - sixteenM_PI2 * (h_0[0] + h_1[0] * q2 + h_2[0] * q2 * q2)));
}

gslpp::complex MVll::H_V_m

(

double

q2

)

The helicity amplitude \( H_V^- \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( H_V^- \)

Definition at line 1383 of file MVll.cpp.

{
    return -(((C_9 + fDeltaC9_m(q2) + Y(q2)) * V_m(q2) - C_9p * V_p(q2)) + MM2 / q2 * (twoMboMM * (C_7 * T_m(q2) - C_7p * T_p(q2)) - sixteenM_PI2 * (h_0[2] + h_1[2] * q2 + h_2[2] * q2 * q2)));
}

gslpp::complex MVll::H_V_p

(

double

q2

)

The helicity amplitude \( H_V^+ \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( H_V^+ \)

Definition at line 1378 of file MVll.cpp.

{
    return -(((C_9 + fDeltaC9_p(q2) + Y(q2)) * V_p(q2) - C_9p * V_m(q2)) + MM2 / q2 * (twoMboMM * (C_7 * T_p(q2) - C_7p * T_m(q2)) - sixteenM_PI2 * (h_0[1] + h_1[1] * q2 + h_2[1] * q2 * q2)));
}

gslpp::complex MVll::I1 ( double  u, double  q2  )

private

The \( I_1 \) function from [40] .

Parameters

[in]	u	dummy variable to be integrated out
[in]	q2	\(q^2\) of the decay

Returns

\( I_1 \)

Definition at line 828 of file MVll.cpp.

{    
    std::pair<double, double > uq2 = std::make_pair(u,q2);

    if (I1Cached[uq2] == 0) {
        ubar = 1. - u;
        xp = .5 + sqrt(0.25 - (Mc2 - gslpp::complex::i()*1.e-10) / (ubar * MM2 + u * q2));
        xm = .5 - sqrt(0.25 - (Mc2 - gslpp::complex::i()*1.e-10) / (ubar * MM2 + u * q2));
        yp = .5 + sqrt(0.25 - (Mc2 - gslpp::complex::i()*1.e-10) / q2);
        ym = .5 - sqrt(0.25 - (Mc2 - gslpp::complex::i()*1.e-10) / q2);
        L1xp = log(1. - 1. / xp) * log(1. - xp) - M_PI2osix + dilog(xp / (xp - 1.));
        L1xm = log(1. - 1. / xm) * log(1. - xm) - M_PI2osix + dilog(xm / (xm - 1.));
        L1yp = log(1. - 1. / yp) * log(1. - yp) - M_PI2osix + dilog(yp / (yp - 1.));
        L1ym = log(1. - 1. / ym) * log(1. - ym) - M_PI2osix + dilog(ym / (ym - 1.));

        cacheI1[uq2] = 1. + twoMc2 / ubar / (MM2 - q2)*(L1xp + L1xm - L1yp - L1ym);
        I1Cached[uq2] = 1;
    }
    
    return cacheI1[uq2];   
}

double MVll::I_1c ( double  q2 )

private

The angular coefficient \( I_{1c} \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( I_{1c} \)

Definition at line 1441 of file MVll.cpp.

{
    return F(q2, b)*((H_V_0(q2).abs2() + H_A_0(q2).abs2()) / 2. + H_P(q2).abs2() + 2. * Mlep2 / q2 * (H_V_0(q2).abs2()
            - H_A_0(q2).abs2()) + beta2(q2) * H_S(q2).abs2());
}

double MVll::I_1s ( double  q2 )

private

The angular coefficient \( I_{1s} \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( I_{1s} \)

Definition at line 1447 of file MVll.cpp.

{
    return F(q2, b)*((beta2(q2) + 2.) / 8. * (H_V_p(q2).abs2() + H_V_m(q2).abs2() + H_A_p(q2).abs2() + H_A_m(q2).abs2()) +
            Mlep2 / q2 * (H_V_p(q2).abs2() + H_V_m(q2).abs2() - H_A_p(q2).abs2() - H_A_m(q2).abs2()));
}

double MVll::I_2c ( double  q2 )

private

The angular coefficient \( I_{2c} \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( I_{2c} \)

Definition at line 1453 of file MVll.cpp.

{
    return -F(q2, b) * beta2(q2) / 2. * (H_V_0(q2).abs2() + H_A_0(q2).abs2());
}

double MVll::I_2s ( double  q2 )

private

The angular coefficient \( I_{2s} \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( I_{2s} \)

Definition at line 1458 of file MVll.cpp.

{
    return F(q2, b) * beta2(q2) / 8. * (H_V_p(q2).abs2() + H_V_m(q2).abs2() + H_A_p(q2).abs2() + H_A_m(q2).abs2());
}

double MVll::I_3 ( double  q2 )

private

The angular coefficient \( I_3 \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( I_3 \)

Definition at line 1463 of file MVll.cpp.

{
    return -F(q2, b) / 2. * ((H_V_p(q2) * H_V_m(q2).conjugate()).real() + (H_A_p(q2) * H_A_m(q2).conjugate()).real());
}

double MVll::I_4 ( double  q2 )

private

The angular coefficient \( I_4 \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( I_4 \)

Definition at line 1468 of file MVll.cpp.

{
    return F(q2, b) * beta2(q2) / 4. * (((H_V_m(q2) + H_V_p(q2)) * H_V_0(q2).conjugate()).real() + ((H_A_m(q2) + H_A_p(q2)) * H_A_0(q2).conjugate()).real());
}

double MVll::I_5 ( double  q2 )

private

The angular coefficient \( I_5 \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( I_5 \)

Definition at line 1473 of file MVll.cpp.

{
    return F(q2, b)*(beta(q2) / 2. * (((H_V_m(q2) - H_V_p(q2)) * H_A_0(q2).conjugate()).real() + ((H_A_m(q2) - H_A_p(q2)) * H_V_0(q2).conjugate()).real()) -
            beta(q2) * Mlep / sqrt(q2)*(H_S(q2).conjugate()*(H_V_p(q2) + H_V_m(q2))).real());
}

double MVll::I_6c ( double  q2 )

private

The angular coefficient \( I_{6c} \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( I_{6c} \)

Definition at line 1484 of file MVll.cpp.

{
    return 2. * F(q2, b) * beta(q2) * Mlep / sqrt(q2)*(H_S(q2).conjugate() * H_V_0(q2)).real();
}

double MVll::I_6s ( double  q2 )

private

The angular coefficient \( I_{6s} \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( I_{6s} \)

Definition at line 1479 of file MVll.cpp.

{
    return F(q2, b) * beta(q2)*(H_V_m(q2)*(H_A_m(q2).conjugate()) - H_V_p(q2)*(H_A_p(q2).conjugate())).real();
}

double MVll::I_7 ( double  q2 )

private

The angular coefficient \( I_7 \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( I_7 \)

Definition at line 1489 of file MVll.cpp.

{
    return F(q2, b)*(beta(q2) / 2. * (((H_V_m(q2) + H_V_p(q2)) * H_A_0(q2).conjugate()).imag() + ((H_A_m(q2) + H_A_p(q2)) * H_V_0(q2).conjugate()).imag()) -
            beta(q2) * Mlep / sqrt(q2)*(H_S(q2).conjugate()*(H_V_m(q2) - H_V_p(q2))).imag());
}

double MVll::I_8 ( double  q2 )

private

The angular coefficient \( I_8 \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( I_8 \)

Definition at line 1495 of file MVll.cpp.

{
    return F(q2, b) * beta2(q2) / 4. * (((H_V_m(q2) - H_V_p(q2)) * H_V_0(q2).conjugate()).imag() + ((H_A_m(q2) - H_A_p(q2)) * H_A_0(q2).conjugate()).imag());
}

double MVll::I_9 ( double  q2 )

private

The angular coefficient \( I_9 \) .

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( I_9 \)

Definition at line 1500 of file MVll.cpp.

{
    return F(q2, b) * beta2(q2) / 2. * ((H_V_p(q2) * H_V_m(q2).conjugate()).imag() + (H_A_p(q2) * H_A_m(q2).conjugate()).imag());
}

double MVll::imDC9fit ( double *  x, double *  p  )

private

The fit function for the imaginary part of the QCDF correction \( \Delta C_9^{\lambda} \).

Parameters

[in]	x	fit variable
[in]	p	fit parameters vector

Returns

\( f_{Im \Delta C_9^{\lambda}} \)

Definition at line 1079 of file MVll.cpp.

{
    return p[0]/x[0] + p[1] + p[2]*x[0] + p[3]*x[0]*x[0] + p[4]*x[0]*x[0]*x[0] + p[5]*x[0]*x[0]*x[0]*x[0] + p[6]*x[0]*x[0]*x[0]*x[0]*x[0] + p[7]*x[0]*x[0]*x[0]*x[0]*x[0]*x[0];
    
    //double thr = 4.*Mc2;
}

double MVll::Integrand_ImTpar_pm ( double  up )

private

The sum of Integrand_ImTparplus() and Integrand_ImTparminus().

Parameters

[in]

up

dummy variable to be integrated out

Returns

\( Im T^{\parallel}_+ \Phi^{\parallel} + Im T^{\parallel}_- \Phi^{\parallel}\)

Definition at line 938 of file MVll.cpp.

                                         {
    return Integrand_ImTparplus(up) + Integrand_ImTparminus(up);
}

double MVll::Integrand_ImTparminus ( double  up )

private

The imaginary part of the integral involving \( T^{\parallel}_- \) at fixed \( q^2 \), according to [40] .

Parameters

[in]

up

dummy variable to be integrated out

Returns

\( Im T^{\parallel}_- \Phi^{\parallel}\)

Definition at line 927 of file MVll.cpp.

{
    double Lambdaplus = mySM.getMesons(meson).getLambdaM();
    gslpp::complex Lambdamin = exp(-tmpq2 / MM / Lambdaplus) / Lambdaplus * (-gsl_sf_expint_Ei(tmpq2 / MM / Lambdaplus) + gslpp::complex::i() * M_PI);
    
    double u = up;
    return ((Tparminus(u, tmpq2)*6. * u * (1. - u)*
            (1 + threeGegen0 * (2. * u - 1)
            + threeGegen1otwo * ((10. * u - 5.)*(2. * u - 1.) - 1.)))/Lambdamin).imag();
}

double MVll::Integrand_ImTparplus ( double  up )

private

The imaginary part of the integral involving \( T^{\parallel}_+ \) at fixed \( q^2 \), according to [40] .

Parameters

[in]

up

dummy variable to be integrated out

Returns

\( Im T^{\parallel}_+ \Phi^{\parallel}\)

Definition at line 908 of file MVll.cpp.

{   
    double u = up;
    return ((Tparplus(u, tmpq2)*6. * u * (1. - u)*
            (1 + threeGegen0 * (2. * u - 1)
            + threeGegen1otwo * ((10. * u - 5.)*(2. * u - 1.) - 1.)))/mySM.getMesons(meson).getLambdaM()).imag();
}

double MVll::Integrand_ImTperpplus ( double  up )

private

The imaginary part of the integral involving \( T^{\perp}_+ \) at fixed \( q^2 \), according to [40] .

Parameters

[in]

up

dummy variable to be integrated out

Returns

\( Im T^{\perp}_+ \Phi^{\perp}\)

Definition at line 892 of file MVll.cpp.

{
    double u = up;
    return (Tperpplus(u, tmpq2)*6. * u * (1. - u)*
            (1 + threeGegen0 * (2. * u - 1)
            + threeGegen1otwo * ((10. * u - 5.)*(2. * u - 1.) - 1.))).imag();
}

double MVll::Integrand_ReTpar_pm ( double  up )

private

The sum of Integrand_ReTparplus() and Integrand_ReTparminus().

Parameters

[in]

up

dummy variable to be integrated out

Returns

\( Re T^{\parallel}_+ \Phi^{\parallel} + Re T^{\parallel}_- \Phi^{\parallel}\)

Definition at line 942 of file MVll.cpp.

                                         {
    return Integrand_ReTparplus(up) + Integrand_ReTparminus(up);
}

double MVll::Integrand_ReTparminus ( double  up )

private

The real part of the integral involving \( T^{\parallel}_- \) at fixed \( q^2 \), according to [40] .

Parameters

[in]

up

dummy variable to be integrated out

Returns

\( Re T^{\parallel}_- \Phi^{\parallel}\)

Definition at line 916 of file MVll.cpp.

{
    double Lambdaplus = mySM.getMesons(meson).getLambdaM();
    gslpp::complex Lambdamin = exp(-tmpq2 / MM / Lambdaplus) / Lambdaplus * (-gsl_sf_expint_Ei(tmpq2 / MM / Lambdaplus) + gslpp::complex::i() * M_PI);
    
    double u = up;
    return ((Tparminus(u, tmpq2)*6. * u * (1. - u)*
            (1 + threeGegen0 * (2. * u - 1)
            + threeGegen1otwo * ((10. * u - 5.)*(2. * u - 1.) - 1.)))/Lambdamin).real();
}

double MVll::Integrand_ReTparplus ( double  up )

private

The real part of the integral involving \( T^{\parallel}_+ \) at fixed \( q^2 \), according to [40] .

Parameters

[in]

up

dummy variable to be integrated out

Returns

\( Re T^{\parallel}_+ \Phi^{\parallel}\)

Definition at line 900 of file MVll.cpp.

{   
    double u = up;
    return ((Tparplus(u, tmpq2)*6. * u * (1. - u)*
            (1 + threeGegen0 * (2. * u - 1)
            + threeGegen1otwo * ((10. * u - 5.)*(2. * u - 1.) - 1.)))/mySM.getMesons(meson).getLambdaM()).real();
}

double MVll::Integrand_ReTperpplus ( double  up )

private

The real part of the integral involving \( T^{\perp}_+ \) at fixed \( q^2 \), according to [40] .

Parameters

[in]

up

dummy variable to be integrated out

Returns

\( Re T^{\perp}_+ \Phi^{\perp}\)

Definition at line 884 of file MVll.cpp.

{
    double u = up;
    return (Tperpplus(u, tmpq2)*6. * u * (1. - u)*
            (1 + threeGegen0 * (2. * u - 1)
            + threeGegen1otwo * ((10. * u - 5.)*(2. * u - 1.) - 1.))).real();
}

double MVll::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 1675 of file MVll.cpp.

{
    updateParameters();

    std::pair<double, double > qbin = std::make_pair(q_min, q_max);

    old_handler = gsl_set_error_handler_off();

    switch (i) {
        case 0:
            if (delta0Cached[qbin] == 0) {
                FD = convertToGslFunction(boost::bind(&MVll::getDelta0, &(*this), _1));
                if (gsl_integration_cquad(&FD, q_min, q_max, 1.e-2, 1.e-1, w_delta, &avaDelta, &errDelta, NULL) != 0) return std::numeric_limits<double>::quiet_NaN();
                cacheDelta0[qbin] = avaDelta;
                delta0Cached[qbin] = 1;
            }
            return cacheDelta0[qbin];
            break;
        case 1:
            if (delta1Cached[qbin] == 0) {
                FD = convertToGslFunction(boost::bind(&MVll::getDelta1, &(*this), _1));
                if (gsl_integration_cquad(&FD, q_min, q_max, 1.e-2, 1.e-1, w_delta, &avaDelta, &errDelta, NULL) != 0) return std::numeric_limits<double>::quiet_NaN();
                cacheDelta1[qbin] = avaDelta;
                delta1Cached[qbin] = 1;
            }
            return cacheDelta1[qbin];
            break;
        case 2:
            if (delta2Cached[qbin] == 0) {
                FD = convertToGslFunction(boost::bind(&MVll::getDelta2, &(*this), _1));
                if (gsl_integration_cquad(&FD, q_min, q_max, 1.e-2, 1.e-1, w_delta, &avaDelta, &errDelta, NULL) != 0) return std::numeric_limits<double>::quiet_NaN();
                cacheDelta2[qbin] = avaDelta;
                delta2Cached[qbin] = 1;
            }
            return cacheDelta2[qbin];
            break;
        case 3:
            if (delta3Cached[qbin] == 0) {
                FD = convertToGslFunction(boost::bind(&MVll::getDelta3, &(*this), _1));
                if (gsl_integration_cquad(&FD, q_min, q_max, 1.e-2, 1.e-1, w_delta, &avaDelta, &errDelta, NULL) != 0) return std::numeric_limits<double>::quiet_NaN();
                cacheDelta3[qbin] = avaDelta;
                delta3Cached[qbin] = 1;
            }
            return cacheDelta3[qbin];
            break;
        case 7:
            if (delta7Cached[qbin] == 0) {
                FD = convertToGslFunction(boost::bind(&MVll::getDelta7, &(*this), _1));
                if (gsl_integration_cquad(&FD, q_min, q_max, 1.e-2, 1.e-1, w_delta, &avaDelta, &errDelta, NULL) != 0) return std::numeric_limits<double>::quiet_NaN();
                cacheDelta7[qbin] = avaDelta;
                delta7Cached[qbin] = 1;
            }
            return cacheDelta7[qbin];
            break;
        case 11:
            if (delta11Cached[qbin] == 0) {
                FD = convertToGslFunction(boost::bind(&MVll::getDelta11, &(*this), _1));
                if (gsl_integration_cquad(&FD, q_min, q_max, 1.e-2, 1.e-1, w_delta, &avaDelta, &errDelta, NULL) != 0) return std::numeric_limits<double>::quiet_NaN();
                cacheDelta11[qbin] = avaDelta;
                delta11Cached[qbin] = 1;
            }
            return cacheDelta11[qbin];
            break;
        default:
            std::stringstream out;
            out << i;
            throw std::runtime_error("integrateDelta: index " + out.str() + " not implemented");
    }

    gsl_set_error_handler(old_handler);

}

double MVll::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 1512 of file MVll.cpp.

{
    updateParameters();

    std::pair<double, double > qbin = std::make_pair(q_min, q_max);

    old_handler = gsl_set_error_handler_off();

    switch (i) {
        case 0:
            if (sigma0Cached[qbin] == 0) {
                FS = convertToGslFunction(boost::bind(&MVll::getSigma1c, &(*this), _1));
                if (gsl_integration_cquad(&FS, q_min, q_max, 1.e-2, 1.e-1, w_sigma, &avaSigma, &errSigma, NULL) != 0) return std::numeric_limits<double>::quiet_NaN();
                cacheSigma0[qbin] = NN*avaSigma;
                sigma0Cached[qbin] = 1;
            }
            return cacheSigma0[qbin];
            break;
        case 1:
            if (sigma1Cached[qbin] == 0) {
                FS = convertToGslFunction(boost::bind(&MVll::getSigma1s, &(*this), _1));
                if (gsl_integration_cquad(&FS, q_min, q_max, 1.e-2, 1.e-1, w_sigma, &avaSigma, &errSigma, NULL) != 0) return std::numeric_limits<double>::quiet_NaN();
                cacheSigma1[qbin] = NN*avaSigma;
                sigma1Cached[qbin] = 1;
            }
            return cacheSigma1[qbin];
            break;
        case 2:
            if (sigma2Cached[qbin] == 0) {
                FS = convertToGslFunction(boost::bind(&MVll::getSigma2c, &(*this), _1));
                if (gsl_integration_cquad(&FS, q_min, q_max, 1.e-2, 1.e-1, w_sigma, &avaSigma, &errSigma, NULL) != 0) return std::numeric_limits<double>::quiet_NaN();
                cacheSigma2[qbin] = NN*avaSigma;
                sigma2Cached[qbin] = 1;
            }
            return cacheSigma2[qbin];
            break;
        case 3:
            if (sigma3Cached[qbin] == 0) {
                FS = convertToGslFunction(boost::bind(&MVll::getSigma2s, &(*this), _1));
                if (gsl_integration_cquad(&FS, q_min, q_max, 1.e-2, 1.e-1, w_sigma, &avaSigma, &errSigma, NULL) != 0) return std::numeric_limits<double>::quiet_NaN();
                cacheSigma3[qbin] = NN*avaSigma;
                sigma3Cached[qbin] = 1;
            }
            return cacheSigma3[qbin];
            break;
        case 4:
            if (sigma4Cached[qbin] == 0) {
                FS = convertToGslFunction(boost::bind(&MVll::getSigma3, &(*this), _1));
                if (gsl_integration_cquad(&FS, q_min, q_max, 1.e-2, 1.e-1, w_sigma, &avaSigma, &errSigma, NULL) != 0) return std::numeric_limits<double>::quiet_NaN();
                cacheSigma4[qbin] = NN*avaSigma;
                sigma4Cached[qbin] = 1;
            }
            return cacheSigma4[qbin];
            break;
        case 5:
            if (sigma5Cached[qbin] == 0) {
                FS = convertToGslFunction(boost::bind(&MVll::getSigma4, &(*this), _1));
                if (gsl_integration_cquad(&FS, q_min, q_max, 1.e-2, 1.e-1, w_sigma, &avaSigma, &errSigma, NULL) != 0) return std::numeric_limits<double>::quiet_NaN();
                cacheSigma5[qbin] = NN*avaSigma;
                sigma5Cached[qbin] = 1;
            }
            return cacheSigma5[qbin];
            break;
        case 6:
            if (sigma6Cached[qbin] == 0) {
                FS = convertToGslFunction(boost::bind(&MVll::getSigma5, &(*this), _1));
                if (gsl_integration_cquad(&FS, q_min, q_max, 1.e-2, 1.e-1, w_sigma, &avaSigma, &errSigma, NULL) != 0) return std::numeric_limits<double>::quiet_NaN();
                cacheSigma6[qbin] = NN*avaSigma;
                sigma6Cached[qbin] = 1;
            }
            return cacheSigma6[qbin];
            break;
        case 7:
            if (sigma7Cached[qbin] == 0) {
                FS = convertToGslFunction(boost::bind(&MVll::getSigma6s, &(*this), _1));
                if (gsl_integration_cquad(&FS, q_min, q_max, 1.e-2, 1.e-1, w_sigma, &avaSigma, &errSigma, NULL) != 0) return std::numeric_limits<double>::quiet_NaN();
                cacheSigma7[qbin] = NN*avaSigma;
                sigma7Cached[qbin] = 1;
            }
            return cacheSigma7[qbin];
            break;
        case 9:
            if (sigma9Cached[qbin] == 0) {
                FS = convertToGslFunction(boost::bind(&MVll::getSigma7, &(*this), _1));
                if (gsl_integration_cquad(&FS, q_min, q_max, 1.e-2, 1.e-1, w_sigma, &avaSigma, &errSigma, NULL) != 0) return std::numeric_limits<double>::quiet_NaN();
                cacheSigma9[qbin] = NN*avaSigma;
                sigma9Cached[qbin] = 1;
            }
            return cacheSigma9[qbin];
            break;
        case 10:
            if (sigma10Cached[qbin] == 0) {
                FS = convertToGslFunction(boost::bind(&MVll::getSigma8, &(*this), _1));
                if (gsl_integration_cquad(&FS, q_min, q_max, 1.e-2, 1.e-1, w_sigma, &avaSigma, &errSigma, NULL) != 0) return std::numeric_limits<double>::quiet_NaN();
                cacheSigma10[qbin] = NN*avaSigma;
                sigma10Cached[qbin] = 1;
            }
            return cacheSigma10[qbin];
            break;
        case 11:
            if (sigma11Cached[qbin] == 0) {
                FS = convertToGslFunction(boost::bind(&MVll::getSigma9, &(*this), _1));
                if (gsl_integration_cquad(&FS, q_min, q_max, 1.e-2, 1.e-1, w_sigma, &avaSigma, &errSigma, NULL) != 0) return std::numeric_limits<double>::quiet_NaN();
                cacheSigma11[qbin] = NN*avaSigma;
                sigma11Cached[qbin] = 1;
            }
            return cacheSigma11[qbin];
            break;
        default:
            std::stringstream out;
            out << i;
            throw std::runtime_error("MVll::integrateSigma: index " + out.str() + " not implemented");
    }

    gsl_set_error_handler(old_handler);

}

double MVll::k2 ( double  q2 )

private

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 1416 of file MVll.cpp.

{
    return (MM4 + q2 * q2 + MV4 - twoMV2 * q2 - twoMM2 * (q2 + MV2)) / fourMM2;
}

double MVll::lambda ( double  q2 )

private

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

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \lambda \)

Definition at line 1431 of file MVll.cpp.

{
    return (MM4 + q2 * q2 + MV4 - twoMV2 * q2 - twoMM2 * (q2 + MV2));
}

double MVll::reDC9fit ( double *  x, double *  p  )

private

The fit function for the real part of the QCDF correction \( \Delta C_9^{\lambda} \).

Parameters

[in]	x	fit variable
[in]	p	fit parameters vector

Returns

\( f_{Re \Delta C_9^{\lambda}} \)

Definition at line 1074 of file MVll.cpp.

{
    return p[0]/x[0] + p[1] + p[2]*x[0] + p[3]*x[0]*x[0] + p[4]*x[0]*x[0]*x[0] + p[5]*x[0]*x[0]*x[0]*x[0] + p[6]*x[0]*x[0]*x[0]*x[0]*x[0]; 
}

double MVll::S_L ( double  q2 )

private

The helicity form factor \( S_L \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( S_L \)

Definition at line 819 of file MVll.cpp.

{
    return -sqrt(lambda(q2)) / twoMM_mbpms * A_0(q2);
}

double MVll::T_0t ( double  q2 )

private

The helicity form factor \( \tilde{T}_0 \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \tilde{T}_0 \)

Definition at line 804 of file MVll.cpp.

{
    return 2 * sqrt(q2) * MV / MM_MMpMV * FF_fit(q2, a_0T23, a_1T23, a_2T23, MRT23_2);
}

double MVll::T_1 ( double  q2 )

private

The transverse form factor \( T_1 \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( T_1 \)

Definition at line 779 of file MVll.cpp.

{
    return FF_fit(q2, a_0T1, a_1T1, a_2T1, MRT1_2);
}

double MVll::T_2 ( double  q2 )

private

The transverse form factor \( T_2 \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( T_2 \)

Definition at line 784 of file MVll.cpp.

{
    return FF_fit(q2, a_0T2_LCSR, a_1T2, a_2T2, MRT2_2);
}

double MVll::T_m ( double  q2 )

private

The helicity form factor \( T_- \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( T_- \)

Definition at line 814 of file MVll.cpp.

{
    return (MM2mMV2 * T_2(q2) + sqrt(lambda(q2)) * T_1(q2)) / twoMM2;
}

double MVll::T_p ( double  q2 )

private

The helicity form factor \( T_+ \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( T_+ \)

Definition at line 809 of file MVll.cpp.

{
    return (MM2mMV2 * T_2(q2) - sqrt(lambda(q2)) * T_1(q2)) / twoMM2;
}

gslpp::complex MVll::Tparminus ( double  u, double  q2  )

private

The \( T^{\parallel}_- \) function from [40] .

Parameters

[in]	u	dummy variable to be integrated out
[in]	q2	\(q^2\) of the decay

Returns

\( T^{\parallel}_- \)

Definition at line 877 of file MVll.cpp.

{
    double ubar = 1. - u;
    return mySM.getQuarks(QCD::DOWN).getCharge()*(8. * C_8Lh / (ubar + u * q2 / MM2)
            + sixMMoMb * H_c(ubar * MM2 + u * q2,mu_h*mu_h) * C_2Lh_bar);
}

gslpp::complex MVll::Tparplus ( double  u, double  q2  )

private

The \( T^{\parallel}_+ \) function from [40] .

Parameters

[in]	u	dummy variable to be integrated out
[in]	q2	\(q^2\) of the decay

Returns

\( T^{\parallel}_+ \)

Definition at line 864 of file MVll.cpp.

{
    Ee = (MM2 - q2) / twoMM;
    ubar = 1. - u;
    arg1 = (fourMc2 - gslpp::complex::i()*1.e-10)/ (ubar * MM2 + u * q2) - 1.;
    B01 = -2. * sqrt(arg1) * arctan(1. / sqrt(arg1));
    arg1 = (fourMc2 - gslpp::complex::i()*1.e-10)/ q2 - 1.;
    B00 = -2. * sqrt(arg1) * arctan(1. / sqrt(arg1));
    
    gslpp::complex tpar = twoMM / Ee / ubar * I1(u, q2) + (ubar * MM2 + u * q2) / Ee / Ee / ubar / ubar * (B01 - B00);
    return MM / Mb * mySM.getQuarks(QCD::UP).getCharge() * tpar*C_2Lh_bar;
}

gslpp::complex MVll::Tperpplus ( double  u, double  q2  )

private

The \( T^{\perp}_+ \) function from [40] .

Parameters

[in]	u	dummy variable to be integrated out
[in]	q2	\(q^2\) of the decay

Returns

\( T^{\perp}_+ \)

Definition at line 850 of file MVll.cpp.

{
    Ee = (MM2 - q2) / twoMM;
    ubar = 1. - u;
    arg1 = (fourMc2 - gslpp::complex::i()*1.e-10)/ (ubar * MM2 + u * q2) - 1.;
    B01 = -2. * sqrt(arg1) * arctan(1. / sqrt(arg1));
    arg1 = (fourMc2 - gslpp::complex::i()*1.e-10)/ q2 - 1.;
    B00 = -2. * sqrt(arg1) * arctan(1. / sqrt(arg1));
    
    gslpp::complex tperp = twoMM / Ee / ubar * I1(u, q2) + q2 / Ee / Ee / ubar / ubar * (B01 - B00);
    return m4downcharge * C_8Lh / (u + ubar * q2 / MM2) + MM / 2. / Mb *
            mySM.getQuarks(QCD::UP).getCharge() * tperp * C_2Lh_bar;
}

void MVll::updateParameters ( )

private

The update parameter method for MVll.

Definition at line 83 of file MVll.cpp.

{
    if (!mySM.getMyFlavour()->getUpdateFlag(meson, vectorM, lep)) return;


    GF = mySM.getGF();
    ale = mySM.getAle();
    Mlep = mySM.getLeptons(lep).getMass();
    MM = mySM.getMesons(meson).getMass();
    MV = mySM.getMesons(vectorM).getMass();
    mu_b = mySM.getMub();
    mu_h = sqrt(mu_b * .5); // From Beneke Neubert
    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();
    width = mySM.getMesons(meson).computeWidth();

    switch (vectorM) {
        case StandardModel::K_star:
            a_0V = mySM.geta_0V();
            a_1V = mySM.geta_1V();
            a_2V = mySM.geta_2V();
            MRV_2 = mySM.getMRV() * mySM.getMRV();

            a_0A0 = mySM.geta_0A0();
            a_1A0 = mySM.geta_1A0();
            a_2A0 = mySM.geta_2A0();
            MRA0_2 = mySM.getMRA0() * mySM.getMRA0();

            a_0A1 = mySM.geta_0A1();
            a_1A1 = mySM.geta_1A1();
            a_2A1 = mySM.geta_2A1();
            MRA1_2 = mySM.getMRA1() * mySM.getMRA1();

            a_0A12 = mySM.geta_0A12();
            a_1A12 = mySM.geta_1A12();
            a_2A12 = mySM.geta_2A12();
            MRA12_2 = mySM.getMRA12() * mySM.getMRA12();

            a_0T1 = mySM.geta_0T1();
            a_1T1 = mySM.geta_1T1();
            a_2T1 = mySM.geta_2T1();
            MRT1_2 = mySM.getMRT1() * mySM.getMRT1();

            a_0T2 = mySM.geta_0T2();
            a_1T2 = mySM.geta_1T2();
            a_2T2 = mySM.geta_2T2();
            MRT2_2 = mySM.getMRT2() * mySM.getMRT2();

            a_0T23 = mySM.geta_0T23();
            a_1T23 = mySM.geta_1T23();
            a_2T23 = mySM.geta_2T23();
            MRT23_2 = mySM.getMRT23() * mySM.getMRT23();

            b = 1;
            break;
        case StandardModel::PHI:
            a_0V = mySM.geta_0Vphi();
            a_1V = mySM.geta_1Vphi();
            a_2V = mySM.geta_2Vphi();
            MRV_2 = mySM.getMRVphi() * mySM.getMRVphi();

            a_0A0 = mySM.geta_0A0phi();
            a_1A0 = mySM.geta_1A0phi();
            a_2A0 = mySM.geta_2A0phi();
            MRA0_2 = mySM.getMRA0phi() * mySM.getMRA0phi();

            a_0A1 = mySM.geta_0A1phi();
            a_1A1 = mySM.geta_1A1phi();
            a_2A1 = mySM.geta_2A1phi();
            MRA1_2 = mySM.getMRA1phi() * mySM.getMRA1phi();

            a_0A12 = mySM.geta_0A12phi();
            a_1A12 = mySM.geta_1A12phi();
            a_2A12 = mySM.geta_2A12phi();
            MRA12_2 = mySM.getMRA12phi() * mySM.getMRA12phi();

            a_0T1 = mySM.geta_0T1phi();
            a_1T1 = mySM.geta_1T1phi();
            a_2T1 = mySM.geta_2T1phi();
            MRT1_2 = mySM.getMRT1phi() * mySM.getMRT1phi();

            a_0T2 = mySM.geta_0T2phi();
            a_1T2 = mySM.geta_1T2phi();
            a_2T2 = mySM.geta_2T2phi();
            MRT2_2 = mySM.getMRT2phi() * mySM.getMRT2phi();

            a_0T23 = mySM.geta_0T23phi();
            a_1T23 = mySM.geta_1T23phi();
            a_2T23 = mySM.geta_2T23phi();
            MRT23_2 = mySM.getMRT23phi() * mySM.getMRT23phi();

            b = 0.489;
            break;
        default:
            std::stringstream out;
            out << vectorM;
            throw std::runtime_error("MVll: vector " + out.str() + " not implemented");
    }


    h_0[0] = mySM.geth_0();
    h_0[1] = mySM.geth_p();
    h_0[2] = mySM.geth_m();

    h_1[0] = mySM.geth_0_1();
    h_1[1] = mySM.geth_p_1();
    h_1[2] = mySM.geth_m_1();

    h_2[0] = mySM.geth_0_2();
    h_2[1] = mySM.geth_p_2();
    h_2[2] = mySM.geth_m_2();

    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_1L_bar = (*(allcoeff[LO]))(0)/2.;
    C_2 = ((*(allcoeff[LO]))(1) + (*(allcoeff[NLO]))(1));
    C_2L_bar = (*(allcoeff[LO]))(1) - (*(allcoeff[LO]))(0)/6.;
    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_8L = (*(allcoeff[LO]))(7);
    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);
    
    allcoeffh = mySM.getMyFlavour()->ComputeCoeffBMll(mu_h); //check the mass scale, scheme fixed to NDR

    C_1Lh_bar = (*(allcoeffh[LO]))(0)/2.;
    C_2Lh_bar = (*(allcoeffh[LO]))(1) - (*(allcoeff[LO]))(0)/6.;
    C_8Lh = (*(allcoeffh[LO]))(7);
    
    checkCache();

    t_p = pow(MM + MV, 2.);
    t_m = pow(MM - MV, 2.);
    t_0 = t_p * (1. - sqrt(1. - t_m / t_p)); /*Modify it for Lattice*/
    z_0 = (sqrt(t_p) - sqrt(t_p - t_0)) / (sqrt(t_p) + sqrt(t_p - t_0));
    //else t_0 = 12.;
    MMpMV = MM + MV;
    MMpMV2 = MMpMV * MMpMV;
    MMmMV = MM - MV;
    MMmMV2 = MMmMV * MMmMV;
    MM2 = MM*MM;
    MM4 = MM2*MM2;
    MV2 = MV*MV;
    MV4 = MV2*MV2;
    MMMV = MM*MV;
    MM2mMV2 = MM2 - MV2;
    fourMV = 4. * MV;
    twoMM2 = 2. * MM2;
    twoMV2 = 2. * MV2;
    onepMMoMV = (1. + MV / MM);
    MM_MMpMV = MM * MMpMV;
    twoMM_mbpms = 2. * MM * (Mb + Ms);
    fourMM2 = 4. * MM2;
    Mlep2 = Mlep*Mlep;
    twoMlepMb = 2. * Mlep*Mb;
    MboMW = Mb / MW;
    MsoMb = Ms / Mb;
    ninetysixM_PI3MM3 = 96. * M_PI * M_PI * M_PI * MM * MM*MM;
    sixteenM_PI2 = 16. * M_PI*M_PI;
    sixteenM_PI2MM2 = sixteenM_PI2 * MM*MM;
    twoMboMM = 2 * Mb / MM;
    H_0_pre = 8. / 27. + 4. / 9. * gslpp::complex::i() * M_PI;
    H_0_WC = (C_3 + 4. / 3. * C_4 + 16. * C_5 + 64. / 3. * C_6);
    H_c_WC = (4. / 3. * C_1 + C_2 + 6. * C_3 + 60. * C_5);
    H_b_WC = (7. * C_3 + 4. / 3. * C_4 + 76. * C_5 + 64. / 3. * C_6);
    mu_b2 = mu_b*mu_b;
    Mc2 = Mc*Mc;
    Mb2 = Mb*Mb;
    fourMc2 = 4. * Mc2;
    fourMb2 = 4. * Mb2;
    logMc = log(Mc2 / mu_b2);
    logMb = log(Mb2 / mu_b2);
    fournineth = 4. / 9.;
    half = 1. / 2.;
    twothird = 2. / 3.;
    ihalfMPI = gslpp::complex::i() * M_PI / 2.;
    twoMM3 = 2. * MM2 * MM;
    C2_inv = 1. / (2. * C_2.real());
    gtilde_1_pre = -16. * pow(MM, 3.)*(MM + MV) * pow(M_PI, 2.);
    gtilde_2_pre = -16. * pow(MM, 3.) * pow(M_PI, 2.) / MMpMV;
    gtilde_3_pre = 64. * pow(MM, 3.) * pow(M_PI, 2.) * MV*MMpMV;
    S_L_pre = (-2. * MM * (Mb + Ms));
    
    M_PI2osix = M_PI * M_PI / 6.;
    twoMM = 2.*MM;
    m4downcharge = -4. * mySM.getQuarks(QCD::DOWN).getCharge();
    threeGegen0 = mySM.getMesons(vectorM).getGegenalpha(0)*3.;
    threeGegen1otwo = mySM.getMesons(vectorM).getGegenalpha(1)*3./2.;
    twoMc2 = 2.*Mc2;

    sixMMoMb = 6. * MM / Mb;
    CF =4./3.;
    
    deltaT_0 = mySM.Als(mu_b) * CF / 4. / M_PI;
    deltaT_1par = mySM.Als(mu_h) * CF / 4. * M_PI / 3. * mySM.getMesons(meson).getDecayconst() *
            mySM.getMesons(vectorM).getDecayconst() / MM; 
    deltaT_1perp = mySM.Als(mu_h) * CF / 4. * M_PI / 3. * mySM.getMesons(meson).getDecayconst() *
            mySM.getFKstarp() / MM; 
            
    F87_0=-32. / 9. * log(mu_b / Mb) + 8. / 27. * M_PI * M_PI - 44. / 9. - 8. / 9. * gslpp::complex::i() * M_PI;
    F87_1 = (4. / 3. * M_PI * M_PI - 40. / 3.);
    F87_2 = (32. / 9. * M_PI * M_PI - 316. / 9.);
    F87_3 = (200. / 27. * M_PI * M_PI - 658. / 9.);

    F89_0 = 104. / 9. - 32. / 27. * M_PI * M_PI ;
    F89_1 = 1184. / 27. - 40. / 9. * M_PI * M_PI;
    F89_2 = (-32. / 3. * M_PI * M_PI + 14212. / 135.);
    F89_3 = (-560. / 27. * M_PI * M_PI + 193444. / 945.);

    F29_0 = (256./243. - 32./81.*gslpp::complex::i()*M_PI - 128./9.*log(Mc/Mb))*log(mu_b/Mb) + 512./81.*log(mu_b/Mb)*log(mu_b/Mb) + 5.4082 - 1.0934 * gslpp::complex::i();
    F29_L1 = 32./81.*log(mu_b/Mb) + (0.48576 + 0.31119 * gslpp::complex::i());
    F29_1 = (-32./405. + 64./45./Mc2*Mb2)*log(mu_b/Mb) + (1.9061 + 0.80843 * gslpp::complex::i());
    F29_2 = (-8./945. + 16./105./Mc2*Mb2/Mc2*Mb2)*log(mu_b/Mb) + (-1.8286 + 2.8428 * gslpp::complex::i());
    F29_3 = (-32./25515. + 64./2835./Mc2*Mb2/Mc2*Mb2/Mc2*Mb2)*log(mu_b/Mb) + (-12.113 + 8.1251 * gslpp::complex::i());
    F29_L1_1 = (0.21951 - 0.14852 * gslpp::complex::i());
    F29_L1_2 = (0.13015 - 0.22155 * gslpp::complex::i());
    F29_L1_3 = (-0.079692 - 0.31214 * gslpp::complex::i());

    F27_0 = 416./81. *log(mu_b/Mb) + 3.8367 + 0.3531 * gslpp::complex::i();
    F27_1 = (1.3098 + 0.60185 * gslpp::complex::i());
    F27_2 = (0.13507 + 0.89014 * gslpp::complex::i());
    F27_3 = (-1.0271 + 0.77168 * gslpp::complex::i());
    F27_L1_1 = (-0.031936 - 0.10981 * gslpp::complex::i());
    F27_L1_2 = (-0.14169 - 0.035553 * gslpp::complex::i());
    F27_L1_3 = (-0.13592 + 0.093 * gslpp::complex::i()); 

    a_0A12_LCSR = a_0A0 * (MM2mMV2) / (8. * MMMV);
    a_0T2_LCSR = a_0T1;

    NN = ((4. * GF * MM * ale * lambda_t) / (sqrt(2.)*4. * M_PI)).abs2();
    
    if (mySM.getMyFlavour()->getUpdateFlag(meson, vectorM, lep)) {
        switch (lep) {
            case StandardModel::MU:
                fit_DeltaC9_p_mumu();
                fit_DeltaC9_m_mumu();
                fit_DeltaC9_0_mumu();
                break;
            case StandardModel::ELECTRON:
                fit_DeltaC9_p_ee();
                fit_DeltaC9_m_ee();
                fit_DeltaC9_0_ee();
                break;
            default:
                std::stringstream out;
                out << lep;
                throw std::runtime_error("MVll: lepton " + out.str() + " not implemented");
        }
    }

    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 (I1_updated == 0) for (it = sigma1Cached.begin(); it != sigma1Cached.end(); ++it) it->second = 0;
    if (I2_updated == 0) for (it = sigma2Cached.begin(); it != sigma2Cached.end(); ++it) it->second = 0;
    if (I3_updated == 0) for (it = sigma3Cached.begin(); it != sigma3Cached.end(); ++it) it->second = 0;
    if (I4_updated == 0) for (it = sigma4Cached.begin(); it != sigma4Cached.end(); ++it) it->second = 0;
    if (I5_updated == 0) for (it = sigma5Cached.begin(); it != sigma5Cached.end(); ++it) it->second = 0;
    if (I6_updated == 0) for (it = sigma6Cached.begin(); it != sigma6Cached.end(); ++it) it->second = 0;
    if (I7_updated == 0) for (it = sigma7Cached.begin(); it != sigma7Cached.end(); ++it) it->second = 0;
    if (I9_updated == 0) for (it = sigma9Cached.begin(); it != sigma9Cached.end(); ++it) it->second = 0;
    if (I10_updated == 0) for (it = sigma10Cached.begin(); it != sigma10Cached.end(); ++it) it->second = 0;
    if (I11_updated == 0) for (it = sigma11Cached.begin(); it != sigma11Cached.end(); ++it) it->second = 0;

    if (I0_updated == 0) for (it = delta0Cached.begin(); it != delta0Cached.end(); ++it) it->second = 0;
    if (I1_updated == 0) for (it = delta1Cached.begin(); it != delta1Cached.end(); ++it) it->second = 0;
    if (I2_updated == 0) for (it = delta2Cached.begin(); it != delta2Cached.end(); ++it) it->second = 0;
    if (I3_updated == 0) for (it = delta3Cached.begin(); it != delta3Cached.end(); ++it) it->second = 0;
    if (I11_updated == 0) for (it = delta11Cached.begin(); it != delta11Cached.end(); ++it) it->second = 0;
    
    std::map<double, unsigned int >::iterator iti;
    if (deltaTparpupdated == 0) for (iti = deltaTparpCached.begin(); iti != deltaTparpCached.end(); ++iti) iti->second = 0;
    if (deltaTparmupdated == 0) for (iti = deltaTparmCached.begin(); iti != deltaTparmCached.end(); ++iti) iti->second = 0;
    if (deltaTperpupdated == 0) for (iti = deltaTparpCached.begin(); iti != deltaTparpCached.end(); ++iti) iti->second = 0;
    
    if (deltaTparpupdated*deltaTparmupdated == 0) for (it = I1Cached.begin(); it != I1Cached.end(); ++it) it->second = 0;

    mySM.getMyFlavour()->setUpdateFlag(meson, vectorM, lep, false);
    
    return;
}

double MVll::V ( double  q2 )

private

The transverse form factor \( V \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( V \)

Definition at line 759 of file MVll.cpp.

{
    return FF_fit(q2, a_0V, a_1V, a_2V, MRV_2);
}

double MVll::V_0t ( double  q2 )

private

The helicity form factor \( \tilde{V}_0 \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( \tilde{V}_0 \)

Definition at line 789 of file MVll.cpp.

{
    return fourMV / sqrt(q2) * FF_fit(q2, a_0A12_LCSR, a_1A12, a_2A12, MRA12_2);
}

double MVll::V_m ( double  q2 )

private

The helicity form factor \( V_- \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( V_- \)

Definition at line 799 of file MVll.cpp.

{
    return half * (onepMMoMV * A_1(q2) + sqrt(lambda(q2)) / (MM_MMpMV) * V(q2));
}

double MVll::V_p ( double  q2 )

private

The helicity form factor \( V_+ \).

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( V_+ \)

Definition at line 794 of file MVll.cpp.

{
    return half * (onepMMoMV * A_1(q2) - sqrt(lambda(q2)) / (MM_MMpMV) * V(q2));
}

gslpp::complex MVll::Y ( double  q2 )

private

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 1368 of file MVll.cpp.

{
    return -half * H_0(q2) * H_0_WC + H_c(q2,mu_b2) * H_c_WC - half * H_b(q2) * H_b_WC;
}

double MVll::z ( double  q2 )

private

The lattice parameter \( z \) from arXiv:1310.3722v3.

Parameters

[in]

q2

\(q^2\) of the decay

Returns

\( z \)

Definition at line 754 of file MVll.cpp.

{
    return ( sqrt(t_p - q2) - sqrt(t_p - t_0)) / (sqrt(t_p - q2) + sqrt(t_p - t_0));
}

Member Data Documentation

gslpp::vector<double> MVll::A0_cache

private

Cache variable

Definition at line 704 of file MVll.h.

unsigned int MVll::A0_updated

private

Cache variable

Definition at line 703 of file MVll.h.

gslpp::vector<double> MVll::A1_cache

private

Cache variable

Definition at line 707 of file MVll.h.

unsigned int MVll::A1_updated

private

Cache variable

Definition at line 706 of file MVll.h.

double MVll::a_0A0

private

LCSR fit parameter

Definition at line 568 of file MVll.h.

double MVll::a_0A1

private

LCSR fit parameter

Definition at line 572 of file MVll.h.

double MVll::a_0A12

private

LCSR fit parameter

Definition at line 576 of file MVll.h.

double MVll::a_0A12_LCSR

private

LCSR fit parameter

Definition at line 561 of file MVll.h.

double MVll::a_0T1

private

LCSR fit parameter

Definition at line 580 of file MVll.h.

double MVll::a_0T2

private

LCSR fit parameter

Definition at line 584 of file MVll.h.

double MVll::a_0T23

private

LCSR fit parameter

Definition at line 588 of file MVll.h.

double MVll::a_0T2_LCSR

private

LCSR fit parameter

Definition at line 562 of file MVll.h.

double MVll::a_0V

private

LCSR fit parameter

Definition at line 564 of file MVll.h.

double MVll::a_1A0

private

LCSR fit parameter

Definition at line 569 of file MVll.h.

double MVll::a_1A1

private

LCSR fit parameter

Definition at line 573 of file MVll.h.

double MVll::a_1A12

private

LCSR fit parameter

Definition at line 577 of file MVll.h.

double MVll::a_1T1

private

LCSR fit parameter

Definition at line 581 of file MVll.h.

double MVll::a_1T2

private

LCSR fit parameter

Definition at line 585 of file MVll.h.

double MVll::a_1T23

private

LCSR fit parameter

Definition at line 589 of file MVll.h.

double MVll::a_1V

private

LCSR fit parameter

Definition at line 565 of file MVll.h.

double MVll::a_2A0

private

LCSR fit parameter

Definition at line 570 of file MVll.h.

double MVll::a_2A1

private

LCSR fit parameter

Definition at line 574 of file MVll.h.

double MVll::a_2A12

private

LCSR fit parameter

Definition at line 578 of file MVll.h.

double MVll::a_2T1

private

LCSR fit parameter

Definition at line 582 of file MVll.h.

double MVll::a_2T2

private

LCSR fit parameter

Definition at line 586 of file MVll.h.

double MVll::a_2T23

private

LCSR fit parameter

Definition at line 590 of file MVll.h.

double MVll::a_2V

private

LCSR fit parameter

Definition at line 566 of file MVll.h.

double MVll::ale

private

Alpha electromagnetic

Definition at line 444 of file MVll.h.

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

private

Vector that contains the Wilson coeffients

Definition at line 593 of file MVll.h.

gslpp::vector<gslpp::complex>** MVll::allcoeffh

private

Vector that contains the Wilson coeffients at scale \(\mu_h\)

Definition at line 594 of file MVll.h.

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

private

Vector that contains the primed Wilson coeffients

Definition at line 595 of file MVll.h.

gslpp::complex MVll::arg1

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 526 of file MVll.h.

double MVll::avaDelta

private

Gsl integral variable

Definition at line 658 of file MVll.h.

double MVll::avaDTPPR

private

Gsl integral variable

Definition at line 659 of file MVll.h.

double MVll::avaSigma

private

Gsl integral variable

Definition at line 657 of file MVll.h.

double MVll::b

private

BF of the decay V -> final states

Definition at line 456 of file MVll.h.

gslpp::complex MVll::B00

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 528 of file MVll.h.

gslpp::complex MVll::B01

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 527 of file MVll.h.

double MVll::beta_cache

private

Cache variable

Definition at line 723 of file MVll.h.

unsigned int MVll::beta_updated

private

Cache variable

Definition at line 722 of file MVll.h.

double MVll::C2_inv

private

Cache variable

Definition at line 515 of file MVll.h.

gslpp::complex MVll::C_1

private

Wilson coeffients \(C_1\)

Definition at line 597 of file MVll.h.

gslpp::complex MVll::C_10

private

Wilson coeffients \(C_{10}\)

Definition at line 611 of file MVll.h.

gslpp::complex MVll::C_10_cache

private

Cache variable

Definition at line 781 of file MVll.h.

unsigned int MVll::C_10_updated

private

Cache variable

Definition at line 780 of file MVll.h.

gslpp::complex MVll::C_10p

private

Wilson coeffients \(C_{10}'\)

Definition at line 617 of file MVll.h.

gslpp::complex MVll::C_10p_cache

private

Cache variable

Definition at line 790 of file MVll.h.

unsigned int MVll::C_10p_updated

private

Cache variable

Definition at line 789 of file MVll.h.

gslpp::complex MVll::C_1_cache

private

Cache variable

Definition at line 757 of file MVll.h.

unsigned int MVll::C_1_updated

private

Cache variable

Definition at line 756 of file MVll.h.

gslpp::complex MVll::C_1L_bar

private

Wilson coeffients \(C_1\)

Definition at line 598 of file MVll.h.

gslpp::complex MVll::C_1Lh_bar

private

Wilson coeffients \(C_1\)

Definition at line 599 of file MVll.h.

gslpp::complex MVll::C_2

private

Wilson coeffients \(C_2\)

Definition at line 600 of file MVll.h.

gslpp::complex MVll::C_2_cache

private

Cache variable

Definition at line 760 of file MVll.h.

unsigned int MVll::C_2_updated

private

Cache variable

Definition at line 759 of file MVll.h.

gslpp::complex MVll::C_2L_bar

private

Leading order Wilson coeffients \(C_2\)

Definition at line 601 of file MVll.h.

gslpp::complex MVll::C_2Lh_bar

private

Leading order Wilson coeffients \(C_2\) at scale \(\mu_h\)

Definition at line 602 of file MVll.h.

gslpp::complex MVll::C_2Lh_cache

private

Cache variable

Definition at line 805 of file MVll.h.

unsigned int MVll::C_2Lh_updated

private

Cache variable

Definition at line 804 of file MVll.h.

gslpp::complex MVll::C_3

private

Wilson coeffients \(C_3\)

Definition at line 603 of file MVll.h.

gslpp::complex MVll::C_3_cache

private

Cache variable

Definition at line 763 of file MVll.h.

unsigned int MVll::C_3_updated

private

Cache variable

Definition at line 762 of file MVll.h.

gslpp::complex MVll::C_4

private

Wilson coeffients \(C_4\)

Definition at line 604 of file MVll.h.

gslpp::complex MVll::C_4_cache

private

Cache variable

Definition at line 766 of file MVll.h.

unsigned int MVll::C_4_updated

private

Cache variable

Definition at line 765 of file MVll.h.

gslpp::complex MVll::C_5

private

Wilson coeffients \(C_5\)

Definition at line 605 of file MVll.h.

gslpp::complex MVll::C_5_cache

private

Cache variable

Definition at line 769 of file MVll.h.

unsigned int MVll::C_5_updated

private

Cache variable

Definition at line 768 of file MVll.h.

gslpp::complex MVll::C_6

private

Wilson coeffients \(C_6\)

Definition at line 606 of file MVll.h.

gslpp::complex MVll::C_6_cache

private

Cache variable

Definition at line 772 of file MVll.h.

unsigned int MVll::C_6_updated

private

Cache variable

Definition at line 771 of file MVll.h.

gslpp::complex MVll::C_7

private

Wilson coeffients \(C_7\)

Definition at line 607 of file MVll.h.

gslpp::complex MVll::C_7_cache

private

Cache variable

Definition at line 775 of file MVll.h.

unsigned int MVll::C_7_updated

private

Cache variable

Definition at line 774 of file MVll.h.

gslpp::complex MVll::C_7p

private

Wilson coeffients \(C_7'\)

Definition at line 615 of file MVll.h.

gslpp::complex MVll::C_7p_cache

private

Cache variable

Definition at line 784 of file MVll.h.

unsigned int MVll::C_7p_updated

private

Cache variable

Definition at line 783 of file MVll.h.

gslpp::complex MVll::C_8L

private

Leading order Wilson coeffients \(C_8\)

Definition at line 608 of file MVll.h.

gslpp::complex MVll::C_8Lh

private

Leading order Wilson coeffients \(C_8\) at scale \(\mu_h\)

Definition at line 609 of file MVll.h.

gslpp::complex MVll::C_8Lh_cache

private

Cache variable

Definition at line 808 of file MVll.h.

unsigned int MVll::C_8Lh_updated

private

Cache variable

Definition at line 807 of file MVll.h.

gslpp::complex MVll::C_9

private

Wilson coeffients \(C_9\)

Definition at line 610 of file MVll.h.

gslpp::complex MVll::C_9_cache

private

Cache variable

Definition at line 778 of file MVll.h.

unsigned int MVll::C_9_updated

private

Cache variable

Definition at line 777 of file MVll.h.

gslpp::complex MVll::C_9p

private

Wilson coeffients \(C_9'\)

Definition at line 616 of file MVll.h.

gslpp::complex MVll::C_9p_cache

private

Cache variable

Definition at line 787 of file MVll.h.

unsigned int MVll::C_9p_updated

private

Cache variable

Definition at line 786 of file MVll.h.

gslpp::complex MVll::C_P

private

Wilson coeffients \(C_P\)

Definition at line 613 of file MVll.h.

gslpp::complex MVll::C_P_cache

private

Cache variable

Definition at line 796 of file MVll.h.

unsigned int MVll::C_P_updated

private

Cache variable

Definition at line 795 of file MVll.h.

gslpp::complex MVll::C_Pp

private

Wilson coeffients \(C_P'\)

Definition at line 619 of file MVll.h.

gslpp::complex MVll::C_Pp_cache

private

Cache variable

Definition at line 802 of file MVll.h.

unsigned int MVll::C_Pp_updated

private

Cache variable

Definition at line 801 of file MVll.h.

gslpp::complex MVll::C_S

private

Wilson coeffients \(C_S\)

Definition at line 612 of file MVll.h.

gslpp::complex MVll::C_S_cache

private

Cache variable

Definition at line 793 of file MVll.h.

unsigned int MVll::C_S_updated

private

Cache variable

Definition at line 792 of file MVll.h.

gslpp::complex MVll::C_Sp

private

Wilson coeffients \(C_S'\)

Definition at line 618 of file MVll.h.

gslpp::complex MVll::C_Sp_cache

private

Cache variable

Definition at line 799 of file MVll.h.

unsigned int MVll::C_Sp_updated

private

Cache variable

Definition at line 798 of file MVll.h.

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

private

Cache variable

Definition at line 689 of file MVll.h.

std::map<std::pair<double, double>, double > MVll::cacheDelta1

private

Cache variable

Definition at line 690 of file MVll.h.

std::map<std::pair<double, double>, double > MVll::cacheDelta11

private

Cache variable

Definition at line 694 of file MVll.h.

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

private

Cache variable

Definition at line 691 of file MVll.h.

std::map<std::pair<double, double>, double > MVll::cacheDelta3

private

Cache variable

Definition at line 692 of file MVll.h.

std::map<std::pair<double, double>, double > MVll::cacheDelta7

private

Cache variable

Definition at line 693 of file MVll.h.

std::map<double, gslpp::complex> MVll::cacheDeltaTparm

private

Cache variable

Definition at line 879 of file MVll.h.

std::map<double, gslpp::complex> MVll::cacheDeltaTparp

private

Cache variable

Definition at line 878 of file MVll.h.

std::map<double, gslpp::complex> MVll::cacheDeltaTperp

private

Cache variable

Definition at line 880 of file MVll.h.

std::map<std::pair<double, double>, gslpp::complex > MVll::cacheI1

private

Cache variable

Definition at line 675 of file MVll.h.

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

private

Cache variable

Definition at line 677 of file MVll.h.

std::map<std::pair<double, double>, double > MVll::cacheSigma1

private

Cache variable

Definition at line 678 of file MVll.h.

std::map<std::pair<double, double>, double > MVll::cacheSigma10

private

Cache variable

Definition at line 686 of file MVll.h.

std::map<std::pair<double, double>, double > MVll::cacheSigma11

private

Cache variable

Definition at line 687 of file MVll.h.

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

private

Cache variable

Definition at line 679 of file MVll.h.

std::map<std::pair<double, double>, double > MVll::cacheSigma3

private

Cache variable

Definition at line 680 of file MVll.h.

std::map<std::pair<double, double>, double > MVll::cacheSigma4

private

Cache variable

Definition at line 681 of file MVll.h.

std::map<std::pair<double, double>, double > MVll::cacheSigma5

private

Cache variable

Definition at line 682 of file MVll.h.

std::map<std::pair<double, double>, double > MVll::cacheSigma6

private

Cache variable

Definition at line 683 of file MVll.h.

std::map<std::pair<double, double>, double > MVll::cacheSigma7

private

Cache variable

Definition at line 684 of file MVll.h.

std::map<std::pair<double, double>, double > MVll::cacheSigma9

private

Cache variable

Definition at line 685 of file MVll.h.

double MVll::CF

private

Cache variable

Definition at line 519 of file MVll.h.

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

private

Cache variable

Definition at line 867 of file MVll.h.

std::map<std::pair<double, double>, unsigned int > MVll::delta11Cached

private

Cache variable

Definition at line 872 of file MVll.h.

std::map<std::pair<double, double>, unsigned int > MVll::delta1Cached

private

Cache variable

Definition at line 868 of file MVll.h.

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

private

Cache variable

Definition at line 869 of file MVll.h.

std::map<std::pair<double, double>, unsigned int > MVll::delta3Cached

private

Cache variable

Definition at line 870 of file MVll.h.

std::map<std::pair<double, double>, unsigned int > MVll::delta7Cached

private

Cache variable

Definition at line 871 of file MVll.h.

double MVll::deltaT_0

private

Cache variable

Definition at line 520 of file MVll.h.

double MVll::deltaT_1par

private

Cache variable

Definition at line 521 of file MVll.h.

double MVll::deltaT_1perp

private

Cache variable

Definition at line 522 of file MVll.h.

std::map<double, unsigned int> MVll::deltaTparmCached

private

Cache variable

Definition at line 875 of file MVll.h.

unsigned int MVll::deltaTparmupdated

private

Cache variable

Definition at line 883 of file MVll.h.

std::map<double, unsigned int> MVll::deltaTparpCached

private

Cache variable

Definition at line 874 of file MVll.h.

unsigned int MVll::deltaTparpupdated

private

Cache variable

Definition at line 882 of file MVll.h.

std::map<double, unsigned int> MVll::deltaTperpCached

private

Cache variable

Definition at line 876 of file MVll.h.

unsigned int MVll::deltaTperpupdated

private

Cache variable

Definition at line 884 of file MVll.h.

gsl_function MVll::DTPPR

private

Gsl integral variable

Definition at line 667 of file MVll.h.

double MVll::Ee

private

Cache variable

Definition at line 523 of file MVll.h.

double MVll::errDelta

private

Gsl integral variable

Definition at line 662 of file MVll.h.

double MVll::errDTPPR

private

Gsl integral variable

Definition at line 663 of file MVll.h.

double MVll::errSigma

private

Gsl integral variable

Definition at line 661 of file MVll.h.

gslpp::complex MVll::F27_0

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 549 of file MVll.h.

gslpp::complex MVll::F27_1

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 550 of file MVll.h.

gslpp::complex MVll::F27_2

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 551 of file MVll.h.

gslpp::complex MVll::F27_3

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 552 of file MVll.h.

gslpp::complex MVll::F27_L1_1

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 553 of file MVll.h.

gslpp::complex MVll::F27_L1_2

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 554 of file MVll.h.

gslpp::complex MVll::F27_L1_3

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 555 of file MVll.h.

gslpp::complex MVll::F29_0

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 541 of file MVll.h.

gslpp::complex MVll::F29_1

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 543 of file MVll.h.

gslpp::complex MVll::F29_2

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 544 of file MVll.h.

gslpp::complex MVll::F29_3

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 545 of file MVll.h.

gslpp::complex MVll::F29_L1

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 542 of file MVll.h.

gslpp::complex MVll::F29_L1_1

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 546 of file MVll.h.

gslpp::complex MVll::F29_L1_2

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 547 of file MVll.h.

gslpp::complex MVll::F29_L1_3

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 548 of file MVll.h.

gslpp::complex MVll::F87_0

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 537 of file MVll.h.

gslpp::complex MVll::F87_1

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 538 of file MVll.h.

gslpp::complex MVll::F87_2

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 539 of file MVll.h.

gslpp::complex MVll::F87_3

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 540 of file MVll.h.

double MVll::F89_0

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 556 of file MVll.h.

double MVll::F89_1

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 557 of file MVll.h.

double MVll::F89_2

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 558 of file MVll.h.

double MVll::F89_3

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 559 of file MVll.h.

unsigned int MVll::F_updated

private

Cache variable

Definition at line 725 of file MVll.h.

gsl_function MVll::FD

private

Gsl integral variable

Definition at line 666 of file MVll.h.

double MVll::fourMb2

private

Cache variable

Definition at line 501 of file MVll.h.

double MVll::fourMc2

private

Cache variable

Definition at line 500 of file MVll.h.

double MVll::fourMM2

private

Cache variable

Definition at line 481 of file MVll.h.

double MVll::fourMV

private

Cache variable

Definition at line 475 of file MVll.h.

double MVll::fournineth

private

Cache variable

Definition at line 507 of file MVll.h.

gsl_function MVll::FS

private

Gsl integral variable

Definition at line 665 of file MVll.h.

double MVll::GF

private

Fermi constant

Definition at line 443 of file MVll.h.

TGraph MVll::gr1

private

Tgraph to be used for fitting the QCDF \(\Delta C_9\)

Definition at line 649 of file MVll.h.

TGraph MVll::gr2

private

Tgraph to be used for fitting the QCDF \(\Delta C_9\)

Definition at line 650 of file MVll.h.

double MVll::gtilde_1_pre

private

Cache variable

Definition at line 512 of file MVll.h.

double MVll::gtilde_2_pre

private

Cache variable

Definition at line 513 of file MVll.h.

double MVll::gtilde_3_pre

private

Cache variable

Definition at line 514 of file MVll.h.

unsigned int MVll::h0_updated

private

Cache variable

Definition at line 817 of file MVll.h.

gslpp::complex MVll::h0Ccache[3]

private

Cache variable

Definition at line 813 of file MVll.h.

unsigned int MVll::h1_updated

private

Cache variable

Definition at line 818 of file MVll.h.

gslpp::complex MVll::h1Ccache[3]

private

Cache variable

Definition at line 814 of file MVll.h.

unsigned int MVll::h2_updated

private

Cache variable

Definition at line 819 of file MVll.h.

gslpp::complex MVll::h2Ccache[3]

private

Cache variable

Definition at line 815 of file MVll.h.

gslpp::complex MVll::h_0[3]

private

Parameter that contains the contribution from the hadronic hamiltonian

Definition at line 457 of file MVll.h.

gslpp::complex MVll::H_0_pre

private

Cache variable

Definition at line 496 of file MVll.h.

gslpp::complex MVll::H_0_WC

private

Cache variable

Definition at line 504 of file MVll.h.

gslpp::complex MVll::h_1[3]

private

Parameter that contains the contribution from the hadronic hamiltonian

Definition at line 458 of file MVll.h.

gslpp::complex MVll::h_2[3]

private

Parameter that contains the contribution from the hadronic hamiltonian

Definition at line 459 of file MVll.h.

unsigned int MVll::H_A0updated

private

Cache variable

Definition at line 830 of file MVll.h.

unsigned int MVll::H_A1updated

private

Cache variable

Definition at line 831 of file MVll.h.

unsigned int MVll::H_A2updated

private

Cache variable

Definition at line 832 of file MVll.h.

gslpp::complex MVll::H_b_WC

private

Cache variable

Definition at line 506 of file MVll.h.

gslpp::complex MVll::H_c_WC

private

Cache variable

Definition at line 505 of file MVll.h.

gslpp::vector<double> MVll::H_Pcache

private

Cache variable

Definition at line 838 of file MVll.h.

unsigned int MVll::H_Pupdated

private

Cache variable

Definition at line 837 of file MVll.h.

gslpp::vector<double> MVll::H_Scache

private

Cache variable

Definition at line 835 of file MVll.h.

unsigned int MVll::H_Supdated

private

Cache variable

Definition at line 834 of file MVll.h.

gslpp::vector<double> MVll::H_V0cache

private

Cache variable

Definition at line 822 of file MVll.h.

unsigned int MVll::H_V0updated

private

Cache variable

Definition at line 821 of file MVll.h.

gslpp::vector<double> MVll::H_V1cache

private

Cache variable

Definition at line 825 of file MVll.h.

unsigned int MVll::H_V1updated

private

Cache variable

Definition at line 824 of file MVll.h.

gslpp::vector<double> MVll::H_V2cache

private

Cache variable

Definition at line 828 of file MVll.h.

unsigned int MVll::H_V2updated

private

Cache variable

Definition at line 827 of file MVll.h.

double MVll::half

private

Cache variable

Definition at line 508 of file MVll.h.

unsigned int MVll::I0_updated

private

Cache variable

Definition at line 840 of file MVll.h.

unsigned int MVll::I10_updated

private

Cache variable

Definition at line 850 of file MVll.h.

unsigned int MVll::I11_updated

private

Cache variable

Definition at line 851 of file MVll.h.

unsigned int MVll::I1_updated

private

Cache variable

Definition at line 841 of file MVll.h.

std::map<std::pair<double, double>, unsigned int > MVll::I1Cached

private

Cache variable

Definition at line 853 of file MVll.h.

unsigned int MVll::I2_updated

private

Cache variable

Definition at line 842 of file MVll.h.

unsigned int MVll::I3_updated

private

Cache variable

Definition at line 843 of file MVll.h.

unsigned int MVll::I4_updated

private

Cache variable

Definition at line 844 of file MVll.h.

unsigned int MVll::I5_updated

private

Cache variable

Definition at line 845 of file MVll.h.

unsigned int MVll::I6_updated

private

Cache variable

Definition at line 846 of file MVll.h.

unsigned int MVll::I7_updated

private

Cache variable

Definition at line 847 of file MVll.h.

unsigned int MVll::I8_updated

private

Cache variable

Definition at line 848 of file MVll.h.

unsigned int MVll::I9_updated

private

Cache variable

Definition at line 849 of file MVll.h.

gslpp::complex MVll::ihalfMPI

private

Cache variable

Definition at line 510 of file MVll.h.

std::vector<double> MVll::ImDeltaC9_0_ee

private

Vector that samples the QCDF \(\Delta C_9\)

Definition at line 632 of file MVll.h.

std::vector<double> MVll::ImDeltaC9_0_mumu

private

Vector that samples the QCDF \(\Delta C_9\)

Definition at line 626 of file MVll.h.

std::vector<double> MVll::ImDeltaC9_m_ee

private

Vector that samples the QCDF \(\Delta C_9\)

Definition at line 630 of file MVll.h.

std::vector<double> MVll::ImDeltaC9_m_mumu

private

Vector that samples the QCDF \(\Delta C_9\)

Definition at line 624 of file MVll.h.

std::vector<double> MVll::ImDeltaC9_p_ee

private

Vector that samples the QCDF \(\Delta C_9\)

Definition at line 628 of file MVll.h.

std::vector<double> MVll::ImDeltaC9_p_mumu

private

Vector that samples the QCDF \(\Delta C_9\)

Definition at line 622 of file MVll.h.

TF1 MVll::imffit

private

TF1 to be used for fitting the QCDF \(\Delta C_9\)

Definition at line 653 of file MVll.h.

TFitResultPtr MVll::imfres_0_ee

private

Vector that contains the fitted QCDF \(\Delta C_9\)

Definition at line 647 of file MVll.h.

TFitResultPtr MVll::imfres_0_mumu

private

Vector that contains the fitted QCDF \(\Delta C_9\)

Definition at line 640 of file MVll.h.

TFitResultPtr MVll::imfres_m_ee

private

Vector that contains the fitted QCDF \(\Delta C_9\)

Definition at line 645 of file MVll.h.

TFitResultPtr MVll::imfres_m_mumu

private

Vector that contains the fitted QCDF \(\Delta C_9\)

Definition at line 638 of file MVll.h.

TFitResultPtr MVll::imfres_p_ee

private

Vector that contains the fitted QCDF \(\Delta C_9\)

Definition at line 643 of file MVll.h.

TFitResultPtr MVll::imfres_p_mumu

private

Vector that contains the fitted QCDF \(\Delta C_9\)

Definition at line 636 of file MVll.h.

gslpp::vector<double> MVll::k2_cache

private

Cache variable

Definition at line 716 of file MVll.h.

unsigned int MVll::k2_updated

private

Cache variable

Definition at line 715 of file MVll.h.

gslpp::complex MVll::L1xm

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 534 of file MVll.h.

gslpp::complex MVll::L1xp

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 533 of file MVll.h.

gslpp::complex MVll::L1ym

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 536 of file MVll.h.

gslpp::complex MVll::L1yp

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 535 of file MVll.h.

gslpp::complex MVll::lambda_t

private

Vckm factor

Definition at line 455 of file MVll.h.

unsigned int MVll::lambda_updated

private

Cache variable

Definition at line 720 of file MVll.h.

StandardModel::lepton MVll::lep

private

Final leptons type

Definition at line 439 of file MVll.h.

double MVll::logMb

private

Cache variable

Definition at line 503 of file MVll.h.

double MVll::logMc

private

Cache variable

Definition at line 502 of file MVll.h.

double MVll::m4downcharge

private

Cache variable

Definition at line 488 of file MVll.h.

double MVll::M_PI2osix

private

Cache variable

Definition at line 486 of file MVll.h.

double MVll::Mb

private

b quark mass

Definition at line 448 of file MVll.h.

double MVll::Mb2

private

Cache variable

Definition at line 499 of file MVll.h.

unsigned int MVll::Mb_Ms_updated

private

Cache variable

Definition at line 749 of file MVll.h.

double MVll::MboMW

private

Cache variable

Definition at line 484 of file MVll.h.

double MVll::Mc

private

c quark mass

Definition at line 451 of file MVll.h.

double MVll::Mc2

private

Cache variable

Definition at line 498 of file MVll.h.

StandardModel::meson MVll::meson

private

Initial meson type

Definition at line 440 of file MVll.h.

double MVll::Mlep

private

Muon mass

Definition at line 445 of file MVll.h.

double MVll::Mlep2

private

Cache variable

Definition at line 482 of file MVll.h.

double MVll::MM

private

Initial meson mass

Definition at line 446 of file MVll.h.

double MVll::MM2

private

Cache variable

Definition at line 469 of file MVll.h.

double MVll::MM2mMV2

private

Cache variable

Definition at line 474 of file MVll.h.

double MVll::MM4

private

Cache variable

Definition at line 470 of file MVll.h.

double MVll::MM_MMpMV

private

Cache variable

Definition at line 477 of file MVll.h.

double MVll::MMmMV

private

Cache variable

Definition at line 467 of file MVll.h.

double MVll::MMmMV2

private

Cache variable

Definition at line 468 of file MVll.h.

double MVll::MMMV

private

Cache variable

Definition at line 473 of file MVll.h.

double MVll::MMpMV

private

Cache variable

Definition at line 465 of file MVll.h.

double MVll::MMpMV2

private

Cache variable

Definition at line 466 of file MVll.h.

double MVll::MRA0_2

private

LCSR fit parameter

Definition at line 571 of file MVll.h.

double MVll::MRA12_2

private

LCSR fit parameter

Definition at line 579 of file MVll.h.

double MVll::MRA1_2

private

LCSR fit parameter

Definition at line 575 of file MVll.h.

double MVll::MRT1_2

private

LCSR fit parameter

Definition at line 583 of file MVll.h.

double MVll::MRT23_2

private

LCSR fit parameter

Definition at line 591 of file MVll.h.

double MVll::MRT2_2

private

LCSR fit parameter

Definition at line 587 of file MVll.h.

double MVll::MRV_2

private

LCSR fit parameter

Definition at line 567 of file MVll.h.

double MVll::Ms

private

s quark mass

Definition at line 452 of file MVll.h.

double MVll::MsoMb

private

Cache variable

Definition at line 485 of file MVll.h.

double MVll::mu_b

private

b mass scale

Definition at line 449 of file MVll.h.

double MVll::mu_b2

private

Cache variable

Definition at line 497 of file MVll.h.

double MVll::mu_h

private

\(\sqrt{\mu_b*\lambda_{QCD}}\)

Definition at line 450 of file MVll.h.

double MVll::MV

private

Final vector meson mass

Definition at line 447 of file MVll.h.

double MVll::MV2

private

Cache variable

Definition at line 471 of file MVll.h.

double MVll::MV4

private

Cache variable

Definition at line 472 of file MVll.h.

double MVll::MW

private

W boson mass

Definition at line 454 of file MVll.h.

std::vector<double> MVll::myq2

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 633 of file MVll.h.

const StandardModel& MVll::mySM

private

Model type

Definition at line 438 of file MVll.h.

gslpp::vector<double> MVll::N_cache

private

Cache variable

Definition at line 697 of file MVll.h.

unsigned int MVll::N_updated

private

Cache variable

Definition at line 696 of file MVll.h.

gslpp::complex MVll::Nc_cache

private

Cache variable

Definition at line 698 of file MVll.h.

double MVll::ninetysixM_PI3MM3

private

Cache variable

Definition at line 492 of file MVll.h.

double MVll::NN

private

Cache variable

Definition at line 517 of file MVll.h.

gsl_error_handler_t* MVll::old_handler

private

GSL error handler store

Definition at line 673 of file MVll.h.

double MVll::onepMMoMV

private

Cache variable

Definition at line 476 of file MVll.h.

std::vector<double> MVll::ReDeltaC9_0_ee

private

Vector that samples the QCDF \(\Delta C_9\)

Definition at line 631 of file MVll.h.

std::vector<double> MVll::ReDeltaC9_0_mumu

private

Vector that samples the QCDF \(\Delta C_9\)

Definition at line 625 of file MVll.h.

std::vector<double> MVll::ReDeltaC9_m_ee

private

Vector that samples the QCDF \(\Delta C_9\)

Definition at line 629 of file MVll.h.

std::vector<double> MVll::ReDeltaC9_m_mumu

private

Vector that samples the QCDF \(\Delta C_9\)

Definition at line 623 of file MVll.h.

std::vector<double> MVll::ReDeltaC9_p_ee

private

Vector that samples the QCDF \(\Delta C_9\)

Definition at line 627 of file MVll.h.

std::vector<double> MVll::ReDeltaC9_p_mumu

private

Vector that samples the QCDF \(\Delta C_9\)

Definition at line 621 of file MVll.h.

TF1 MVll::reffit

private

TF1 to be used for fitting the QCDF \(\Delta C_9\)

Definition at line 652 of file MVll.h.

TFitResultPtr MVll::refres_0_ee

private

Vector that contains the fitted QCDF \(\Delta C_9\)

Definition at line 646 of file MVll.h.

TFitResultPtr MVll::refres_0_mumu

private

Vector that contains the fitted QCDF \(\Delta C_9\)

Definition at line 639 of file MVll.h.

TFitResultPtr MVll::refres_m_ee

private

Vector that contains the fitted QCDF \(\Delta C_9\)

Definition at line 644 of file MVll.h.

TFitResultPtr MVll::refres_m_mumu

private

Vector that contains the fitted QCDF \(\Delta C_9\)

Definition at line 637 of file MVll.h.

TFitResultPtr MVll::refres_p_ee

private

Vector that contains the fitted QCDF \(\Delta C_9\)

Definition at line 642 of file MVll.h.

TFitResultPtr MVll::refres_p_mumu

private

Vector that contains the fitted QCDF \(\Delta C_9\)

Definition at line 635 of file MVll.h.

double MVll::S_L_pre

private

Cache variable

Definition at line 516 of file MVll.h.

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

private

Cache variable

Definition at line 855 of file MVll.h.

std::map<std::pair<double, double>, unsigned int > MVll::sigma10Cached

private

Cache variable

Definition at line 864 of file MVll.h.

std::map<std::pair<double, double>, unsigned int > MVll::sigma11Cached

private

Cache variable

Definition at line 865 of file MVll.h.

std::map<std::pair<double, double>, unsigned int > MVll::sigma1Cached

private

Cache variable

Definition at line 856 of file MVll.h.

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

private

Cache variable

Definition at line 857 of file MVll.h.

std::map<std::pair<double, double>, unsigned int > MVll::sigma3Cached

private

Cache variable

Definition at line 858 of file MVll.h.

std::map<std::pair<double, double>, unsigned int > MVll::sigma4Cached

private

Cache variable

Definition at line 859 of file MVll.h.

std::map<std::pair<double, double>, unsigned int > MVll::sigma5Cached

private

Cache variable

Definition at line 860 of file MVll.h.

std::map<std::pair<double, double>, unsigned int > MVll::sigma6Cached

private

Cache variable

Definition at line 861 of file MVll.h.

std::map<std::pair<double, double>, unsigned int > MVll::sigma7Cached

private

Cache variable

Definition at line 862 of file MVll.h.

std::map<std::pair<double, double>, unsigned int > MVll::sigma9Cached

private

Cache variable

Definition at line 863 of file MVll.h.

double MVll::sixMMoMb

private

Cache variable

Definition at line 518 of file MVll.h.

double MVll::sixteenM_PI2

private

Cache variable

Definition at line 493 of file MVll.h.

double MVll::sixteenM_PI2MM2

private

Cache variable

Definition at line 494 of file MVll.h.

gslpp::vector<double> MVll::SL_cache

private

Cache variable

Definition at line 752 of file MVll.h.

unsigned int MVll::SL_updated

private

Cache variable

Definition at line 751 of file MVll.h.

unsigned int MVll::SR_updated

private

Cache variable

Definition at line 754 of file MVll.h.

gslpp::vector<double> MVll::T1_cache

private

Cache variable

Definition at line 710 of file MVll.h.

unsigned int MVll::T1_updated

private

Cache variable

Definition at line 709 of file MVll.h.

gslpp::vector<double> MVll::T2_cache

private

Cache variable

Definition at line 713 of file MVll.h.

unsigned int MVll::T2_updated

private

Cache variable

Definition at line 712 of file MVll.h.

double MVll::t_0

private

Cache variable

Definition at line 463 of file MVll.h.

gslpp::vector<double> MVll::T_cache

private

Cache variable

Definition at line 887 of file MVll.h.

double MVll::t_m

private

Cache variable

Definition at line 462 of file MVll.h.

double MVll::t_p

private

Cache variable

Definition at line 461 of file MVll.h.

unsigned int MVll::T_updated

private

Cache variable

Definition at line 886 of file MVll.h.

double MVll::threeGegen0

private

Cache variable

Definition at line 489 of file MVll.h.

double MVll::threeGegen1otwo

private

Cache variable

Definition at line 490 of file MVll.h.

gslpp::vector<double> MVll::TL0_cache

private

Cache variable

Definition at line 743 of file MVll.h.

unsigned int MVll::TL0_updated

private

Cache variable

Definition at line 742 of file MVll.h.

unsigned int MVll::TL1_updated

private

Cache variable

Definition at line 730 of file MVll.h.

unsigned int MVll::TL2_updated

private

Cache variable

Definition at line 731 of file MVll.h.

double MVll::tmpq2

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 655 of file MVll.h.

unsigned int MVll::TR0_updated

private

Cache variable

Definition at line 747 of file MVll.h.

unsigned int MVll::TR1_updated

private

Cache variable

Definition at line 736 of file MVll.h.

unsigned int MVll::TR2_updated

private

Cache variable

Definition at line 737 of file MVll.h.

double MVll::twoMboMM

private

Cache variable

Definition at line 495 of file MVll.h.

double MVll::twoMc2

private

Cache variable

Definition at line 491 of file MVll.h.

double MVll::twoMlepMb

private

Cache variable

Definition at line 483 of file MVll.h.

double MVll::twoMM

private

Cache variable

Definition at line 487 of file MVll.h.

double MVll::twoMM2

private

Cache variable

Definition at line 478 of file MVll.h.

double MVll::twoMM3

private

Cache variable

Definition at line 511 of file MVll.h.

double MVll::twoMM_mbpms

private

Cache variable

Definition at line 480 of file MVll.h.

double MVll::twoMV2

private

Cache variable

Definition at line 479 of file MVll.h.

double MVll::twothird

private

Cache variable

Definition at line 509 of file MVll.h.

gslpp::complex MVll::ubar

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 525 of file MVll.h.

gslpp::vector<double> MVll::V_cache

private

Cache variable

Definition at line 701 of file MVll.h.

unsigned int MVll::V_updated

private

Cache variable

Definition at line 700 of file MVll.h.

StandardModel::meson MVll::vectorM

private

Final vector meson type

Definition at line 441 of file MVll.h.

gslpp::vector<double> MVll::VL0_cache

private

Cache variable

Definition at line 740 of file MVll.h.

unsigned int MVll::VL0_updated

private

Cache variable

Definition at line 739 of file MVll.h.

unsigned int MVll::VL1_updated

private

Cache variable

Definition at line 727 of file MVll.h.

unsigned int MVll::VL2_updated

private

Cache variable

Definition at line 728 of file MVll.h.

unsigned int MVll::VR0_updated

private

Cache variable

Definition at line 745 of file MVll.h.

unsigned int MVll::VR1_updated

private

Cache variable

Definition at line 733 of file MVll.h.

unsigned int MVll::VR2_updated

private

Cache variable

Definition at line 734 of file MVll.h.

gsl_integration_cquad_workspace* MVll::w_delta

private

Gsl integral variable

Definition at line 671 of file MVll.h.

gsl_integration_cquad_workspace* MVll::w_DTPPR

private

Gsl integral variable

Definition at line 669 of file MVll.h.

gsl_integration_cquad_workspace* MVll::w_sigma

private

Gsl integral variable

Definition at line 670 of file MVll.h.

double MVll::width

private

Initial meson width

Definition at line 453 of file MVll.h.

gslpp::complex MVll::xm

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 530 of file MVll.h.

gslpp::complex MVll::xp

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 529 of file MVll.h.

gslpp::vector<double> MVll::Ycache

private

Cache variable

Definition at line 811 of file MVll.h.

gslpp::complex MVll::ym

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 532 of file MVll.h.

gslpp::complex MVll::yp

private

Variable used to compute the QCDF \(\Delta C_9\)

Definition at line 531 of file MVll.h.

unsigned int MVll::Yupdated

private

Cache variable

Definition at line 810 of file MVll.h.

double MVll::z_0

private

Cache variable

Definition at line 464 of file MVll.h.

unsigned int MVll::z_updated

private

Cache variable

Definition at line 718 of file MVll.h.

---

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

• MVll.h
• MVll.cpp