 |
Belle II Software
release-05-01-25
|
12 #include "EvtGenBase/EvtComplex.hh"
30 class EvtBSemiTauonicHelicityAmplitudeCalculator {
60 const double AS1,
const double AR3,
61 const double bottomMass,
const double charmMass,
62 const EvtComplex& CV1,
const EvtComplex& CV2,
const EvtComplex& CS1,
const EvtComplex& CS2,
const EvtComplex& CT,
63 const double parentMass,
const double DMass,
const double DstarMass);
76 EvtComplex
helAmp(
double mtau,
int tauhel,
int Dhel,
double w,
double costau)
const;
93 EvtComplex
helAmp(
const EvtComplex& CV1,
const EvtComplex& CV2,
const EvtComplex& CS1,
const EvtComplex& CS2,
const EvtComplex& CT,
94 double mtau,
int tauhel,
int Dhel,
double w,
double costau)
const;
104 double Lep(
const double mtau,
int tauhel,
int whel,
double q2,
double costau)
const;
113 double Lep(
const double mtau,
int tauhel,
double q2,
double costau)
const;
125 double Lep(
const double mtau,
int tauhel,
int whel1,
int whel2,
double q2,
double costau)
const;
133 double HadV1(
int Dhel,
int whel,
double w)
const;
141 double HadV2(
int Dhel,
int whel,
double w)
const;
149 double HadS1(
int Dhel,
double w)
const;
157 double HadS2(
int Dhel,
double w)
const;
167 double HadT(
int Dhel,
int whel1,
int whel2,
double w)
const;
178 double helampSM(
double mtau,
int tauhel,
int Dhel,
double w,
double costau)
const;
190 double helampV1(
double mtau,
int tauhel,
int Dhel,
double w,
double costau)
const;
202 double helampV2(
double mtau,
int tauhel,
int Dhel,
double w,
double costau)
const;
214 double helampS1(
double mtau,
int tauhel,
int Dhel,
double w,
double costau)
const;
226 double helampS2(
double mtau,
int tauhel,
int Dhel,
double w,
double costau)
const;
238 double helampT(
double mtau,
int tauhel,
int Dhel,
double w,
double costau)
const;
244 double eta(
int whel)
const {
return (whel == 2) ? -1 : 1;}
252 double hp(
double w)
const;
255 double hm(
double w)
const;
258 double hA1(
double w)
const;
261 double hV(
double w)
const;
264 double hA2(
double w)
const;
267 double hA3(
double w)
const;
272 double hS(
double w)
const;
275 double hP(
double w)
const;
280 double hT(
double w)
const;
283 double hT1(
double w)
const;
286 double hT2(
double w)
const;
289 double hT3(
double w)
const;
295 double z(
double w)
const;
301 double ffV1(
double w)
const;
307 double ffS1(
double w)
const;
313 double ffA1(
double w)
const;
319 double ffR1(
double w)
const;
325 double ffR2(
double w)
const;
331 double ffR3(
double w)
const;
337 double ffV11()
const {
return 1.;}
343 double ffA11()
const {
return 1.;}
349 double dS1(
double w)
const;
355 double dR3(
double w)
const;
371 double mD(
int Dhel)
const;
377 double r(
int Dhel)
const {
return mD(Dhel) /
m_mB;}
389 double v(
double mtau,
double q2)
const;
396 double q2(
int Dhel,
double w)
const;
403 double qh2(
int Dhel,
double w)
const;
411 double q2min(
double mtau)
const {
return mtau * mtau;}
424 double wfunc(
int Dhel,
double q2)
const {
return (1. +
r(Dhel) *
r(Dhel) -
q2 /
m_mB /
m_mB) / 2. /
r(Dhel);}
429 double wmin()
const {
return 1.;};
436 double wmax(
double mtau,
int Dhel)
const {
return wfunc(Dhel,
q2min(mtau));}
double hT2(double w) const
D* tensor form factor h_{T2}(w).
double dS1(double w) const
HQET correction factor for the scalar form factor for B->Dtaunu.
double HadV1(int Dhel, int whel, double w) const
The function to calculate the Hadronic Amplitudes of left handed (V-A) type contribution.
double aR3() const
HQET correction factor for the uncertainty of 1/m_Q correction.
void setRhoA12(double v)
Sets the form factor parameter rho_A1^2.
double hA1(double w) const
HQET D* axial vector form factor h_{A1}(w).
double z(double w) const
CLN form factor z.
double hT(double w) const
D tensor form factor h_T(w) in terms of vector form factors.
double m_rhoA12
Form factor slope parameters rho_A1^2.
double HadV2(int Dhel, int whel, double w) const
The function to calculate the Hadronic Amplitudes of right handed (V+A) type contribution.
double helampT(double mtau, int tauhel, int Dhel, double w, double costau) const
Helicity Amplitudes of tensor type contribution.
void setCT(const EvtComplex &v)
Sets the Wilson coeffcient CT.
EvtBSemiTauonicHelicityAmplitudeCalculator()
The default constructor.
virtual ~EvtBSemiTauonicHelicityAmplitudeCalculator()
The destructor.
double m_mDst
daughter vector (D*) meson mass.
void setAS1(double v)
Sets the form factor 1/m_Q correction parameter a_S1.
double m_ffR11
Form factor parameter R_1(1).
double m_mB
parent (B) meson mass.
void setCV2(const EvtComplex &v)
Sets the Wilson coeffcient CV2.
double ffS1(double w) const
CLN form factor S1.
double wfunc(int Dhel, double q2) const
Calculate the velocity transfer variable w.
double hT1(double w) const
D* tensor form factor h_{T1}(w) in terms of axial vector form factors.
double getRho12() const
parameter accessor
EvtComplex m_CV1
Wilson coefficient CV1.
double ffR3(double w) const
CLN form factor R3.
double qh2(int Dhel, double w) const
Function to calculate the q^2 divided by the square of parent mass (m_B^2).
double q2(int Dhel, double w) const
Function to calculate the q^2 of the decay (square of l+nu invariant mass).
double v(double mtau, double q2) const
Function to calculate the tau velocity.
void setMBottom(double m)
Returns the bottom quark mass.
void setMDst(double m)
Sets the daughter vector (D) meson mass.
double getMDst() const
Returns the daughter vector (D*) meson mass.
double getMD() const
Returns the daughter scalar (D) meson mass.
double getMBottom() const
Returns the bottom quark mass.
void setMB(double m)
Sets the parent (B) meson mass.
double hA3(double w) const
HQET D* axial vector form factor h_{A3}(w).
void setCV1(const EvtComplex &v)
Sets the Wilson coeffcient CV1.
void setRho12(double v)
Sets the form factor parameter rho_1^2.
double getAS1() const
Returns form factor 1/m_Q correction factor a_S1.
void setR21(double v)
Sets the form factor parameter R_2(1).
double ffV1(double w) const
CLN form factor V1.
double rq() const
Ratio of the charm quark mass to the charm quark mass.
double r(int Dhel) const
Ratio of the daughter meson mass to the parent meson.
EvtComplex getCS1() const
Returns the Wilson coeffcient CS1.
double dR3(double w) const
HQET correction factor for the scalar form factor for B->D*taunu.
double m_ffR21
Form factor parameter R_2(1).
double wmax(double mtau, int Dhel) const
Maximum value of the velocity transfer variable w.
double hm(double w) const
HQET D vector form factor h_-(w).
double getMB() const
Returns the parent (B) meson mass.
double aS1() const
HQET correction factor for the uncertainty of 1/m_Q correction.
double hV(double w) const
HQET D* axial vector form factor h_V(w).
double getR21() const
Returns form factor parameter R_2(1).
void setAR3(double v)
Sets the form factor 1/m_Q correction parameter a_R3.
double helampS2(double mtau, int tauhel, int Dhel, double w, double costau) const
Helicity Amplitudes of scalar (S-P) type contribution.
double q2max(int Dhel) const
Maximum value of the q^2.
double hS(double w) const
D scalar form factor h_S(w) in terms of vector form factors.
double ffV11() const
Form factor normalization factor for B->Dlnu.
double m_mD
daughter scalar (D) meson mass.
double eta(int whel) const
The metric factor.
double Lep(const double mtau, int tauhel, int whel, double q2, double costau) const
The function to calculate the Leptonic Amplitudes for B->D*taunu decay of the vector type contributio...
Abstract base class for different kinds of events.
EvtComplex getCV1() const
Returns the Wilson coeffcient CV1.
double m_aR3
1/mQ correcion factor a_R3.
double getR11() const
Returns form factor parameter R_1(1).
double getMCharm() const
Returns the charm quark mass.
double wmin() const
Minimum value of the velocity transfer variable w.
EvtComplex getCS2() const
Returns the Wilson coeffcient CS2.
double HadS2(int Dhel, double w) const
The function to calculate the Hadronic Amplitudes of scalar (S-P) type contribution.
void setCS1(const EvtComplex &v)
Sets the Wilson coeffcient CS1.
double m_aS1
1/mQ correcion factor a_S1.
double ffR1(double w) const
CLN form factor R1.
EvtComplex m_CS1
Wilson coefficient CS1.
void setCS2(const EvtComplex &v)
Sets the Wilson coeffcient CS2.
EvtComplex getCV2() const
Returns the Wilson coeffcient CV2.
double helampS1(double mtau, int tauhel, int Dhel, double w, double costau) const
Helicity Amplitudes of scalar (S+P) type contribution.
EvtComplex getCT() const
Returns the Wilson coeffcient CT.
double HadT(int Dhel, int whel1, int whel2, double w) const
The function to calculate the Hadronic Amplitudes of tensor type contribution.
double q2min(double mtau) const
Minimum value of the q^2.
double m_mBottom
b quark mass (running mass at m_b scale), used for scalar form factor term
double ffA1(double w) const
CLN form factor A1.
double HadS1(int Dhel, double w) const
The function to calculate the Hadronic Amplitudes of scalar (S+P) type contribution.
double hT3(double w) const
D* tensor form factor h_{T3}(w).
double helampSM(double mtau, int tauhel, int Dhel, double w, double costau) const
Helicity Amplitudes of SM (left handed) contribution.
void setR11(double v)
Sets the form factor parameter R_1(1).
bool chkDhel(int Dhel) const
sanity checkers
void setMD(double m)
Sets the daughter scalar (D) meson mass.
double helampV2(double mtau, int tauhel, int Dhel, double w, double costau) const
Helicity Amplitudes of right handed (V+A) contribution.
double getAR3() const
Returns form factor 1/m_Q correction factor a_R3.
double mD(int Dhel) const
Daughter D(*) meson mass.
EvtComplex m_CS2
Wilson coefficient CS2.
double ffA11() const
Form factor normalization factor for B->D*lnu.
double m_mCharm
c quark mass (running mass at m_b scale), used for scalar form factor term )
EvtComplex helAmp(double mtau, int tauhel, int Dhel, double w, double costau) const
The function calculates the helicity amplitude.
double helampV1(double mtau, int tauhel, int Dhel, double w, double costau) const
Helicity Amplitudes of left handed (V-A) contribution.
double getRhoA12() const
Returns form factor parameter rho_A1^2.
bool chktauhel(int tauhel) const
Function to check if tauhel is in the valid range.
double hp(double w) const
HQET D vector form factor h_+(w).
bool chkwhel(int whel) const
Function to check if whel is in the valid range.
void setMCharm(double m)
Returns the charm quark mass.
double hP(double w) const
D* pseudo scalar form factor h_P(w) in terms of axial vector form factors.
double hA2(double w) const
HQET D* axial vector form factor h_{A2}(w).
EvtComplex m_CT
Wilson coefficient CT.
double ffR2(double w) const
CLN form factor R2.
EvtComplex m_CV2
Wilson coefficient CV2.
double m_rho12
Form factor slope parameters rho_1^2.