Changeset View
Changeset View
Standalone View
Standalone View
src/EvtGenModels/EvtBBScalar.cpp
Show All 36 Lines | |||||
const EvtComplex EvtBBScalar::V_us_star = EvtComplex( 0.22, 0 ); | const EvtComplex EvtBBScalar::V_us_star = EvtComplex( 0.22, 0 ); | ||||
const EvtComplex EvtBBScalar::a1 = EvtComplex( 1.05, 0 ); | const EvtComplex EvtBBScalar::a1 = EvtComplex( 1.05, 0 ); | ||||
const EvtComplex EvtBBScalar::V_tb = EvtComplex( 0.99915, 0 ); | const EvtComplex EvtBBScalar::V_tb = EvtComplex( 0.99915, 0 ); | ||||
const EvtComplex EvtBBScalar::V_ts_star = | const EvtComplex EvtBBScalar::V_ts_star = | ||||
EvtComplex( -0.04029 - 0.000813 * cos( 60 / 180 * pi ), | EvtComplex( -0.04029 - 0.000813 * cos( 60 / 180 * pi ), | ||||
-0.000813 * cos( 60 / 180 * pi ) ); | -0.000813 * cos( 60 / 180 * pi ) ); | ||||
const EvtComplex EvtBBScalar::a4 = EvtComplex( -387.3e-4, -121e-4 ); | const EvtComplex EvtBBScalar::a4 = EvtComplex( -387.3e-4, -121e-4 ); | ||||
const EvtComplex EvtBBScalar::a6 = EvtComplex( -555.3e-4, -121e-4 ); | const EvtComplex EvtBBScalar::a6 = EvtComplex( -555.3e-4, -121e-4 ); | ||||
const double EvtBBScalar::x[] = {420.96, -10485.50, 100639.97, -433916.61, | const double EvtBBScalar::x[] = { 420.96, -10485.50, 100639.97, -433916.61, | ||||
613780.15}; | 613780.15 }; | ||||
const double EvtBBScalar::y[] = {292.62, -735.73}; | const double EvtBBScalar::y[] = { 292.62, -735.73 }; | ||||
const double EvtBBScalar::m_s = 0.120; | const double EvtBBScalar::m_s = 0.120; | ||||
const double EvtBBScalar::m_u = 0.029 * 0.120; | const double EvtBBScalar::m_u = 0.029 * 0.120; | ||||
const double EvtBBScalar::m_b = 4.88; | const double EvtBBScalar::m_b = 4.88; | ||||
EvtBBScalar::EvtBBScalar() : EvtDecayAmp(), _massRatio( 0 ), _baryonMassSum( 0 ) | EvtBBScalar::EvtBBScalar() : EvtDecayAmp(), _massRatio( 0 ), _baryonMassSum( 0 ) | ||||
{ | { | ||||
FormFactor dummy; | FormFactor dummy; | ||||
dummy.value = 0.36; | dummy.value = 0.36; | ||||
▲ Show 20 Lines • Show All 183 Lines • ▼ Show 20 Lines | if ( _baryonCombination.test( Lambda ) and _baryonCombination.test( Proton ) ) { | ||||
EvtGenReport( EVTGEN_ERROR, "EvtGen" ) | EvtGenReport( EVTGEN_ERROR, "EvtGen" ) | ||||
<< "EvtBBScalar::baryon_fS: Don't know what to do with this type as the first or second baryon\n"; | << "EvtBBScalar::baryon_fS: Don't know what to do with this type as the first or second baryon\n"; | ||||
exit( 2 ); | exit( 2 ); | ||||
} | } | ||||
} | } | ||||
double EvtBBScalar::D_A( double t ) const | double EvtBBScalar::D_A( double t ) const | ||||
{ | { | ||||
const double d_tilde[] = {x[0] - 1.5 * y[0], -478}; | const double d_tilde[] = { x[0] - 1.5 * y[0], -478 }; | ||||
const vector<double> v_d_tilde( d_tilde, d_tilde + 2 ); | const vector<double> v_d_tilde( d_tilde, d_tilde + 2 ); | ||||
return formFactorFit( t, v_d_tilde ); | return formFactorFit( t, v_d_tilde ); | ||||
} | } | ||||
double EvtBBScalar::F_A( double t ) const | double EvtBBScalar::F_A( double t ) const | ||||
{ | { | ||||
const double f_tilde[] = {2. / 3 * x[0] + 0.5 * y[0], -478}; | const double f_tilde[] = { 2. / 3 * x[0] + 0.5 * y[0], -478 }; | ||||
const vector<double> v_f_tilde( f_tilde, f_tilde + 2 ); | const vector<double> v_f_tilde( f_tilde, f_tilde + 2 ); | ||||
return formFactorFit( t, v_f_tilde ); | return formFactorFit( t, v_f_tilde ); | ||||
} | } | ||||
double EvtBBScalar::D_P( double t ) const | double EvtBBScalar::D_P( double t ) const | ||||
{ | { | ||||
const double d_bar[] = {1.5 * y[0] * _massRatio, /*-952*/ 0}; | const double d_bar[] = { 1.5 * y[0] * _massRatio, /*-952*/ 0 }; | ||||
const vector<double> v_d_bar( d_bar, d_bar + 2 ); | const vector<double> v_d_bar( d_bar, d_bar + 2 ); | ||||
return formFactorFit( t, v_d_bar ); | return formFactorFit( t, v_d_bar ); | ||||
} | } | ||||
double EvtBBScalar::F_P( double t ) const | double EvtBBScalar::F_P( double t ) const | ||||
{ | { | ||||
const double f_bar[] = {( x[0] - 0.5 * y[0] ) * _massRatio, /*-952*/ 0}; | const double f_bar[] = { ( x[0] - 0.5 * y[0] ) * _massRatio, /*-952*/ 0 }; | ||||
const vector<double> v_f_bar( f_bar, f_bar + 2 ); | const vector<double> v_f_bar( f_bar, f_bar + 2 ); | ||||
return formFactorFit( t, v_f_bar ); | return formFactorFit( t, v_f_bar ); | ||||
} | } | ||||
double EvtBBScalar::D_S( double t ) const | double EvtBBScalar::D_S( double t ) const | ||||
{ | { | ||||
return -1.5 * _massRatio * G_n( t ); | return -1.5 * _massRatio * G_n( t ); | ||||
} | } | ||||
▲ Show 20 Lines • Show All 221 Lines • Show Last 20 Lines |