Page MenuHomeHEPForge
Differential D90 Diff 388 src/EvtGenModels/Evtbs2llGammaFFMNT.cpp

Changeset View

Standalone View

View Options

src/EvtGenModels/Evtbs2llGammaFFMNT.cpp

Show First 20 Lines • Show All 52 Lines • ▼ Show 20 Linesvoid Evtbs2llGammaFFMNT::getPhotonFF( int decay_id, double fb, EvtId parent,
EvtComplex& Fv, EvtComplex& Fa, EvtComplex& Fv, EvtComplex& Fa,
EvtComplex& Ftv, EvtComplex& Fta ) EvtComplex& Ftv, EvtComplex& Fta )
{ {
EvtComplex unit1( 1.0, 0.0 ); // real unit EvtComplex unit1( 1.0, 0.0 ); // real unit
EvtComplex uniti( 0.0, 1.0 ); // imaginary unit EvtComplex uniti( 0.0, 1.0 ); // imaginary unit
EvtComplex unit0( 0.0, 0.0 ); // complex zero unit EvtComplex unit0( 0.0, 0.0 ); // complex zero unit
// characteristics of resonances rho, omega, phi // characteristics of resonances rho, omega, phi
double M_res[] = {0.7758, 0.78259, 1.019456}; // particle masses, Gev double M_res[] = { 0.7758, 0.78259, 1.019456 }; // particle masses, Gev
double Gamma[] = {0.1503, 0.00849, 0.00426}; // particle widthes, Gev double Gamma[] = { 0.1503, 0.00849, 0.00426 }; // particle widthes, Gev
double f_lept[] = {5.04, 17.1, -13.2}; // decay constants f_i double f_lept[] = { 5.04, 17.1, -13.2 }; // decay constants f_i
double g_plus[] = {0.27, -0.27, -0.38}; // and form-factors g+(0) double g_plus[] = { 0.27, -0.27, -0.38 }; // and form-factors g+(0)
g_plus[0] = g_plus[0] / sqrt( 2.0 ); // by D.Melikhov, N.Nikitin, K.Toms, g_plus[0] = g_plus[0] / sqrt( 2.0 ); // by D.Melikhov, N.Nikitin, K.Toms,
g_plus[1] = g_plus[1] / sqrt( 2.0 ); // Phys.At.Nucl. 68, p.1842 (2005) g_plus[1] = g_plus[1] / sqrt( 2.0 ); // Phys.At.Nucl. 68, p.1842 (2005)
double hatq2 = q2 / pow( M1, 2 ); double hatq2 = q2 / pow( M1, 2 );
// E - photon energy in the B-meson rest frame // E - photon energy in the B-meson rest frame
double E = 0.5 * M1 * ( 1 - hatq2 ); double E = 0.5 * M1 * ( 1 - hatq2 );
// parametrs for form-factors Fv, Ftv, Fa, Fta // parametrs for form-factors Fv, Ftv, Fa, Fta
//(by D.Melikhov, N.Nikitin, K.Toms, Yad. Fiz. 62, No 11) //(by D.Melikhov, N.Nikitin, K.Toms, Yad. Fiz. 62, No 11)
double beta[] = {0.28, 0.30, 0.26, 0.33}; // beta, Gev^(-1) double beta[] = { 0.28, 0.30, 0.26, 0.33 }; // beta, Gev^(-1)
double Delta[] = {0.04, 0.04, 0.30, 0.30}; // Delta, Gev double Delta[] = { 0.04, 0.04, 0.30, 0.30 }; // Delta, Gev
// form-factors // form-factors
EvtComplex Ftvq0, Ftaq0, Ftv00, Fta00; EvtComplex Ftvq0, Ftaq0, Ftv00, Fta00;
Fv = unit1 * beta[0] * fb * M1 / ( Delta[0] + E ); // Fv(q^2) Fv = unit1 * beta[0] * fb * M1 / ( Delta[0] + E ); // Fv(q^2)
Ftvq0 = unit1 * beta[1] * fb * M1 / ( Delta[1] + E ); // Ftv(q^2,0) Ftvq0 = unit1 * beta[1] * fb * M1 / ( Delta[1] + E ); // Ftv(q^2,0)
Fa = unit1 * beta[2] * fb * M1 / ( Delta[2] + E ); // Fa(q^2) Fa = unit1 * beta[2] * fb * M1 / ( Delta[2] + E ); // Fa(q^2)
Ftaq0 = unit1 * beta[3] * fb * M1 / ( Delta[3] + E ); // Fta(q^2,0) Ftaq0 = unit1 * beta[3] * fb * M1 / ( Delta[3] + E ); // Fta(q^2,0)
Ftv00 = unit1 * beta[1] * fb * M1 / ( Delta[1] + 0.5 * M1 ); // Ftv(0,0) Ftv00 = unit1 * beta[1] * fb * M1 / ( Delta[1] + 0.5 * M1 ); // Ftv(0,0)
▲ Show 20 Lines • Show All 108 Lines • Show Last 20 Lines