Changeset View
Changeset View
Standalone View
Standalone View
src/EvtGenModels/EvtbTosllMSFF.cpp
Show First 20 Lines • Show All 72 Lines • ▼ Show 20 Lines | if ( ( parent == EvtPDL::getId( std::string( "B+" ) ) && | ||||
( parent == EvtPDL::getId( std::string( "B0" ) ) && | ( parent == EvtPDL::getId( std::string( "B0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "K_S0" ) ) ) || | daught == EvtPDL::getId( std::string( "K_S0" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "anti-B0" ) ) && | ( parent == EvtPDL::getId( std::string( "anti-B0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "K_S0" ) ) ) || | daught == EvtPDL::getId( std::string( "K_S0" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "B0" ) ) && | ( parent == EvtPDL::getId( std::string( "B0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "K_L0" ) ) ) || | daught == EvtPDL::getId( std::string( "K_L0" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "anti-B0" ) ) && | ( parent == EvtPDL::getId( std::string( "anti-B0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "K_L0" ) ) ) ) { | daught == EvtPDL::getId( std::string( "K_L0" ) ) ) ) { | ||||
double ff0[] = {0.36, 0.36, 0.35}; | double ff0[] = { 0.36, 0.36, 0.35 }; | ||||
double sigma1[] = {0.43, 0.70, 0.43}; | double sigma1[] = { 0.43, 0.70, 0.43 }; | ||||
double sigma2[] = {0.00, 0.27, 0.00}; | double sigma2[] = { 0.00, 0.27, 0.00 }; | ||||
int eq_num[] = {9, 10, 9}; | int eq_num[] = { 9, 10, 9 }; | ||||
double M_P2 = 5.37 * 5.37; // GeV^2 for B^0_s - meson | double M_P2 = 5.37 * 5.37; // GeV^2 for B^0_s - meson | ||||
double M_V2 = 5.42 * 5.42; // GeV^2 for B^*_s - meson | double M_V2 = 5.42 * 5.42; // GeV^2 for B^*_s - meson | ||||
fp = equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | fp = equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | ||||
f0 = equation9_10( ff0[1], M_V2, t, sigma1[1], sigma2[1], eq_num[1] ); | f0 = equation9_10( ff0[1], M_V2, t, sigma1[1], sigma2[1], eq_num[1] ); | ||||
ft = equation9_10( ff0[2], M_P2, t, sigma1[2], sigma2[2], eq_num[2] ); | ft = equation9_10( ff0[2], M_P2, t, sigma1[2], sigma2[2], eq_num[2] ); | ||||
models_counter = models_counter + 1; | models_counter = models_counter + 1; | ||||
// EvtGenReport(EVTGEN_NOTICE,"EvtGen") <<"\n The function EvtbTosllMSFF::getVectorFF(...) passed." | // EvtGenReport(EVTGEN_NOTICE,"EvtGen") <<"\n The function EvtbTosllMSFF::getVectorFF(...) passed." | ||||
// << "\n B -> K transition form factors" | // << "\n B -> K transition form factors" | ||||
// << std::endl; | // << std::endl; | ||||
} | } | ||||
// B -> \pi transition form factors | // B -> \pi transition form factors | ||||
if ( ( parent == EvtPDL::getId( std::string( "B+" ) ) && | if ( ( parent == EvtPDL::getId( std::string( "B+" ) ) && | ||||
daught == EvtPDL::getId( std::string( "pi+" ) ) ) || | daught == EvtPDL::getId( std::string( "pi+" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "B-" ) ) && | ( parent == EvtPDL::getId( std::string( "B-" ) ) && | ||||
daught == EvtPDL::getId( std::string( "pi-" ) ) ) || | daught == EvtPDL::getId( std::string( "pi-" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "B0" ) ) && | ( parent == EvtPDL::getId( std::string( "B0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "pi0" ) ) ) || | daught == EvtPDL::getId( std::string( "pi0" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "anti-B0" ) ) && | ( parent == EvtPDL::getId( std::string( "anti-B0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "pi0" ) ) ) ) { | daught == EvtPDL::getId( std::string( "pi0" ) ) ) ) { | ||||
double ff0[] = {0.29, 0.29, 0.28}; | double ff0[] = { 0.29, 0.29, 0.28 }; | ||||
double sigma1[] = {0.48, 0.76, 0.48}; | double sigma1[] = { 0.48, 0.76, 0.48 }; | ||||
double sigma2[] = {0.00, 0.28, 0.00}; | double sigma2[] = { 0.00, 0.28, 0.00 }; | ||||
int eq_num[] = {9, 10, 9}; | int eq_num[] = { 9, 10, 9 }; | ||||
double M_P2 = 5.27 * 5.27; // GeV^2 for B^0 - meson | double M_P2 = 5.27 * 5.27; // GeV^2 for B^0 - meson | ||||
double M_V2 = 5.32 * 5.32; // GeV^2 for B^* - meson | double M_V2 = 5.32 * 5.32; // GeV^2 for B^* - meson | ||||
fp = equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | fp = equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | ||||
f0 = equation9_10( ff0[1], M_V2, t, sigma1[1], sigma2[1], eq_num[1] ); | f0 = equation9_10( ff0[1], M_V2, t, sigma1[1], sigma2[1], eq_num[1] ); | ||||
ft = equation9_10( ff0[2], M_P2, t, sigma1[2], sigma2[2], eq_num[2] ); | ft = equation9_10( ff0[2], M_P2, t, sigma1[2], sigma2[2], eq_num[2] ); | ||||
models_counter = models_counter + 1; | models_counter = models_counter + 1; | ||||
// EvtGenReport(EVTGEN_NOTICE,"EvtGen") <<"\n The function EvtbTosllMSFF::getVectorFF(...) passed." | // EvtGenReport(EVTGEN_NOTICE,"EvtGen") <<"\n The function EvtbTosllMSFF::getVectorFF(...) passed." | ||||
// << "\n B -> pi transition form factors" | // << "\n B -> pi transition form factors" | ||||
// << std::endl; | // << std::endl; | ||||
} | } | ||||
// B_d -> \eta transition form factors | // B_d -> \eta transition form factors | ||||
if ( ( parent == EvtPDL::getId( std::string( "B0" ) ) && | if ( ( parent == EvtPDL::getId( std::string( "B0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "eta" ) ) ) || | daught == EvtPDL::getId( std::string( "eta" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "anti-B0" ) ) && | ( parent == EvtPDL::getId( std::string( "anti-B0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "eta" ) ) ) ) { | daught == EvtPDL::getId( std::string( "eta" ) ) ) ) { | ||||
double ff0[] = {0.36, 0.36, 0.36}; | double ff0[] = { 0.36, 0.36, 0.36 }; | ||||
double sigma1[] = {0.60, 0.80, 0.58}; | double sigma1[] = { 0.60, 0.80, 0.58 }; | ||||
double sigma2[] = {0.20, 0.40, 0.18}; | double sigma2[] = { 0.20, 0.40, 0.18 }; | ||||
int eq_num[] = {9, 10, 9}; | int eq_num[] = { 9, 10, 9 }; | ||||
double M_P2 = 5.27 * 5.27; // GeV^2 for B_d^0 - meson | double M_P2 = 5.27 * 5.27; // GeV^2 for B_d^0 - meson | ||||
double M_V2 = 5.32 * 5.32; // GeV^2 for B_d^* - meson | double M_V2 = 5.32 * 5.32; // GeV^2 for B_d^* - meson | ||||
fp = equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | fp = equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | ||||
fp = -0.5 * fp; | fp = -0.5 * fp; | ||||
f0 = equation9_10( ff0[1], M_V2, t, sigma1[1], sigma2[1], eq_num[1] ); | f0 = equation9_10( ff0[1], M_V2, t, sigma1[1], sigma2[1], eq_num[1] ); | ||||
f0 = -0.5 * f0; | f0 = -0.5 * f0; | ||||
ft = equation9_10( ff0[2], M_P2, t, sigma1[2], sigma2[2], eq_num[2] ); | ft = equation9_10( ff0[2], M_P2, t, sigma1[2], sigma2[2], eq_num[2] ); | ||||
ft = -0.5 * ft; | ft = -0.5 * ft; | ||||
models_counter = models_counter + 1; | models_counter = models_counter + 1; | ||||
// EvtGenReport(EVTGEN_NOTICE,"EvtGen") <<"\n The function EvtbTosllMSFF::getVectorFF(...) passed." | // EvtGenReport(EVTGEN_NOTICE,"EvtGen") <<"\n The function EvtbTosllMSFF::getVectorFF(...) passed." | ||||
// << "\n Bd -> eta transition form factors" | // << "\n Bd -> eta transition form factors" | ||||
// << std::endl; | // << std::endl; | ||||
} | } | ||||
// B_d -> \eta' transition form factors | // B_d -> \eta' transition form factors | ||||
if ( ( parent == EvtPDL::getId( std::string( "B0" ) ) && | if ( ( parent == EvtPDL::getId( std::string( "B0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "eta'" ) ) ) || | daught == EvtPDL::getId( std::string( "eta'" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "anti-B0" ) ) && | ( parent == EvtPDL::getId( std::string( "anti-B0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "eta'" ) ) ) ) { | daught == EvtPDL::getId( std::string( "eta'" ) ) ) ) { | ||||
double ff0[] = {0.36, 0.36, 0.39}; | double ff0[] = { 0.36, 0.36, 0.39 }; | ||||
double sigma1[] = {0.60, 0.80, 0.58}; | double sigma1[] = { 0.60, 0.80, 0.58 }; | ||||
double sigma2[] = {0.20, 0.45, 0.18}; | double sigma2[] = { 0.20, 0.45, 0.18 }; | ||||
int eq_num[] = {9, 10, 9}; | int eq_num[] = { 9, 10, 9 }; | ||||
double M_P2 = 5.27 * 5.27; // GeV^2 for B_d^0 - meson | double M_P2 = 5.27 * 5.27; // GeV^2 for B_d^0 - meson | ||||
double M_V2 = 5.32 * 5.32; // GeV^2 for B_d^* - meson | double M_V2 = 5.32 * 5.32; // GeV^2 for B_d^* - meson | ||||
fp = equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | fp = equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | ||||
f0 = equation9_10( ff0[1], M_V2, t, sigma1[1], sigma2[1], eq_num[1] ); | f0 = equation9_10( ff0[1], M_V2, t, sigma1[1], sigma2[1], eq_num[1] ); | ||||
ft = equation9_10( ff0[2], M_P2, t, sigma1[2], sigma2[2], eq_num[2] ); | ft = equation9_10( ff0[2], M_P2, t, sigma1[2], sigma2[2], eq_num[2] ); | ||||
models_counter = models_counter + 1; | models_counter = models_counter + 1; | ||||
// EvtGenReport(EVTGEN_NOTICE,"EvtGen") <<"\n The function EvtbTosllMSFF::getVectorFF(...) passed." | // EvtGenReport(EVTGEN_NOTICE,"EvtGen") <<"\n The function EvtbTosllMSFF::getVectorFF(...) passed." | ||||
// << "\n Bd -> eta transition form factors" | // << "\n Bd -> eta transition form factors" | ||||
// << std::endl; | // << std::endl; | ||||
} | } | ||||
// B_s -> \eta transition form factors | // B_s -> \eta transition form factors | ||||
if ( ( parent == EvtPDL::getId( std::string( "B_s0" ) ) && | if ( ( parent == EvtPDL::getId( std::string( "B_s0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "eta" ) ) ) || | daught == EvtPDL::getId( std::string( "eta" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "anti-B_s0" ) ) && | ( parent == EvtPDL::getId( std::string( "anti-B_s0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "eta" ) ) ) ) { | daught == EvtPDL::getId( std::string( "eta" ) ) ) ) { | ||||
double ff0[] = {0.36, 0.36, 0.36}; | double ff0[] = { 0.36, 0.36, 0.36 }; | ||||
double sigma1[] = {0.60, 0.80, 0.58}; | double sigma1[] = { 0.60, 0.80, 0.58 }; | ||||
double sigma2[] = {0.20, 0.40, 0.18}; | double sigma2[] = { 0.20, 0.40, 0.18 }; | ||||
int eq_num[] = {9, 10, 9}; | int eq_num[] = { 9, 10, 9 }; | ||||
double M_P2 = 5.37 * 5.37; // GeV^2 for B_s^0 - meson | double M_P2 = 5.37 * 5.37; // GeV^2 for B_s^0 - meson | ||||
double M_V2 = 5.42 * 5.42; // GeV^2 for B_s^* - meson | double M_V2 = 5.42 * 5.42; // GeV^2 for B_s^* - meson | ||||
fp = equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | fp = equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | ||||
f0 = equation9_10( ff0[1], M_V2, t, sigma1[1], sigma2[1], eq_num[1] ); | f0 = equation9_10( ff0[1], M_V2, t, sigma1[1], sigma2[1], eq_num[1] ); | ||||
ft = equation9_10( ff0[2], M_P2, t, sigma1[2], sigma2[2], eq_num[2] ); | ft = equation9_10( ff0[2], M_P2, t, sigma1[2], sigma2[2], eq_num[2] ); | ||||
models_counter = models_counter + 1; | models_counter = models_counter + 1; | ||||
// EvtGenReport(EVTGEN_NOTICE,"EvtGen") <<"\n The function EvtbTosllMSFF::getVectorFF(...) passed." | // EvtGenReport(EVTGEN_NOTICE,"EvtGen") <<"\n The function EvtbTosllMSFF::getVectorFF(...) passed." | ||||
// << "\n Bs -> eta transition form factors" | // << "\n Bs -> eta transition form factors" | ||||
// << std::endl; | // << std::endl; | ||||
} | } | ||||
// B_s -> \eta' transition form factors | // B_s -> \eta' transition form factors | ||||
if ( ( parent == EvtPDL::getId( std::string( "B_s0" ) ) && | if ( ( parent == EvtPDL::getId( std::string( "B_s0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "eta'" ) ) ) || | daught == EvtPDL::getId( std::string( "eta'" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "anti-B_s0" ) ) && | ( parent == EvtPDL::getId( std::string( "anti-B_s0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "eta'" ) ) ) ) { | daught == EvtPDL::getId( std::string( "eta'" ) ) ) ) { | ||||
double ff0[] = {0.36, 0.36, 0.39}; | double ff0[] = { 0.36, 0.36, 0.39 }; | ||||
double sigma1[] = {0.60, 0.80, 0.58}; | double sigma1[] = { 0.60, 0.80, 0.58 }; | ||||
double sigma2[] = {0.20, 0.45, 0.18}; | double sigma2[] = { 0.20, 0.45, 0.18 }; | ||||
int eq_num[] = {9, 10, 9}; | int eq_num[] = { 9, 10, 9 }; | ||||
double M_P2 = 5.37 * 5.37; // GeV^2 for B_s^0 - meson | double M_P2 = 5.37 * 5.37; // GeV^2 for B_s^0 - meson | ||||
double M_V2 = 5.42 * 5.42; // GeV^2 for B_s^* - meson | double M_V2 = 5.42 * 5.42; // GeV^2 for B_s^* - meson | ||||
fp = equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | fp = equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | ||||
f0 = equation9_10( ff0[1], M_V2, t, sigma1[1], sigma2[1], eq_num[1] ); | f0 = equation9_10( ff0[1], M_V2, t, sigma1[1], sigma2[1], eq_num[1] ); | ||||
ft = equation9_10( ff0[2], M_P2, t, sigma1[2], sigma2[2], eq_num[2] ); | ft = equation9_10( ff0[2], M_P2, t, sigma1[2], sigma2[2], eq_num[2] ); | ||||
models_counter = models_counter + 1; | models_counter = models_counter + 1; | ||||
// EvtGenReport(EVTGEN_NOTICE,"EvtGen") <<"\n The function EvtbTosllMSFF::getVectorFF(...) passed." | // EvtGenReport(EVTGEN_NOTICE,"EvtGen") <<"\n The function EvtbTosllMSFF::getVectorFF(...) passed." | ||||
// << "\n Bs -> eta transition form factors" | // << "\n Bs -> eta transition form factors" | ||||
// << std::endl; | // << std::endl; | ||||
} | } | ||||
// B_s -> f_0(980) transition form factors | // B_s -> f_0(980) transition form factors | ||||
if ( ( parent == EvtPDL::getId( std::string( "B_s0" ) ) && | if ( ( parent == EvtPDL::getId( std::string( "B_s0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "f_0" ) ) ) || | daught == EvtPDL::getId( std::string( "f_0" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "anti-B_s0" ) ) && | ( parent == EvtPDL::getId( std::string( "anti-B_s0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "f_0" ) ) ) ) { | daught == EvtPDL::getId( std::string( "f_0" ) ) ) ) { | ||||
double ff0[] = {0.238, 0.238, 0.308}; | double ff0[] = { 0.238, 0.238, 0.308 }; | ||||
double sigma1[] = {1.50, 0.53, 1.46}; | double sigma1[] = { 1.50, 0.53, 1.46 }; | ||||
double sigma2[] = {0.58, -0.36, 0.58}; | double sigma2[] = { 0.58, -0.36, 0.58 }; | ||||
int eq_num[] = {10, 10, 10}; | int eq_num[] = { 10, 10, 10 }; | ||||
double M_P2 = 5.366 * 5.366; // GeV^2 for B_s^0 - meson | double M_P2 = 5.366 * 5.366; // GeV^2 for B_s^0 - meson | ||||
fp = 0.0 - | fp = 0.0 - | ||||
equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | ||||
f0 = 0.0 - | f0 = 0.0 - | ||||
equation9_10( ff0[1], M_P2, t, sigma1[1], sigma2[1], eq_num[1] ); | equation9_10( ff0[1], M_P2, t, sigma1[1], sigma2[1], eq_num[1] ); | ||||
ft = equation9_10( ff0[2], M_P2, t, sigma1[2], sigma2[2], eq_num[2] ); | ft = equation9_10( ff0[2], M_P2, t, sigma1[2], sigma2[2], eq_num[2] ); | ||||
Show All 31 Lines | void EvtbTosllMSFF::getVectorFF( EvtId parent, EvtId daught, double t, | ||||
if ( ( parent == EvtPDL::getId( std::string( "B+" ) ) && | if ( ( parent == EvtPDL::getId( std::string( "B+" ) ) && | ||||
daught == EvtPDL::getId( std::string( "K*+" ) ) ) || | daught == EvtPDL::getId( std::string( "K*+" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "B-" ) ) && | ( parent == EvtPDL::getId( std::string( "B-" ) ) && | ||||
daught == EvtPDL::getId( std::string( "K*-" ) ) ) || | daught == EvtPDL::getId( std::string( "K*-" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "B0" ) ) && | ( parent == EvtPDL::getId( std::string( "B0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "K*0" ) ) ) || | daught == EvtPDL::getId( std::string( "K*0" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "anti-B0" ) ) && | ( parent == EvtPDL::getId( std::string( "anti-B0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "anti-K*0" ) ) ) ) { | daught == EvtPDL::getId( std::string( "anti-K*0" ) ) ) ) { | ||||
double ff0[] = {0.44, 0.45, 0.36, 0.32, 0.39, 0.39, 0.27}; | double ff0[] = { 0.44, 0.45, 0.36, 0.32, 0.39, 0.39, 0.27 }; | ||||
double sigma1[] = {0.45, 0.46, 0.64, 1.23, 0.45, 0.72, 1.31}; | double sigma1[] = { 0.45, 0.46, 0.64, 1.23, 0.45, 0.72, 1.31 }; | ||||
double sigma2[] = {0.00, 0.00, 0.36, 0.38, 0.00, 0.62, 0.41}; | double sigma2[] = { 0.00, 0.00, 0.36, 0.38, 0.00, 0.62, 0.41 }; | ||||
int eq_num[] = {9, 9, 10, 10, 9, 10, 10}; | int eq_num[] = { 9, 9, 10, 10, 9, 10, 10 }; | ||||
double M_P2 = 5.37 * 5.37; // GeV^2 for B^0_s - meson | double M_P2 = 5.37 * 5.37; // GeV^2 for B^0_s - meson | ||||
double M_V2 = 5.42 * 5.42; // GeV^2 for B^*_s - meson | double M_V2 = 5.42 * 5.42; // GeV^2 for B^*_s - meson | ||||
v = equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | v = equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | ||||
a0 = equation9_10( ff0[1], M_P2, t, sigma1[1], sigma2[1], eq_num[1] ); | a0 = equation9_10( ff0[1], M_P2, t, sigma1[1], sigma2[1], eq_num[1] ); | ||||
a1 = equation9_10( ff0[2], M_V2, t, sigma1[2], sigma2[2], eq_num[2] ); | a1 = equation9_10( ff0[2], M_V2, t, sigma1[2], sigma2[2], eq_num[2] ); | ||||
a2 = equation9_10( ff0[3], M_V2, t, sigma1[3], sigma2[3], eq_num[3] ); | a2 = equation9_10( ff0[3], M_V2, t, sigma1[3], sigma2[3], eq_num[3] ); | ||||
Show All 12 Lines | void EvtbTosllMSFF::getVectorFF( EvtId parent, EvtId daught, double t, | ||||
if ( ( parent == EvtPDL::getId( std::string( "B+" ) ) && | if ( ( parent == EvtPDL::getId( std::string( "B+" ) ) && | ||||
daught == EvtPDL::getId( std::string( "rho+" ) ) ) || | daught == EvtPDL::getId( std::string( "rho+" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "B-" ) ) && | ( parent == EvtPDL::getId( std::string( "B-" ) ) && | ||||
daught == EvtPDL::getId( std::string( "rho-" ) ) ) || | daught == EvtPDL::getId( std::string( "rho-" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "B0" ) ) && | ( parent == EvtPDL::getId( std::string( "B0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "rho0" ) ) ) || | daught == EvtPDL::getId( std::string( "rho0" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "anti-B0" ) ) && | ( parent == EvtPDL::getId( std::string( "anti-B0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "rho0" ) ) ) ) { | daught == EvtPDL::getId( std::string( "rho0" ) ) ) ) { | ||||
double ff0[] = {0.31, 0.30, 0.26, 0.24, 0.27, 0.27, 0.19}; | double ff0[] = { 0.31, 0.30, 0.26, 0.24, 0.27, 0.27, 0.19 }; | ||||
double sigma1[] = {0.59, 0.54, 0.73, 1.40, 0.60, 0.74, 1.42}; | double sigma1[] = { 0.59, 0.54, 0.73, 1.40, 0.60, 0.74, 1.42 }; | ||||
double sigma2[] = {0.00, 0.00, 0.10, 0.50, 0.00, 0.19, 0.51}; | double sigma2[] = { 0.00, 0.00, 0.10, 0.50, 0.00, 0.19, 0.51 }; | ||||
int eq_num[] = {9, 9, 10, 10, 9, 10, 10}; | int eq_num[] = { 9, 9, 10, 10, 9, 10, 10 }; | ||||
double M_P2 = 5.27 * 5.27; // GeV^2 for B - meson | double M_P2 = 5.27 * 5.27; // GeV^2 for B - meson | ||||
double M_V2 = 5.32 * 5.32; // GeV^2 for B^* - meson | double M_V2 = 5.32 * 5.32; // GeV^2 for B^* - meson | ||||
v = equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | v = equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | ||||
a0 = equation9_10( ff0[1], M_P2, t, sigma1[1], sigma2[1], eq_num[1] ); | a0 = equation9_10( ff0[1], M_P2, t, sigma1[1], sigma2[1], eq_num[1] ); | ||||
a1 = equation9_10( ff0[2], M_V2, t, sigma1[2], sigma2[2], eq_num[2] ); | a1 = equation9_10( ff0[2], M_V2, t, sigma1[2], sigma2[2], eq_num[2] ); | ||||
a2 = equation9_10( ff0[3], M_V2, t, sigma1[3], sigma2[3], eq_num[3] ); | a2 = equation9_10( ff0[3], M_V2, t, sigma1[3], sigma2[3], eq_num[3] ); | ||||
t1 = equation9_10( ff0[4], M_P2, t, sigma1[4], sigma2[4], eq_num[4] ); | t1 = equation9_10( ff0[4], M_P2, t, sigma1[4], sigma2[4], eq_num[4] ); | ||||
t2 = equation9_10( ff0[5], M_V2, t, sigma1[5], sigma2[5], eq_num[5] ); | t2 = equation9_10( ff0[5], M_V2, t, sigma1[5], sigma2[5], eq_num[5] ); | ||||
t3 = equation9_10( ff0[6], M_V2, t, sigma1[6], sigma2[6], eq_num[6] ); | t3 = equation9_10( ff0[6], M_V2, t, sigma1[6], sigma2[6], eq_num[6] ); | ||||
models_counter = models_counter + 1; | models_counter = models_counter + 1; | ||||
// EvtGenReport(EVTGEN_NOTICE,"EvtGen") <<"\n The function EvtbTosllMSFF::getVectorFF(...) passed." | // EvtGenReport(EVTGEN_NOTICE,"EvtGen") <<"\n The function EvtbTosllMSFF::getVectorFF(...) passed." | ||||
// << "\n barB -> bar rho transition form factors" | // << "\n barB -> bar rho transition form factors" | ||||
// << std::endl; | // << std::endl; | ||||
} | } | ||||
// \bar B -> \omega transition form factors (exactly as for \bar B -> \rho^0 ff!) | // \bar B -> \omega transition form factors (exactly as for \bar B -> \rho^0 ff!) | ||||
if ( ( parent == EvtPDL::getId( std::string( "B0" ) ) && | if ( ( parent == EvtPDL::getId( std::string( "B0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "omega" ) ) ) || | daught == EvtPDL::getId( std::string( "omega" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "anti-B0" ) ) && | ( parent == EvtPDL::getId( std::string( "anti-B0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "omega" ) ) ) ) { | daught == EvtPDL::getId( std::string( "omega" ) ) ) ) { | ||||
double ff0[] = {0.31, 0.30, 0.26, 0.24, 0.27, 0.27, 0.19}; | double ff0[] = { 0.31, 0.30, 0.26, 0.24, 0.27, 0.27, 0.19 }; | ||||
double sigma1[] = {0.59, 0.54, 0.73, 1.40, 0.60, 0.74, 1.42}; | double sigma1[] = { 0.59, 0.54, 0.73, 1.40, 0.60, 0.74, 1.42 }; | ||||
double sigma2[] = {0.00, 0.00, 0.10, 0.50, 0.00, 0.19, 0.51}; | double sigma2[] = { 0.00, 0.00, 0.10, 0.50, 0.00, 0.19, 0.51 }; | ||||
int eq_num[] = {9, 9, 10, 10, 9, 10, 10}; | int eq_num[] = { 9, 9, 10, 10, 9, 10, 10 }; | ||||
double M_P2 = 5.27 * 5.27; // GeV^2 for B - meson | double M_P2 = 5.27 * 5.27; // GeV^2 for B - meson | ||||
double M_V2 = 5.32 * 5.32; // GeV^2 for B^* - meson | double M_V2 = 5.32 * 5.32; // GeV^2 for B^* - meson | ||||
v = equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | v = equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | ||||
a0 = equation9_10( ff0[1], M_P2, t, sigma1[1], sigma2[1], eq_num[1] ); | a0 = equation9_10( ff0[1], M_P2, t, sigma1[1], sigma2[1], eq_num[1] ); | ||||
a1 = equation9_10( ff0[2], M_V2, t, sigma1[2], sigma2[2], eq_num[2] ); | a1 = equation9_10( ff0[2], M_V2, t, sigma1[2], sigma2[2], eq_num[2] ); | ||||
a2 = equation9_10( ff0[3], M_V2, t, sigma1[3], sigma2[3], eq_num[3] ); | a2 = equation9_10( ff0[3], M_V2, t, sigma1[3], sigma2[3], eq_num[3] ); | ||||
t1 = equation9_10( ff0[4], M_P2, t, sigma1[4], sigma2[4], eq_num[4] ); | t1 = equation9_10( ff0[4], M_P2, t, sigma1[4], sigma2[4], eq_num[4] ); | ||||
t2 = equation9_10( ff0[5], M_V2, t, sigma1[5], sigma2[5], eq_num[5] ); | t2 = equation9_10( ff0[5], M_V2, t, sigma1[5], sigma2[5], eq_num[5] ); | ||||
t3 = equation9_10( ff0[6], M_V2, t, sigma1[6], sigma2[6], eq_num[6] ); | t3 = equation9_10( ff0[6], M_V2, t, sigma1[6], sigma2[6], eq_num[6] ); | ||||
models_counter = models_counter + 1; | models_counter = models_counter + 1; | ||||
// EvtGenReport(EVTGEN_NOTICE,"EvtGen") <<"\n The function EvtbTosllMSFF::getVectorFF(...) passed." | // EvtGenReport(EVTGEN_NOTICE,"EvtGen") <<"\n The function EvtbTosllMSFF::getVectorFF(...) passed." | ||||
// << "\n barB -> omega transition form factors" | // << "\n barB -> omega transition form factors" | ||||
// << std::endl; | // << std::endl; | ||||
} | } | ||||
// \bar Bs -> phi transition form factors | // \bar Bs -> phi transition form factors | ||||
if ( ( parent == EvtPDL::getId( std::string( "B_s0" ) ) && | if ( ( parent == EvtPDL::getId( std::string( "B_s0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "phi" ) ) ) || | daught == EvtPDL::getId( std::string( "phi" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "anti-B_s0" ) ) && | ( parent == EvtPDL::getId( std::string( "anti-B_s0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "phi" ) ) ) ) { | daught == EvtPDL::getId( std::string( "phi" ) ) ) ) { | ||||
double ff0[] = {0.44, 0.42, 0.34, 0.31, 0.38, 0.38, 0.26}; | double ff0[] = { 0.44, 0.42, 0.34, 0.31, 0.38, 0.38, 0.26 }; | ||||
double sigma1[] = {0.62, 0.55, 0.73, 1.30, 0.62, 0.83, 1.41}; | double sigma1[] = { 0.62, 0.55, 0.73, 1.30, 0.62, 0.83, 1.41 }; | ||||
double sigma2[] = {0.20, 0.12, 0.42, 0.52, 0.20, 0.71, 0.57}; | double sigma2[] = { 0.20, 0.12, 0.42, 0.52, 0.20, 0.71, 0.57 }; | ||||
int eq_num[] = {9, 9, 10, 10, 9, 10, 10}; | int eq_num[] = { 9, 9, 10, 10, 9, 10, 10 }; | ||||
double M_P2 = 5.37 * 5.37; // GeV^2 for B^0_s - meson | double M_P2 = 5.37 * 5.37; // GeV^2 for B^0_s - meson | ||||
double M_V2 = 5.42 * 5.42; // GeV^2 for B^*_s - meson | double M_V2 = 5.42 * 5.42; // GeV^2 for B^*_s - meson | ||||
v = equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | v = equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | ||||
a0 = equation9_10( ff0[1], M_P2, t, sigma1[1], sigma2[1], eq_num[1] ); | a0 = equation9_10( ff0[1], M_P2, t, sigma1[1], sigma2[1], eq_num[1] ); | ||||
a1 = equation9_10( ff0[2], M_V2, t, sigma1[2], sigma2[2], eq_num[2] ); | a1 = equation9_10( ff0[2], M_V2, t, sigma1[2], sigma2[2], eq_num[2] ); | ||||
a2 = equation9_10( ff0[3], M_V2, t, sigma1[3], sigma2[3], eq_num[3] ); | a2 = equation9_10( ff0[3], M_V2, t, sigma1[3], sigma2[3], eq_num[3] ); | ||||
t1 = equation9_10( ff0[4], M_P2, t, sigma1[4], sigma2[4], eq_num[4] ); | t1 = equation9_10( ff0[4], M_P2, t, sigma1[4], sigma2[4], eq_num[4] ); | ||||
t2 = equation9_10( ff0[5], M_V2, t, sigma1[5], sigma2[5], eq_num[5] ); | t2 = equation9_10( ff0[5], M_V2, t, sigma1[5], sigma2[5], eq_num[5] ); | ||||
t3 = equation9_10( ff0[6], M_V2, t, sigma1[6], sigma2[6], eq_num[6] ); | t3 = equation9_10( ff0[6], M_V2, t, sigma1[6], sigma2[6], eq_num[6] ); | ||||
models_counter = models_counter + 1; | models_counter = models_counter + 1; | ||||
// EvtGenReport(EVTGEN_NOTICE,"EvtGen") <<"\n The function EvtbTosllMSFF::getVectorFF(...) passed." | // EvtGenReport(EVTGEN_NOTICE,"EvtGen") <<"\n The function EvtbTosllMSFF::getVectorFF(...) passed." | ||||
// << "\n barBs -> phi transition form factors" | // << "\n barBs -> phi transition form factors" | ||||
// << std::endl; | // << std::endl; | ||||
} | } | ||||
// \bar Bs -> K* (without \bar !) transition form factors | // \bar Bs -> K* (without \bar !) transition form factors | ||||
if ( ( parent == EvtPDL::getId( std::string( "B_s0" ) ) && | if ( ( parent == EvtPDL::getId( std::string( "B_s0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "anti-K*0" ) ) ) || | daught == EvtPDL::getId( std::string( "anti-K*0" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "anti-B_s0" ) ) && | ( parent == EvtPDL::getId( std::string( "anti-B_s0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "K*0" ) ) ) ) { | daught == EvtPDL::getId( std::string( "K*0" ) ) ) ) { | ||||
double ff0[] = {0.38, 0.37, 0.29, 0.26, 0.32, 0.32, 0.23}; | double ff0[] = { 0.38, 0.37, 0.29, 0.26, 0.32, 0.32, 0.23 }; | ||||
double sigma1[] = {0.66, 0.60, 0.86, 1.32, 0.66, 0.98, 1.42}; | double sigma1[] = { 0.66, 0.60, 0.86, 1.32, 0.66, 0.98, 1.42 }; | ||||
double sigma2[] = {0.30, 0.16, 0.60, 0.54, 0.31, 0.90, 0.62}; | double sigma2[] = { 0.30, 0.16, 0.60, 0.54, 0.31, 0.90, 0.62 }; | ||||
int eq_num[] = {9, 9, 10, 10, 9, 10, 10}; | int eq_num[] = { 9, 9, 10, 10, 9, 10, 10 }; | ||||
double M_P2 = 5.27 * 5.27; // GeV^2 for B - meson | double M_P2 = 5.27 * 5.27; // GeV^2 for B - meson | ||||
double M_V2 = 5.32 * 5.32; // GeV^2 for B^* - meson | double M_V2 = 5.32 * 5.32; // GeV^2 for B^* - meson | ||||
v = equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | v = equation9_10( ff0[0], M_P2, t, sigma1[0], sigma2[0], eq_num[0] ); | ||||
a0 = equation9_10( ff0[1], M_P2, t, sigma1[1], sigma2[1], eq_num[1] ); | a0 = equation9_10( ff0[1], M_P2, t, sigma1[1], sigma2[1], eq_num[1] ); | ||||
a1 = equation9_10( ff0[2], M_V2, t, sigma1[2], sigma2[2], eq_num[2] ); | a1 = equation9_10( ff0[2], M_V2, t, sigma1[2], sigma2[2], eq_num[2] ); | ||||
a2 = equation9_10( ff0[3], M_V2, t, sigma1[3], sigma2[3], eq_num[3] ); | a2 = equation9_10( ff0[3], M_V2, t, sigma1[3], sigma2[3], eq_num[3] ); | ||||
Show All 13 Lines | void EvtbTosllMSFF::getVectorFF( EvtId parent, EvtId daught, double t, | ||||
if ( ( parent == EvtPDL::getId( std::string( "B+" ) ) && | if ( ( parent == EvtPDL::getId( std::string( "B+" ) ) && | ||||
daught == EvtPDL::getId( std::string( "K_1+" ) ) ) || | daught == EvtPDL::getId( std::string( "K_1+" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "B-" ) ) && | ( parent == EvtPDL::getId( std::string( "B-" ) ) && | ||||
daught == EvtPDL::getId( std::string( "K_1-" ) ) ) || | daught == EvtPDL::getId( std::string( "K_1-" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "B0" ) ) && | ( parent == EvtPDL::getId( std::string( "B0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "K_10" ) ) ) || | daught == EvtPDL::getId( std::string( "K_10" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "anti-B0" ) ) && | ( parent == EvtPDL::getId( std::string( "anti-B0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "anti-K_10" ) ) ) ) { | daught == EvtPDL::getId( std::string( "anti-K_10" ) ) ) ) { | ||||
double ff0A[] = {0.450, 0.340, 0.41, 0.22, 0.31, 0.310, 0.28}; | double ff0A[] = { 0.450, 0.340, 0.41, 0.22, 0.31, 0.310, 0.28 }; | ||||
double sigma1A[] = {1.600, 0.635, 1.51, 2.40, 2.01, 0.629, 1.36}; | double sigma1A[] = { 1.600, 0.635, 1.51, 2.40, 2.01, 0.629, 1.36 }; | ||||
double sigma2A[] = {0.974, 0.211, 1.18, 1.78, 1.50, 0.387, 0.72}; | double sigma2A[] = { 0.974, 0.211, 1.18, 1.78, 1.50, 0.387, 0.72 }; | ||||
double ff0B[] = {-0.37, -0.29, -0.17, -0.45, -0.25, -0.250, -0.11}; | double ff0B[] = { -0.37, -0.29, -0.17, -0.45, -0.25, -0.250, -0.11 }; | ||||
double sigma1B[] = {1.72, 0.729, 0.919, 1.34, 1.59, 0.378, -1.61}; | double sigma1B[] = { 1.72, 0.729, 0.919, 1.34, 1.59, 0.378, -1.61 }; | ||||
double sigma2B[] = {0.912, 0.074, 0.855, 0.69, 0.79, -0.755, 10.2}; | double sigma2B[] = { 0.912, 0.074, 0.855, 0.69, 0.79, -0.755, 10.2 }; | ||||
int eq_num[] = {10, 10, 10, 10, 10, 10, 10}; | int eq_num[] = { 10, 10, 10, 10, 10, 10, 10 }; | ||||
double MM2 = 5.279 * 5.279; // GeV^2 | double MM2 = 5.279 * 5.279; // GeV^2 | ||||
double MB = 5.279; // GeV | double MB = 5.279; // GeV | ||||
double MK1 = 1.272; // GeV | double MK1 = 1.272; // GeV | ||||
double MK1A = 1.31; // GeV | double MK1A = 1.31; // GeV | ||||
double MK1B = 1.34; // GeV | double MK1B = 1.34; // GeV | ||||
double sinK = sin( thetaK ); // sin(-34^o) | double sinK = sin( thetaK ); // sin(-34^o) | ||||
▲ Show 20 Lines • Show All 68 Lines • ▼ Show 20 Lines | void EvtbTosllMSFF::getVectorFF( EvtId parent, EvtId daught, double t, | ||||
if ( ( parent == EvtPDL::getId( std::string( "B+" ) ) && | if ( ( parent == EvtPDL::getId( std::string( "B+" ) ) && | ||||
daught == EvtPDL::getId( std::string( "K'_1+" ) ) ) || | daught == EvtPDL::getId( std::string( "K'_1+" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "B-" ) ) && | ( parent == EvtPDL::getId( std::string( "B-" ) ) && | ||||
daught == EvtPDL::getId( std::string( "K'_1-" ) ) ) || | daught == EvtPDL::getId( std::string( "K'_1-" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "B0" ) ) && | ( parent == EvtPDL::getId( std::string( "B0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "K'_10" ) ) ) || | daught == EvtPDL::getId( std::string( "K'_10" ) ) ) || | ||||
( parent == EvtPDL::getId( std::string( "anti-B0" ) ) && | ( parent == EvtPDL::getId( std::string( "anti-B0" ) ) && | ||||
daught == EvtPDL::getId( std::string( "anti-K'_10" ) ) ) ) { | daught == EvtPDL::getId( std::string( "anti-K'_10" ) ) ) ) { | ||||
double ff0A[] = {0.450, 0.340, 0.41, 0.22, 0.31, 0.310, 0.28}; | double ff0A[] = { 0.450, 0.340, 0.41, 0.22, 0.31, 0.310, 0.28 }; | ||||
double sigma1A[] = {1.600, 0.635, 1.51, 2.40, 2.01, 0.629, 1.36}; | double sigma1A[] = { 1.600, 0.635, 1.51, 2.40, 2.01, 0.629, 1.36 }; | ||||
double sigma2A[] = {0.974, 0.211, 1.18, 1.78, 1.50, 0.387, 0.72}; | double sigma2A[] = { 0.974, 0.211, 1.18, 1.78, 1.50, 0.387, 0.72 }; | ||||
double ff0B[] = {-0.37, -0.29, -0.17, -0.45, -0.25, -0.250, -0.11}; | double ff0B[] = { -0.37, -0.29, -0.17, -0.45, -0.25, -0.250, -0.11 }; | ||||
double sigma1B[] = {1.72, 0.729, 0.919, 1.34, 1.59, 0.378, -1.61}; | double sigma1B[] = { 1.72, 0.729, 0.919, 1.34, 1.59, 0.378, -1.61 }; | ||||
double sigma2B[] = {0.912, 0.074, 0.855, 0.69, 0.79, -0.755, 10.2}; | double sigma2B[] = { 0.912, 0.074, 0.855, 0.69, 0.79, -0.755, 10.2 }; | ||||
int eq_num[] = {10, 10, 10, 10, 10, 10, 10}; | int eq_num[] = { 10, 10, 10, 10, 10, 10, 10 }; | ||||
double MM2 = 5.279 * 5.279; // GeV^2 | double MM2 = 5.279 * 5.279; // GeV^2 | ||||
double MB = 5.279; // GeV | double MB = 5.279; // GeV | ||||
double MK1 = 1.403; // GeV | double MK1 = 1.403; // GeV | ||||
double MK1A = 1.31; // GeV | double MK1A = 1.31; // GeV | ||||
double MK1B = 1.34; // GeV | double MK1B = 1.34; // GeV | ||||
double sinK = sin( thetaK ); // sin(-34^o) | double sinK = sin( thetaK ); // sin(-34^o) | ||||
▲ Show 20 Lines • Show All 124 Lines • Show Last 20 Lines |