1 //--------------------------------------------------------------------------
4 // This software is part of the EvtGen package developed jointly
5 // for the BaBar and CLEO collaborations. If you use all or part
6 // of it, please give an appropriate acknowledgement.
8 // Copyright Information: See EvtGen/COPYRIGHT
9 // Copyright (C) 1998 Caltech, UCSB
11 // Module: EvtSLPole.cc
13 // Description: Routine to implement semileptonic decays according
14 // to light cone sum rules
16 // Modification history:
18 // DJL April 23, 1998 Module created
20 //------------------------------------------------------------------------
22 #include "EvtGenBase/EvtPatches.hh"
24 #include "EvtGenBase/EvtParticle.hh"
25 #include "EvtGenBase/EvtGenKine.hh"
26 #include "EvtGenBase/EvtPDL.hh"
27 #include "EvtGenBase/EvtReport.hh"
28 #include "EvtGenModels/EvtSLPole.hh"
29 #include "EvtGenModels/EvtSLPoleFF.hh"
30 #include "EvtGenBase/EvtSemiLeptonicScalarAmp.hh"
31 #include "EvtGenBase/EvtSemiLeptonicVectorAmp.hh"
32 #include "EvtGenBase/EvtSemiLeptonicTensorAmp.hh"
35 EvtSLPole::~EvtSLPole() {}
37 std::string EvtSLPole::getName(){
44 EvtDecayBase* EvtSLPole::clone(){
50 void EvtSLPole::decay( EvtParticle *p ){
52 p->initializePhaseSpace(getNDaug(),getDaugs(),_resetDaughterTree);
53 calcamp->CalcAmp(p,_amp2,SLPoleffmodel);
57 void EvtSLPole::initProbMax(){
59 EvtId parnum,mesnum,lnum,nunum;
61 parnum = getParentId();
66 double mymaxprob = calcamp->CalcMaxProb(parnum,mesnum,
67 lnum,nunum,SLPoleffmodel);
69 setProbMax(mymaxprob);
74 void EvtSLPole::init(){
78 //We expect the parent to be a scalar
79 //and the daughters to be X lepton neutrino
81 checkSpinParent(EvtSpinType::SCALAR);
82 checkSpinDaughter(1,EvtSpinType::DIRAC);
83 checkSpinDaughter(2,EvtSpinType::NEUTRINO);
85 EvtSpinType::spintype mesontype=EvtPDL::getSpinType(getDaug(0));
87 SLPoleffmodel = new EvtSLPoleFF(getNArg(),getArgs());
89 if ( mesontype==EvtSpinType::SCALAR ) {
90 calcamp = new EvtSemiLeptonicScalarAmp;
92 if ( mesontype==EvtSpinType::VECTOR ) {
93 calcamp = new EvtSemiLeptonicVectorAmp;
95 if ( mesontype==EvtSpinType::TENSOR ) {
96 calcamp = new EvtSemiLeptonicTensorAmp;
99 _resetDaughterTree=false;
100 if ( getArgStr(getNArg()-1) == "true") _resetDaughterTree=true;