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 tensor-> scalar scalar
15 // Modification history:
17 // RYD November 24, 1996 Module created
19 //------------------------------------------------------------------------
22 #include "EvtGenBase/EvtPatches.hh"
24 #include "EvtGenBase/EvtVector4C.hh"
25 #include "EvtGenBase/EvtParticle.hh"
26 #include "EvtGenBase/EvtGenKine.hh"
27 #include "EvtGenBase/EvtPDL.hh"
28 #include "EvtGenBase/EvtTensor4C.hh"
29 #include "EvtGenBase/EvtReport.hh"
30 #include "EvtGenModels/EvtTSS.hh"
35 std::string EvtTSS::getName(){
42 EvtDecayBase* EvtTSS::clone(){
50 // check that there are 0 arguments
55 checkSpinParent(EvtSpinType::TENSOR);
57 checkSpinDaughter(0,EvtSpinType::SCALAR);
58 checkSpinDaughter(1,EvtSpinType::SCALAR);
63 void EvtTSS::initProbMax() {
69 void EvtTSS::decay( EvtParticle *p){
71 p->initializePhaseSpace(getNDaug(),getDaugs());
73 EvtVector4R moms1 = p->getDaug(0)->getP4();
75 double norm = 1.0/(moms1.d3mag()*moms1.d3mag());
77 vertex(0,norm*(p->epsTensor(0).cont1(EvtVector4C(moms1))*(moms1)));
78 vertex(1,norm*(p->epsTensor(1).cont1(EvtVector4C(moms1))*(moms1)));
79 vertex(2,norm*(p->epsTensor(2).cont1(EvtVector4C(moms1))*(moms1)));
80 vertex(3,norm*(p->epsTensor(3).cont1(EvtVector4C(moms1))*(moms1)));
81 vertex(4,norm*(p->epsTensor(4).cont1(EvtVector4C(moms1))*(moms1)));