69 EvtId parnum, mesnum, lnum, nunum;
71 parnum = getParentId();
76 double mymaxprob =
calcamp->CalcMaxProb(parnum, mesnum,
80 static EvtId EM = EvtPDL::getId(
"e-");
81 static EvtId EP = EvtPDL::getId(
"e+");
82 static EvtId MUM = EvtPDL::getId(
"mu-");
83 static EvtId MUP = EvtPDL::getId(
"mu+");
84 static EvtId TAUM = EvtPDL::getId(
"tau-");
85 static EvtId TAUP = EvtPDL::getId(
"tau+");
87 if (lnum == EP || lnum == EM || lnum == MUP || lnum == MUM) {
88 setProbMax(mymaxprob);
91 if (lnum == TAUP || lnum == TAUM) {
108 checkSpinParent(EvtSpinType::SCALAR);
110 checkSpinDaughter(1, EvtSpinType::DIRAC);
111 checkSpinDaughter(2, EvtSpinType::NEUTRINO);
113 EvtSpinType::spintype d1type = EvtPDL::getSpinType(getDaug(0));
114 if (d1type == EvtSpinType::SCALAR) {
115 if (getNArg() == 3) {
117 calcamp =
new EvtSemiLeptonicScalarAmp;
119 EvtGenReport(EVTGEN_ERROR,
"EvtGen") <<
"HQET3 model for scalar meson daughters needs 2 arguments. Sorry." << endl;
122 }
else if (d1type == EvtSpinType::VECTOR) {
123 if (getNArg() == 5) {
125 calcamp =
new EvtSemiLeptonicVectorAmp;
127 EvtGenReport(EVTGEN_ERROR,
"EvtGen") <<
"HQET3 model for vector meson daughtersneeds 4 arguments. Sorry." << endl;
131 EvtGenReport(EVTGEN_ERROR,
"EvtGen") <<
"HQET3 model handles only scalar and vector meson daughters. Sorry." << endl;