Page MenuHomeHEPForge
Differential D87 Diff 389 src/EvtGenModels/EvtbTosllBallFF.cpp

Changeset View

Standalone View

View Options

src/EvtGenModels/EvtbTosllBallFF.cpp

Show All 24 Lines
#include <cmath> #include <cmath>
EvtbTosllBallFF::EvtbTosllBallFF( int ffmodel ) EvtbTosllBallFF::EvtbTosllBallFF( int ffmodel )
{ {
_theFFModel = ffmodel; _theFFModel = ffmodel;
} }
static EvtId IdKp, IdKm, IdKS, IdK0, IdaK0, IdKL, IdPip, IdPim, IdPi0, IdEta,
IdEtap;
static bool sfirst = true;
void EvtbTosllBallFF::getScalarFF( EvtId parent, EvtId daught, double t, void EvtbTosllBallFF::getScalarFF( EvtId parent, EvtId daught, double t,
double /*mass*/, double& fp, double& f0, double /*mass*/, double& fp, double& f0,
double& ft ) double& ft )
{ {
if ( sfirst ) {
sfirst = false;
IdKp = EvtPDL::getId( "K+" );
IdKm = EvtPDL::getId( "K-" );
IdKS = EvtPDL::getId( "K_S0" );
IdK0 = EvtPDL::getId( "K0" );
IdaK0 = EvtPDL::getId( "anti-K0" );
IdKL = EvtPDL::getId( "K_L0" );
IdPip = EvtPDL::getId( "pi+" );
IdPim = EvtPDL::getId( "pi-" );
IdPi0 = EvtPDL::getId( "pi0" );
IdEta = EvtPDL::getId( "eta" );
IdEtap = EvtPDL::getId( "eta'" );
}
int model = _theFFModel; int model = _theFFModel;
double m = EvtPDL::getMeanMass( parent ); double m = EvtPDL::getMeanMass( parent );
double md = EvtPDL::getMeanMass( daught ); double md = EvtPDL::getMeanMass( daught );
double shat = t / ( m * m ); double shat = t / ( m * m );
double shat2 = shat * shat; double shat2 = shat * shat;
double shat3 = shat2 * shat; double shat3 = shat2 * shat;
if ( daught == EvtPDL::getId( std::string( "K+" ) ) || if ( daught == IdKp || daught == IdKm || daught == IdKS || daught == IdK0 ||
daught == EvtPDL::getId( std::string( "K-" ) ) || daught == IdaK0 || daught == IdKL ) {
daught == EvtPDL::getId( std::string( "K_S0" ) ) ||
daught == EvtPDL::getId( std::string( "K0" ) ) ||
daught == EvtPDL::getId( std::string( "anti-K0" ) ) ||
daught == EvtPDL::getId( std::string( "K_L0" ) ) ) {
// B --> K form factors // B --> K form factors
if ( model == 1 ) { if ( model == 1 ) {
//this is Ali-Ball '01 (or really Ali-Ball'99 minimum allowed) //this is Ali-Ball '01 (or really Ali-Ball'99 minimum allowed)
fp = 0.278 * exp( 1.568 * shat + 0.470 * shat2 + 0.885 * shat3 ); fp = 0.278 * exp( 1.568 * shat + 0.470 * shat2 + 0.885 * shat3 );
f0 = 0.278 * exp( 0.740 * shat + 0.080 * shat2 + 0.425 * shat3 ); f0 = 0.278 * exp( 0.740 * shat + 0.080 * shat2 + 0.425 * shat3 );
ft = 0.300 * exp( 1.600 * shat + 0.501 * shat2 + 0.796 * shat3 ); ft = 0.300 * exp( 1.600 * shat + 0.501 * shat2 + 0.796 * shat3 );
} }
if ( model == 2 ) { if ( model == 2 ) {
Show All 32 Lines if ( daught == IdKp || daught == IdKm || daught == IdKS || daught == IdK0 ||
if ( model == 7 ) { if ( model == 7 ) {
// Ball-Zwicky LCSR '05 "set 4" (mb = 4.6) // Ball-Zwicky LCSR '05 "set 4" (mb = 4.6)
fp = ( 0.1903 / ( 1. - ( t / 29.30 ) ) ) + fp = ( 0.1903 / ( 1. - ( t / 29.30 ) ) ) +
( 0.1478 / ( 1. - ( t / 29.30 ) ) / ( 1. - ( t / 29.30 ) ) ); ( 0.1478 / ( 1. - ( t / 29.30 ) ) / ( 1. - ( t / 29.30 ) ) );
f0 = ( 0.3338 / ( 1. - ( t / 38.98 / 38.98 ) ) ); f0 = ( 0.3338 / ( 1. - ( t / 38.98 / 38.98 ) ) );
ft = ( 0.1851 / ( 1. - ( t / 29.30 ) ) ) + ft = ( 0.1851 / ( 1. - ( t / 29.30 ) ) ) +
( 0.1905 / ( 1. - ( t / 29.30 ) ) / ( 1. - ( t / 29.30 ) ) ); ( 0.1905 / ( 1. - ( t / 29.30 ) ) / ( 1. - ( t / 29.30 ) ) );
} }
} else if ( daught == EvtPDL::getId( std::string( "pi+" ) ) || } else if ( daught == IdPip || daught == IdPim || daught == IdPi0 ) {
daught == EvtPDL::getId( std::string( "pi-" ) ) ||
daught == EvtPDL::getId( std::string( "pi0" ) ) ) {
if ( model == 1 ) { if ( model == 1 ) {
// B --> pi form factors from Ball-Zwicky'01 (tabulated in hep-ph/0306251) // B --> pi form factors from Ball-Zwicky'01 (tabulated in hep-ph/0306251)
fp = 0.261 / ( 1. - 2.03 * shat + 1.293 * shat * shat ); fp = 0.261 / ( 1. - 2.03 * shat + 1.293 * shat * shat );
f0 = 0.261 / ( 1. - 0.27 * shat - 0.752 * shat * shat ); f0 = 0.261 / ( 1. - 0.27 * shat - 0.752 * shat * shat );
ft = 0.296 / ( 1. - 1.28 * shat + 0.193 * shat * shat ); ft = 0.296 / ( 1. - 1.28 * shat + 0.193 * shat * shat );
} }
// The following two (2) and (3) should preferably be replaced with // The following two (2) and (3) should preferably be replaced with
Show All 39 Lines if ( daught == IdKp || daught == IdKm || daught == IdKS || daught == IdK0 ||
if ( model == 7 ) { if ( model == 7 ) {
// Ball-Zwicky LCSR '05 "set 4" (mb = 4.6) // Ball-Zwicky LCSR '05 "set 4" (mb = 4.6)
fp = ( 0.944 / ( 1. - ( t / ( 5.32 * 5.32 ) ) ) ) + fp = ( 0.944 / ( 1. - ( t / ( 5.32 * 5.32 ) ) ) ) +
( -0.669 / ( 1. - ( t / 34.27 ) ) ); ( -0.669 / ( 1. - ( t / 34.27 ) ) );
f0 = ( 0.270 / ( 1. - ( t / 33.63 ) ) ); f0 = ( 0.270 / ( 1. - ( t / 33.63 ) ) );
ft = ( 0.152 / ( 1. - ( t / ( 5.32 * 5.32 ) ) ) ) + ft = ( 0.152 / ( 1. - ( t / ( 5.32 * 5.32 ) ) ) ) +
( 0.122 / ( 1. - ( t / 28.40 ) ) / ( 1. - ( t / 28.40 ) ) ); ( 0.122 / ( 1. - ( t / 28.40 ) ) / ( 1. - ( t / 28.40 ) ) );
} }
} else if ( daught == EvtPDL::getId( std::string( "eta" ) ) || } else if ( daught == IdEta || daught == IdEtap ) {
daught == EvtPDL::getId( std::string( "eta'" ) ) ) {
if ( model == 1 ) { if ( model == 1 ) {
// B --> eta form factors // B --> eta form factors
fp = 0.261 / ( 1. - 2.03 * shat + 1.293 * shat * shat ); fp = 0.261 / ( 1. - 2.03 * shat + 1.293 * shat * shat );
f0 = 0.261 / ( 1. - 0.27 * shat - 0.752 * shat * shat ); f0 = 0.261 / ( 1. - 0.27 * shat - 0.752 * shat * shat );
ft = 0.296 / ( 1. - 1.28 * shat + 0.193 * shat * shat ); ft = 0.296 / ( 1. - 1.28 * shat + 0.193 * shat * shat );
} }
if ( model == 6 ) { if ( model == 6 ) {
// Ball-Zwicky LCSR '05 "set 2" (mb = 4.8) // Ball-Zwicky LCSR '05 "set 2" (mb = 4.8)
Show All 12 Lines if ( daught == IdKp || daught == IdKm || daught == IdKS || daught == IdK0 ||
( 0.1841 / ( 1. - ( t / 28.40 ) ) / ( 1. - ( t / 28.40 ) ) ); ( 0.1841 / ( 1. - ( t / 28.40 ) ) / ( 1. - ( t / 28.40 ) ) );
} }
} }
// cout << "shat "<<shat<<"\t"<<"fp "<<fp<<"\t"<<"f0 "<<f0<<"\t" // cout << "shat "<<shat<<"\t"<<"fp "<<fp<<"\t"<<"f0 "<<f0<<"\t"
// <<"ft "<<ft<<endl; // <<"ft "<<ft<<endl;
} }
static EvtId IdBs, IdaBs, IdRhop, IdRhom, IdRho0, IdOmega, IdKst0, IdaKst0,
IdKstp, IdKstm, IdPhi, IdK1p, IdK1m;
static bool vfirst = true;
void EvtbTosllBallFF::getVectorFF( EvtId parent, EvtId daught, double t, void EvtbTosllBallFF::getVectorFF( EvtId parent, EvtId daught, double t,
double /*mass*/, double& a1, double& a2, double /*mass*/, double& a1, double& a2,
double& a0, double& v, double& t1, double& a0, double& v, double& t1,
double& t2, double& t3 ) double& t2, double& t3 )
{ {
if ( vfirst ) {
vfirst = false;
IdBs = EvtPDL::getId( "B_s0" );
IdaBs = EvtPDL::getId( "anti-B_s0" );
IdRhop = EvtPDL::getId( "rho+" );
IdRhom = EvtPDL::getId( "rho-" );
IdRho0 = EvtPDL::getId( "rho0" );
IdOmega = EvtPDL::getId( "omega" );
IdKst0 = EvtPDL::getId( "K*0" );
IdaKst0 = EvtPDL::getId( "anti-K*0" );
IdKstp = EvtPDL::getId( "K*+" );
IdKstm = EvtPDL::getId( "K*-" );
IdPhi = EvtPDL::getId( "phi" );
IdK1p = EvtPDL::getId( "K_1+" );
IdK1m = EvtPDL::getId( "K_1-" );
}
int model = _theFFModel; int model = _theFFModel;
double m = EvtPDL::getMeanMass( parent ); double m = EvtPDL::getMeanMass( parent );
double md = EvtPDL::getMeanMass( daught ); double md = EvtPDL::getMeanMass( daught );
double shat = t / ( m * m ); double shat = t / ( m * m );
double shat2 = shat * shat; double shat2 = shat * shat;
if ( if ( parent == IdBs || parent == IdaBs ) { // B_s decay form factors
// B_s decay form factors if ( daught == IdKst0 || daught == IdaKst0 ) { // B_s --> K* form factors
parent == EvtPDL::getId( std::string( "B_s0" ) ) ||
parent == EvtPDL::getId( std::string( "anti-B_s0" ) ) ) {
if (
// B_s --> K* form factors
daught == EvtPDL::getId( std::string( "K*0" ) ) ||
daught == EvtPDL::getId( std::string( "anti-K*0" ) ) ) {
if ( model == 6 ) { if ( model == 6 ) {
// Ball-Zwicky LCSR '05 (mb = 4.8) // Ball-Zwicky LCSR '05 (mb = 4.8)
a1 = 0.231 / ( 1. - ( t / 32.94 ) ); a1 = 0.231 / ( 1. - ( t / 32.94 ) );
a2 = ( -0.011 / ( 1. - ( t / 40.14 ) ) ) + a2 = ( -0.011 / ( 1. - ( t / 40.14 ) ) ) +
( 0.192 / ( 1. - ( t / 40.14 ) ) / ( 1. - ( t / 40.14 ) ) ); ( 0.192 / ( 1. - ( t / 40.14 ) ) / ( 1. - ( t / 40.14 ) ) );
a0 = ( 2.813 / ( 1. - ( t / ( 5.37 * 5.37 ) ) ) ) + a0 = ( 2.813 / ( 1. - ( t / ( 5.37 * 5.37 ) ) ) ) +
( -2.509 / ( 1. - ( t / 31.58 ) ) ); ( -2.509 / ( 1. - ( t / 31.58 ) ) );
v = ( 2.351 / ( 1. - ( t / ( 5.42 * 5.42 ) ) ) ) + v = ( 2.351 / ( 1. - ( t / ( 5.42 * 5.42 ) ) ) ) +
( -2.039 / ( 1. - ( t / 33.10 ) ) ); ( -2.039 / ( 1. - ( t / 33.10 ) ) );
t1 = ( 2.047 / ( 1. - ( t / ( 5.42 * 5.42 ) ) ) ) + t1 = ( 2.047 / ( 1. - ( t / ( 5.42 * 5.42 ) ) ) ) +
( -1.787 / ( 1 - ( t / 32.83 ) ) ); ( -1.787 / ( 1 - ( t / 32.83 ) ) );
t2 = 0.260 / ( 1. - ( t / 33.01 ) ); t2 = 0.260 / ( 1. - ( t / 33.01 ) );
double t3tilde = ( 0.043 / ( 1. - ( t / 39.38 ) ) ) + double t3tilde = ( 0.043 / ( 1. - ( t / 39.38 ) ) ) +
( 0.217 / ( 1. - ( t / 39.38 ) ) / ( 0.217 / ( 1. - ( t / 39.38 ) ) /
( 1. - ( t / 39.38 ) ) ); ( 1. - ( t / 39.38 ) ) );
t3 = 0.0; t3 = 0.0;
if ( fabs( t ) > 1e-10 ) { if ( fabs( t ) > 1e-10 ) {
t3 = ( m * m - md * md ) * ( t3tilde - t2 ) / t; t3 = ( m * m - md * md ) * ( t3tilde - t2 ) / t;
} }
} }
} }
else if ( else if ( daught == IdPhi ) { // B_s --> phi form factors
// B_s --> phi form factors
daught == EvtPDL::getId( std::string( "phi" ) ) ) {
if ( model == 6 ) { if ( model == 6 ) {
// Ball-Zwicky LCSR '05 (mb = 4.8) // Ball-Zwicky LCSR '05 (mb = 4.8)
a1 = 0.308 / ( 1. - ( t / 36.54 ) ); a1 = 0.308 / ( 1. - ( t / 36.54 ) );
a2 = ( -0.054 / ( 1. - ( t / 48.94 ) ) ) + a2 = ( -0.054 / ( 1. - ( t / 48.94 ) ) ) +
( 0.288 / ( 1. - ( t / 48.94 ) ) / ( 1. - ( t / 48.94 ) ) ); ( 0.288 / ( 1. - ( t / 48.94 ) ) / ( 1. - ( t / 48.94 ) ) );
a0 = ( 3.310 / ( 1. - ( t / ( 5.28 * 5.28 ) ) ) ) + a0 = ( 3.310 / ( 1. - ( t / ( 5.28 * 5.28 ) ) ) ) +
( -2.835 / ( 1. - ( t / 31.57 ) ) ); ( -2.835 / ( 1. - ( t / 31.57 ) ) );
v = ( 1.484 / ( 1. - ( t / ( 5.32 * 5.32 ) ) ) ) + v = ( 1.484 / ( 1. - ( t / ( 5.32 * 5.32 ) ) ) ) +
( -1.049 / ( 1. - ( t / 39.52 ) ) ); ( -1.049 / ( 1. - ( t / 39.52 ) ) );
t1 = ( 1.303 / ( 1. - ( t / ( 5.32 * 5.32 ) ) ) ) + t1 = ( 1.303 / ( 1. - ( t / ( 5.32 * 5.32 ) ) ) ) +
( -0.954 / ( 1 - ( t / 38.28 ) ) ); ( -0.954 / ( 1 - ( t / 38.28 ) ) );
t2 = 0.348 / t2 = 0.348 /
( 1. - ( t / 37.21 ) ); //0.349 caused a pole in T3 for q2 -> 0. 0.348 has correct behaviour ( 1. - ( t / 37.21 ) ); //0.349 caused a pole in T3 for q2 -> 0. 0.348 has correct behaviour
double t3tilde = ( 0.027 / ( 1. - ( t / 45.56 ) ) ) + double t3tilde = ( 0.027 / ( 1. - ( t / 45.56 ) ) ) +
( 0.321 / ( 1. - ( t / 45.56 ) ) / ( 0.321 / ( 1. - ( t / 45.56 ) ) /
( 1. - ( t / 45.56 ) ) ); ( 1. - ( t / 45.56 ) ) );
t3 = 0.0; t3 = 0.0;
if ( fabs( t ) > 1e-10 ) { if ( fabs( t ) > 1e-10 ) {
t3 = ( m * m - md * md ) * ( t3tilde - t2 ) / t; t3 = ( m * m - md * md ) * ( t3tilde - t2 ) / t;
} }
} }
} }
} else if ( daught == EvtPDL::getId( std::string( "K*+" ) ) || } else if ( daught == IdKstp || daught == IdKstm || daught == IdKst0 ||
daught == EvtPDL::getId( std::string( "K*-" ) ) || daught == IdaKst0 || daught == IdK1p || daught == IdK1m ) {
daught == EvtPDL::getId( std::string( "K*0" ) ) ||
daught == EvtPDL::getId( std::string( "anti-K*0" ) ) ||
daught == EvtPDL::getId( std::string( "K_1+" ) ) ||
daught == EvtPDL::getId( std::string( "K_1-" ) )
) {
if ( model == 1 ) { if ( model == 1 ) {
//this is Ali-Ball '01 //this is Ali-Ball '01
a1 = 0.294 * exp( 0.656 * shat + 0.456 * shat2 ); a1 = 0.294 * exp( 0.656 * shat + 0.456 * shat2 );
a2 = 0.246 * exp( 1.237 * shat + 0.822 * shat2 ); a2 = 0.246 * exp( 1.237 * shat + 0.822 * shat2 );
a0 = 0.412 * exp( 1.543 * shat + 0.954 * shat2 ); a0 = 0.412 * exp( 1.543 * shat + 0.954 * shat2 );
v = 0.399 * exp( 1.537 * shat + 1.123 * shat2 ); v = 0.399 * exp( 1.537 * shat + 1.123 * shat2 );
t1 = 0.334 * exp( 1.575 * shat + 1.140 * shat2 ); t1 = 0.334 * exp( 1.575 * shat + 1.140 * shat2 );
▲ Show 20 Lines • Show All 95 Lines • ▼ Show 20 Lines if ( parent == IdBs || parent == IdaBs ) { // B_s decay form factors
t1 = ( 0.31 / ( 1 - ( t / ( 5.412 * 5.412 ) ) ) ) * t1 = ( 0.31 / ( 1 - ( t / ( 5.412 * 5.412 ) ) ) ) *
( 1 + ( -4.6 ) * ( z - z0 + 0.5 * ( ( z * z ) - ( z0 * z0 ) ) ) ); ( 1 + ( -4.6 ) * ( z - z0 + 0.5 * ( ( z * z ) - ( z0 * z0 ) ) ) );
t2 = ( 0.31 / ( 1 - ( t / ( 5.829 * 5.829 ) ) ) ) * t2 = ( 0.31 / ( 1 - ( t / ( 5.829 * 5.829 ) ) ) ) *
( 1 + ( -3.2 ) * ( z - z0 + 0.5 * ( ( z * z ) - ( z0 * z0 ) ) ) ); ( 1 + ( -3.2 ) * ( z - z0 + 0.5 * ( ( z * z ) - ( z0 * z0 ) ) ) );
t3 = ( 0.22 / ( 1 - ( t / ( 5.829 * 5.829 ) ) ) ) * t3 = ( 0.22 / ( 1 - ( t / ( 5.829 * 5.829 ) ) ) ) *
( 1 + ( 1 +
( -10.3 ) * ( z - z0 + 0.5 * ( ( z * z ) - ( z0 * z0 ) ) ) ); ( -10.3 ) * ( z - z0 + 0.5 * ( ( z * z ) - ( z0 * z0 ) ) ) );
} }
} else if ( daught == EvtPDL::getId( std::string( "rho+" ) ) || } else if ( daught == IdRhop || daught == IdRhom || daught == IdRho0 ) {
daught == EvtPDL::getId( std::string( "rho-" ) ) ||
daught == EvtPDL::getId( std::string( "rho0" ) ) ) {
if ( model == 1 ) { if ( model == 1 ) {
// B --> rho form factors // B --> rho form factors
a1 = 0.261 / ( 1. - 0.29 * shat - 0.415 * shat * shat ); a1 = 0.261 / ( 1. - 0.29 * shat - 0.415 * shat * shat );
a2 = 0.223 / ( 1. - 0.93 * shat - 0.092 * shat * shat ); a2 = 0.223 / ( 1. - 0.93 * shat - 0.092 * shat * shat );
a0 = 0.372 / ( 1. - 1.40 * shat + 0.437 * shat * shat ); a0 = 0.372 / ( 1. - 1.40 * shat + 0.437 * shat * shat );
v = 0.338 / ( 1. - 1.37 * shat + 0.315 * shat * shat ); v = 0.338 / ( 1. - 1.37 * shat + 0.315 * shat * shat );
t1 = 0.285 / ( 1. - 1.41 * shat + 0.361 * shat * shat ); t1 = 0.285 / ( 1. - 1.41 * shat + 0.361 * shat * shat );
t2 = 0.285 / ( 1. - 0.28 * shat - 0.500 * shat * shat ); t2 = 0.285 / ( 1. - 0.28 * shat - 0.500 * shat * shat );
Show All 15 Lines if ( parent == IdBs || parent == IdaBs ) { // B_s decay form factors
double t3tilde = ( 0.022 / ( 1. - ( t / 40.88 ) ) ) + double t3tilde = ( 0.022 / ( 1. - ( t / 40.88 ) ) ) +
( 0.246 / ( 1. - ( t / 40.88 ) ) / ( 0.246 / ( 1. - ( t / 40.88 ) ) /
( 1. - ( t / 40.88 ) ) ); ( 1. - ( t / 40.88 ) ) );
t3 = 0.0; t3 = 0.0;
if ( fabs( t ) > 1e-10 ) { if ( fabs( t ) > 1e-10 ) {
t3 = ( m * m - md * md ) * ( t3tilde - t2 ) / t; t3 = ( m * m - md * md ) * ( t3tilde - t2 ) / t;
} }
} }
} else if ( daught == EvtPDL::getId( std::string( "omega" ) ) ) { } else if ( daught == IdOmega ) {
if ( model == 1 ) { if ( model == 1 ) {
// B --> rho form factors // B --> rho form factors
a1 = 0.261 / ( 1. - 0.29 * shat - 0.415 * shat * shat ); a1 = 0.261 / ( 1. - 0.29 * shat - 0.415 * shat * shat );
a2 = 0.223 / ( 1. - 0.93 * shat - 0.092 * shat * shat ); a2 = 0.223 / ( 1. - 0.93 * shat - 0.092 * shat * shat );
a0 = 0.372 / ( 1. - 1.40 * shat + 0.437 * shat * shat ); a0 = 0.372 / ( 1. - 1.40 * shat + 0.437 * shat * shat );
v = 0.338 / ( 1. - 1.37 * shat + 0.315 * shat * shat ); v = 0.338 / ( 1. - 1.37 * shat + 0.315 * shat * shat );
t1 = 0.285 / ( 1. - 1.41 * shat + 0.361 * shat * shat ); t1 = 0.285 / ( 1. - 1.41 * shat + 0.361 * shat * shat );
t2 = 0.285 / ( 1. - 0.28 * shat - 0.500 * shat * shat ); t2 = 0.285 / ( 1. - 0.28 * shat - 0.500 * shat * shat );
Show All 25 Lines