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

#include <EvtTensorParticle.hh>

Inheritance diagram for EvtTensorParticle:

Public Member Functions

 EvtTensorParticle ()
virtual ~EvtTensorParticle ()
void init (EvtId part_n, double e, double px, double py, double pz)
void init (EvtId part_n, const EvtVector4R &p4)
EvtTensor4C epsTensorParent (int i) const
EvtTensor4C epsTensor (int i) const
EvtSpinDensity rotateToHelicityBasis () const
EvtSpinDensity rotateToHelicityBasis (double alpha, double beta, double gamma) const
Public Member Functions inherited from EvtParticle
 EvtParticle ()
virtual ~EvtParticle ()
virtual EvtVector4C epsParent (int i) const
virtual EvtVector4C eps (int i) const
virtual EvtVector4C epsParentPhoton (int i)
virtual EvtVector4C epsPhoton (int i)
virtual EvtDiracSpinor spParent (int) const
virtual EvtDiracSpinor sp (int) const
virtual EvtDiracSpinor spParentNeutrino () const
virtual EvtDiracSpinor spNeutrino () const
void addDaug (EvtParticle *node)
void decay ()
void deleteTree ()
void deleteDaughters (bool keepChannel=false)
void setChannel (int i)
void setGeneratorFlag (int flag)
int getGeneratorFlag ()
void setIntFlag (std::vector< int > vi)
std::vector< int > getIntFlag ()
void makeDaughters (int ndaug, EvtId *id)
double initializePhaseSpace (int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
EvtParticlegetDaug (int i)
EvtParticlenextIter (EvtParticle *rootOfTree=0)
void makeStdHep (EvtStdHep &stdhep, EvtSecondary &secondary, EvtId *stable_parent_ihep)
void makeStdHep (EvtStdHep &stdhep)
EvtVector4R getP4Lab ()
EvtVector4R getP4Restframe ()
EvtVector4R get4Pos ()
EvtParticlegetParent ()
void insertDaugPtr (int idaug, EvtParticle *partptr)
double mass () const
int firstornot () const
void setFirstOrNot ()
void resetFirstOrNot ()
EvtId getId () const
EvtSpinType::spintype getSpinType () const
int getSpinStates () const
const EvtVector4RgetP4 () const
void setP4 (const EvtVector4R &p4)
int getChannel () const
int getNDaug () const
void resetNDaug ()
void printTree () const
void printTreeRec (int level) const
std::string writeTreeRec (std::string) const
void dumpTree () const
void dumpTreeRec (int level, int dj) const
std::string treeStr () const
std::string treeStrRec (int level) const
void printParticle () const
void setLifetime (double tau)
void setLifetime ()
double getLifetime ()
void setDiagonalSpinDensity ()
void setVectorSpinDensity ()
void setPolarizedSpinDensity (double r00, double r11, double r22)
void setSpinDensityForward (const EvtSpinDensity &rho)
void setSpinDensityForwardHelicityBasis (const EvtSpinDensity &rho)
void setSpinDensityForwardHelicityBasis (const EvtSpinDensity &rho, double alpha, double beta, double gamma)
EvtSpinDensity getSpinDensityForward ()
void setSpinDensityBackward (const EvtSpinDensity &rho)
EvtSpinDensity getSpinDensityBackward ()
void noLifeTime ()
void setId (EvtId id)
void initDecay (bool useMinMass=false)
void generateMassTree ()
double compMassProb ()
void setMass (double m)
bool isInitialized ()
bool hasValidP4 ()
bool isDecayed ()
double * decayProb ()
void setDecayProb (double p)
void setInclusiveMode (int im)
int getInclusiveMode ()

Additional Inherited Members

Static Public Attributes inherited from EvtParticle
static int _NextLevelDauNum = 0
static EvtId _NextLevelId [20]
static EvtVector4R _NextLevelP4 [20]
Protected Member Functions inherited from EvtParticle
void setp (double e, double px, double py, double pz)
void setp (const EvtVector4R &p4)
void setpart_num (EvtId particle_number)
Protected Attributes inherited from EvtParticle
bool _validP4

Detailed Description

Definition at line 27 of file EvtTensorParticle.hh.

Constructor & Destructor Documentation

◆ EvtTensorParticle()

EvtTensorParticle::EvtTensorParticle ( )
inline

Definition at line 30 of file EvtTensorParticle.hh.

30{}

◆ ~EvtTensorParticle()

EvtTensorParticle::~EvtTensorParticle ( )
virtual

Definition at line 32 of file EvtTensorParticle.cc.

32{}

Member Function Documentation

◆ epsTensor()

EvtTensor4C EvtTensorParticle::epsTensor ( int i) const
virtual

Returns tensor in the particles own restframe for a spin 2 particle.

Reimplemented from EvtParticle.

Definition at line 72 of file EvtTensorParticle.cc.

72{ return eps[i]; } // eps

◆ epsTensorParent()

EvtTensor4C EvtTensorParticle::epsTensorParent ( int i) const
virtual

Returns tensor in the parents restframe for a spin 2 particle.

Reimplemented from EvtParticle.

Definition at line 63 of file EvtTensorParticle.cc.

63 {
64
65 EvtTensor4C temp = eps[i];
66
67 temp.applyBoostTo( this->getP4() );
68 return temp;
69
70} // epsParent
const EvtVector4R & getP4() const
void applyBoostTo(const EvtVector4R &p4)

◆ init() [1/2]

void EvtTensorParticle::init ( EvtId part_n,
const EvtVector4R & p4 )
virtual

Initialiaze particle with id and 4momentum.

Implements EvtParticle.

Definition at line 34 of file EvtTensorParticle.cc.

34 {
35
36 init( part_n, p4.get( 0 ), p4.get( 1 ), p4.get( 2 ), p4.get( 3 ) );
37
39}
void setLifetime()
void init(EvtId part_n, double e, double px, double py, double pz)
double get(int i) const

◆ init() [2/2]

void EvtTensorParticle::init ( EvtId part_n,
double e,
double px,
double py,
double pz )

Definition at line 41 of file EvtTensorParticle.cc.

41 {
42
43 _validP4 = true;
44 setp( e, px, py, pz );
45 setpart_num( part_n );
46
47 eps[0].setdiag( 0.0, -1.0 / sqrt( 6.0 ), -1.0 / sqrt( 6.0 ), 2.0 / sqrt( 6.0 ) );
48 eps[1].setdiag( 0.0, 1.0 / sqrt( 2.0 ), -1.0 / sqrt( 2.0 ), 0.0 );
49 eps[2].setdiag( 0.0, 0.0, 0.0, 0.0 );
50 eps[3].setdiag( 0.0, 0.0, 0.0, 0.0 );
51 eps[4].setdiag( 0.0, 0.0, 0.0, 0.0 );
52
53 eps[2].set( 1, 2, EvtComplex( 1.0 / sqrt( 2.0 ), 0.0 ) );
54 eps[2].set( 2, 1, EvtComplex( 1.0 / sqrt( 2.0 ), 0.0 ) );
55 eps[3].set( 1, 3, EvtComplex( 1.0 / sqrt( 2.0 ), 0.0 ) );
56 eps[3].set( 3, 1, EvtComplex( 1.0 / sqrt( 2.0 ), 0.0 ) );
57 eps[4].set( 2, 3, EvtComplex( 1.0 / sqrt( 2.0 ), 0.0 ) );
58 eps[4].set( 3, 2, EvtComplex( 1.0 / sqrt( 2.0 ), 0.0 ) );
59
61}
void setp(double e, double px, double py, double pz)
void setpart_num(EvtId particle_number)

Referenced by EvtKstarstargamma::decay(), init(), and EvtParticleFactory::particleFactory().

◆ rotateToHelicityBasis() [1/2]

EvtSpinDensity EvtTensorParticle::rotateToHelicityBasis ( ) const
virtual

Returns a rotation matrix need to rotate the basis state to the helicity basis. The EvtSpinDensity matrix is just use as a matrix here. This function is to be implemented in each derived class.

Implements EvtParticle.

Definition at line 74 of file EvtTensorParticle.cc.

74 {
75
76 static EvtVector4C eplus( 0.0, -1.0 / sqrt( 2.0 ), EvtComplex( 0.0, -1.0 / sqrt( 2.0 ) ),
77 0.0 );
78 static EvtVector4C ezero( 0.0, 0.0, 0.0, 1.0 );
79 static EvtVector4C eminus( 0.0, 1.0 / sqrt( 2.0 ), EvtComplex( 0.0, -1.0 / sqrt( 2.0 ) ),
80 0.0 );
81
82 static EvtTensor4C dPpp( directProd( eplus, eplus ) );
83 static EvtTensor4C dPp0( directProd( eplus, ezero ) );
84 static EvtTensor4C dP0p( directProd( ezero, eplus ) );
85 static EvtTensor4C dPpm( directProd( eplus, eminus ) );
86 static EvtTensor4C dP00( directProd( ezero, ezero ) );
87 static EvtTensor4C dPmp( directProd( eminus, eplus ) );
88 static EvtTensor4C dPmm( directProd( eminus, eminus ) );
89 static EvtTensor4C dPm0( directProd( eminus, ezero ) );
90 static EvtTensor4C dP0m( directProd( ezero, eminus ) );
91
92 static EvtTensor4C es0( conj( dPpp ) );
93 static EvtTensor4C es1( conj( ( 1 / sqrt( 2.0 ) ) * dPp0 + ( 1 / sqrt( 2.0 ) ) * dP0p ) );
94 static EvtTensor4C es2( conj( ( 1 / sqrt( 6.0 ) ) * dPpm + ( 2 / sqrt( 6.0 ) ) * dP00 +
95 ( 1 / sqrt( 6.0 ) ) * dPmp ) );
96 static EvtTensor4C es3( conj( ( 1 / sqrt( 2.0 ) ) * dPm0 + ( 1 / sqrt( 2.0 ) ) * dP0m ) );
97 static EvtTensor4C es4( conj( dPmm ) );
98
99 EvtSpinDensity R;
100 R.SetDim( 5 );
101
102 for ( int j = 0; j < 5; j++ )
103 {
104 R.Set( 0, j, cont( es0, eps[j] ) );
105 R.Set( 1, j, cont( es1, eps[j] ) );
106 R.Set( 2, j, cont( es2, eps[j] ) );
107 R.Set( 3, j, cont( es3, eps[j] ) );
108 R.Set( 4, j, cont( es4, eps[j] ) );
109 }
110 return R;
111}
Evt3Rank3C conj(const Evt3Rank3C &t2)
Evt3Rank3C directProd(const EvtVector3C &c1, const EvtVector3C &c2, const EvtVector3C &c3)
EvtComplex cont(const EvtTensor4C &t1, const EvtTensor4C &t2)
complex_t R(double Q2, double M2, double G, double Mp2, double Mm2)
Definition TUtil.h:22

◆ rotateToHelicityBasis() [2/2]

EvtSpinDensity EvtTensorParticle::rotateToHelicityBasis ( double alpha,
double beta,
double gamma ) const
virtual

Implements EvtParticle.

Definition at line 113 of file EvtTensorParticle.cc.

114 {
115
116 EvtTensor4C es[5];
117
118 static EvtVector4C eplus( 0.0, -1.0 / sqrt( 2.0 ), EvtComplex( 0.0, -1.0 / sqrt( 2.0 ) ),
119 0.0 );
120 static EvtVector4C ezero( 0.0, 0.0, 0.0, 1.0 );
121 static EvtVector4C eminus( 0.0, 1.0 / sqrt( 2.0 ), EvtComplex( 0.0, -1.0 / sqrt( 2.0 ) ),
122 0.0 );
123
124 eplus.applyRotateEuler( alpha, beta, gamma );
125 ezero.applyRotateEuler( alpha, beta, gamma );
126 eminus.applyRotateEuler( alpha, beta, gamma );
127
128 for ( int i = 0; i < 5; i++ ) es[i].zero();
129
130 es[0] = directProd( eplus, eplus );
131 es[1] = ( 1 / sqrt( 2.0 ) ) * directProd( eplus, ezero ) +
132 ( 1 / sqrt( 2.0 ) ) * directProd( ezero, eplus );
133 es[2] = ( 1 / sqrt( 6.0 ) ) * directProd( eplus, eminus ) +
134 ( 2 / sqrt( 6.0 ) ) * directProd( ezero, ezero ) +
135 ( 1 / sqrt( 6.0 ) ) * directProd( eminus, eplus );
136 es[3] = ( 1 / sqrt( 2.0 ) ) * directProd( eminus, ezero ) +
137 ( 1 / sqrt( 2.0 ) ) * directProd( ezero, eminus );
138 es[4] = directProd( eminus, eminus );
139
140 for ( int i = 0; i < 5; i++ ) es[i] = conj( es[i] );
141
142 EvtSpinDensity R;
143 R.SetDim( 5 );
144
145 for ( int i = 0; i < 5; i++ )
146 for ( int j = 0; j < 5; j++ ) R.Set( i, j, cont( es[i], eps[j] ) );
147
148 return R;
149}
double alpha
const DifNumber zero

The documentation for this class was generated from the following files: