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) 2004 Cornell
11 // Module: EvtVPHOtoV.cc
13 // Description: Routine to decay vpho -> vector
15 // Modification history:
17 // Ryd March 9, 2004 Module created
19 //------------------------------------------------------------------------
22 #include "EvtGenBase/EvtParticle.hh"
23 #include "EvtGenBase/EvtGenKine.hh"
24 #include "EvtGenBase/EvtPDL.hh"
25 #include "EvtGenBase/EvtVector4C.hh"
26 #include "EvtGenBase/EvtVector4R.hh"
27 #include "EvtGenBase/EvtReport.hh"
28 #include "EvtGenModels/EvtVPHOtoV.hh"
31 EvtVPHOtoV::~EvtVPHOtoV() {}
33 std::string EvtVPHOtoV::getName(){
40 EvtDecayBase* EvtVPHOtoV::clone(){
42 return new EvtVPHOtoV;
46 void EvtVPHOtoV::init(){
48 // check that there are 0 arguments
51 // check that there are 1 daughters
54 // check the parent and daughter spins
55 checkSpinParent(EvtSpinType::VECTOR);
56 checkSpinDaughter(0,EvtSpinType::VECTOR);
59 void EvtVPHOtoV::initProbMax() {
65 void EvtVPHOtoV::decay( EvtParticle *p){
67 p->initializePhaseSpace(getNDaug(),getDaugs());
69 EvtParticle* d=p->getDaug(0);
71 d->setP4(p->getP4Restframe());
73 vertex(0,0,p->eps(0)*p->epsParent(0).conj());
74 vertex(1,0,p->eps(1)*p->epsParent(0).conj());
75 vertex(2,0,p->eps(2)*p->epsParent(0).conj());
77 vertex(0,1,p->eps(0)*p->epsParent(1).conj());
78 vertex(1,1,p->eps(1)*p->epsParent(1).conj());
79 vertex(2,1,p->eps(2)*p->epsParent(1).conj());
81 vertex(0,2,p->eps(0)*p->epsParent(2).conj());
82 vertex(1,2,p->eps(1)*p->epsParent(2).conj());
83 vertex(2,2,p->eps(2)*p->epsParent(2).conj());