BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtVVSPwave Class Reference

#include <EvtVVSPwave.hh>

Inheritance diagram for EvtVVSPwave:

Public Member Functions

 EvtVVSPwave ()
virtual ~EvtVVSPwave ()
void getName (std::string &name)
EvtDecayBase * clone ()
void decay (EvtParticle *p)
void init ()
void initProbMax ()
Public Member Functions inherited from EvtDecayAmp
void makeDecay (EvtParticle *p)
void setWeight (double weight)
void vertex (const EvtComplex &amp)
void vertex (int i1, const EvtComplex &amp)
void vertex (int i1, int i2, const EvtComplex &amp)
void vertex (int i1, int i2, int i3, const EvtComplex &amp)
void vertex (int *i1, const EvtComplex &amp)
virtual ~EvtDecayAmp ()
Public Member Functions inherited from EvtDecayBase
virtual std::string commandName ()
virtual void command (std::string cmd)
double getProbMax (double prob)
double resetProbMax (double prob)
 EvtDecayBase ()
virtual ~EvtDecayBase ()
virtual bool matchingDecay (const EvtDecayBase &other) const
EvtId getParentId ()
double getBranchingFraction ()
void disableCheckQ ()
void checkQ ()
int getNDaug ()
EvtId * getDaugs ()
EvtId getDaug (int i)
int getNArg ()
int getPHOTOS ()
void setPHOTOS ()
void setVerbose ()
void setSummary ()
double * getArgs ()
std::string * getArgsStr ()
double getArg (int j)
std::string getArgStr (int j)
std::string getModelName ()
int getDSum ()
int summary ()
int verbose ()
void saveDecayInfo (EvtId ipar, int ndaug, EvtId *daug, int narg, std::vector< std::string > &args, std::string name, double brfr)
void printSummary ()
void setProbMax (double prbmx)
void noProbMax ()
void checkNArg (int a1, int a2=-1, int a3=-1, int a4=-1)
void checkNDaug (int d1, int d2=-1)
void checkSpinParent (EvtSpinType::spintype sp)
void checkSpinDaughter (int d1, EvtSpinType::spintype sp)
virtual int nRealDaughters ()

Additional Inherited Members

Static Public Member Functions inherited from EvtDecayBase
static void findMasses (EvtParticle *p, int ndaugs, EvtId daugs[10], double masses[10])
static void findMass (EvtParticle *p)
static double findMaxMass (EvtParticle *p)
Protected Member Functions inherited from EvtDecayBase
bool daugsDecayedByParentModel ()
Protected Attributes inherited from EvtDecayAmp
EvtAmp _amp2
Protected Attributes inherited from EvtDecayBase
bool _daugsDecayedByParentModel

Detailed Description

Definition at line 28 of file EvtVVSPwave.hh.

Constructor & Destructor Documentation

◆ EvtVVSPwave()

EvtVVSPwave::EvtVVSPwave ( )
inline

Definition at line 31 of file EvtVVSPwave.hh.

31{}

Referenced by clone().

◆ ~EvtVVSPwave()

EvtVVSPwave::~EvtVVSPwave ( )
virtual

Definition at line 36 of file EvtVVSPwave.cc.

36{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtVVSPwave::clone ( )
virtual

Implements EvtDecayBase.

Definition at line 40 of file EvtVVSPwave.cc.

40{ return new EvtVVSPwave; }

◆ decay()

void EvtVVSPwave::decay ( EvtParticle * p)
virtual

Implements EvtDecayBase.

Definition at line 60 of file EvtVVSPwave.cc.

60 {
61
62 p->initializePhaseSpace( getNDaug(), getDaugs() );
63
64 EvtComplex as( getArg( 0 ) * cos( getArg( 1 ) ), getArg( 0 ) * sin( getArg( 1 ) ) );
65 EvtComplex ap( getArg( 2 ) * cos( getArg( 3 ) ), getArg( 2 ) * sin( getArg( 3 ) ) );
66 EvtComplex ad( getArg( 4 ) * cos( getArg( 5 ) ), getArg( 4 ) * sin( getArg( 5 ) ) );
67
68 EvtParticle *v, *s;
69 v = p->getDaug( 0 );
70 s = p->getDaug( 1 );
71
72 if ( ap != EvtComplex( 0.0, 0.0 ) )
73 { // implement P-wave pingrg,2008-11-24
74 EvtVector4C epi[3], epf[3];
75
76 epi[0] = p->eps( 0 ); // for parent polarization
77 epi[1] = p->eps( 1 );
78 epi[2] = p->eps( 2 );
79
80 epf[0] = v->eps( 0 ); // for sun vector polarization
81 epf[1] = v->eps( 1 );
82 epf[2] = v->eps( 2 );
83
84 EvtVector4R momv = p->getDaug( 0 )->getP4();
85 EvtVector3C p1( momv.get( 1 ), momv.get( 2 ), momv.get( 3 ) );
86
87 EvtVector3C ei0( epi[0].get( 1 ), epi[0].get( 2 ), epi[0].get( 3 ) );
88 EvtVector3C ei1( epi[1].get( 1 ), epi[1].get( 2 ), epi[1].get( 3 ) );
89 EvtVector3C ei2( epi[2].get( 1 ), epi[2].get( 2 ), epi[2].get( 3 ) );
90
91 EvtVector3C ef0( epf[0].get( 1 ), epf[0].get( 2 ), epf[0].get( 3 ) );
92 EvtVector3C ef1( epf[1].get( 1 ), epf[1].get( 2 ), epf[1].get( 3 ) );
93 EvtVector3C ef2( epf[2].get( 1 ), epf[2].get( 2 ), epf[2].get( 3 ) );
94
95 double norm = 1 / momv.d3mag();
96 EvtVector3C q0 = p1.cross( ef0 );
97 EvtVector3C q1 = p1.cross( ef1 );
98 EvtVector3C q2 = p1.cross( ef2 );
99
100 vertex( 0, 0, norm * ei0.conj() * q0 );
101 vertex( 0, 1, norm * ei0.conj() * q1 );
102 vertex( 0, 2, norm * ei0.conj() * q2 );
103
104 vertex( 1, 0, norm * ei1.conj() * q0 );
105 vertex( 1, 1, norm * ei1.conj() * q1 );
106 vertex( 1, 2, norm * ei1.conj() * q2 );
107
108 vertex( 2, 0, norm * ei2.conj() * q0 );
109 vertex( 2, 1, norm * ei2.conj() * q1 );
110 vertex( 2, 2, norm * ei2.conj() * q2 );
111 }
112 else if ( as != EvtComplex( 0.0, 0.0 ) || ad != EvtComplex( 0.0, 0.0 ) )
113 { // implemet s- and d-wave
114
115 // EvtParticle *v,*s;
116 // v=p->getDaug(0);
117 // s=p->getDaug(1);
118
119 EvtTensor4C d, g;
120
121 g.setdiag( 1.0, -1.0, -1.0, -1.0 );
122
123 d = ad * ( ( 1.0 / ( v->getP4().d3mag() * v->getP4().d3mag() ) ) *
124 directProd( v->getP4(), v->getP4() ) +
125 ( 1 / 3.0 ) * g ) +
126 as * g;
127
128 EvtVector4C ep0, ep1, ep2;
129
130 ep0 = d.cont1( p->eps( 0 ) );
131 ep1 = d.cont1( p->eps( 1 ) );
132 ep2 = d.cont1( p->eps( 2 ) );
133
134 vertex( 0, 0, ep0.cont( v->eps( 0 ).conj() ) );
135 vertex( 0, 1, ep0.cont( v->eps( 1 ).conj() ) );
136 vertex( 0, 2, ep0.cont( v->eps( 2 ).conj() ) );
137
138 vertex( 1, 0, ep1.cont( v->eps( 0 ).conj() ) );
139 vertex( 1, 1, ep1.cont( v->eps( 1 ).conj() ) );
140 vertex( 1, 2, ep1.cont( v->eps( 2 ).conj() ) );
141
142 vertex( 2, 0, ep2.cont( v->eps( 0 ).conj() ) );
143 vertex( 2, 1, ep2.cont( v->eps( 1 ).conj() ) );
144 vertex( 2, 2, ep2.cont( v->eps( 2 ).conj() ) );
145 }
146 else
147 {
148 report( ERROR, "EvtGen" ) << "In EvtVectorToVectorScalar.cc" << endl;
149 report( ERROR, "EvtGen" ) << "Not all zero value for parameters are need!!" << endl;
150 ::abort();
151 }
152 return;
153}
p1
double p1[4]
Definition BabayagaNLO.cxx:41
directProd
Evt3Rank3C directProd(const EvtVector3C &c1, const EvtVector3C &c2, const EvtVector3C &c3)
report
ostream & report(Severity severity, const char *facility)
Definition EvtReport.cc:34
ERROR
@ ERROR
Definition EvtReport.hh:49
s
XmlRpcServer s
Definition HelloServer.cpp:11
sin
double sin(const BesAngle a)
Definition InstallArea/x86_64-el9-gcc13-dbg/include/MdcGeom/BesAngle.h:185
cos
double cos(const BesAngle a)
Definition InstallArea/x86_64-el9-gcc13-dbg/include/MdcGeom/BesAngle.h:187
v
**********Class see also m_nmax DOUBLE PRECISION m_amel DOUBLE PRECISION m_x2 DOUBLE PRECISION m_alfinv DOUBLE PRECISION m_Xenph INTEGER m_KeyWtm INTEGER m_idyfs DOUBLE PRECISION m_zini DOUBLE PRECISION m_q2 DOUBLE PRECISION m_Wt_KF DOUBLE PRECISION m_WtCut INTEGER m_KFfin *COMMON c_KarLud $ !Input CMS energy[GeV] $ !CMS energy after beam spread beam strahlung[GeV] $ !Beam energy spread[GeV] $ !z boost due to beam spread $ !electron beam mass *ff pair spectrum $ !minimum v
Definition KarLud.h:35
EvtDecayAmp::vertex
void vertex(const EvtComplex &amp)
Definition EvtDecayAmp.hh:36
EvtDecayBase::getArg
double getArg(int j)
Definition EvtDecayBase.cc:520
EvtDecayBase::getDaugs
EvtId * getDaugs()
Definition EvtDecayBase.hh:68
EvtParticle::getP4
const EvtVector4R & getP4() const
Definition EvtParticle.cc:108
EvtParticle::getDaug
EvtParticle * getDaug(int i)
Definition EvtParticle.cc:80
EvtParticle::eps
virtual EvtVector4C eps(int i) const
Definition EvtParticle.cc:592
EvtParticle::initializePhaseSpace
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition EvtParticle.cc:1042
EvtTensor4C::setdiag
void setdiag(double t00, double t11, double t22, double t33)
Definition EvtTensor4C.cc:190
EvtTensor4C::cont1
EvtVector4C cont1(const EvtVector4C &v4) const
Definition EvtTensor4C.cc:410
EvtVector4C::cont
EvtComplex cont(const EvtVector4C &v4) const
Definition EvtVector4C.hh:137
EvtVector4R::get
double get(int i) const
Definition EvtVector4R.hh:171
EvtVector4R::d3mag
double d3mag() const
Definition EvtVector4R.cc:168

◆ getName()

void EvtVVSPwave::getName ( std::string & name)
virtual

Implements EvtDecayBase.

Definition at line 38 of file EvtVVSPwave.cc.

38{ model_name = "VVS_PWAVE"; }

◆ init()

void EvtVVSPwave::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 42 of file EvtVVSPwave.cc.

42 {
43
44 // check that there are 6 arguments
45 checkNArg( 6 );
46 checkNDaug( 2 );
47
48 checkSpinParent( EvtSpinType::VECTOR );
49 checkSpinDaughter( 0, EvtSpinType::VECTOR );
50 checkSpinDaughter( 1, EvtSpinType::SCALAR );
51}
EvtDecayBase::checkSpinDaughter
void checkSpinDaughter(int d1, EvtSpinType::spintype sp)
Definition EvtDecayBase.cc:492
EvtDecayBase::checkSpinParent
void checkSpinParent(EvtSpinType::spintype sp)
Definition EvtDecayBase.cc:479
EvtDecayBase::checkNDaug
void checkNDaug(int d1, int d2=-1)
Definition EvtDecayBase.cc:465
EvtDecayBase::checkNArg
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
Definition EvtDecayBase.cc:449

◆ initProbMax()

void EvtVVSPwave::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 53 of file EvtVVSPwave.cc.

53 {
54
55 // probmax is 1.0 for all possible decays I think!
56
57 setProbMax( 1.0 );
58}
EvtDecayBase::setProbMax
void setProbMax(double prbmx)
Definition EvtDecayBase.cc:261
The documentation for this class was generated from the following files: