a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
Public Member Functions | Private Attributes | List of all members
MVlnu Class Reference

#include <MVlnu.h>

Detailed Description

Definition at line 21 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 BGLflag
 
bool btocNPpmflag
 
bool CLNflag
 
bool DMflag
 
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_ia reference to an object of type StandardModel
[in]meson_iinitial meson of the decay
[in]vector_ifinal vector meson of the decay
[in]lep_ifinal leptons of the decay

Definition at line 19 of file MVlnu.cpp.

20: mySM(SM_i), ig(ROOT::Math::Integration::kADAPTIVESINGULAR, ROOT::Math::Integration::kGAUSS51), wf()
21{
22 lep = lep_i;
23 meson = meson_i;
24 vectorM = vector_i;
25 CLNflag = false;
26 BGLflag = false;
27 DMflag = false;
28 btocNPpmflag = false;
29 NPanalysis = false;
30
31 checkcache_int_tau = false;
32 checkcache_int_mu = false;
33 checkcache_int_el = false;
34
35 double max_double = std::numeric_limits<double>::max();
36
37 hA1w1_cache = max_double;
38 rho2_cache = max_double;
39 R1w1_cache = max_double;
40 R2w1_cache = max_double;
41 N_A_cache = max_double;
42 N_1_cache =max_double;
43 N_2_cache = max_double;
44 j_A_cache= max_double;
45 j_0_cache = max_double;
46 j_1_cache = max_double;
47 j_2_cache = max_double;
48 k_A_cache = max_double;
49 k_0_cache = max_double;
50 k_1_cache = max_double;
51 k_2_cache = max_double;
52 l_A_cache = max_double;
53
54 af0_cache = max_double;
55 af1_cache = max_double;
56 af2_cache = max_double;
57 ag0_cache = max_double;
58 ag1_cache = max_double;
59 ag2_cache = max_double;
60 aF11_cache = max_double;
61 aF12_cache = max_double;
62 aF13_cache = max_double;
63 aF21_cache = max_double;
64 aF22_cache = max_double;
65 aF23_cache = max_double;
66
67 af_1_cache = max_double;
68 ag_1_cache = max_double;
69 aF1_1_cache = max_double;
70 aF2_1_cache = max_double;
71 af_2_cache = max_double;
72 ag_2_cache = max_double;
73 aF1_2_cache = max_double;
74 aF2_2_cache = max_double;
75 af_3_cache = max_double;
76 ag_3_cache = max_double;
77 aF1_3_cache = max_double;
78 aF2_3_cache = max_double;
79 af_4_cache = max_double;
80 ag_4_cache = max_double;
81 aF1_4_cache = max_double;
82 aF2_4_cache = max_double;
83 af_5_cache = max_double;
84 ag_5_cache = max_double;
85 aF1_5_cache = max_double;
86 aF2_5_cache = max_double;
87 af_6_cache = max_double;
88 ag_6_cache = max_double;
89 aF1_6_cache = max_double;
90 aF2_6_cache = max_double;
91 af_7_cache = max_double;
92 ag_7_cache = max_double;
93 aF1_7_cache = max_double;
94 aF2_7_cache = max_double;
95 af_8_cache = max_double;
96 ag_8_cache = max_double;
97 aF1_8_cache = max_double;
98 aF2_8_cache = max_double;
99 af_9_cache = max_double;
100 ag_9_cache = max_double;
101 aF1_9_cache = max_double;
102 aF2_9_cache = max_double;
103 af_10_cache = max_double;
104 ag_10_cache = max_double;
105 aF1_10_cache = max_double;
106 aF2_10_cache = max_double;
107 bf_1_cache = max_double;
108 bg_1_cache = max_double;
109 bF1_1_cache = max_double;
110 bF2_1_cache = max_double;
111 bf_2_cache = max_double;
112 bg_2_cache = max_double;
113 bF1_2_cache = max_double;
114 bF2_2_cache = max_double;
115 bf_3_cache = max_double;
116 bg_3_cache = max_double;
117 bF1_3_cache = max_double;
118 bF2_3_cache = max_double;
119 bf_4_cache = max_double;
120 bg_4_cache = max_double;
121 bF1_4_cache = max_double;
122 bF2_4_cache = max_double;
123 bf_5_cache = max_double;
124 bg_5_cache = max_double;
125 bF1_5_cache = max_double;
126 bF2_5_cache = max_double;
127 bf_6_cache = max_double;
128 bg_6_cache = max_double;
129 bF1_6_cache = max_double;
130 bF2_6_cache = max_double;
131 bf_7_cache = max_double;
132 bg_7_cache = max_double;
133 bF1_7_cache = max_double;
134 bF2_7_cache = max_double;
135 bf_8_cache = max_double;
136 bg_8_cache = max_double;
137 bF1_8_cache = max_double;
138 bF2_8_cache = max_double;
139 bf_9_cache = max_double;
140 bg_9_cache = max_double;
141 bF1_9_cache = max_double;
142 bF2_9_cache = max_double;
143 bf_10_cache = max_double;
144 bg_10_cache = max_double;
145 bF1_10_cache = max_double;
146 bF2_10_cache = max_double;
147
148 CS_cache = max_double;
149 CSp_cache = max_double;
150 CP_cache = max_double;
151 CPp_cache = max_double;
152 CV_cache = max_double;
153 CVp_cache = max_double;
154 CA_cache = max_double;
155 CAp_cache = max_double;
156 CT_cache = max_double;
157 CTp_cache = max_double;
158
159 ig.SetRelTolerance(1.E-6); // set relative tolerance
160 ig.SetAbsTolerance(1.E-6); // set absoulte tolerance
161
162}
MVlnu::btocNPpmflag
bool btocNPpmflag
Definition: MVlnu.h:166
MVlnu::meson
QCD::meson meson
Definition: MVlnu.h:160
MVlnu::BGLflag
bool BGLflag
Definition: MVlnu.h:164
MVlnu::mySM
const StandardModel & mySM
Definition: MVlnu.h:158
MVlnu::NPanalysis
bool NPanalysis
Definition: MVlnu.h:167
MVlnu::vectorM
QCD::meson vectorM
Definition: MVlnu.h:161
MVlnu::DMflag
bool DMflag
Definition: MVlnu.h:165
MVlnu::CLNflag
bool CLNflag
Definition: MVlnu.h:163
MVlnu::lep
QCD::lepton lep
Definition: MVlnu.h:159

◆ ~MVlnu()

MVlnu::~MVlnu ( )
virtual

Destructor.

Definition at line 164 of file MVlnu.cpp.

164 {
165}

Member Function Documentation

◆ get_hA1()

double MVlnu::get_hA1 ( double  w)

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

Returns
\( h_{A1}(w) \)

Definition at line 1518 of file MVlnu.cpp.

1519{
1520 updateParameters();
1521 double q2 = (2. * MM * MV)*(w0 - w);
1522
1523 return hA1(q2);
1524}
MVlnu::w0
double w0
Definition: MVlnu.h:177
MVlnu::MM
double MM
Definition: MVlnu.h:175
MVlnu::MV
double MV
Definition: MVlnu.h:176

◆ get_hA1w1()

double MVlnu::get_hA1w1 ( )

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

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

Definition at line 1511 of file MVlnu.cpp.

1512{
1513 updateParameters();
1514
1515 return hA1(q2max);
1516}

◆ get_hA2()

double MVlnu::get_hA2 ( double  w)

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

Returns
\( h_{A2}(w) \)

Definition at line 1526 of file MVlnu.cpp.

1527{
1528 updateParameters();
1529 double q2 = (2. * MM * MV)*(w0 - w);
1530
1531 return (hA1(q2) * (1. + w) - RV * (A0(q2) * (1. + MV_o_MM) - A2(q2) * (MV_o_MM - w))) / (1. + MV_o_MM*MV_o_MM - 2.* MV_o_MM * w);
1532 // return (hA1(q2) * (1. + w) - A0(q2) * RV * (1. + MV_o_MM) - A2(q2) * RV * (w - MV_o_MM)) / (1. + MV_o_MM*MV_o_MM - 2.*MV_o_MM*w) ;
1533}
MVlnu::MV_o_MM
double MV_o_MM
Definition: MVlnu.h:169
MVlnu::RV
double RV
Definition: MVlnu.h:178

◆ get_hA3()

double MVlnu::get_hA3 ( double  w)

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

Returns
\( h_{A3}(w) \)

Definition at line 1535 of file MVlnu.cpp.

1536{
1537 updateParameters();
1538 double q2 = (2. * MM * MV)*(w0 - w);
1539
1540 return A2(q2) * RV - MV_o_MM * get_hA2(w) ;
1541}
MVlnu::get_hA2
double get_hA2(double w)
return at
Definition: MVlnu.cpp:1526

◆ get_hV()

double MVlnu::get_hV ( double  w)

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

Returns
\( h_{V}(w) \)

Definition at line 1543 of file MVlnu.cpp.

1544{
1545 updateParameters();
1546 double q2 = (2. * MM * MV)*(w0 - w);
1547
1548 return V(q2) * RV;
1549}

◆ get_R0()

double MVlnu::get_R0 ( double  w)

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

Returns
\( R_{0}(w) \)

Definition at line 1571 of file MVlnu.cpp.

1572{
1573 updateParameters();
1574 double q2 = (2. * MM * MV)*(w0 - w);
1575
1576 if (CLNflag) return R0(q2);
1577 else if (BGLflag) return A0(q2) / A1(q2);
1578 return 0.;
1579}

◆ get_R1()

double MVlnu::get_R1 ( double  w)

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

Returns
\( R_{1}(w) \)

Definition at line 1551 of file MVlnu.cpp.

1552{
1553 updateParameters();
1554 double q2 = (2. * MM * MV)*(w0 - w);
1555
1556 if (CLNflag) return R1(q2);
1557 else if (BGLflag) return V(q2) * RV / hA1(q2);
1558 return 0.;
1559}

◆ get_R2()

double MVlnu::get_R2 ( double  w)

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

Returns
\( R_{2}(w) \)

Definition at line 1561 of file MVlnu.cpp.

1562{
1563 updateParameters();
1564 double q2 = (2. * MM * MV)*(w0 - w);
1565
1566 if (CLNflag) return R2(q2);
1567 else if (BGLflag) return A2(q2) * RV / hA1(q2);
1568 return 0.;
1569}

◆ 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 1492 of file MVlnu.cpp.

1493{
1494 updateParameters();
1495
1496 double aF10 = (MM - MV)*(phi_F1(0.) / phi_f(0.)) * af0;
1497 return af0 * af0 + af1 * af1 + af2 * af2 + aF10 * aF10 + aF11 * aF11 + aF12*aF12 + aF13*aF13;
1498}

◆ 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 1500 of file MVlnu.cpp.

1501{
1502 updateParameters();
1503
1504 double z0 = (sqrt(w0 + 1.) - M_SQRT2) / (sqrt(w0 + 1.) + M_SQRT2);
1505 double PfacF2z0 = (z0 - zP1) / (1. - z0 * zP1)*(z0 - zP2) / (1. - z0 * zP2)*(z0 - zP3) / (1. - z0 * zP3);
1506 double phiF2z0 = phi_F2(z0);
1507 double aF20 = PfacF2z0 * phiF2z0 * 2. * F1_BGL(0.) / (MM * MM - MV * MV) - aF21 * z0 - aF22 * z0*z0 - aF23 * z0*z0*z0;
1508 return aF20 * aF20 + aF21 * aF21 + aF22*aF22 + aF23*aF23;
1509}

◆ 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 1484 of file MVlnu.cpp.

1485{
1486 updateParameters();
1487
1488 return ag0 * ag0 + ag1 * ag1 + ag2*ag2;
1489
1490}

◆ 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_minminmum cos angle value of the integral
[in]chi_maxmaximum cos angle value of the integral
Returns
\( <d\Gamma/dcos\chi> \)

Definition at line 1453 of file MVlnu.cpp.

1454{
1455 updateParameters();
1456
1457 double intJ1s = integrateJ(1, q2min, q2max);
1458 double intJ1c = integrateJ(2, q2min, q2max);
1459 double intJ2s = integrateJ(3, q2min, q2max);
1460 double intJ2c = integrateJ(4, q2min, q2max);
1461 double intJ3 = integrateJ(5, q2min, q2max);
1462 double intJ9 = integrateJ(12, q2min, q2max);
1463
1464 return ((chi_max - chi_min)*(3. * intJ1c + 6. * intJ1s - intJ2c - 2. * intJ2s) / 4. +
1465 (sin(2. * chi_max) - sin(2. * chi_min)) / 2. * intJ3 -
1466 (cos(2. * chi_max) - cos(2. * chi_min)) / 2. * intJ9) / (2. * M_PI);
1467}

◆ 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_minminmum cos angle value of the integral
[in]cl_maxmaximum cos angle value of the integral
Returns
\( <d\Gamma/dcos\theta_{l}> \)

Definition at line 1382 of file MVlnu.cpp.

1383{
1384 updateParameters();
1385
1386 double intJ1s = integrateJ(1, q2min, q2max);
1387 double intJ1c = integrateJ(2, q2min, q2max);
1388 double intJ2s = integrateJ(3, q2min, q2max);
1389 double intJ2c = integrateJ(4, q2min, q2max);
1390 double intJ6s = integrateJ(8, q2min, q2max);
1391 double intJ6c = integrateJ(9, q2min, q2max);
1392
1393 return 3. / 8. * ((cl_max - cl_min)*(intJ1c + 2. * intJ1s)+
1394 (cl_max * cl_max - cl_min * cl_min) / 2. * (intJ6c + 2. * intJ6s)+
1395 (2. / 3. * (cl_max * cl_max * cl_max - cl_min * cl_min * cl_min)-(cl_max - cl_min))*(intJ2c + 2. * intJ2s));
1396}

◆ 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_minminmum cos angle value of the integral
[in]cV_maxmaximum cos angle value of the integral
Returns
\( <d\Gamma/dcos\theta_{V}> \)

Definition at line 1417 of file MVlnu.cpp.

1418{
1419 updateParameters();
1420
1421 double intJ1s = integrateJ(1, q2min, q2max);
1422 double intJ1c = integrateJ(2, q2min, q2max);
1423 double intJ2s = integrateJ(3, q2min, q2max);
1424 double intJ2c = integrateJ(4, q2min, q2max);
1425
1426 return 3. / 8. * ((cV_max * cV_max * cV_max - cV_min * cV_min * cV_min) / 3. * (3. * intJ1c - intJ2c)+
1427 ((cV_max - cV_min)-(cV_max * cV_max * cV_max - cV_min * cV_min * cV_min) / 3.)*(3. * intJ1s - intJ2s));
1428}

◆ 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_minmaximum q^2 of the integral
[in]w_maxminimum q^2 of the integral
Returns
\( <d\Gamma/dw> \)

Definition at line 1346 of file MVlnu.cpp.

1347{
1348 updateParameters();
1349
1350 double q2_min = (2. * MM * MV)*(w0 - w_max); // min is Mlep*Mlep;
1351 double q2_max = (2. * MM * MV)*(w0 - w_min); // max is (MM-MV)*(MM-MV);
1352
1353 double intJ1s = integrateJ(1, q2_min, q2_max);
1354 double intJ1c = integrateJ(2, q2_min, q2_max);
1355 double intJ2s = integrateJ(3, q2_min, q2_max);
1356 double intJ2c = integrateJ(4, q2_min, q2_max);
1357
1358 return 3. / 4. * (2. * intJ1s + intJ1c) - 1. / 4. * (2. * intJ2s + intJ2c);
1359}

◆ getFL()

double MVlnu::getFL ( )

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

Returns
\( <F_{L}> \)

Definition at line 1469 of file MVlnu.cpp.

1470{
1471 updateParameters();
1472
1473 double intJ1s = integrateJ(1, q2min, q2max);
1474 double intJ1c = integrateJ(2, q2min, q2max);
1475 double intJ2s = integrateJ(3, q2min, q2max);
1476 double intJ2c = integrateJ(4, q2min, q2max);
1477
1478 double DeltaJL = (3. * intJ1c - intJ2c) / 4.;
1479 double DeltaJ = 3. / 4. * (2. * intJ1s + intJ1c) - 1. / 4. * (2. * intJ2s + intJ2c);
1480 return DeltaJL/DeltaJ;
1481
1482}

◆ getMwidth()

double MVlnu::getMwidth ( )
inline

The width of the meson M.

Returns
\( \Gamma_M \)

Definition at line 152 of file MVlnu.h.

152 {
153 updateParameters();
154 return width;
155 }
MVlnu::width
double width
Definition: MVlnu.h:182

◆ getPlep()

double MVlnu::getPlep ( )

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

Returns
\( <P_{\ell}> \)

Definition at line 1647 of file MVlnu.cpp.

1648{
1649 updateParameters();
1650
1651 double DeltaGammaPlus = integrateGpm(1,q2min,q2max);
1652 double DeltaGammaMinus = integrateGpm(2,q2min,q2max);
1653
1654 return (DeltaGammaPlus-DeltaGammaMinus)/(DeltaGammaPlus+DeltaGammaMinus);
1655
1656}

◆ initializeMVlnuParameters()

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

Definition at line 167 of file MVlnu.cpp.

168{
169 CLNflag = mySM.getFlavour().getFlagCLN();
170 BGLflag = mySM.getFlavour().getFlagBGL();
171 DMflag = mySM.getFlavour().getFlagDM();
172 btocNPpmflag = (mySM.getModelName().compare("RealWeakEFTCCPM") == 0);
173 NPanalysis = (mySM.getModelName().compare("RealWeakEFTCCPM") == 0 || mySM.getModelName().compare("RealWeakEFTCC") == 0);
174
175 if (CLNflag + BGLflag + DMflag != true) throw std::runtime_error("MVlnu: Set only one among CLNflag, BGLflag, DMflag to true");
176 else mvlnuParameters = make_vector<std::string>();
177 if (CLNflag) {
178 mvlnuParameters.clear();
179 if (vectorM == StandardModel::D_star_P) mvlnuParameters = make_vector<std::string>()
180 << "hA1w1" << "rho2" << "R1w1" << "R2w1"
181 << "N_A" << "N_1" << "N_2" << "j_A" << "j_0" << "j_1" << "j_2"
182 << "k_A" << "k_0" << "k_1" << "k_2" << "l_A";
183 }
184 else if (BGLflag) {
185 mvlnuParameters.clear();
186 if (vectorM == StandardModel::D_star_P) mvlnuParameters = make_vector<std::string>()
187 << "af0" << "af1" << "af2" << "ag0" << "ag1" << "ag2"
188 << "aF11" << "aF12" << "aF13" << "aF21" << "aF22" << "aF23"
189 << "mBcstV1" << "mBcstV2" << "mBcstV3" << "mBcstV4"
190 << "mBcstA1" << "mBcstA2" << "mBcstA3" << "mBcstA4"
191 << "mBcstP1" << "mBcstP2" << "mBcstP3"
192 << "chiTV" << "chiTA" << "chiTP" << "nI";
193 }
194 else if (DMflag){
195 mvlnuParameters.clear();
196 if (vectorM == StandardModel::D_star_P) mvlnuParameters = make_vector<std::string>()
197 << "af_1" << "ag_1" << "aF1_1" << "aF2_1"
198 << "af_2" << "ag_2" << "aF1_2" << "aF2_2"
199 << "af_3" << "ag_3" << "aF1_3" << "aF2_3"
200 << "af_4" << "ag_4" << "aF1_4" << "aF2_4"
201 << "af_5" << "ag_5" << "aF1_5" << "aF2_5"
202 << "af_6" << "ag_6" << "aF1_6" << "aF2_6"
203 << "af_7" << "ag_7" << "aF1_7" << "aF2_7"
204 << "af_8" << "ag_8" << "aF1_8" << "aF2_8"
205 << "af_9" << "ag_9" << "aF1_9" << "aF2_9"
206 << "af_10" << "ag_10" << "aF1_10" << "aF2_10"
207 << "bf_1" << "bg_1" << "bF1_1" << "bF2_1"
208 << "bf_2" << "bg_2" << "bF1_2" << "bF2_2"
209 << "bf_3" << "bg_3" << "bF1_3" << "bF2_3"
210 << "bf_4" << "bg_4" << "bF1_4" << "bF2_4"
211 << "bf_5" << "bg_5" << "bF1_5" << "bF2_5"
212 << "bf_6" << "bg_6" << "bF1_6" << "bF2_6"
213 << "bf_7" << "bg_7" << "bF1_7" << "bF2_7"
214 << "bf_8" << "bg_8" << "bF1_8" << "bF2_8"
215 << "bf_9" << "bg_9" << "bF1_9" << "bF2_9"
216 << "bf_10" << "bg_10" << "bF1_10" << "bF2_10";
217 }
218 else {
219 std::stringstream out;
220 out << vectorM;
221 throw std::runtime_error("MVlnu: vector " + out.str() + " not implemented");
222 }
223
224 mySM.initializeMeson(meson);
225 mySM.initializeMeson(vectorM);
226 return mvlnuParameters;
227}
Flavour::getFlagBGL
bool getFlagBGL() const
Definition: Flavour.h:344
Flavour::getFlagCLN
bool getFlagCLN() const
Definition: Flavour.h:340
Flavour::getFlagDM
bool getFlagDM() const
Definition: Flavour.h:348
MVlnu::mvlnuParameters
std::vector< std::string > mvlnuParameters
Definition: MVlnu.h:162
Model::getModelName
std::string getModelName() const
A method to fetch the name of the model.
Definition: Model.h:59
QCD::D_star_P
@ D_star_P
Definition: QCD.h:352
QCD::initializeMeson
void initializeMeson(QCD::meson meson_i) const
A method to initialize a meson.
Definition: QCD.cpp:280
StandardModel::getFlavour
const Flavour & getFlavour() const
Definition: StandardModel.h:1033

Member Data Documentation

◆ BGLflag

bool MVlnu::BGLflag
private

A flag for switching to BGL parameterization

Definition at line 164 of file MVlnu.h.

◆ btocNPpmflag

bool MVlnu::btocNPpmflag
private

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

Definition at line 166 of file MVlnu.h.

◆ CLNflag

bool MVlnu::CLNflag
private

A flag for switching to CLN parameterization

Definition at line 163 of file MVlnu.h.

◆ DMflag

bool MVlnu::DMflag
private

A flag for switching to Dispersive Matrix method

Definition at line 165 of file MVlnu.h.

◆ GF

double MVlnu::GF
private

Fermi constant

Definition at line 172 of file MVlnu.h.

◆ lep

QCD::lepton MVlnu::lep
private

Final leptons type

Definition at line 159 of file MVlnu.h.

◆ Mb

double MVlnu::Mb
private

b quark mass

Definition at line 180 of file MVlnu.h.

◆ Mc

double MVlnu::Mc
private

charm quark mass

Definition at line 181 of file MVlnu.h.

◆ meson

QCD::meson MVlnu::meson
private

Initial meson type

Definition at line 160 of file MVlnu.h.

◆ Mlep

double MVlnu::Mlep
private

Lepton mass

Definition at line 173 of file MVlnu.h.

◆ MM

double MVlnu::MM
private

Initial meson mass

Definition at line 175 of file MVlnu.h.

◆ Mnu

double MVlnu::Mnu
private

Neutrino mass

Definition at line 174 of file MVlnu.h.

◆ mu_b

double MVlnu::mu_b
private

b mass scale

Definition at line 179 of file MVlnu.h.

◆ MV

double MVlnu::MV
private

Final vector meson mass

Definition at line 176 of file MVlnu.h.

◆ MV_o_MM

double MVlnu::MV_o_MM
private

Definition at line 169 of file MVlnu.h.

◆ mvlnuParameters

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

The string of mandatory MVlnu parameters

Definition at line 162 of file MVlnu.h.

◆ mySM

const StandardModel& MVlnu::mySM
private

Model type

Definition at line 158 of file MVlnu.h.

◆ NPanalysis

bool MVlnu::NPanalysis
private

A flag to switch to BSM analysis

Definition at line 167 of file MVlnu.h.

◆ RV

double MVlnu::RV
private

Dimensionless meson - vector mass ratio

Definition at line 178 of file MVlnu.h.

◆ sqrtMV_o_MM

double MVlnu::sqrtMV_o_MM
private

Definition at line 170 of file MVlnu.h.

◆ vectorM

QCD::meson MVlnu::vectorM
private

Final vector meson type

Definition at line 161 of file MVlnu.h.

◆ w0

double MVlnu::w0
private

Kinematic variable w at q2=0

Definition at line 177 of file MVlnu.h.

◆ width

double MVlnu::width
private

Initial meson width

Definition at line 182 of file MVlnu.h.

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