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
13 // Description: Routine to decay scalar -> vector + scalar
15 // Modification history:
17 // RYD November 24, 1996 Module created
19 //------------------------------------------------------------------------
21 #include "EvtGenBase/EvtPatches.hh"
23 #include "EvtGenBase/EvtParticle.hh"
24 #include "EvtGenBase/EvtGenKine.hh"
25 #include "EvtGenBase/EvtReport.hh"
26 #include "EvtGenBase/EvtPDL.hh"
27 #include "EvtGenBase/EvtVector4C.hh"
28 #include "EvtGenModels/EvtSVS.hh"
33 std::string EvtSVS::getName(){
40 EvtDecayBase* EvtSVS::clone(){
48 // check that there are 0 arguments
52 checkSpinParent(EvtSpinType::SCALAR);
54 checkSpinDaughter(0,EvtSpinType::VECTOR);
55 checkSpinDaughter(1,EvtSpinType::SCALAR);
62 void EvtSVS::initProbMax(){
68 void EvtSVS::decay( EvtParticle *p){
70 p->initializePhaseSpace(getNDaug(),getDaugs());
74 double massv = v->mass();
75 EvtVector4R momv = v->getP4();
76 EvtVector4R moms = p->getDaug(1)->getP4();
77 double m_parent = p->mass();
78 EvtVector4R p4_parent = momv+moms;
80 double norm=massv/(momv.d3mag()*m_parent);
81 p4_parent = norm*p4_parent;
82 vertex(0,p4_parent*(v->epsParent(0)));
83 vertex(1,p4_parent*(v->epsParent(1)));
84 vertex(2,p4_parent*(v->epsParent(2)));