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: EvtVSPPwave.cc
13 // Description: Routine to decay vector-> scalar photon in P-wave
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/EvtTensor4C.hh"
26 #include "EvtGenBase/EvtPDL.hh"
27 #include "EvtGenModels/EvtVSPPwave.hh"
28 #include "EvtGenBase/EvtReport.hh"
30 #include "EvtGenBase/EvtVector4C.hh"
32 EvtVSPPwave::~EvtVSPPwave() {}
34 std::string EvtVSPPwave::getName(){
41 EvtDecayBase* EvtVSPPwave::clone(){
43 return new EvtVSPPwave;
47 void EvtVSPPwave::init(){
49 // check that there are 0 arguments
53 checkSpinParent(EvtSpinType::VECTOR);
55 checkSpinDaughter(0,EvtSpinType::SCALAR);
56 checkSpinDaughter(1,EvtSpinType::PHOTON);
60 void EvtVSPPwave::initProbMax(){
65 void EvtVSPPwave::decay( EvtParticle *p){
67 p->initializePhaseSpace(getNDaug(),getDaugs());
70 gamma = p->getDaug(1);
73 EvtVector4R momgamma = gamma->getP4();
75 //work in the parent ,p, rest frame.
77 p4_p.set(m_p, 0.0, 0.0, 0.0 );
79 // Put phase space results into the daughters.
83 double norm=1/(m_p*momgamma.d3mag());
85 tds = dual(EvtGenFunctions::directProd(norm*p4_p,momgamma));
87 vertex(0,0,(tds.cont1( p->eps(0))).cont(
88 gamma->epsParentPhoton(0).conj() ) );
89 vertex(0,1,(tds.cont1( p->eps(0))).cont(
90 gamma->epsParentPhoton(1).conj() ) );
92 vertex(1,0,(tds.cont1( p->eps(1))).cont(
93 gamma->epsParentPhoton(0).conj() ) );
94 vertex(1,1,(tds.cont1( p->eps(1))).cont(
95 gamma->epsParentPhoton(1).conj() ) );
97 vertex(2,0,(tds.cont1( p->eps(2))).cont(
98 gamma->epsParentPhoton(0).conj() ) );
99 vertex(2,1,(tds.cont1( p->eps(2))).cont(
100 gamma->epsParentPhoton(1).conj() ) );