Commit 4873b1b0 authored by Frank Siegert's avatar Frank Siegert

Merge Ds+ decay improvements by Katharina Danziger from tmp-bachelor2 branch:

c30347
c30330
c30303
parent cb82fcb2
......@@ -69,8 +69,7 @@ void VA_P_P::SetModelParameters( struct GeneralModel model )
Vxx=Tools::Vcd;
}
else if (kf0==kf_D_s_plus) {
if (kf1==kf_eta_prime_958 || kf1==kf_eta)
if (kf1==kf_eta_prime_958 || kf1==kf_eta || kf1==kf_f_0_980)
Vxx=Tools::Vcs;
else if (kf1==kf_K || kf1==kf_K_S || kf1==kf_K_L)
Vxx=Tools::Vcd;
......
......@@ -39,6 +39,20 @@ BallZwicky::BallZwicky(GeneralModel model, double* masses,
m_fplus_r = 0.963;
m_fplus_m2 = 2.45*2.45;
}
else if (kf0==kf_D_s_plus) {
if (kf1==kf_eta) {
m_fplus_0 = 0.432;
m_fplus_alpha = 0.252;
m_fplus_r = 0.284;
m_fplus_m2 = 1.9685*1.9685;
}
else if (kf1==kf_eta_prime_958) {
m_fplus_0 = 0.520;
m_fplus_alpha = 0.252;
m_fplus_r = 0.284;
m_fplus_m2 = 1.9685*1.9685;
}
}
m_fplus_0 = model("fplus_0",m_fplus_0);
m_fplus_r = model("fplus_r",m_fplus_r);
......
......@@ -119,14 +119,139 @@ PoleFit::PoleFit(GeneralModel model,double* masses, const Flavour_Vector& flavs,
}
else if (kf0==kf_D_s_plus) {
if (kf1==kf_f_0_980) {
// hep-ph 0701108
m_fplus_F0 = 1.7*0.27;
m_fplus_m0 = m_m0;
m_fplus_a0 = -0.87;
m_fplus_b0 = -0.17;
m_fplus_c0 = -0.37;
m_fplus_d0 = 1.46;
switch( int(model("FORM_FACTOR_TUNE", 3)+0.5)) {
case 2:
// hep-ph 0304193 QCD sum rules
m_fplus_F0 = 0.44;
m_fplus_m0 = 1.6;
m_fplus_a0 = 1.0;
break;
case 3:
// hep-ph 0701108 3-point QCD sum rules (2007)
m_fplus_F0 = 1.7*0.27;
m_fplus_m0 = m_m0;
m_fplus_a0 = -0.87;
m_fplus_b0 = -0.17;
m_fplus_c0 = -0.37;
m_fplus_d0 = 1.46;
break;
case 5:
// hep-ph 1002.2466 light-front quark model
m_fplus_F0 = 0.46;
m_fplus_m0 = m_m0;
m_fplus_a0 = 0.29;
m_fplus_b0 = 0.07;
break;
case 6:
// hep-ph 0907.3201 experimental fit
m_fplus_F0 = 1.0;
m_fplus_m0 = 1.7;
m_fplus_a0 = 1.0;
break;
case 7:
// doi:10.1209/0295-5075/90/61001 3-point QCD sum rules (2010)
m_fplus_F0 = 0.48;
m_fplus_m0 = m_m0;
m_fplus_a0 = -0.81;
m_fplus_b0 = -0.18;
m_fplus_c0 = -0.19;
m_fplus_d0 = 0.86;
break;
default:
msg_Error()<<METHOD<<": You chose a form factor tune which does not "
<<"exist. Aborting."<<std::endl;
abort();
}
}
else if (kf1==kf_eta) {
switch( int(model("FORM_FACTOR_TUNE", 3)+0.5)) {
case 1:
// hep-ph 1209.3386 lattice QCD
m_fplus_F0 = 0.75;
m_fplus_m0 = 1.9;
m_fplus_a0 = -(1.0+0.21);
m_fplus_b0 = 0.21;
break;
case 4:
// hep-ph 1011.6046 light-cone sum rules (2011)
m_fplus_F0 = 0.45;
m_fplus_m0 = 1.9685;
m_fplus_a0 = -1.96;
m_fplus_b0 = 1.12;
break;
case 5:
// hep-ph 1011.6046 light-front quark model 1
m_fplus_F0 = 0.50;
m_fplus_m0 = 1.9685;
m_fplus_a0 = -1.17;
m_fplus_b0 = 0.34;
break;
case 8:
// hep-ph 1011.6046 light-front quark model 2
m_fplus_F0 = 0.48;
m_fplus_m0 = 1.9685;
m_fplus_a0 = -1.11;
m_fplus_b0 = 0.25;
m_fmin_F0 = -0.44;
m_fmin_m0 = 1.9685;
m_fmin_a0 = -2.05;
m_fmin_b0 = 1.08;
break;
default:
msg_Error()<<METHOD<<": You chose a form factor tune which does not "
<<"exist. Aborting."<<std::endl;
abort();
}
}
else if (kf1==kf_eta_prime_958) {
switch( int(model("FORM_FACTOR_TUNE", 3)+0.5)) {
case 1:
// hep-ph 1209.3386 lattice QCD
m_fplus_F0 = 0.86;
m_fplus_m0 = 1.9;
m_fplus_a0 = -(1.0+0.21);
m_fplus_b0 = 0.21;
break;
case 4:
// hep-ph 1011.6046 light-cone sum rules (2011)
m_fplus_F0 = 0.55;
m_fplus_m0 = 1.9685;
m_fplus_a0 = -1.96;
m_fplus_b0 = 1.12;
break;
case 5:
// hep-ph 1011.6046 light-front quark model 1
m_fplus_F0 = 0.62;
m_fplus_m0 = 1.9685;
m_fplus_a0 = -1.17;
m_fplus_b0 = 0.34;
break;
case 8:
// hep-ph 1011.6046 light-front quark model 2
m_fplus_F0 = 0.60;
m_fplus_m0 = 1.9685;
m_fplus_a0 = -1.11;
m_fplus_b0 = 0.25;
break;
default:
msg_Error()<<METHOD<<": You chose a form factor tune which does not "
<<"exist. Aborting."<<std::endl;
abort();
}
}
else if (kf1==kf_K) {
// PhysRevD.78.054002 light-front quark model
m_fplus_F0 = 0.67;
m_fplus_m0 = 1.9685;
m_fplus_a0 = -1.17;
m_fplus_b0 = 0.29;
m_f0_F0 = 0.67;
m_f0_m0 = 1.9685;
m_f0_a0 = -0.50;
m_f0_b0 = -0.005;
}
}
else if (kf0==kf_B_c) {
// hep-ph/0007169
......
......@@ -48,6 +48,31 @@ Polynomial::Polynomial(GeneralModel model, double* masses,
}
}
else if (kf0==kf_D_s_plus) {
if (kf1==kf_eta) {
// hep-ph/0107137
// QCD sum rules
m_fplus_0 = 0.50;
m_fplus_lambda = 1.0108;
m_fplus_m2 = 1.9*1.9;
m_f0_0 = 0.0;
m_f0_lambda = 0.0;
m_f0_m2 = 1.9*1.9;
}
else if (kf1==kf_eta_prime_958) {
// hep-ph/0107137
// QCD sum rules
m_fplus_0 = 0.61745;
m_fplus_lambda = 1.0108;
m_fplus_m2 = 1.9*1.9;
m_f0_0 = 0.0;
m_f0_lambda = 0.0;
m_f0_m2 = 1.9*1.9;
}
}
m_fplus_0 = model("fplus_0",m_fplus_0);
m_fplus_lambda = model("fplus_lambda",m_fplus_lambda);
m_fplus_m2 = model("fplus_m2",m_fplus_m2);
......
......@@ -131,26 +131,115 @@ PoleFit::PoleFit(GeneralModel model, double* masses,
}
else if (kf0==kf_D_s_plus) {
if (kf1==kf_phi_1020) {
m_V_F0 = 0.9;
switch( int(model("FORM_FACTOR_TUNE", 1)+0.5)) {
case 1:
// hep-lat 1311.6669v1 lattice QCD
m_V_F0 = 1.059;
m_V_m0 = 2.112;
m_V_a0 = -1.0;
m_A0_F0 = 0.706;
m_A0_m0 = 2.459;
m_A0_a0 = -1.0;
m_A1_F0 = 0.615;
m_A1_m0 = 2.459;
m_A1_a0 = -1.0;
m_A2_F0 = 0.457;
m_A2_m0 = 2.459;
m_A2_a0 = -1.0;
break;
case 2:
// hep-ph 0308259 QCD sum rules
m_V_F0 = 1.21;
m_V_m0 = 2.08;
m_V_a0 = -1.0;
m_A0_F0 = 0.42;
m_A0_m0 = 1.9;
m_A0_a0 = -1.0;
m_A1_F0 = 0.55;
m_A1_m0 = 1.9;
m_A1_a0 = -1.0;
m_A2_F0 = 0.59;
m_A2_m0 = 1.9;
m_A2_a0 = -1.0;
break;
case 4:
// hep-ph/0401181 light-cone sum rules
m_V_F0 = 0.9;
m_V_m0 = 1.9685;
m_V_a0 = -2.82;
m_V_b0 = 1.51;
m_A0_F0 = 0.56;
m_A0_m0 = 1.9685;
m_A0_a0 = -0.13;
m_A0_b0 = -0.46;
m_A1_F0 = 0.65;
m_A1_m0 = 1.9685;
m_A1_a0 = -1.36;
m_A1_b0 = -0.31;
m_A2_F0 = 0.85;
m_A2_m0 = 1.9685;
m_A2_a0 = -4.5;
m_A2_b0 = 5.55;
break;
case 5:
// PhysRevD.78.054002 light-front quark model
m_V_F0 = 0.91;
m_V_m0 = 1.9685;
m_V_a0 = -1.15;
m_V_b0 = 0.39;
m_A0_F0 = 0.62;
m_A0_m0 = 1.9685;
m_A0_a0 = -1.02;
m_A0_b0 = 0.45;
m_A1_F0 = 0.61;
m_A1_m0 = 1.9685;
m_A1_a0 = -0.69;
m_A1_b0 = 0.02;
m_A2_F0 = 0.58;
m_A2_m0 = 1.9685;
m_A2_a0 = -0.98;
m_A2_b0 = 0.27;
break;
default:
msg_Error()<<METHOD<<": You chose a form factor tune which does not "
<<"exist. Aborting."<<std::endl;
abort();
}
}
else if (kf1==kf_K_star_892) {
// doi:10.1103/PhysRevD.78.054002 light-front quark model
m_V_F0 = 0.79;
m_V_m0 = 1.9685;
m_V_a0 = -2.82;
m_V_b0 = 1.51;
m_V_a0 = -1.24;
m_V_b0 = 0.50;
m_A0_F0 = 0.56;
m_A0_F0 = 0.54;
m_A0_m0 = 1.9685;
m_A0_a0 = -0.13;
m_A0_b0 = 0.46;
m_A0_a0 = -1.03;
m_A0_b0 = 0.62;
m_A1_F0 = 0.65;
m_A1_F0 = 0.53;
m_A1_m0 = 1.9685;
m_A1_a0 = -1.36;
m_A1_b0 = -0.31;
m_A1_a0 = -0.53;
m_A1_b0 = 0.03;
m_A2_F0 = 0.85;
m_A2_F0 = 0.49;
m_A2_m0 = 1.9685;
m_A2_a0 = -4.5;
m_A2_b0 = 5.55;
}
m_A2_a0 = -1.03;
m_A2_b0 = 0.34;
}
}
m_V_F0 = model("V_F0",m_V_F0); m_V_F1 = model("V_F1",m_V_F1);
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment