BOSS 7.0.4
BESIII Offline Software System
Loading...
Searching...
No Matches
Public Member Functions | List of all members
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 void getName (std::string &name)=0
 
virtual void decay (EvtParticle *p)=0
 
virtual void makeDecay (EvtParticle *p)=0
 
virtual EvtDecayBase * clone ()=0
 
virtual void init ()
 
virtual void initProbMax ()
 
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 32 of file EvtVVSPwave.hh.

32{}

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 45 of file EvtVVSPwave.cc.

45 {
46
47 return new EvtVVSPwave;
48
49}

◆ decay()

void EvtVVSPwave::decay ( EvtParticle *  p)
virtual

Implements EvtDecayBase.

Definition at line 70 of file EvtVVSPwave.cc.

70 {
71
72 p->initializePhaseSpace(getNDaug(),getDaugs());
73
74 EvtComplex as(getArg(0)*cos(getArg(1)),getArg(0)*sin(getArg(1)));
75 EvtComplex ap(getArg(2)*cos(getArg(3)),getArg(2)*sin(getArg(3)));
76 EvtComplex ad(getArg(4)*cos(getArg(5)),getArg(4)*sin(getArg(5)));
77
78 EvtParticle *v,*s;
79 v=p->getDaug(0);
80 s=p->getDaug(1);
81
82 if (ap!=EvtComplex(0.0,0.0)) { // implement P-wave pingrg,2008-11-24
83 EvtVector4C epi[3],epf[3];
84
85 epi[0] = p->eps(0); //for parent polarization
86 epi[1] = p->eps(1);
87 epi[2] = p->eps(2);
88
89 epf[0] = v->eps(0); //for sun vector polarization
90 epf[1] = v->eps(1);
91 epf[2] = v->eps(2);
92
93 EvtVector4R momv = p->getDaug(0)->getP4();
94 EvtVector3C p1(momv.get(1),momv.get(2),momv.get(3));
95
96 EvtVector3C ei0(epi[0].get(1),epi[0].get(2),epi[0].get(3));
97 EvtVector3C ei1(epi[1].get(1),epi[1].get(2),epi[1].get(3));
98 EvtVector3C ei2(epi[2].get(1),epi[2].get(2),epi[2].get(3));
99
100 EvtVector3C ef0(epf[0].get(1),epf[0].get(2),epf[0].get(3));
101 EvtVector3C ef1(epf[1].get(1),epf[1].get(2),epf[1].get(3));
102 EvtVector3C ef2(epf[2].get(1),epf[2].get(2),epf[2].get(3));
103
104 double norm = 1/momv.d3mag();
105 EvtVector3C q0 = p1.cross(ef0);
106 EvtVector3C q1 = p1.cross(ef1);
107 EvtVector3C q2 = p1.cross(ef2);
108
109 vertex(0,0,norm*ei0.conj()*q0);
110 vertex(0,1,norm*ei0.conj()*q1);
111 vertex(0,2,norm*ei0.conj()*q2);
112
113 vertex(1,0,norm*ei1.conj()*q0);
114 vertex(1,1,norm*ei1.conj()*q1);
115 vertex(1,2,norm*ei1.conj()*q2);
116
117 vertex(2,0,norm*ei2.conj()*q0);
118 vertex(2,1,norm*ei2.conj()*q1);
119 vertex(2,2,norm*ei2.conj()*q2);
120
121 } else if (as!=EvtComplex(0.0,0.0) || ad!=EvtComplex(0.0,0.0) ){// implemet s- and d-wave
122
123 // EvtParticle *v,*s;
124 // v=p->getDaug(0);
125 // s=p->getDaug(1);
126
127 EvtTensor4C d,g;
128
129 g.setdiag(1.0,-1.0,-1.0,-1.0);
130
131 d=ad*((1.0/(v->getP4().d3mag()*v->getP4().d3mag()))*directProd(v->getP4(),v->getP4())+(1/3.0)*g)+
132 as*g;
133
134 EvtVector4C ep0,ep1,ep2;
135
136 ep0=d.cont1(p->eps(0));
137 ep1=d.cont1(p->eps(1));
138 ep2=d.cont1(p->eps(2));
139
140 vertex(0,0,ep0.cont(v->eps(0).conj()));
141 vertex(0,1,ep0.cont(v->eps(1).conj()));
142 vertex(0,2,ep0.cont(v->eps(2).conj()));
143
144 vertex(1,0,ep1.cont(v->eps(0).conj()));
145 vertex(1,1,ep1.cont(v->eps(1).conj()));
146 vertex(1,2,ep1.cont(v->eps(2).conj()));
147
148 vertex(2,0,ep2.cont(v->eps(0).conj()));
149 vertex(2,1,ep2.cont(v->eps(1).conj()));
150 vertex(2,2,ep2.cont(v->eps(2).conj()));
151 } else
152 {
153 report(ERROR,"EvtGen") << "In EvtVectorToVectorScalar.cc"<<endl;
154 report(ERROR,"EvtGen") << "Not all zero value for parameters are need!!"<<endl;
155 ::abort();
156 }
157 return ;
158
159}
directProd
Evt3Rank3C directProd(const EvtVector3C &c1, const EvtVector3C &c2, const EvtVector3C &c3)
report
ostream & report(Severity severity, const char *facility)
Definition: EvtReport.cc:36
ERROR
@ ERROR
Definition: EvtReport.hh:49
s
XmlRpcServer s
Definition: HelloServer.cpp:11
sin
double sin(const BesAngle a)
Definition: InstallArea/include/MdcGeom/MdcGeom/BesAngle.h:210
cos
double cos(const BesAngle a)
Definition: InstallArea/include/MdcGeom/MdcGeom/BesAngle.h:213
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:37
EvtDecayBase::getArg
double getArg(int j)
Definition: EvtDecayBase.cc:564
EvtDecayBase::getDaugs
EvtId * getDaugs()
Definition: EvtDecayBase.hh:65
EvtParticle::getP4
const EvtVector4R & getP4() const
Definition: EvtParticle.cc:121
EvtParticle::getDaug
EvtParticle * getDaug(int i)
Definition: EvtParticle.cc:85
EvtParticle::eps
virtual EvtVector4C eps(int i) const
Definition: EvtParticle.cc:576
EvtParticle::initializePhaseSpace
double initializePhaseSpace(int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
Definition: EvtParticle.cc:1071
EvtTensor4C::setdiag
void setdiag(double t00, double t11, double t22, double t33)
Definition: EvtTensor4C.cc:218
EvtTensor4C::cont1
EvtVector4C cont1(const EvtVector4C &v4) const
Definition: EvtTensor4C.cc:461
EvtVector3C::cross
EvtVector3C cross(const EvtVector3C &v2)
Definition: EvtVector3C.cc:48
EvtVector4C::cont
EvtComplex cont(const EvtVector4C &v4) const
Definition: EvtVector4C.hh:146
EvtVector4R::get
double get(int i) const
Definition: EvtVector4R.hh:179
EvtVector4R::d3mag
double d3mag() const
Definition: EvtVector4R.cc:186

◆ getName()

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

Implements EvtDecayBase.

Definition at line 38 of file EvtVVSPwave.cc.

38 {
39
40 model_name="VVS_PWAVE";
41
42}

◆ init()

void EvtVVSPwave::init ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 51 of file EvtVVSPwave.cc.

51 {
52
53 // check that there are 6 arguments
54 checkNArg(6);
55 checkNDaug(2);
56
57 checkSpinParent(EvtSpinType::VECTOR);
58 checkSpinDaughter(0,EvtSpinType::VECTOR);
59 checkSpinDaughter(1,EvtSpinType::SCALAR);
60}
EvtDecayBase::checkSpinDaughter
void checkSpinDaughter(int d1, EvtSpinType::spintype sp)
Definition: EvtDecayBase.cc:533
EvtDecayBase::checkSpinParent
void checkSpinParent(EvtSpinType::spintype sp)
Definition: EvtDecayBase.cc:520
EvtDecayBase::checkNDaug
void checkNDaug(int d1, int d2=-1)
Definition: EvtDecayBase.cc:504
EvtDecayBase::checkNArg
void checkNArg(int a1, int a2=-1, int a3=-1, int a4=-1)
Definition: EvtDecayBase.cc:482

◆ initProbMax()

void EvtVVSPwave::initProbMax ( )
virtual

Reimplemented from EvtDecayBase.

Definition at line 62 of file EvtVVSPwave.cc.

62 {
63
64 //probmax is 1.0 for all possible decays I think!
65
66 setProbMax(1.0);
67
68}
EvtDecayBase::setProbMax
void setProbMax(double prbmx)
Definition: EvtDecayBase.cc:297
The documentation for this class was generated from the following files: