9 #ifndef ThePEG_LorentzRSSpinor_H
10 #define ThePEG_LorentzRSSpinor_H
12 #include "ThePEG/Config/ThePEG.h"
13 #include "ThePEG/Vectors/ThreeVector.h"
15 #include "LorentzRSSpinor.fh"
16 #include "LorentzRSSpinorBar.h"
17 #include "LorentzSpinorBar.h"
18 #include "LorentzSpinor.h"
19 #include "LorentzPolarizationVector.h"
72 template<
typename Value>
83 for(
unsigned int ix=0;ix<4;++ix)
84 for(
unsigned int iy=0;iy<4;++iy)
85 _spin[ix][iy]=Value();
93 complex<Value> c1, complex<Value> d1,
94 complex<Value> a2, complex<Value> b2,
95 complex<Value> c2, complex<Value> d2,
96 complex<Value> a3, complex<Value> b3,
97 complex<Value> c3, complex<Value> d3,
98 complex<Value> a4, complex<Value> b4,
99 complex<Value> c4, complex<Value> d4,
115 assert( i >= 0 && i <= 3 && j>=0 && j<=3);
123 assert( i >= 0 && i <= 3 && j>=0 && j<=3);
297 for(ix=0;ix<4;++ix) {
298 temp =
_spin[3][ix]*vec.t();
299 temp -=
_spin[0][ix]*vec.x();
300 temp -=
_spin[1][ix]*vec.y();
301 temp -=
_spin[2][ix]*vec.z();
315 for(ix=0;ix<4;++ix) {
316 temp = - (
_spin[0][ix]*vec.x() +
_spin[1][ix]*vec.y()+
317 _spin[2][ix]*vec.z() ) +
_spin[3][ix]*vec.t();
362 template <
typename ValueB>
363 complex<typename BinaryOpTraits<Value,ValueB>::MulT>
365 complex<typename BinaryOpTraits<Value,ValueB>::MulT> output;
368 left*(fbar(3,0)*
_spin[3][0]+fbar(3,1)*
_spin[3][1])
369 +right*(fbar(3,2)*
_spin[3][2]+fbar(3,3)*
_spin[3][3]);
370 for(iz=0;iz<3;++iz) {
372 left*(fbar(iz,0)*
_spin[iz][0]+fbar(iz,1)*
_spin[iz][1])
373 +right*(fbar(iz,2)*
_spin[iz][2]+fbar(iz,3)*
_spin[iz][3]);
384 template <
typename ValueB>
387 typedef complex<typename BinaryOpTraits<Value,ValueB>::MulT> ResultT;
389 for(
size_t iz=0;iz<4;++iz)
410 #ifndef ThePEG_TEMPLATES_IN_CC_FILE
411 #include "LorentzRSSpinor.tcc"