#include <EvtPhotonParticle.hh>

Inheritance diagram for EvtPhotonParticle:

Public Member Functions
	EvtPhotonParticle ()

virtual	~EvtPhotonParticle ()

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

void	init (EvtId part_n, const EvtVector4R &p4)

EvtVector4C	epsParentPhoton (int i)

EvtVector4C	epsPhoton (int i)

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

virtual EvtTensor4C	epsTensorParent (int i) const

virtual EvtTensor4C	epsTensor (int i) const

virtual void	init (EvtId part_n, const EvtVector4R &p4)=0

void	addDaug (EvtParticle *node)

void	decay ()

void	deleteTree ()

void	deleteDaughters (bool keepChannel=false)

void	setChannel (int i)

void	setGeneratorFlag (int flag)

int	getGeneratorFlag ()

void	makeDaughters (int ndaug, EvtId *id)

double	initializePhaseSpace (int numdaughter, EvtId *daughters, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)

EvtParticle *	getDaug (int i)

EvtParticle *	nextIter (EvtParticle *rootOfTree=0)

void	makeStdHep (EvtStdHep &stdhep, EvtSecondary &secondary, EvtId *stable_parent_ihep)

void	makeStdHep (EvtStdHep &stdhep)

EvtVector4R	getP4Lab ()

EvtVector4R	getP4Restframe ()

EvtVector4R	get4Pos ()

EvtParticle *	getParent ()

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 EvtVector4R &	getP4 () 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)

virtual EvtSpinDensity	rotateToHelicityBasis () const =0

virtual EvtSpinDensity	rotateToHelicityBasis (double alpha, double beta, double gamma) const =0

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
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 31 of file EvtPhotonParticle.hh.

Constructor & Destructor Documentation

◆ EvtPhotonParticle()

EvtPhotonParticle::EvtPhotonParticle ( )

inline

Definition at line 35 of file EvtPhotonParticle.hh.

35{}

◆ ~EvtPhotonParticle()

EvtPhotonParticle::~EvtPhotonParticle ( )

virtual

Definition at line 31 of file EvtPhotonParticle.cc.

31{}

Member Function Documentation

◆ epsParentPhoton()

EvtVector4C EvtPhotonParticle::epsParentPhoton ( int i )

virtual

Returns polarization vector in the parents restframe for a photon.

Reimplemented from EvtParticle.

Definition at line 53 of file EvtPhotonParticle.cc.

                                                   {
 
  if (!_evalBasis){
 
    _evalBasis=1;
    eps1.set(EvtComplex(0.0,0.0), EvtComplex(-1.0/sqrt(2.0),0.0),
         EvtComplex(0.0,-1.0/sqrt(2.0)), EvtComplex(0.0,0.0));
    eps2.set(EvtComplex(0.0,0.0), EvtComplex(1.0/sqrt(2.0),0.0),
         EvtComplex(0.0,-1.0/sqrt(2.0)), EvtComplex(0.0,0.0));
    
    // These are for photon along z axis.  Rotate to get 
    // correct direction...
    
    double phi,theta;
 
    EvtVector4R p=this->getP4();
 
    double px=p.get(1);
    double py=p.get(2);
    double pz=p.get(3);
 
    phi = atan2(py,px);
    theta = acos(pz/sqrt(px*px+py*py+pz*pz));
    eps1.applyRotateEuler(phi,theta,-phi);
    eps2.applyRotateEuler(phi,theta,-phi);
    
  }
 
 
  EvtVector4C temp;
  
  switch(i) {
 
  case 0:
    temp=eps1;
    break;
  case 1:
    temp=eps2;
    break;
  default:
    report(ERROR,"EvtGen") << "EvtPhotonParticle.cc: Asked "
               << "for state:"<<i<<endl;
    ::abort();
    break;
  }
 
  return temp;
}

Referenced by EvtLambdaP_BarGamma::decay().

◆ epsPhoton()

EvtVector4C EvtPhotonParticle::epsPhoton ( int i )

virtual

Returns polarization vector in the particles own restframe for a photon.

Reimplemented from EvtParticle.

Definition at line 102 of file EvtPhotonParticle.cc.

                                             {
 
  report(ERROR,"EvtGen") << "EvtPhotonParticle.cc: Can not get "
             << "state in photons restframe."<<endl;;
  ::abort();
  return EvtVector4C();
 
}

◆ init() [1/2]

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

virtual

Initialiaze particle with id and 4momentum.

Implements EvtParticle.

Definition at line 33 of file EvtPhotonParticle.cc.

                                                              {
 
  init(part_n,p4.get(0),p4.get(1),
       p4.get(2),p4.get(3));
}

◆ init() [2/2]

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

Definition at line 39 of file EvtPhotonParticle.cc.

                                                                               {
 
  _validP4=true;
  setp(e,px,py,pz);
  setpart_num(part_n);
 
  setLifetime();
  
  //defere calculation of basis vectors untill they are needed!
  _evalBasis=0;
 
}

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

◆ rotateToHelicityBasis() [1/2]

EvtSpinDensity EvtPhotonParticle::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 112 of file EvtPhotonParticle.cc.

                                                              {
 
  EvtVector4C eplus(0.0,-1.0/sqrt(2.0),
            EvtComplex(0.0,-1.0/sqrt(2.0)),0.0);
  EvtVector4C eminus(0.0,1.0/sqrt(2.0),
             EvtComplex(0.0,-1.0/sqrt(2.0)),0.0);
 
  //Really uggly have to cast away constness because the
  //function epsParentPhoton caches the state vectors...
  EvtVector4C e1=((EvtParticle*)this)->epsParentPhoton(0);
  EvtVector4C e2=((EvtParticle*)this)->epsParentPhoton(1);
 
 
  EvtSpinDensity R;
  R.SetDim(2);
 
  R.Set(0,0,(eplus.conj())*e1);
  R.Set(0,1,(eplus.conj())*e2);
  
  R.Set(1,0,(eminus.conj())*e1);
  R.Set(1,1,(eminus.conj())*e2);
  
  return R;
    
}

◆ rotateToHelicityBasis() [2/2]

EvtSpinDensity EvtPhotonParticle::rotateToHelicityBasis	(	double	alpha,
		double	beta,
		double	gamma
	)		const

virtual

Implements EvtParticle.

Definition at line 139 of file EvtPhotonParticle.cc.

                                               {
 
  EvtVector4C eplus(0.0,-1.0/sqrt(2.0),
            EvtComplex(0.0,-1.0/sqrt(2.0)),0.0);
  EvtVector4C eminus(0.0,1.0/sqrt(2.0),
             EvtComplex(0.0,-1.0/sqrt(2.0)),0.0);
  
  eplus.applyRotateEuler(alpha,beta,gamma);
  eminus.applyRotateEuler(alpha,beta,gamma);
  
  
  //Really uggly have to cast away constness because the
  //function epsParentPhoton caches the state vectors...
  EvtVector4C e1=((EvtParticle*)this)->epsParentPhoton(0);
  EvtVector4C e2=((EvtParticle*)this)->epsParentPhoton(1);
  
  EvtSpinDensity R;
  R.SetDim(2);
  
  R.Set(0,0,(eplus.conj())*e1);
  R.Set(0,1,(eplus.conj())*e2);
  
  R.Set(1,0,(eminus.conj())*e1);
  R.Set(1,1,(eminus.conj())*e2);
  
  return R;
 
}

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

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ EvtPhotonParticle()

◆ ~EvtPhotonParticle()

Member Function Documentation

◆ epsParentPhoton()

◆ epsPhoton()

◆ init() [1/2]

◆ init() [2/2]

◆ rotateToHelicityBasis() [1/2]

◆ rotateToHelicityBasis() [2/2]