#include <MPlnu.h>

Detailed Description

Definition at line 21 of file MPlnu.h.

Public Member Functions
double	get_f0 (double q2)
	return f0 form factor at \( q^2 \) More...

double	get_fplus (double q2)
	return fplus form factor at \( q^2 \) More...

double	get_fT (double q2)
	return fT form factor at \( q^2 \) More...

double	get_strong_unitarity_BGL ()
	Strong Unitarity constraint for BGL parameters using HQET. More...

double	get_unitarity_0plus_BGL ()
	Weak Unitarity constraint for BGL parameters related to 0+ resonances. More...

double	get_unitarity_1min_BGL ()
	Weak Unitarity constraint for BGL parameters related to 1- resonances. More...

double	getaf0_0 ()
	The BGL parameter \( a_0^{f_0}\). More...

double	getDeltaGammaDeltaw (double w_min, double w_max)
	The integral of \( d\Gamma/dw \) from \(w_{min}\) to \(w_{max}\). More...

double	getMwidth ()
	The width of the meson M. More...

std::vector< std::string >	initializeMPlnuParameters ()

	MPlnu (const StandardModel &SM_i, QCD::meson meson_i, QCD::meson pseudoscalar_i, QCD::lepton lep_i)
	Constructor. More...

virtual	~MPlnu ()
	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	MP

std::vector< std::string >	mplnuParameters

double	mu_b

const StandardModel &	mySM

bool	NPanalysis

QCD::meson	pseudoscalarM

double	RV

double	w0

double	width

double	z0

Constructor & Destructor Documentation

◆ MPlnu()

MPlnu::MPlnu	(	const StandardModel &	SM_i,
		QCD::meson	meson_i,
		QCD::meson	pseudoscalar_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 20 of file MPlnu.cpp.

 : mySM(SM_i)
 {
     lep = lep_i;
     meson = meson_i;
     pseudoscalarM = pseudoscalar_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();
     
     fplusz0_cache = max_double;
     rho1to2_cache = max_double;
     N_0_cache = max_double;
     alpha_0_cache = max_double;
     alpha_p_cache = max_double;
     beta_0_cache = max_double;
     beta_p_cache = max_double;
     gamma_0_cache = max_double;
     gamma_p_cache = max_double;
     af0_1_cache = max_double;
     af0_2_cache = max_double;
     afplus_0_cache = max_double;
     afplus_1_cache = max_double;
     afplus_2_cache = max_double;
     
 #if NBGL == 3
     af0_3_cache = max_double;
     afplus_3_cache = max_double;
 #endif 
     
     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;
 }

◆ ~MPlnu()

MPlnu::~MPlnu ( )

virtual

Destructor.

Definition at line 69 of file MPlnu.cpp.

70 {}

Member Function Documentation

◆ get_f0()

double MPlnu::get_f0 ( double q2 )

return f0 form factor at \( q^2 \)

Returns: \( f_{0}(q^2) \)

Definition at line 629 of file MPlnu.cpp.

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

◆ get_fplus()

double MPlnu::get_fplus ( double q2 )

return fplus form factor at \( q^2 \)

Returns: \( f_{+}(q^2) \)

Definition at line 622 of file MPlnu.cpp.

◆ get_fT()

double MPlnu::get_fT ( double q2 )

return fT form factor at \( q^2 \)

Returns: \( q^2 \)

Definition at line 636 of file MPlnu.cpp.

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

◆ get_strong_unitarity_BGL()

double MPlnu::get_strong_unitarity_BGL ( )

Strong Unitarity constraint for BGL parameters using HQET.

Returns: \( \Sum_i (af0_i^2) \)

Definition at line 609 of file MPlnu.cpp.

 {
     updateParameters();
  
 #if NBGL == 3
     return 1707.54 * afplus_0 * afplus_0 + 1299.57 * afplus_0 * afplus_1 + 442.82 * afplus_1 * afplus_1 - 356.01 * afplus_0 * afplus_2 

◆ get_unitarity_0plus_BGL()

double MPlnu::get_unitarity_0plus_BGL ( )

Weak Unitarity constraint for BGL parameters related to 0+ resonances.

Returns: \( \Sum_i (af0_i^2) \)

Definition at line 598 of file MPlnu.cpp.

 {
     updateParameters();
  

◆ get_unitarity_1min_BGL()

double MPlnu::get_unitarity_1min_BGL ( )

Weak Unitarity constraint for BGL parameters related to 1- resonances.

Returns: \( \Sum_i (af+_i^2) \)

Definition at line 587 of file MPlnu.cpp.

 {
     updateParameters();
  

◆ getaf0_0()

double MPlnu::getaf0_0 ( )

inline

The BGL parameter \( a_0^{f_0}\).

Returns: \( a_0^{f_0}\)

Definition at line 96 of file MPlnu.h.

     {
         updateParameters();
         return af0_0;
     }

◆ getDeltaGammaDeltaw()

double MPlnu::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 571 of file MPlnu.cpp.

 {
     updateParameters();
     
     double q2_min = (2. * MM * MP)*(w0 - w_max); // min is Mlep*Mlep;
     double q2_max = (2. * MM * MP)*(w0 - w_min); // max is (MM-MP)*(MM-MP);
     
     double intJ1 = integrateJ(1, q2_min, q2_max);
     double intJ3 = integrateJ(3, q2_min, q2_max);

◆ getMwidth()

double MPlnu::getMwidth ( )

inline

The width of the meson M.

Returns: \( \Gamma_M \)

Definition at line 86 of file MPlnu.h.

     {
         updateParameters();
         return width;
     }

◆ initializeMPlnuParameters()

std::vector< std::string > MPlnu::initializeMPlnuParameters ( )

Definition at line 72 of file MPlnu.cpp.

 {
     CLNflag = mySM.getFlavour().getFlagCLN();
     btocNPpmflag = (mySM.getModelName().compare("RealWeakEFTCCPM") == 0);
     NPanalysis = (mySM.getModelName().compare("RealWeakEFTCCPM") == 0 || mySM.getModelName().compare("RealWeakEFTCC") == 0);
  
     if (pseudoscalarM == StandardModel::D_P) mplnuParameters = make_vector<std::string>()
         << "af0_1" << "af0_2" << "afplus_0" << "afplus_1" << "afplus_2"
 #if NBGL == 3
         << "af0_3" << "afplus_3"
 #endif                
         << "mBc1m_1" << "mBc1m_2" << "mBc1m_3" << "mBc1m_4"
         << "mBc0p_1" << "mBc0p_2" << "chitildeT" << "chiL" << "nI";
     else {
         std::stringstream out;
         out << pseudoscalarM;
         throw std::runtime_error("MPlnu: vector " + out.str() + " not implemented");
     }
  
     if (CLNflag) {
         mplnuParameters.clear();
         if (pseudoscalarM == StandardModel::D_P) mplnuParameters = make_vector<std::string>()
             << "fplusz0" << "rho1to2"
             << "N_0" << "alpha_0" << "alpha_p" << "beta_0" << "beta_p" << "gamma_0" << "gamma_p";
     }
  
     mySM.initializeMeson(meson);
     mySM.initializeMeson(pseudoscalarM);
     return mplnuParameters;
 }

Member Data Documentation

◆ btocNPpmflag

bool MPlnu::btocNPpmflag

private

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

Definition at line 109 of file MPlnu.h.

◆ CLNflag

bool MPlnu::CLNflag

private

A flag for switching to CLN parameterization

Definition at line 108 of file MPlnu.h.

◆ GF

double MPlnu::GF

private

Fermi constant

Definition at line 112 of file MPlnu.h.

◆ lep

QCD::lepton MPlnu::lep

private

Final leptons type

Definition at line 104 of file MPlnu.h.

◆ Mb

double MPlnu::Mb

private

b quark mass

Definition at line 121 of file MPlnu.h.

◆ Mc

double MPlnu::Mc

private

charm quark mass

Definition at line 122 of file MPlnu.h.

◆ meson

QCD::meson MPlnu::meson

private

Initial meson type

Definition at line 105 of file MPlnu.h.

◆ Mlep

double MPlnu::Mlep

private

Lepton mass

Definition at line 113 of file MPlnu.h.

◆ MM

double MPlnu::MM

private

Initial meson mass

Definition at line 115 of file MPlnu.h.

◆ Mnu

double MPlnu::Mnu

private

Neutrino mass

Definition at line 114 of file MPlnu.h.

◆ MP

double MPlnu::MP

private

Final pseudoscalar meson mass

Definition at line 116 of file MPlnu.h.

◆ mplnuParameters

std::vector<std::string> MPlnu::mplnuParameters

private

The string of mandatory MPlnu parameters

Definition at line 107 of file MPlnu.h.

◆ mu_b

double MPlnu::mu_b

private

b mass scale

Definition at line 120 of file MPlnu.h.

◆ mySM

const StandardModel& MPlnu::mySM

private

Model type

Definition at line 103 of file MPlnu.h.

◆ NPanalysis

bool MPlnu::NPanalysis

private

A flag to switch to BSM analysis

Definition at line 110 of file MPlnu.h.

◆ pseudoscalarM

QCD::meson MPlnu::pseudoscalarM

private

Final vector meson type

Definition at line 106 of file MPlnu.h.

◆ RV

double MPlnu::RV

private

Dimensionless meson - vector mass ratio

Definition at line 119 of file MPlnu.h.

◆ w0

double MPlnu::w0

private

Kinematic variable w at q2=0

Definition at line 117 of file MPlnu.h.

◆ width

double MPlnu::width

private

Initial meson width

Definition at line 123 of file MPlnu.h.

◆ z0

double MPlnu::z0

private

Kinematic variable z at q2=0

Definition at line 118 of file MPlnu.h.

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

Detailed Description

Public Member Functions

Private Attributes

Constructor & Destructor Documentation

◆ MPlnu()

◆ ~MPlnu()

Member Function Documentation

◆ get_f0()

◆ get_fplus()

◆ get_fT()

◆ get_strong_unitarity_BGL()

◆ get_unitarity_0plus_BGL()

◆ get_unitarity_1min_BGL()

◆ getaf0_0()

◆ getDeltaGammaDeltaw()

◆ getMwidth()

◆ initializeMPlnuParameters()

Member Data Documentation

◆ btocNPpmflag

◆ CLNflag

◆ GF

◆ lep

◆ Mb

◆ Mc

◆ meson

◆ Mlep

◆ MM

◆ Mnu

◆ MP

◆ mplnuParameters

◆ mu_b

◆ mySM

◆ NPanalysis

◆ pseudoscalarM

◆ RV

◆ w0

◆ width

◆ z0