EvtDiracSpinor.hh File Reference

#include "EvtGenBase/EvtComplex.hh"
#include "EvtGenBase/EvtVector4R.hh"
#include "EvtGenBase/EvtVector3R.hh"

Go to the source code of this file.

Classes
class	EvtDiracSpinor

Functions
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)
EvtDiracSpinor	operator+ (const EvtDiracSpinor &u1, const EvtDiracSpinor &u2)
EvtDiracSpinor	operator- (const EvtDiracSpinor &u1, const EvtDiracSpinor &u2)
std::ostream &	operator<< (std::ostream &s, const EvtDiracSpinor &c)
EvtDiracSpinor	operator* (const EvtComplex &c, const EvtDiracSpinor &d)

Function Documentation

boostTo() [1/2]

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

Definition at line 57 of file EvtDiracSpinor.cc.

                                       {
 
  EvtDiracSpinor tmp(sp);
  tmp.applyBoostTo(boost);
  return tmp;
 
}

boostTo() [2/2]

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

Definition at line 55 of file EvtDiracSpinor.cc.

                                    {
 
  EvtDiracSpinor tmp(sp);
  tmp.applyBoostTo(p4);
  return tmp;
 
}

Referenced by EvtbTosllAmp::CalcMaxProb(), EvtSemiLeptonicAmp::CalcMaxProb(), EvtBtoXsll::decay(), EvtConExc::decay(), EvtJetSet::decay(), EvtKstarstargamma::decay(), EvtLunda::decay(), EvtLundCharm::decay(), EvtPythia::decay(), EvtTauola::decay(), EvtVectorParticle::epsParent(), EvtAmpA1(), EvtAmpA2(), EvtDecayAngleChi(), EvtParticle::get4Pos(), EvtParticle::getP4Lab(), EvtRaritaSchwingerParticle::init(), and EvtDiracParticle::init().

EvtLeptonACurrent()

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

Definition at line 63 of file EvtDiracSpinor.cc.

                                                                               {
 
  EvtVector4C temp;
 
  EvtGammaMatrix mat;
 
  // no conjugate here; done in the multiplication
  // yes this is stupid and fooled me to for a long time (ryd)
 
  mat = EvtGammaMatrix::v0()-EvtGammaMatrix::va0();
  temp.set(0,d*(mat*dp));
 
  mat = EvtGammaMatrix::v1()-EvtGammaMatrix::va1();
  temp.set(1,d*(mat*dp));
 
  mat = EvtGammaMatrix::v2()-EvtGammaMatrix::va2();
  temp.set(2,d*(mat*dp));
 
  mat = EvtGammaMatrix::v3()-EvtGammaMatrix::va3();
  temp.set(3,d*(mat*dp));
  
  return temp;
}

Referenced by EvtbTosllScalarAmp::CalcAmp(), and EvtbTosllVectorAmp::CalcAmp().

EvtLeptonPCurrent()

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

Definition at line 67 of file EvtDiracSpinor.cc.

                                                                              {
 
  EvtComplex temp;
 
  // no conjugate here; done in the multiplication
  // yes this is stupid and fooled me to for a long time (ryd)
  static EvtGammaMatrix m=EvtGammaMatrix::g0()*EvtGammaMatrix::g5();
  temp=d*(m*dp);
  
  return temp;
}

Referenced by EvtHypNonLepton::calcAmp().

EvtLeptonSCurrent()

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

Definition at line 65 of file EvtDiracSpinor.cc.

                                                                              {
 
  EvtComplex temp;
 
  // no conjugate here; done in the multiplication
  // yes this is stupid and fooled me to for a long time (ryd)
 
  temp=d*(EvtGammaMatrix::g0()*dp);
  
  return temp;
}

Referenced by EvtHypNonLepton::calcAmp().

EvtLeptonTCurrent()

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

Definition at line 69 of file EvtDiracSpinor.cc.

                                                                               {
 
  EvtTensor4C temp;
  temp.zero();
  EvtComplex i2(0,0.5);
 
  static EvtGammaMatrix mat01=EvtGammaMatrix::g0()*
    (EvtGammaMatrix::g0()*EvtGammaMatrix::g1()-
     EvtGammaMatrix::g1()*EvtGammaMatrix::g0());
  static EvtGammaMatrix mat02=EvtGammaMatrix::g0()*
    (EvtGammaMatrix::g0()*EvtGammaMatrix::g2()-
     EvtGammaMatrix::g2()*EvtGammaMatrix::g0());
  static EvtGammaMatrix mat03=EvtGammaMatrix::g0()*
    (EvtGammaMatrix::g0()*EvtGammaMatrix::g3()-
    EvtGammaMatrix::g3()*EvtGammaMatrix::g0());
  static EvtGammaMatrix mat12=EvtGammaMatrix::g0()*
    (EvtGammaMatrix::g1()*EvtGammaMatrix::g2()-
    EvtGammaMatrix::g2()*EvtGammaMatrix::g1());
  static EvtGammaMatrix mat13=EvtGammaMatrix::g0()*
    (EvtGammaMatrix::g1()*EvtGammaMatrix::g3()-
     EvtGammaMatrix::g3()*EvtGammaMatrix::g1());
  static EvtGammaMatrix mat23=EvtGammaMatrix::g0()*
    (EvtGammaMatrix::g2()*EvtGammaMatrix::g3()-
     EvtGammaMatrix::g3()*EvtGammaMatrix::g2());
 
 
  temp.set(0,1,i2*(d*(mat01*dp)));
  temp.set(1,0,-temp.get(0,1));
 
  temp.set(0,2,i2*(d*(mat02*dp)));
  temp.set(2,0,-temp.get(0,2));
 
  temp.set(0,3,i2*(d*(mat03*dp)));
  temp.set(3,0,-temp.get(0,3));
 
  temp.set(1,2,i2*(d*(mat12*dp)));
  temp.set(2,1,-temp.get(1,2));
 
  temp.set(1,3,i2*(d*(mat13*dp)));
  temp.set(3,1,-temp.get(1,3));
 
  temp.set(2,3,i2*(d*(mat23*dp)));
  temp.set(3,2,-temp.get(2,3));
  
  return temp;
}

EvtLeptonVACurrent()

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

Definition at line 59 of file EvtDiracSpinor.cc.

                                                                                {
 
  //Old code; below is a new specialized code that does it more efficiently.
  //EvtGammaMatrix mat;
  //EvtVector4C temp;
  //mat.va0();
  //temp.set(0,d*(mat*dp));
  //mat.va1();
  //temp.set(1,d*(mat*dp));
  //mat.va2();
  //temp.set(2,d*(mat*dp));
  //mat.va3();
  //temp.set(3,d*(mat*dp));
  //return temp;
 
 
  EvtComplex u02=::conj(d.spinor[0]-d.spinor[2]);  
  EvtComplex u13=::conj(d.spinor[1]-d.spinor[3]);  
 
  EvtComplex v02=dp.spinor[0]-dp.spinor[2];
  EvtComplex v13=dp.spinor[1]-dp.spinor[3];
 
  EvtComplex a=u02*v02;
  EvtComplex b=u13*v13;
 
  EvtComplex c=u02*v13;
  EvtComplex e=u13*v02;
 
  return EvtVector4C(a+b,-(c+e),EvtComplex(0,1)*(c-e),b-a);
 
  
}

Referenced by EvtSemiLeptonicBaryonAmp::CalcAmp(), EvtSemiLeptonicTensorAmp::CalcAmp(), EvtSemiLeptonicVectorAmp::CalcAmp(), EvtKstarnunu::decay(), EvtLNuGamma::decay(), EvtSll::decay(), EvtSLN::decay(), EvtTauHadnu::decay(), EvtTaulnunu::decay(), EvtTauScalarnu::decay(), and EvtTauVectornu::decay().

EvtLeptonVCurrent()

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

Definition at line 61 of file EvtDiracSpinor.cc.

                                                                               {
 
  EvtVector4C temp;
 
  // no conjugate here; done in the multiplication
  // yes this is stupid and fooled me to for a long time (ryd)
 
  temp.set(0,d*(EvtGammaMatrix::v0()*dp));
  temp.set(1,d*(EvtGammaMatrix::v1()*dp));
  temp.set(2,d*(EvtGammaMatrix::v2()*dp));
  temp.set(3,d*(EvtGammaMatrix::v3()*dp));
  
  return temp;
}

Referenced by EvtbTosllScalarAmp::CalcAmp(), EvtbTosllVectorAmp::CalcAmp(), and EvtVll::decay().

operator*()

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

Definition at line 75 of file EvtDiracSpinor.cc.

                                                                       {
     EvtDiracSpinor result;
     result.spinor[0] = c*d.spinor[0];
     result.spinor[1] = c*d.spinor[1];
     result.spinor[2] = c*d.spinor[2];
     result.spinor[3] = c*d.spinor[3];
     return result;
 }

operator+()

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

inline

Definition at line 145 of file EvtDiracSpinor.hh.

                                         {
  
  return EvtDiracSpinor(u1)+=u2;
 
}

operator-()

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

inline

Definition at line 162 of file EvtDiracSpinor.hh.

                                         {
  
  return EvtDiracSpinor(u1)-=u2;
 
}

operator<<()

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

rotateEuler()

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

Definition at line 53 of file EvtDiracSpinor.cc.

                                                     {
 
  EvtDiracSpinor tmp(sp);
  tmp.applyRotateEuler(alpha,beta,gamma);
  return tmp;
 
}

Referenced by EvtJTO3P::decay().