84 checkSpinParent(EvtSpinType::SCALAR);
86 checkSpinDaughter(1, EvtSpinType::DIRAC);
87 checkSpinDaughter(2, EvtSpinType::NEUTRINO);
89 EvtSpinType::spintype d1type = EvtPDL::getSpinType(getDaug(0));
91 if (d1type == EvtSpinType::VECTOR) {
92 B2INFO(
"EvtBSemiTauonic::init()>> Initializing for decays to a vector type meson ");
94 const double rhoa12 = getArg(0);
95 const double R11 = getArg(1);
96 const double R21 = getArg(2);
97 const double aR3 = getArg(3);
99 const double rho12 = 0;
100 const double aS1 = 0;
102 const double m_b = getArg(4);
103 const double m_c = getArg(5);
105 EvtComplex Coeffs[5];
106 Coeffs[0] = EvtComplex(getArg(6) * cos(getArg(7)), getArg(6) * sin(getArg(7)));
107 Coeffs[1] = EvtComplex(getArg(8) * cos(getArg(9)), getArg(8) * sin(getArg(9)));
108 Coeffs[2] = EvtComplex(getArg(10) * cos(getArg(11)), getArg(10) * sin(getArg(11)));
109 Coeffs[3] = EvtComplex(getArg(12) * cos(getArg(13)), getArg(12) * sin(getArg(13)));
110 Coeffs[4] = EvtComplex(getArg(14) * cos(getArg(15)), getArg(14) * sin(getArg(15)));
113 Coeffs[0], Coeffs[1], Coeffs[2], Coeffs[3], Coeffs[4],
114 EvtPDL::getMeanMass(getParentId()),
116 EvtPDL::getMeanMass(getDaug(0)));
118 }
else if (d1type == EvtSpinType::SCALAR) {
119 B2INFO(
"EvtBSemiTauonic::init()>> Initializing for decays to a scalar type meson ");
121 const double rho12 = getArg(0);
122 const double aS1 = getArg(1);
124 const double rhoa12 = 0;
125 const double R11 = 0;
126 const double R21 = 0;
127 const double aR3 = 0;
129 const double m_b = getArg(2);
130 const double m_c = getArg(3);
132 EvtComplex Coeffs[5];
133 Coeffs[0] = EvtComplex(getArg(4) * cos(getArg(5)), getArg(4) * sin(getArg(5)));
134 Coeffs[1] = EvtComplex(getArg(6) * cos(getArg(7)), getArg(6) * sin(getArg(7)));
135 Coeffs[2] = EvtComplex(getArg(8) * cos(getArg(9)), getArg(8) * sin(getArg(9)));
136 Coeffs[3] = EvtComplex(getArg(10) * cos(getArg(11)), getArg(10) * sin(getArg(11)));
137 Coeffs[4] = EvtComplex(getArg(12) * cos(getArg(13)), getArg(12) * sin(getArg(13)));
140 Coeffs[0], Coeffs[1], Coeffs[2], Coeffs[3], Coeffs[4],
141 EvtPDL::getMeanMass(getParentId()),
142 EvtPDL::getMeanMass(getDaug(0)),
146 B2ERROR(
"BSemiTauonic model handles only scalar and vector meson daughters.");