12 #include "EvtGenBase/EvtReport.hh"
13 #include "EvtGenBase/EvtId.hh"
14 #include "EvtGenBase/EvtPDL.hh"
18 #include "generators/evtgen/models/EvtBGLFF.h"
21 EvtBGLFF::EvtBGLFF(
double bglap_0,
double bglap_1,
double bglap_2,
double bglap_3,
double bgla0_0,
double bgla0_1,
double bgla0_2,
40 EvtBGLFF::EvtBGLFF(
double bgla_0,
double bgla_1,
double bglb_0,
double bglb_1,
double bglc_1,
double bglc_2)
62 double t,
double mass,
double* fp,
double* f0)
65 double mb = EvtPDL::getMeanMass(parent);
67 double w = ((mb * mb) + (mass * mass) - t) / (2.0 * mb * mass);
68 const double z = (sqrt(w + 1) - sqrt(2.)) / (sqrt(w + 1) + sqrt(2.));
69 const double p_i = 1.0;
71 double phi_sub = (1 + r) * (1 - z) + 2.0 * sqrt(r) * (1 + z);
73 double g_sub = (4.0 * r) / pow(1 + r, 2);
75 double phi_p = 1.1213 * pow(1 + z, 2.) * sqrt(1 - z) * pow(phi_sub, -5);
77 double phi_0 = 0.5299 * (1 + z) * pow(1 - z, 1.5) * pow(phi_sub, -4);
80 *fp = g_sub * (
ap_0 +
ap_1 * z +
ap_2 * z * z +
ap_3 * z * z * z) / (p_i * phi_p);
82 *f0 = g_sub * (
a0_0 +
a0_1 * z +
a0_2 * z * z +
a0_3 * z * z * z) / (p_i * phi_0);
93 double* a2f,
double* vf,
double* a0f)
97 double mb = EvtPDL::getMeanMass(parent);
98 double w = ((mb * mb) + (mass * mass) - t) / (2. * mb * mass);
103 const double r = mass / mb;
104 double z = (sqrt(w + 1.) - sqrt(2.)) / (sqrt(w + 1) + sqrt(2.));
105 const double rstar = (2. * sqrt(mb * mass)) / (mb + mass);
106 const double chiT_plus33 = 5.28e-4;
107 const double chiT_minus33 = 3.07e-4;
108 const double n_i = 2.6;
109 const double axialvector_poles[4] = {6.730, 6.736, 7.135, 7.142};
110 const double vector_poles[4] = {6.337, 6.899, 7.012, 7.280};
112 const double c_0 = (mb - mass) / mb * sqrt(0.5) / (1 + r + 2. * sqrt(r)) *
b_0;
114 double phi_g = sqrt(256. * n_i / (3. * M_PI * chiT_plus33)) * r * r * (1. + z) * (1. + z) / sqrt(1. - z) / pow((1. + r) *
115 (1. - z) + 2. * sqrt(r) * (1. + z), 4.);
116 double phi_f = 1. / (mb * mb) * sqrt(16. * n_i / (3. * M_PI * chiT_minus33)) * r * (1. + z) * pow(1. - z ,
117 1.5) / pow((1. + r) * (1. - z) + 2. * sqrt(r) * (1. + z), 4.);
118 double phi_F1 = 1. / (mb * mb * mb) * sqrt(8. * n_i / (3. * M_PI * chiT_minus33)) * r * (1. + z) * pow(1. - z ,
119 2.5) / pow((1. + r) * (1. - z) + 2 * sqrt(r) * (1. + z), 5.);
123 const double term3 = sqrt((mb + mass) * (mb + mass) - (mb - mass) * (mb - mass));
124 for (
int i = 0; i < 4; i++) {
125 double term1 = sqrt((mb + mass) * (mb + mass) - vector_poles[i] * vector_poles[i]);
126 double term2 = sqrt((mb + mass) * (mb + mass) - axialvector_poles[i] * axialvector_poles[i]);
127 double z_p1 = (term1 - term3) / (term1 + term3);
128 p_g = p_g * (z - z_p1) / (1 - z * z_p1);
129 double z_p2 = (term2 - term3) / (term2 + term3);
130 p_f = p_f * (z - z_p2) / (1 - z * z_p2);
133 double g = 1. / p_g / phi_g * (
a_0 +
a_1 * z);
134 double f = 1. / p_f / phi_f * (
b_0 +
b_1 * z);
135 double F1 = 1. / p_f / phi_F1 * (c_0 +
c_1 * z +
c_2 * z * z);
137 double ha1 = f / sqrt(mb * mass) / (1. + w);
138 double r1 = (w + 1.) * mb * mass * g / f;
139 double r2 = (w - r) / (w - 1) - F1 / mb / (w - 1) / f;
141 *a1f = (w + 1.) / 2. * rstar * ha1;
142 *a2f = (r2 / rstar) * ha1;
143 *vf = (r1 / rstar) * ha1;
151 double*,
double*,
double*)
154 EvtGenReport(EVTGEN_ERROR,
"EvtGen") <<
"Not implemented :gettensorff in EvtBGLFF.\n";
162 double*,
double*,
double*)
165 EvtGenReport(EVTGEN_ERROR,
"EvtGen") <<
"Not implemented :getbaryonff in EvtBGLFF.\n";
171 double*,
double*,
double*,
double*)
174 EvtGenReport(EVTGEN_ERROR,
"EvtGen") <<
"Not implemented :getdiracff in EvtBGLFF.\n";
180 double*,
double*,
double*,
double*,
double*,
double*)
183 EvtGenReport(EVTGEN_ERROR,
"EvtGen") <<
"Not implemented :getraritaff in EvtBGLFF.\n";