BOSS 7.0.9
BESIII Offline Software System
Loading...
Searching...
No Matches
Public Member Functions | Friends | List of all members
EvtDiracSpinor Class Reference

#include <EvtDiracSpinor.hh>

Public Member Functions

 EvtDiracSpinor ()
 
 EvtDiracSpinor (const EvtComplex &sp0, const EvtComplex &sp1, const EvtComplex &sp2, const EvtComplex &sp3)
 
virtual ~EvtDiracSpinor ()
 
 EvtDiracSpinor (const EvtDiracSpinor &dspinor)
 
EvtDiracSpinor & operator= (const EvtDiracSpinor &dspinor)
 
EvtDiracSpinor & operator+= (const EvtDiracSpinor &u2)
 
EvtDiracSpinor & operator-= (const EvtDiracSpinor &u2)
 
void set (const EvtComplex &sp0, const EvtComplex &sp1, const EvtComplex &sp2, const EvtComplex &sp3)
 
void set_spinor (int i, const EvtComplex &sp)
 
const EvtComplex & get_spinor (int i) const
 
EvtDiracSpinor conj () const
 
void applyRotateEuler (double alpha, double beta, double gamma)
 
void applyBoostTo (const EvtVector4R &p4)
 
void applyBoostTo (const EvtVector3R &boost)
 
EvtDiracSpinor adjoint () const
 
 EvtDiracSpinor ()
 
 EvtDiracSpinor (const EvtComplex &sp0, const EvtComplex &sp1, const EvtComplex &sp2, const EvtComplex &sp3)
 
virtual ~EvtDiracSpinor ()
 
 EvtDiracSpinor (const EvtDiracSpinor &dspinor)
 
EvtDiracSpinor & operator= (const EvtDiracSpinor &dspinor)
 
EvtDiracSpinor & operator+= (const EvtDiracSpinor &u2)
 
EvtDiracSpinor & operator-= (const EvtDiracSpinor &u2)
 
void set (const EvtComplex &sp0, const EvtComplex &sp1, const EvtComplex &sp2, const EvtComplex &sp3)
 
void set_spinor (int i, const EvtComplex &sp)
 
const EvtComplex & get_spinor (int i) const
 
EvtDiracSpinor conj () const
 
void applyRotateEuler (double alpha, double beta, double gamma)
 
void applyBoostTo (const EvtVector4R &p4)
 
void applyBoostTo (const EvtVector3R &boost)
 
EvtDiracSpinor adjoint () const
 

Friends

EvtDiracSpinor rotateEuler (const EvtDiracSpinor &sp, double alpha, double beta, double gamma)
 
EvtDiracSpinor boostTo (const EvtDiracSpinor &sp, const EvtVector4R p4)
 
EvtDiracSpinor boostTo (const EvtDiracSpinor &sp, const EvtVector3R boost)
 
EvtVector4C EvtLeptonVACurrent (const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
 
EvtVector4C EvtLeptonVCurrent (const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
 
EvtVector4C EvtLeptonACurrent (const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
 
EvtComplex EvtLeptonSCurrent (const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
 
EvtComplex EvtLeptonPCurrent (const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
 
EvtTensor4C EvtLeptonTCurrent (const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
 
EvtTensor4C EvtLeptonTg5Current (const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
 
EvtDiracSpinor operator+ (const EvtDiracSpinor &u1, const EvtDiracSpinor &u2)
 
EvtDiracSpinor operator- (const EvtDiracSpinor &u1, const EvtDiracSpinor &u2)
 
EvtDiracSpinor operator* (const EvtComplex &c, const EvtDiracSpinor &d)
 
std::ostream & operator<< (std::ostream &s, const EvtDiracSpinor &c)
 
EvtDiracSpinor rotateEuler (const EvtDiracSpinor &sp, double alpha, double beta, double gamma)
 
EvtDiracSpinor boostTo (const EvtDiracSpinor &sp, const EvtVector4R p4)
 
EvtDiracSpinor boostTo (const EvtDiracSpinor &sp, const EvtVector3R boost)
 
EvtVector4C EvtLeptonVACurrent (const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
 
EvtVector4C EvtLeptonVCurrent (const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
 
EvtVector4C EvtLeptonACurrent (const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
 
EvtComplex EvtLeptonSCurrent (const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
 
EvtComplex EvtLeptonPCurrent (const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
 
EvtTensor4C EvtLeptonTCurrent (const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
 
EvtTensor4C EvtLeptonTg5Current (const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
 
EvtDiracSpinor operator+ (const EvtDiracSpinor &u1, const EvtDiracSpinor &u2)
 
EvtDiracSpinor operator- (const EvtDiracSpinor &u1, const EvtDiracSpinor &u2)
 
EvtDiracSpinor operator* (const EvtComplex &c, const EvtDiracSpinor &d)
 
std::ostream & operator<< (std::ostream &s, const EvtDiracSpinor &c)
 

Detailed Description

Definition at line 52 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.hh.

Constructor & Destructor Documentation

◆ EvtDiracSpinor() [1/6]

EvtDiracSpinor::EvtDiracSpinor ( )
inline

Definition at line 111 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.hh.

111 {
112
113 spinor[0]=EvtComplex(); spinor[1]=EvtComplex();
114 spinor[2]=EvtComplex(); spinor[3]=EvtComplex();
115
116}

◆ EvtDiracSpinor() [2/6]

EvtDiracSpinor::EvtDiracSpinor ( const EvtComplex &  sp0,
const EvtComplex &  sp1,
const EvtComplex &  sp2,
const EvtComplex &  sp3 
)

Definition at line 35 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

36 {
37 set(sp0,sp1,sp2,sp3);
38}
EvtDiracSpinor::set
void set(const EvtComplex &sp0, const EvtComplex &sp1, const EvtComplex &sp2, const EvtComplex &sp3)
Definition: bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc:40

◆ ~EvtDiracSpinor() [1/2]

EvtDiracSpinor::~EvtDiracSpinor ( )
virtual

Definition at line 33 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

33{}

◆ EvtDiracSpinor() [3/6]

EvtDiracSpinor::EvtDiracSpinor ( const EvtDiracSpinor &  dspinor)
inline

Definition at line 118 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.hh.

118 {
119
120 spinor[0]=dspinor.spinor[0];
121 spinor[1]=dspinor.spinor[1];
122 spinor[2]=dspinor.spinor[2];
123 spinor[3]=dspinor.spinor[3];
124
125}

◆ EvtDiracSpinor() [4/6]

EvtDiracSpinor::EvtDiracSpinor ( )
inline

◆ EvtDiracSpinor() [5/6]

EvtDiracSpinor::EvtDiracSpinor ( const EvtComplex &  sp0,
const EvtComplex &  sp1,
const EvtComplex &  sp2,
const EvtComplex &  sp3 
)

◆ ~EvtDiracSpinor() [2/2]

virtual EvtDiracSpinor::~EvtDiracSpinor ( )
virtual

◆ EvtDiracSpinor() [6/6]

EvtDiracSpinor::EvtDiracSpinor ( const EvtDiracSpinor &  dspinor)
inline

Member Function Documentation

◆ adjoint() [1/2]

EvtDiracSpinor EvtDiracSpinor::adjoint ( ) const

Definition at line 387 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

388{
389 EvtDiracSpinor d = this->conj(); // first conjugate, then multiply with gamma0
390 EvtGammaMatrix g0 = EvtGammaMatrix::g0();
391 EvtDiracSpinor result; // automatically initialized to 0
392
393 for (int i=0; i<4; ++i)
394 for (int j=0; j<4; ++j)
395 result.spinor[i] += d.spinor[j] * g0.gamma[i][j];
396
397 return result;
398}

Referenced by EvtLambdaP_BarGamma::decay().

◆ adjoint() [2/2]

EvtDiracSpinor EvtDiracSpinor::adjoint ( ) const

◆ applyBoostTo() [1/4]

void EvtDiracSpinor::applyBoostTo ( const EvtVector3R &  boost)

Definition at line 107 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

107 {
108
109 double bx,by,bz,gamma,b2,f1,f2;
110 EvtComplex spinorp[4];
111
112 bx=boost.get(0);
113 by=boost.get(1);
114 bz=boost.get(2);
115 b2=bx*bx+by*by+bz*bz;
116
117 if (b2==0.0){
118 return;
119 }
120
121 assert(b2<1.0);
122
123 gamma=1.0/sqrt(1-b2);
124
125 f1=sqrt((gamma+1.0)/2.0);
126 f2=f1*gamma/(gamma+1.0);
127
128 spinorp[0]=f1*spinor[0]+f2*bz*spinor[2]+
129 f2*EvtComplex(bx,-by)*spinor[3];
130 spinorp[1]=f1*spinor[1]+f2*EvtComplex(bx,by)*spinor[2]-
131 f2*bz*spinor[3];
132 spinorp[2]=f2*bz*spinor[0]+f2*EvtComplex(bx,-by)*spinor[1]+
133 f1*spinor[2];
134 spinorp[3]=f2*EvtComplex(bx,by)*spinor[0]-
135 f2*bz*spinor[1]+f1*spinor[3];
136
137 spinor[0]=spinorp[0];
138 spinor[1]=spinorp[1];
139 spinor[2]=spinorp[2];
140 spinor[3]=spinorp[3];
141
142 return;
143}
f1
TFile * f1
Definition: DataBase/tau_mode.c:4

◆ applyBoostTo() [2/4]

void EvtDiracSpinor::applyBoostTo ( const EvtVector3R &  boost)

◆ applyBoostTo() [3/4]

void EvtDiracSpinor::applyBoostTo ( const EvtVector4R &  p4)

Definition at line 93 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

93 {
94
95 double e=p4.get(0);
96
97 EvtVector3R boost(p4.get(1)/e,p4.get(2)/e,p4.get(3)/e);
98
99 applyBoostTo(boost);
100
101 return;
102
103}

Referenced by EvtRaritaSchwinger::applyBoostTo(), applyBoostTo(), and boostTo().

◆ applyBoostTo() [4/4]

void EvtDiracSpinor::applyBoostTo ( const EvtVector4R &  p4)

◆ applyRotateEuler() [1/2]

void EvtDiracSpinor::applyRotateEuler ( double  alpha,
double  beta,
double  gamma 
)

Definition at line 145 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

146 {
147
148 EvtComplex retVal[4];
149
150 double cb2=cos(0.5*beta);
151 double sb2=sin(0.5*beta);
152 double capg2=cos(0.5*(alpha+gamma));
153 double camg2=cos(0.5*(alpha-gamma));
154 double sapg2=sin(0.5*(alpha+gamma));
155 double samg2=sin(0.5*(alpha-gamma));
156
157 EvtComplex m11(cb2*capg2,-cb2*sapg2);
158 EvtComplex m12(-sb2*camg2,sb2*samg2);
159 EvtComplex m21(sb2*camg2,sb2*samg2);
160 EvtComplex m22(cb2*capg2,cb2*sapg2);
161
162 retVal[0]=m11*spinor[0]+m12*spinor[1];
163 retVal[1]=m21*spinor[0]+m22*spinor[1];
164 retVal[2]=m11*spinor[2]+m12*spinor[3];
165 retVal[3]=m21*spinor[2]+m22*spinor[3];
166
167 spinor[0]=retVal[0];
168 spinor[1]=retVal[1];
169 spinor[2]=retVal[2];
170 spinor[3]=retVal[3];
171
172 return;
173}
sin
double sin(const BesAngle a)
Definition: BesAngle.h:210
cos
double cos(const BesAngle a)
Definition: BesAngle.h:213
alpha
const double alpha
Definition: FastVertexFit.cxx:4

Referenced by operator*(), rotateEuler(), and EvtDiracParticle::rotateToHelicityBasis().

◆ applyRotateEuler() [2/2]

void EvtDiracSpinor::applyRotateEuler ( double  alpha,
double  beta,
double  gamma 
)

◆ conj() [1/2]

EvtDiracSpinor EvtDiracSpinor::conj ( ) const

Definition at line 177 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

177 {
178
179 EvtDiracSpinor sp;
180
181 for ( int i=0; i<4; i++)
182 sp.set_spinor(i,::conj(spinor[i]));
183
184 return sp;
185}

Referenced by adjoint(), and conj().

◆ conj() [2/2]

EvtDiracSpinor EvtDiracSpinor::conj ( ) const

◆ get_spinor() [1/2]

const EvtComplex & EvtDiracSpinor::get_spinor ( int  i) const

Definition at line 60 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

60 {
61
62 return spinor[i];
63
64}

Referenced by dirProd(), operator*(), and EvtRaritaSchwinger::setSpinor().

◆ get_spinor() [2/2]

const EvtComplex & EvtDiracSpinor::get_spinor ( int  i) const

◆ operator+=() [1/2]

EvtDiracSpinor & EvtDiracSpinor::operator+= ( const EvtDiracSpinor &  u2)
inline

Definition at line 138 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.hh.

138 {
139
140 spinor[0]+=u2.spinor[0];
141 spinor[1]+=u2.spinor[1];
142 spinor[2]+=u2.spinor[2];
143 spinor[3]+=u2.spinor[3];
144
145 return *this;
146}

◆ operator+=() [2/2]

EvtDiracSpinor & EvtDiracSpinor::operator+= ( const EvtDiracSpinor &  u2)
inline

◆ operator-=() [1/2]

EvtDiracSpinor & EvtDiracSpinor::operator-= ( const EvtDiracSpinor &  u2)
inline

Definition at line 155 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.hh.

155 {
156
157 spinor[0]-=u2.spinor[0];
158 spinor[1]-=u2.spinor[1];
159 spinor[2]-=u2.spinor[2];
160 spinor[3]-=u2.spinor[3];
161
162 return *this;
163}

◆ operator-=() [2/2]

EvtDiracSpinor & EvtDiracSpinor::operator-= ( const EvtDiracSpinor &  u2)
inline

◆ operator=() [1/2]

EvtDiracSpinor & EvtDiracSpinor::operator= ( const EvtDiracSpinor &  dspinor)
inline

Definition at line 127 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.hh.

127 {
128
129 spinor[0]=dspinor.spinor[0];
130 spinor[1]=dspinor.spinor[1];
131 spinor[2]=dspinor.spinor[2];
132 spinor[3]=dspinor.spinor[3];
133
134 return *this;
135
136}

◆ operator=() [2/2]

EvtDiracSpinor & EvtDiracSpinor::operator= ( const EvtDiracSpinor &  dspinor)
inline

◆ set() [1/2]

void EvtDiracSpinor::set ( const EvtComplex &  sp0,
const EvtComplex &  sp1,
const EvtComplex &  sp2,
const EvtComplex &  sp3 
)

Definition at line 40 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

41 {
42
43 spinor[0]=sp0;spinor[1]=sp1;spinor[2]=sp2;spinor[3]=sp3;
44}

Referenced by EvtDiracSpinor(), EvtRaritaSchwinger::getSpinor(), EvtRaritaSchwingerParticle::init(), EvtDiracParticle::init(), EvtNeutrinoParticle::init(), EvtDiracParticle::rotateToHelicityBasis(), and EvtRaritaSchwingerParticle::rotateToHelicityBasis().

◆ set() [2/2]

void EvtDiracSpinor::set ( const EvtComplex &  sp0,
const EvtComplex &  sp1,
const EvtComplex &  sp2,
const EvtComplex &  sp3 
)

◆ set_spinor() [1/2]

void EvtDiracSpinor::set_spinor ( int  i,
const EvtComplex &  sp 
)

Definition at line 46 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

46 {
47
48 spinor[i]=sp;
49}

Referenced by conj(), and operator*().

◆ set_spinor() [2/2]

void EvtDiracSpinor::set_spinor ( int  i,
const EvtComplex &  sp 
)

Friends And Related Function Documentation

◆ boostTo [1/4]

EvtDiracSpinor boostTo ( const EvtDiracSpinor &  sp,
const EvtVector3R  boost 
)
friend

Definition at line 84 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

85 {
86
87 EvtDiracSpinor tmp(sp);
88 tmp.applyBoostTo(boost);
89 return tmp;
90
91}

◆ boostTo [2/4]

EvtDiracSpinor boostTo ( const EvtDiracSpinor &  sp,
const EvtVector3R  boost 
)
friend

Definition at line 84 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

85 {
86
87 EvtDiracSpinor tmp(sp);
88 tmp.applyBoostTo(boost);
89 return tmp;
90
91}

◆ boostTo [3/4]

EvtDiracSpinor boostTo ( const EvtDiracSpinor &  sp,
const EvtVector4R  p4 
)
friend

Definition at line 75 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

76 {
77
78 EvtDiracSpinor tmp(sp);
79 tmp.applyBoostTo(p4);
80 return tmp;
81
82}

◆ boostTo [4/4]

EvtDiracSpinor boostTo ( const EvtDiracSpinor &  sp,
const EvtVector4R  p4 
)
friend

Definition at line 75 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

76 {
77
78 EvtDiracSpinor tmp(sp);
79 tmp.applyBoostTo(p4);
80 return tmp;
81
82}

◆ EvtLeptonACurrent [1/2]

EvtVector4C EvtLeptonACurrent ( const EvtDiracSpinor &  d,
const EvtDiracSpinor &  dp 
)
friend

Definition at line 236 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

236 {
237
238 EvtVector4C temp;
239
240 EvtGammaMatrix mat;
241
242 // no conjugate here; done in the multiplication
243 // yes this is stupid and fooled me to for a long time (ryd)
244
245 mat = EvtGammaMatrix::v0()-EvtGammaMatrix::va0();
246 temp.set(0,d*(mat*dp));
247
248 mat = EvtGammaMatrix::v1()-EvtGammaMatrix::va1();
249 temp.set(1,d*(mat*dp));
250
251 mat = EvtGammaMatrix::v2()-EvtGammaMatrix::va2();
252 temp.set(2,d*(mat*dp));
253
254 mat = EvtGammaMatrix::v3()-EvtGammaMatrix::va3();
255 temp.set(3,d*(mat*dp));
256
257 return temp;
258}

◆ EvtLeptonACurrent [2/2]

EvtVector4C EvtLeptonACurrent ( const EvtDiracSpinor &  d,
const EvtDiracSpinor &  dp 
)
friend

Definition at line 236 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

236 {
237
238 EvtVector4C temp;
239
240 EvtGammaMatrix mat;
241
242 // no conjugate here; done in the multiplication
243 // yes this is stupid and fooled me to for a long time (ryd)
244
245 mat = EvtGammaMatrix::v0()-EvtGammaMatrix::va0();
246 temp.set(0,d*(mat*dp));
247
248 mat = EvtGammaMatrix::v1()-EvtGammaMatrix::va1();
249 temp.set(1,d*(mat*dp));
250
251 mat = EvtGammaMatrix::v2()-EvtGammaMatrix::va2();
252 temp.set(2,d*(mat*dp));
253
254 mat = EvtGammaMatrix::v3()-EvtGammaMatrix::va3();
255 temp.set(3,d*(mat*dp));
256
257 return temp;
258}

◆ EvtLeptonPCurrent [1/2]

EvtComplex EvtLeptonPCurrent ( const EvtDiracSpinor &  d,
const EvtDiracSpinor &  dp 
)
friend

Definition at line 272 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

272 {
273
274 EvtComplex temp;
275
276 // no conjugate here; done in the multiplication
277 // yes this is stupid and fooled me to for a long time (ryd)
278 static EvtGammaMatrix m=EvtGammaMatrix::g0()*EvtGammaMatrix::g5();
279 temp=d*(m*dp);
280
281 return temp;
282}

◆ EvtLeptonPCurrent [2/2]

EvtComplex EvtLeptonPCurrent ( const EvtDiracSpinor &  d,
const EvtDiracSpinor &  dp 
)
friend

Definition at line 272 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

272 {
273
274 EvtComplex temp;
275
276 // no conjugate here; done in the multiplication
277 // yes this is stupid and fooled me to for a long time (ryd)
278 static EvtGammaMatrix m=EvtGammaMatrix::g0()*EvtGammaMatrix::g5();
279 temp=d*(m*dp);
280
281 return temp;
282}

◆ EvtLeptonSCurrent [1/2]

EvtComplex EvtLeptonSCurrent ( const EvtDiracSpinor &  d,
const EvtDiracSpinor &  dp 
)
friend

Definition at line 260 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

260 {
261
262 EvtComplex temp;
263
264 // no conjugate here; done in the multiplication
265 // yes this is stupid and fooled me to for a long time (ryd)
266
267 temp=d*(EvtGammaMatrix::g0()*dp);
268
269 return temp;
270}

◆ EvtLeptonSCurrent [2/2]

EvtComplex EvtLeptonSCurrent ( const EvtDiracSpinor &  d,
const EvtDiracSpinor &  dp 
)
friend

Definition at line 260 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

260 {
261
262 EvtComplex temp;
263
264 // no conjugate here; done in the multiplication
265 // yes this is stupid and fooled me to for a long time (ryd)
266
267 temp=d*(EvtGammaMatrix::g0()*dp);
268
269 return temp;
270}

◆ EvtLeptonTCurrent [1/2]

EvtTensor4C EvtLeptonTCurrent ( const EvtDiracSpinor &  d,
const EvtDiracSpinor &  dp 
)
friend

Definition at line 284 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

284 {
285
286 EvtTensor4C temp;
287 temp.zero();
288 EvtComplex i2(0,0.5);
289
290 static EvtGammaMatrix mat01=EvtGammaMatrix::g0()*
291 (EvtGammaMatrix::g0()*EvtGammaMatrix::g1()-
292 EvtGammaMatrix::g1()*EvtGammaMatrix::g0());
293 static EvtGammaMatrix mat02=EvtGammaMatrix::g0()*
294 (EvtGammaMatrix::g0()*EvtGammaMatrix::g2()-
295 EvtGammaMatrix::g2()*EvtGammaMatrix::g0());
296 static EvtGammaMatrix mat03=EvtGammaMatrix::g0()*
297 (EvtGammaMatrix::g0()*EvtGammaMatrix::g3()-
298 EvtGammaMatrix::g3()*EvtGammaMatrix::g0());
299 static EvtGammaMatrix mat12=EvtGammaMatrix::g0()*
300 (EvtGammaMatrix::g1()*EvtGammaMatrix::g2()-
301 EvtGammaMatrix::g2()*EvtGammaMatrix::g1());
302 static EvtGammaMatrix mat13=EvtGammaMatrix::g0()*
303 (EvtGammaMatrix::g1()*EvtGammaMatrix::g3()-
304 EvtGammaMatrix::g3()*EvtGammaMatrix::g1());
305 static EvtGammaMatrix mat23=EvtGammaMatrix::g0()*
306 (EvtGammaMatrix::g2()*EvtGammaMatrix::g3()-
307 EvtGammaMatrix::g3()*EvtGammaMatrix::g2());
308
309
310 temp.set(0,1,i2*(d*(mat01*dp)));
311 temp.set(1,0,-temp.get(0,1));
312
313 temp.set(0,2,i2*(d*(mat02*dp)));
314 temp.set(2,0,-temp.get(0,2));
315
316 temp.set(0,3,i2*(d*(mat03*dp)));
317 temp.set(3,0,-temp.get(0,3));
318
319 temp.set(1,2,i2*(d*(mat12*dp)));
320 temp.set(2,1,-temp.get(1,2));
321
322 temp.set(1,3,i2*(d*(mat13*dp)));
323 temp.set(3,1,-temp.get(1,3));
324
325 temp.set(2,3,i2*(d*(mat23*dp)));
326 temp.set(3,2,-temp.get(2,3));
327
328 return temp;
329}
EvtTensor4C::set
void set(int i, int j, const EvtComplex &c)
Definition: bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtTensor4C.hh:122
EvtTensor4C::get
const EvtComplex & get(int i, int j) const
Definition: bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtTensor4C.hh:126

◆ EvtLeptonTCurrent [2/2]

EvtTensor4C EvtLeptonTCurrent ( const EvtDiracSpinor &  d,
const EvtDiracSpinor &  dp 
)
friend

Definition at line 284 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

284 {
285
286 EvtTensor4C temp;
287 temp.zero();
288 EvtComplex i2(0,0.5);
289
290 static EvtGammaMatrix mat01=EvtGammaMatrix::g0()*
291 (EvtGammaMatrix::g0()*EvtGammaMatrix::g1()-
292 EvtGammaMatrix::g1()*EvtGammaMatrix::g0());
293 static EvtGammaMatrix mat02=EvtGammaMatrix::g0()*
294 (EvtGammaMatrix::g0()*EvtGammaMatrix::g2()-
295 EvtGammaMatrix::g2()*EvtGammaMatrix::g0());
296 static EvtGammaMatrix mat03=EvtGammaMatrix::g0()*
297 (EvtGammaMatrix::g0()*EvtGammaMatrix::g3()-
298 EvtGammaMatrix::g3()*EvtGammaMatrix::g0());
299 static EvtGammaMatrix mat12=EvtGammaMatrix::g0()*
300 (EvtGammaMatrix::g1()*EvtGammaMatrix::g2()-
301 EvtGammaMatrix::g2()*EvtGammaMatrix::g1());
302 static EvtGammaMatrix mat13=EvtGammaMatrix::g0()*
303 (EvtGammaMatrix::g1()*EvtGammaMatrix::g3()-
304 EvtGammaMatrix::g3()*EvtGammaMatrix::g1());
305 static EvtGammaMatrix mat23=EvtGammaMatrix::g0()*
306 (EvtGammaMatrix::g2()*EvtGammaMatrix::g3()-
307 EvtGammaMatrix::g3()*EvtGammaMatrix::g2());
308
309
310 temp.set(0,1,i2*(d*(mat01*dp)));
311 temp.set(1,0,-temp.get(0,1));
312
313 temp.set(0,2,i2*(d*(mat02*dp)));
314 temp.set(2,0,-temp.get(0,2));
315
316 temp.set(0,3,i2*(d*(mat03*dp)));
317 temp.set(3,0,-temp.get(0,3));
318
319 temp.set(1,2,i2*(d*(mat12*dp)));
320 temp.set(2,1,-temp.get(1,2));
321
322 temp.set(1,3,i2*(d*(mat13*dp)));
323 temp.set(3,1,-temp.get(1,3));
324
325 temp.set(2,3,i2*(d*(mat23*dp)));
326 temp.set(3,2,-temp.get(2,3));
327
328 return temp;
329}

◆ EvtLeptonTg5Current [1/2]

EvtTensor4C EvtLeptonTg5Current ( const EvtDiracSpinor &  d,
const EvtDiracSpinor &  dp 
)
friend

Definition at line 331 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

331 {//sigma^{mu,nu} gamma_5
332
333 EvtTensor4C temp;
334 temp.zero();
335 EvtComplex i2(0,0.5);
336
337 static EvtGammaMatrix mat01=EvtGammaMatrix::g0()*
338 (EvtGammaMatrix::g0()*EvtGammaMatrix::g1()-
339 EvtGammaMatrix::g1()*EvtGammaMatrix::g0())*EvtGammaMatrix::g5();
340 static EvtGammaMatrix mat02=EvtGammaMatrix::g0()*
341 (EvtGammaMatrix::g0()*EvtGammaMatrix::g2()-
342 EvtGammaMatrix::g2()*EvtGammaMatrix::g0())*EvtGammaMatrix::g5();
343 static EvtGammaMatrix mat03=EvtGammaMatrix::g0()*
344 (EvtGammaMatrix::g0()*EvtGammaMatrix::g3()-
345 EvtGammaMatrix::g3()*EvtGammaMatrix::g0())*EvtGammaMatrix::g5();
346 static EvtGammaMatrix mat12=EvtGammaMatrix::g0()*
347 (EvtGammaMatrix::g1()*EvtGammaMatrix::g2()-
348 EvtGammaMatrix::g2()*EvtGammaMatrix::g1())*EvtGammaMatrix::g5();
349 static EvtGammaMatrix mat13=EvtGammaMatrix::g0()*
350 (EvtGammaMatrix::g1()*EvtGammaMatrix::g3()-
351 EvtGammaMatrix::g3()*EvtGammaMatrix::g1())*EvtGammaMatrix::g5();
352 static EvtGammaMatrix mat23=EvtGammaMatrix::g0()*
353 (EvtGammaMatrix::g2()*EvtGammaMatrix::g3()-
354 EvtGammaMatrix::g3()*EvtGammaMatrix::g2())*EvtGammaMatrix::g5();
355
356
357 temp.set(0,1,i2*(d*(mat01*dp)));
358 temp.set(1,0,-temp.get(0,1));
359
360 temp.set(0,2,i2*(d*(mat02*dp)));
361 temp.set(2,0,-temp.get(0,2));
362
363 temp.set(0,3,i2*(d*(mat03*dp)));
364 temp.set(3,0,-temp.get(0,3));
365
366 temp.set(1,2,i2*(d*(mat12*dp)));
367 temp.set(2,1,-temp.get(1,2));
368
369 temp.set(1,3,i2*(d*(mat13*dp)));
370 temp.set(3,1,-temp.get(1,3));
371
372 temp.set(2,3,i2*(d*(mat23*dp)));
373 temp.set(3,2,-temp.get(2,3));
374
375 return temp;
376}

◆ EvtLeptonTg5Current [2/2]

EvtTensor4C EvtLeptonTg5Current ( const EvtDiracSpinor &  d,
const EvtDiracSpinor &  dp 
)
friend

Definition at line 331 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

331 {//sigma^{mu,nu} gamma_5
332
333 EvtTensor4C temp;
334 temp.zero();
335 EvtComplex i2(0,0.5);
336
337 static EvtGammaMatrix mat01=EvtGammaMatrix::g0()*
338 (EvtGammaMatrix::g0()*EvtGammaMatrix::g1()-
339 EvtGammaMatrix::g1()*EvtGammaMatrix::g0())*EvtGammaMatrix::g5();
340 static EvtGammaMatrix mat02=EvtGammaMatrix::g0()*
341 (EvtGammaMatrix::g0()*EvtGammaMatrix::g2()-
342 EvtGammaMatrix::g2()*EvtGammaMatrix::g0())*EvtGammaMatrix::g5();
343 static EvtGammaMatrix mat03=EvtGammaMatrix::g0()*
344 (EvtGammaMatrix::g0()*EvtGammaMatrix::g3()-
345 EvtGammaMatrix::g3()*EvtGammaMatrix::g0())*EvtGammaMatrix::g5();
346 static EvtGammaMatrix mat12=EvtGammaMatrix::g0()*
347 (EvtGammaMatrix::g1()*EvtGammaMatrix::g2()-
348 EvtGammaMatrix::g2()*EvtGammaMatrix::g1())*EvtGammaMatrix::g5();
349 static EvtGammaMatrix mat13=EvtGammaMatrix::g0()*
350 (EvtGammaMatrix::g1()*EvtGammaMatrix::g3()-
351 EvtGammaMatrix::g3()*EvtGammaMatrix::g1())*EvtGammaMatrix::g5();
352 static EvtGammaMatrix mat23=EvtGammaMatrix::g0()*
353 (EvtGammaMatrix::g2()*EvtGammaMatrix::g3()-
354 EvtGammaMatrix::g3()*EvtGammaMatrix::g2())*EvtGammaMatrix::g5();
355
356
357 temp.set(0,1,i2*(d*(mat01*dp)));
358 temp.set(1,0,-temp.get(0,1));
359
360 temp.set(0,2,i2*(d*(mat02*dp)));
361 temp.set(2,0,-temp.get(0,2));
362
363 temp.set(0,3,i2*(d*(mat03*dp)));
364 temp.set(3,0,-temp.get(0,3));
365
366 temp.set(1,2,i2*(d*(mat12*dp)));
367 temp.set(2,1,-temp.get(1,2));
368
369 temp.set(1,3,i2*(d*(mat13*dp)));
370 temp.set(3,1,-temp.get(1,3));
371
372 temp.set(2,3,i2*(d*(mat23*dp)));
373 temp.set(3,2,-temp.get(2,3));
374
375 return temp;
376}

◆ EvtLeptonVACurrent [1/2]

EvtVector4C EvtLeptonVACurrent ( const EvtDiracSpinor &  d,
const EvtDiracSpinor &  dp 
)
friend

Definition at line 187 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

187 {
188
189 //Old code; below is a new specialized code that does it more efficiently.
190 //EvtGammaMatrix mat;
191 //EvtVector4C temp;
192 //mat.va0();
193 //temp.set(0,d*(mat*dp));
194 //mat.va1();
195 //temp.set(1,d*(mat*dp));
196 //mat.va2();
197 //temp.set(2,d*(mat*dp));
198 //mat.va3();
199 //temp.set(3,d*(mat*dp));
200 //return temp;
201
202
203 EvtComplex u02=::conj(d.spinor[0]-d.spinor[2]);
204 EvtComplex u13=::conj(d.spinor[1]-d.spinor[3]);
205
206 EvtComplex v02=dp.spinor[0]-dp.spinor[2];
207 EvtComplex v13=dp.spinor[1]-dp.spinor[3];
208
209 EvtComplex a=u02*v02;
210 EvtComplex b=u13*v13;
211
212 EvtComplex c=u02*v13;
213 EvtComplex e=u13*v02;
214
215 return EvtVector4C(a+b,-(c+e),EvtComplex(0,1)*(c-e),b-a);
216
217
218}
b
const double b
Definition: slope.cxx:9

◆ EvtLeptonVACurrent [2/2]

EvtVector4C EvtLeptonVACurrent ( const EvtDiracSpinor &  d,
const EvtDiracSpinor &  dp 
)
friend

Definition at line 187 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

187 {
188
189 //Old code; below is a new specialized code that does it more efficiently.
190 //EvtGammaMatrix mat;
191 //EvtVector4C temp;
192 //mat.va0();
193 //temp.set(0,d*(mat*dp));
194 //mat.va1();
195 //temp.set(1,d*(mat*dp));
196 //mat.va2();
197 //temp.set(2,d*(mat*dp));
198 //mat.va3();
199 //temp.set(3,d*(mat*dp));
200 //return temp;
201
202
203 EvtComplex u02=::conj(d.spinor[0]-d.spinor[2]);
204 EvtComplex u13=::conj(d.spinor[1]-d.spinor[3]);
205
206 EvtComplex v02=dp.spinor[0]-dp.spinor[2];
207 EvtComplex v13=dp.spinor[1]-dp.spinor[3];
208
209 EvtComplex a=u02*v02;
210 EvtComplex b=u13*v13;
211
212 EvtComplex c=u02*v13;
213 EvtComplex e=u13*v02;
214
215 return EvtVector4C(a+b,-(c+e),EvtComplex(0,1)*(c-e),b-a);
216
217
218}

◆ EvtLeptonVCurrent [1/2]

EvtVector4C EvtLeptonVCurrent ( const EvtDiracSpinor &  d,
const EvtDiracSpinor &  dp 
)
friend

Definition at line 220 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

220 {
221
222 EvtVector4C temp;
223
224 // no conjugate here; done in the multiplication
225 // yes this is stupid and fooled me to for a long time (ryd)
226
227 temp.set(0,d*(EvtGammaMatrix::v0()*dp));
228 temp.set(1,d*(EvtGammaMatrix::v1()*dp));
229 temp.set(2,d*(EvtGammaMatrix::v2()*dp));
230 temp.set(3,d*(EvtGammaMatrix::v3()*dp));
231
232 return temp;
233}

◆ EvtLeptonVCurrent [2/2]

EvtVector4C EvtLeptonVCurrent ( const EvtDiracSpinor &  d,
const EvtDiracSpinor &  dp 
)
friend

Definition at line 220 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

220 {
221
222 EvtVector4C temp;
223
224 // no conjugate here; done in the multiplication
225 // yes this is stupid and fooled me to for a long time (ryd)
226
227 temp.set(0,d*(EvtGammaMatrix::v0()*dp));
228 temp.set(1,d*(EvtGammaMatrix::v1()*dp));
229 temp.set(2,d*(EvtGammaMatrix::v2()*dp));
230 temp.set(3,d*(EvtGammaMatrix::v3()*dp));
231
232 return temp;
233}

◆ operator* [1/2]

EvtDiracSpinor operator* ( const EvtComplex &  c,
const EvtDiracSpinor &  d 
)
friend

Definition at line 378 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

378 {
379 EvtDiracSpinor result;
380 result.spinor[0] = c*d.spinor[0];
381 result.spinor[1] = c*d.spinor[1];
382 result.spinor[2] = c*d.spinor[2];
383 result.spinor[3] = c*d.spinor[3];
384 return result;
385 }

◆ operator* [2/2]

EvtDiracSpinor operator* ( const EvtComplex &  c,
const EvtDiracSpinor &  d 
)
friend

Definition at line 378 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

378 {
379 EvtDiracSpinor result;
380 result.spinor[0] = c*d.spinor[0];
381 result.spinor[1] = c*d.spinor[1];
382 result.spinor[2] = c*d.spinor[2];
383 result.spinor[3] = c*d.spinor[3];
384 return result;
385 }

◆ operator+ [1/2]

EvtDiracSpinor operator+ ( const EvtDiracSpinor &  u1,
const EvtDiracSpinor &  u2 
)
friend

Definition at line 148 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.hh.

149 {
150
151 return EvtDiracSpinor(u1)+=u2;
152
153}

◆ operator+ [2/2]

EvtDiracSpinor operator+ ( const EvtDiracSpinor &  u1,
const EvtDiracSpinor &  u2 
)
friend

Definition at line 148 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.hh.

149 {
150
151 return EvtDiracSpinor(u1)+=u2;
152
153}

◆ operator- [1/2]

EvtDiracSpinor operator- ( const EvtDiracSpinor &  u1,
const EvtDiracSpinor &  u2 
)
friend

Definition at line 165 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.hh.

166 {
167
168 return EvtDiracSpinor(u1)-=u2;
169
170}

◆ operator- [2/2]

EvtDiracSpinor operator- ( const EvtDiracSpinor &  u1,
const EvtDiracSpinor &  u2 
)
friend

Definition at line 165 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.hh.

166 {
167
168 return EvtDiracSpinor(u1)-=u2;
169
170}

◆ operator<< [1/2]

std::ostream & operator<< ( std::ostream &  s,
const EvtDiracSpinor &  c 
)
friend

◆ operator<< [2/2]

std::ostream & operator<< ( std::ostream &  s,
const EvtDiracSpinor &  c 
)
friend

◆ rotateEuler [1/2]

EvtDiracSpinor rotateEuler ( const EvtDiracSpinor &  sp,
double  alpha,
double  beta,
double  gamma 
)
friend

Definition at line 66 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

67 {
68
69 EvtDiracSpinor tmp(sp);
70 tmp.applyRotateEuler(alpha,beta,gamma);
71 return tmp;
72
73}

◆ rotateEuler [2/2]

EvtDiracSpinor rotateEuler ( const EvtDiracSpinor &  sp,
double  alpha,
double  beta,
double  gamma 
)
friend

Definition at line 66 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDiracSpinor.cc.

67 {
68
69 EvtDiracSpinor tmp(sp);
70 tmp.applyRotateEuler(alpha,beta,gamma);
71 return tmp;
72
73}
The documentation for this class was generated from the following files: