MVlnu Class Reference

#include <MVlnu.h>

Detailed Description

Definition at line 16 of file MVlnu.h.

Public Member Functions
double	get_hA1 (double w)
	return \( h_{A1} \) at \( w = w \) More...
double	get_hA1w1 ()
	return A1 form factor at \( w = 1 \) More...
double	get_hA2 (double w)
	return \( h_{A2} \) at \( w = w \) More...
double	get_hA3 (double w)
	return \( h_{A3} \) at \( w = w \) More...
double	get_hV (double w)
	return \( h_{V} \) at \( w = w \) More...
double	get_R0 (double w)
	return \( R_0 \) at \( w = w \) More...
double	get_R1 (double w)
	return \( R_1 \) at \( w = w \) More...
double	get_R2 (double w)
	return \( R_2 \) at \( w = w \) More...
double	get_unitarity_A_BGL ()
	Axial unitarity constraint for BGL parameters. More...
double	get_unitarity_P_BGL ()
	Pseudoscalar unitarity constraint for BGL parameters. More...
double	get_unitarity_V_BGL ()
	Vector unitarity constraint for BGL parameters. More...
double	getDeltaGammaDeltachi (double chi_min, double chi_max)
	The integral of \( d\Gamma/dcos\chi \) from \(cos\chi\) to \(cos\chi\). More...
double	getDeltaGammaDeltacl (double cl_min, double cl_max)
	The integral of \( d\Gamma/dcos\theta_{l} \) from \(cos\theta_{l,min}\) to \(cos\theta_{l,max}\). More...
double	getDeltaGammaDeltacV (double cV_min, double cV_max)
	The integral of \( d\Gamma/dcos\theta_{V} \) from \(cos\theta_{V,min}\) to \(cos\theta_{V,max}\). More...
double	getDeltaGammaDeltaw (double w_min, double w_max)
	The integral of \( d\Gamma/dw \) from \(w_{min}\) to \(w_{max}\). More...
double	getFL ()
	Binned D* polarization fraction \( F_{L} \). More...
double	getMwidth ()
	The width of the meson M. More...
double	getPlep ()
	Binned lepton helicity asymmetry \( P_{\ell} = (<G_{+}> - <G_{-}>) / (<G_{+}> + <G_{-}>) \). More...
std::vector< std::string >	initializeMVlnuParameters ()
	MVlnu (const StandardModel &SM_i, QCD::meson meson_i, QCD::meson vector_i, QCD::lepton lep_i)
	Constructor. More...
virtual	~MVlnu ()
	Destructor. More...

Private Attributes
bool	btocNPpmflag
bool	CLNflag
double	GF
QCD::lepton	lep
double	Mb
double	Mc
QCD::meson	meson
double	Mlep
double	MM
double	Mnu
double	mu_b
double	MV
double	MV_o_MM
std::vector< std::string >	mvlnuParameters
const StandardModel &	mySM
bool	NPanalysis
double	RV
double	sqrtMV_o_MM
QCD::meson	vectorM
double	w0
double	width

Constructor & Destructor Documentation

MVlnu()

MVlnu::MVlnu	(	const StandardModel &	SM_i,
		QCD::meson	meson_i,
		QCD::meson	vector_i,
		QCD::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 18 of file MVlnu.cpp.

: mySM(SM_i)
{
    lep = lep_i;
    meson = meson_i;
    vectorM = vector_i;
    CLNflag = false;
    btocNPpmflag = false;
    
    w_J = gsl_integration_cquad_workspace_alloc (100); 
    
    checkcache_int_tau = 0;
    checkcache_int_mu = 0;
    checkcache_int_el = 0;
    
    double max_double = std::numeric_limits<double>::max();
    
    hA1w1_cache = max_double;
    rho2_cache = max_double;
    R1w1_cache = max_double;
    R2w1_cache = max_double;
    N_A_cache = max_double;
    N_1_cache =max_double;
    N_2_cache = max_double;
    j_A_cache= max_double;
    j_0_cache = max_double;
    j_1_cache = max_double;
    j_2_cache = max_double;
    k_A_cache = max_double;
    k_0_cache = max_double;
    k_1_cache = max_double;
    k_2_cache = max_double;
    l_A_cache = max_double;
 
    af0_cache = max_double;
    af1_cache = max_double;
    af2_cache = max_double;
    ag0_cache = max_double;
    ag1_cache = max_double;
    ag2_cache = max_double;
    aF11_cache = max_double;
    aF12_cache = max_double;
    aF21_cache = max_double;
    aF22_cache = max_double;
 
    CS_cache = max_double;
    CSp_cache = max_double;
    CP_cache = max_double;
    CPp_cache = max_double;
    CV_cache = max_double;
    CVp_cache = max_double;
    CA_cache = max_double;
    CAp_cache = max_double;
    CT_cache = max_double;
    CTp_cache = max_double;
}

~MVlnu()

MVlnu::~MVlnu ( )

virtual

Destructor.

Definition at line 75 of file MVlnu.cpp.

              {
}

Member Function Documentation

get_hA1()

double MVlnu::get_hA1

(

double

w

)

return \( h_{A1} \) at \( w = w \)

Returns

\( h_{A1}(w) \)

Definition at line 1178 of file MVlnu.cpp.

{
    updateParameters();
 
    return hA1(q2max);
}
 
double MVlnu::get_hA1(double w)

get_hA1w1()

double MVlnu::get_hA1w1

(

)

return A1 form factor at \( w = 1 \)

Returns

\( h_{A1}((MM-MV)^2) \)

Definition at line 1171 of file MVlnu.cpp.

get_hA2()

double MVlnu::get_hA2

(

double

w

)

return \( h_{A2} \) at \( w = w \)

Returns

\( h_{A2}(w) \)

Definition at line 1186 of file MVlnu.cpp.

{

get_hA3()

double MVlnu::get_hA3

(

double

w

)

return \( h_{A3} \) at \( w = w \)

Returns

\( h_{A3}(w) \)

Definition at line 1195 of file MVlnu.cpp.

{

get_hV()

double MVlnu::get_hV

(

double

w

)

return \( h_{V} \) at \( w = w \)

Returns

\( h_{V}(w) \)

Definition at line 1203 of file MVlnu.cpp.

{

get_R0()

double MVlnu::get_R0

(

double

w

)

return \( R_0 \) at \( w = w \)

Returns

\( R_{0}(w) \)

Definition at line 1229 of file MVlnu.cpp.

{

get_R1()

double MVlnu::get_R1

(

double

w

)

return \( R_1 \) at \( w = w \)

Returns

\( R_{1}(w) \)

Definition at line 1211 of file MVlnu.cpp.

{

get_R2()

double MVlnu::get_R2

(

double

w

)

return \( R_2 \) at \( w = w \)

Returns

\( R_{2}(w) \)

Definition at line 1220 of file MVlnu.cpp.

{

get_unitarity_A_BGL()

double MVlnu::get_unitarity_A_BGL

(

)

Axial unitarity constraint for BGL parameters.

Returns

\( \sum_i (af_i^2 + aF1_i^2) \)

Definition at line 1152 of file MVlnu.cpp.

{

get_unitarity_P_BGL()

double MVlnu::get_unitarity_P_BGL

(

)

Pseudoscalar unitarity constraint for BGL parameters.

Returns

\( \sum_i aF2_i^2 \)

Definition at line 1160 of file MVlnu.cpp.

{
    updateParameters();
 

get_unitarity_V_BGL()

double MVlnu::get_unitarity_V_BGL

(

)

Vector unitarity constraint for BGL parameters.

Returns

\( \sum_i ag_i^2 \)

Definition at line 1144 of file MVlnu.cpp.

{

getDeltaGammaDeltachi()

double MVlnu::getDeltaGammaDeltachi	(	double	chi_min,
		double	chi_max
	)

The integral of \( d\Gamma/dcos\chi \) from \(cos\chi\) to \(cos\chi\).

Parameters

[in]	chi_min	minmum cos angle value of the integral
[in]	chi_max	maximum cos angle value of the integral

Returns

\( <d\Gamma/dcos\chi> \)

Definition at line 1113 of file MVlnu.cpp.

{
    updateParameters();
 
    double intJ1s = integrateJ(1, q2min, q2max);
    double intJ1c = integrateJ(2, q2min, q2max);
    double intJ2s = integrateJ(3, q2min, q2max);
    double intJ2c = integrateJ(4, q2min, q2max);
    double intJ3 = integrateJ(5, q2min, q2max);

getDeltaGammaDeltacl()

double MVlnu::getDeltaGammaDeltacl	(	double	cl_min,
		double	cl_max
	)

The integral of \( d\Gamma/dcos\theta_{l} \) from \(cos\theta_{l,min}\) to \(cos\theta_{l,max}\).

Parameters

[in]	cl_min	minmum cos angle value of the integral
[in]	cl_max	maximum cos angle value of the integral

Returns

\( <d\Gamma/dcos\theta_{l}> \)

Definition at line 1042 of file MVlnu.cpp.

{
    updateParameters();
 
    double intJ1s = integrateJ(1, q2min, q2max);
    double intJ1c = integrateJ(2, q2min, q2max);
    double intJ2s = integrateJ(3, q2min, q2max);
    double intJ2c = integrateJ(4, q2min, q2max);
    double intJ6s = integrateJ(8, q2min, q2max);

getDeltaGammaDeltacV()

double MVlnu::getDeltaGammaDeltacV	(	double	cV_min,
		double	cV_max
	)

The integral of \( d\Gamma/dcos\theta_{V} \) from \(cos\theta_{V,min}\) to \(cos\theta_{V,max}\).

Parameters

[in]	cV_min	minmum cos angle value of the integral
[in]	cV_max	maximum cos angle value of the integral

Returns

\( <d\Gamma/dcos\theta_{V}> \)

Definition at line 1077 of file MVlnu.cpp.

{
    updateParameters();
 
    double intJ1s = integrateJ(1, q2min, q2max);
    double intJ1c = integrateJ(2, q2min, q2max);

getDeltaGammaDeltaw()

double MVlnu::getDeltaGammaDeltaw	(	double	w_min,
		double	w_max
	)

The integral of \( d\Gamma/dw \) from \(w_{min}\) to \(w_{max}\).

Parameters

[in]	w_min	maximum q^2 of the integral
[in]	w_max	minimum q^2 of the integral

Returns

\( <d\Gamma/dw> \)

Definition at line 1006 of file MVlnu.cpp.

{
    updateParameters();
 
    return 3. / 4. * (2. * J1s(q2) + J1c(q2)) - 1. / 4. * (2. * J2s(q2) + J2c(q2));
}
 
double MVlnu::getDeltaGammaDeltaw(double w_min, double w_max)
{
    updateParameters();
 
    double q2_min = (2. * MM * MV)*(w0 - w_max); // min is Mlep*Mlep;
    double q2_max = (2. * MM * MV)*(w0 - w_min); // max is (MM-MV)*(MM-MV);
 
    double intJ1s = integrateJ(1, q2_min, q2_max);

getFL()

double MVlnu::getFL

(

)

Binned D* polarization fraction \( F_{L} \).

Returns

\( <F_{L}> \)

Definition at line 1129 of file MVlnu.cpp.

{
    updateParameters();
 
    double intJ1s = integrateJ(1, q2min, q2max);
    double intJ1c = integrateJ(2, q2min, q2max);
    double intJ2s = integrateJ(3, q2min, q2max);
    double intJ2c = integrateJ(4, q2min, q2max);

getMwidth()

double MVlnu::getMwidth ( )

inline

The width of the meson M.

Returns

\( \Gamma_M \)

Definition at line 147 of file MVlnu.h.

                      {
        updateParameters();
        return width;
    }

getPlep()

double MVlnu::getPlep

(

)

Binned lepton helicity asymmetry \( P_{\ell} = (<G_{+}> - <G_{-}>) / (<G_{+}> + <G_{-}>) \).

Returns

\( <P_{\ell}> \)

Definition at line 1306 of file MVlnu.cpp.

                                                     : index " + out.str() + " not implemented");
    }
}
 
double MVlnu::getPlep()
{
    updateParameters();

initializeMVlnuParameters()

std::vector< std::string > MVlnu::initializeMVlnuParameters

(

)

Definition at line 78 of file MVlnu.cpp.

{
    CLNflag = mySM.getFlavour().getFlagCLN();
    btocNPpmflag = (mySM.getModelName().compare("RealWeakEFTCCPM") == 0);
    NPanalysis = (mySM.getModelName().compare("RealWeakEFTCCPM") == 0 || mySM.getModelName().compare("RealWeakEFTCC") == 0);
 
    if (vectorM == StandardModel::D_star_P) mvlnuParameters = make_vector<std::string>()
        << "af0" << "af1" << "af2" << "ag0" << "ag1" << "ag2"
        << "aF11" << "aF12" << "aF21" << "aF22"
        << "mBcstV1" << "mBcstV2" << "mBcstV3" << "mBcstV4"
        << "mBcstA1" << "mBcstA2" << "mBcstA3" << "mBcstA4"
        << "mBcstP1" << "mBcstP2" << "mBcstP3"
        << "chiTV" << "chiTA" << "chiTP" << "nI";
    else {
        std::stringstream out;
        out << vectorM;
        throw std::runtime_error("MVlnu: vector " + out.str() + " not implemented");
    }
 
    if (CLNflag) {
        mvlnuParameters.clear();
        if (vectorM == StandardModel::D_star_P) mvlnuParameters = make_vector<std::string>()
            << "hA1w1" << "rho2" << "R1w1" << "R2w1"
            << "N_A" << "N_1" << "N_2" << "j_A" << "j_0" << "j_1" << "j_2"
            << "k_A" << "k_0" << "k_1" << "k_2" << "l_A";
    }    
 
    mySM.initializeMeson(meson);
    mySM.initializeMeson(vectorM);
    return mvlnuParameters;
}

Member Data Documentation

btocNPpmflag

bool MVlnu::btocNPpmflag

private

A flag for switching to the +/- basis for NP Wilson coefficients

Definition at line 159 of file MVlnu.h.

CLNflag

bool MVlnu::CLNflag

private

A flag for switching to CLN parameterization

Definition at line 158 of file MVlnu.h.

GF

double MVlnu::GF

private

Fermi constant

Definition at line 165 of file MVlnu.h.

lep

QCD::lepton MVlnu::lep

private

Final leptons type

Definition at line 154 of file MVlnu.h.

Mb

double MVlnu::Mb

private

b quark mass

Definition at line 173 of file MVlnu.h.

Mc

double MVlnu::Mc

private

charm quark mass

Definition at line 174 of file MVlnu.h.

meson

QCD::meson MVlnu::meson

private

Initial meson type

Definition at line 155 of file MVlnu.h.

Mlep

double MVlnu::Mlep

private

Lepton mass

Definition at line 166 of file MVlnu.h.

MM

double MVlnu::MM

private

Initial meson mass

Definition at line 168 of file MVlnu.h.

Mnu

double MVlnu::Mnu

private

Neutrino mass

Definition at line 167 of file MVlnu.h.

mu_b

double MVlnu::mu_b

private

b mass scale

Definition at line 172 of file MVlnu.h.

MV

double MVlnu::MV

private

Final vector meson mass

Definition at line 169 of file MVlnu.h.

MV_o_MM

double MVlnu::MV_o_MM

private

Definition at line 162 of file MVlnu.h.

mvlnuParameters

std::vector<std::string> MVlnu::mvlnuParameters

private

The string of mandatory MVlnu parameters

Definition at line 157 of file MVlnu.h.

mySM

const StandardModel& MVlnu::mySM

private

Model type

Definition at line 153 of file MVlnu.h.

NPanalysis

bool MVlnu::NPanalysis

private

A flag to switch to BSM analysis

Definition at line 160 of file MVlnu.h.

RV

double MVlnu::RV

private

Dimensionless meson - vector mass ratio

Definition at line 171 of file MVlnu.h.

sqrtMV_o_MM

double MVlnu::sqrtMV_o_MM

private

Definition at line 163 of file MVlnu.h.

vectorM

QCD::meson MVlnu::vectorM

private

Final vector meson type

Definition at line 156 of file MVlnu.h.

w0

double MVlnu::w0

private

Kinematic variable w at q2=0

Definition at line 170 of file MVlnu.h.

width

double MVlnu::width

private

Initial meson width

Definition at line 175 of file MVlnu.h.

---

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

• MVlnu.h
• MVlnu.cpp