#include <PhotosHEPEVTParticle.h>

Inheritance diagram for Photospp::PhotosHEPEVTParticle:

Public Member Functions
	~PhotosHEPEVTParticle ()

	PhotosHEPEVTParticle (int pdgid, int status, double px, double py, double pz, double e, double m, int ms, int me, int ds, int de)

void	addDaughter (PhotosParticle *daughter)

void	setMothers (std::vector< PhotosParticle * > mothers)

void	setDaughters (std::vector< PhotosParticle * > daughters)

std::vector< PhotosParticle * >	getMothers ()

std::vector< PhotosParticle * >	getDaughters ()

std::vector< PhotosParticle * >	getAllDecayProducts ()

bool	checkMomentumConservation ()

PhotosHEPEVTParticle *	createNewParticle (int pdg_id, int status, double mass, double px, double py, double pz, double e)

void	createHistoryEntry ()

void	createSelfDecayVertex (PhotosParticle *out)

bool	isDaughterOf (PhotosHEPEVTParticle *p)

bool	isMotherOf (PhotosHEPEVTParticle *p)

void	print ()

void	setPdgID (int pdg_id)

void	setStatus (int statu)

void	setMass (double mass)

int	getPdgID ()

int	getStatus ()

double	getMass ()

double	getPx ()

double	getPy ()

double	getPz ()

double	getE ()

void	setPx (double px)

void	setPy (double py)

void	setPz (double pz)

void	setE (double e)

int	getBarcode ()

void	setBarcode (int barcode)

void	setEvent (PhotosHEPEVTEvent *event)

int	getFirstMotherIndex ()

int	getSecondMotherIndex ()

int	getDaughterRangeStart ()

int	getDaughterRangeEnd ()

	~PhotosHEPEVTParticle ()

	PhotosHEPEVTParticle (int pdgid, int status, double px, double py, double pz, double e, double m, int ms, int me, int ds, int de)

void	addDaughter (PhotosParticle *daughter)

void	setMothers (std::vector< PhotosParticle * > mothers)

void	setDaughters (std::vector< PhotosParticle * > daughters)

std::vector< PhotosParticle * >	getMothers ()

std::vector< PhotosParticle * >	getDaughters ()

std::vector< PhotosParticle * >	getAllDecayProducts ()

bool	checkMomentumConservation ()

PhotosHEPEVTParticle *	createNewParticle (int pdg_id, int status, double mass, double px, double py, double pz, double e)

void	createHistoryEntry ()

void	createSelfDecayVertex (PhotosParticle *out)

bool	isDaughterOf (PhotosHEPEVTParticle *p)

bool	isMotherOf (PhotosHEPEVTParticle *p)

void	print ()

void	setPdgID (int pdg_id)

void	setStatus (int statu)

void	setMass (double mass)

int	getPdgID ()

int	getStatus ()

double	getMass ()

double	getPx ()

double	getPy ()

double	getPz ()

double	getE ()

void	setPx (double px)

void	setPy (double py)

void	setPz (double pz)

void	setE (double e)

int	getBarcode ()

void	setBarcode (int barcode)

void	setEvent (PhotosHEPEVTEvent *event)

int	getFirstMotherIndex ()

int	getSecondMotherIndex ()

int	getDaughterRangeStart ()

int	getDaughterRangeEnd ()

bool	hasDaughters ()

bool	hasDaughters ()

PhotosParticle *	findLastSelf ()

PhotosParticle *	findLastSelf ()

std::vector< PhotosParticle * >	findProductionMothers ()

std::vector< PhotosParticle * >	findProductionMothers ()

std::vector< PhotosParticle * >	getDecayTree ()

std::vector< PhotosParticle * >	getDecayTree ()

void	boostToRestFrame (PhotosParticle *boost)

void	boostToRestFrame (PhotosParticle *boost)

void	boostDaughtersToRestFrame (PhotosParticle *boost)

void	boostDaughtersToRestFrame (PhotosParticle *boost)

void	boostFromRestFrame (PhotosParticle *boost)

void	boostFromRestFrame (PhotosParticle *boost)

void	boostDaughtersFromRestFrame (PhotosParticle *boost)

void	boostDaughtersFromRestFrame (PhotosParticle *boost)

void	boostAlongZ (double pz, double e)

void	boostAlongZ (double pz, double e)

void	rotate (int axis, double phi, int second_axis=Z_AXIS)

void	rotate (int axis, double phi, int second_axis=Z_AXIS)

void	rotateDaughters (int axis, double phi, int second_axis=Z_AXIS)

void	rotateDaughters (int axis, double phi, int second_axis=Z_AXIS)

double	getRotationAngle (int axis, int second_axis=Z_AXIS)

double	getRotationAngle (int axis, int second_axis=Z_AXIS)

double	getP ()

double	getP (int axis)

double	getP ()

double	getP (int axis)

void	setP (int axis, double p_component)

void	setP (int axis, double p_component)

virtual double	getVirtuality ()

virtual double	getVirtuality ()

Static Public Attributes
static const int	STABLE =1

static const int	DECAYED =2

static const int	HISTORY =3

static const int	X_AXIS =1

static const int	Y_AXIS =2

static const int	Z_AXIS =3

static const int	Z0 = 23

static const int	HIGGS = 25

static const int	HIGGS_H = 35

static const int	HIGGS_A = 36

static const int	HIGGS_PLUS = 37

static const int	HIGGS_MINUS = -37

static const int	W_PLUS = 24

static const int	W_MINUS = -24

static const int	GAMMA = 22

static const int	TAU_PLUS = -15

static const int	TAU_MINUS = 15

static const int	TAU_NEUTRINO = 16

static const int	TAU_ANTINEUTRINO = -16

static const int	MUON_PLUS = -13

static const int	MUON_MINUS = 13

static const int	MUON_NEUTRINO = 14

static const int	MUON_ANTINEUTRINO = -14

static const int	POSITRON = -11

static const int	ELECTRON = 11

static const int	ELECTRON_NEUTRINO = 12

static const int	ELECTRON_ANTINEUTRINO = -12

static const int	UP = 2

static const int	ANTIUP = -2

static const int	DOWN = 1

static const int	ANTIDOWN = -1

static const int	OTHER = 0

Private Member Functions
void	setDaughterRangeStart (int i)

void	setDaughterRangeEnd (int i)

void	setDaughterRangeStart (int i)

void	setDaughterRangeEnd (int i)

Private Attributes
PhotosHEPEVTEvent *	m_event

int	m_barcode

int	m_first_mother

int	m_second_mother

int	m_daughter_start

int	m_daughter_end

int	m_pdgid

int	m_status

double	m_px

double	m_py

double	m_pz

double	m_e

double	m_generated_mass

vector< PhotosHEPEVTParticle * >	cache

Detailed Description

Definition at line 32 of file PhotosHEPEVTParticle.h.

Constructor & Destructor Documentation

PhotosHEPEVTParticle::~PhotosHEPEVTParticle ( )

Default destructor

Definition at line 7 of file PhotosHEPEVTParticle.cxx.

References cache, i, and m_barcode.

 {
   // Cleanup particles that do not belong to event
   for(unsigned int i=0;i<cache.size();i++)
     if(cache[i]->m_barcode<0)
       delete cache[i];
 }

PhotosHEPEVTParticle::PhotosHEPEVTParticle	(	int	pdgid,
		int	status,
		double	px,
		double	py,
		double	pz,
		double	e,
		double	m,
		int	ms,
		int	me,
		int	ds,
		int	de
	)

Default constructor

Definition at line 15 of file PhotosHEPEVTParticle.cxx.

References de, m, m_barcode, m_daughter_end, m_daughter_start, m_e, m_event, m_first_mother, m_generated_mass, m_pdgid, m_px, m_py, m_pz, m_second_mother, m_status, me, pz, and status.

Referenced by createNewParticle().

                                                                                                                                                     {
   m_px = px;
   m_py = py;
   m_pz = pz;
   m_e  = e;
   m_generated_mass = m;
 
   m_pdgid  = pdgid;
   m_status = status;
 
   m_first_mother   = ms;
   m_second_mother  = me;
   m_daughter_start = ds;
   m_daughter_end   = de;
 
   m_barcode = -1;
   m_event = NULL;
 }

Photospp::PhotosHEPEVTParticle::~PhotosHEPEVTParticle ( )

Default destructor

Photospp::PhotosHEPEVTParticle::PhotosHEPEVTParticle	(	int	pdgid,
		int	status,
		double	px,
		double	py,
		double	pz,
		double	e,
		double	m,
		int	ms,
		int	me,
		int	ds,
		int	de
	)

Default constructor

Member Function Documentation

void PhotosHEPEVTParticle::addDaughter ( PhotosParticle * daughter )

virtual

Add a new daughter to this particle

Implements Photospp::PhotosParticle.

Definition at line 35 of file PhotosHEPEVTParticle.cxx.

References Photospp::Log::Fatal(), Photospp::PhotosParticle::getBarcode(), getDaughterRangeEnd(), getDaughterRangeStart(), Photospp::PhotosParticle::getMothers(), Photospp::PhotosHEPEVTEvent::getParticle(), i, m, m_daughter_end, m_daughter_start, m_event, move, setBarcode(), setDaughterRangeEnd(), setDaughterRangeStart(), Photospp::PhotosParticle::setMothers(), and Photospp::PhotosHEPEVTEvent::setParticle().

 {
   if(!m_event) Log::Fatal("PhotosHEPEVTParticle::addDaughter - particle not in event record");
 
   std::vector<PhotosParticle*> mothers = daughter->getMothers();
   
   mothers.push_back( (PhotosParticle*)this );
 
   daughter->setMothers(mothers);
   
   int bc = daughter->getBarcode();
 
   if(m_daughter_end < 0)
   {
     m_daughter_start = bc;
     m_daughter_end   = bc;
   }
   // if it's in the middle of the event record
   else if(m_daughter_end != bc-1)
   {
     PhotosHEPEVTParticle *newPart = m_event->getParticle(bc);
     
     // Move all particles one spot down the list, to make place for new particle
     for(int i=bc-1;i>m_daughter_end;i--)
     {
       PhotosHEPEVTParticle *move = m_event->getParticle(i);
       move->setBarcode(i+1);
       m_event->setParticle(i+1,move);
     }
 
     m_daughter_end++;
     newPart->setBarcode(m_daughter_end);
     m_event->setParticle(m_daughter_end,newPart);
     
     // Now: correct all pointers before new particle
     for(int i=0;i<m_daughter_end;i++)
     {
       PhotosHEPEVTParticle *check = m_event->getParticle(i);
       int m = check->getDaughterRangeEnd();
       if(m!=-1 && m>m_daughter_end)
       {
         check->setDaughterRangeEnd(m+1);
         check->setDaughterRangeStart(check->getDaughterRangeStart()+1);
       }
     }
   }
   else m_daughter_end = bc;
 }

void Photospp::PhotosHEPEVTParticle::addDaughter ( PhotosParticle * daughter )

virtual

Add a new daughter to this particle

Implements Photospp::PhotosParticle.

void PhotosParticle::boostAlongZ	(	double	boost_pz,
		double	boost_e
	)

inherited

Do a Lorenz transformation along the Z axis.

Boost this vector along the Z direction. Assume no momentum components in the X or Y directions.

Definition at line 155 of file PhotosParticle.cxx.

References Photospp::PhotosParticle::getE(), Photospp::PhotosParticle::getPz(), p, Photospp::PhotosParticle::setE(), Photospp::PhotosParticle::setPz(), and sqrt().

Referenced by Photospp::PhotosParticle::boostFromRestFrame(), and Photospp::PhotosParticle::boostToRestFrame().

 {
         // Boost along the Z axis
         double m_tau=sqrt(boost_e*boost_e-boost_pz*boost_pz);
 
         double p=getPz();
         double e=getE();
 
         setPz((boost_e*p + boost_pz*e)/m_tau);
         setE((boost_pz*p + boost_e*e )/m_tau);
 }

void Photospp::PhotosParticle::boostAlongZ	(	double	pz,
		double	e
	)

inherited

Do a Lorenz transformation along the Z axis.

void PhotosParticle::boostDaughtersFromRestFrame ( PhotosParticle * boost )

inherited

Transform this particles four momentum from the lab frame to the rest frame of the parameter PhotosParticle.

Definition at line 63 of file PhotosParticle.cxx.

References Photospp::PhotosParticle::getAllDecayProducts(), and Photospp::PhotosParticle::hasDaughters().

Referenced by Photospp::PH_HEPEVT_Interface::get().

 {
         if(!hasDaughters()) //if there are no daughters
         return;
 
         // get all daughters, granddaughters, etc. then rotate and boost them
         vector<PhotosParticle*> list = getAllDecayProducts();
         vector<PhotosParticle*>::iterator pcl_itr = list.begin();
 
         for(;pcl_itr != list.end();pcl_itr++)
         {
                 (*pcl_itr)->boostFromRestFrame(tau_momentum);
         }
 
         //checkMomentumConservation();
 }

void Photospp::PhotosParticle::boostDaughtersFromRestFrame ( PhotosParticle * boost )

inherited

Transform this particles four momentum from the lab frame to the rest frame of the parameter PhotosParticle.

void PhotosParticle::boostDaughtersToRestFrame ( PhotosParticle * boost )

inherited

Transform the four momentum of all the daughters recursively into the frame of the "particle" PhotosParticle.

Definition at line 80 of file PhotosParticle.cxx.

References Photospp::PhotosParticle::getAllDecayProducts(), and Photospp::PhotosParticle::hasDaughters().

Referenced by Photospp::PH_HEPEVT_Interface::get().

 {
         if(!hasDaughters()) //if there are no daughters
         return;
 
         // get all daughters, granddaughters, etc. then rotate and boost them
         vector<PhotosParticle*> list = getAllDecayProducts();
         vector<PhotosParticle*>::iterator pcl_itr = list.begin();
 
         for(;pcl_itr != list.end();pcl_itr++)
         {
                 (*pcl_itr)->boostToRestFrame(tau_momentum);
         }
 
         //checkMomentumConservation();
 }

void Photospp::PhotosParticle::boostDaughtersToRestFrame ( PhotosParticle * boost )

inherited

Transform the four momentum of all the daughters recursively into the frame of the "particle" PhotosParticle.

void PhotosParticle::boostFromRestFrame ( PhotosParticle * boost )

inherited

Transform this particles four momentum from the rest frame of the paramter PhotosParticle, back into the lab frame.

Definition at line 114 of file PhotosParticle.cxx.

References Photospp::PhotosParticle::boostAlongZ(), Photospp::PhotosParticle::getE(), Photospp::PhotosParticle::getP(), Photospp::PhotosParticle::getRotationAngle(), phi, Photospp::PhotosParticle::rotate(), theta, Photospp::PhotosParticle::X_AXIS, and Photospp::PhotosParticle::Y_AXIS.

Referenced by Photospp::PH_HEPEVT_Interface::get().

 {
         //get the rotation angles
         //and boost z
 
         double theta = tau_momentum->getRotationAngle(Y_AXIS);
         tau_momentum->rotate(Y_AXIS,theta);
 
         double phi = tau_momentum->getRotationAngle(X_AXIS);
         tau_momentum->rotate(Y_AXIS,-theta);
 
         //Now rotate coordinates to get boost in Z direction.
         rotate(Y_AXIS,theta);
         rotate(X_AXIS,phi);
         boostAlongZ(tau_momentum->getP(),tau_momentum->getE());
         rotate(X_AXIS,-phi);
         rotate(Y_AXIS,-theta);
 }

void Photospp::PhotosParticle::boostFromRestFrame ( PhotosParticle * boost )

inherited

Transform this particles four momentum from the rest frame of the paramter PhotosParticle, back into the lab frame.

void Photospp::PhotosParticle::boostToRestFrame ( PhotosParticle * boost )

inherited

Transform this particles four momentum from the lab frome into the rest frame of the paramter PhotosParticle.

void PhotosParticle::boostToRestFrame ( PhotosParticle * boost )

inherited

Transform this particles four momentum from the lab frome into the rest frame of the paramter PhotosParticle.

Definition at line 98 of file PhotosParticle.cxx.

References Photospp::PhotosParticle::boostAlongZ(), Photospp::PhotosParticle::getE(), Photospp::PhotosParticle::getP(), Photospp::PhotosParticle::getRotationAngle(), phi, Photospp::PhotosParticle::rotate(), theta, Photospp::PhotosParticle::X_AXIS, and Photospp::PhotosParticle::Y_AXIS.

Referenced by Photospp::PH_HEPEVT_Interface::get().

 {
         double theta = tau_momentum->getRotationAngle(Y_AXIS);
         tau_momentum->rotate(Y_AXIS,theta);
 
         double phi = tau_momentum->getRotationAngle(X_AXIS);
         tau_momentum->rotate(Y_AXIS,-theta);
 
         //Now rotate coordinates to get boost in Z direction.
         rotate(Y_AXIS,theta);
         rotate(X_AXIS,phi);
         boostAlongZ(-1*tau_momentum->getP(),tau_momentum->getE());
         rotate(X_AXIS,-phi);
         rotate(Y_AXIS,-theta);
 }

bool Photospp::PhotosHEPEVTParticle::checkMomentumConservation ( )

virtual

Check that the 4 momentum in conserved in the decay of this particle

Implements Photospp::PhotosParticle.

bool PhotosHEPEVTParticle::checkMomentumConservation ( )

virtual

Check that the 4 momentum in conserved in the decay of this particle

Implements Photospp::PhotosParticle.

Definition at line 225 of file PhotosHEPEVTParticle.cxx.

References fabs(), getE(), getFirstMotherIndex(), Photospp::PhotosHEPEVTEvent::getParticle(), getPx(), getPy(), getPz(), getSecondMotherIndex(), i, m2(), m_daughter_end, m_daughter_start, m_event, print(), pz, Photospp::Log::RedirectOutput(), Photospp::Log::RevertOutput(), sqrt(), and Photospp::Log::Warning().

                                                     {
 
   if(!m_event)           return true;
   if(m_daughter_end < 0) return true;
 
   PhotosHEPEVTParticle *buf = m_event->getParticle(m_daughter_start);
 
   int first_mother_idx  = buf->getFirstMotherIndex();
   int second_mother_idx = buf->getSecondMotherIndex();
 
   double px =0.0, py =0.0, pz =0.0, e =0.0;
   double px2=0.0, py2=0.0, pz2=0.0, e2=0.0;
 
   for(int i=m_daughter_start;i<=m_daughter_end;i++)
   {
     buf = m_event->getParticle(i);
     px += buf->getPx();
     py += buf->getPy();
     pz += buf->getPz();
     e  += buf->getE ();
   }
 
   if(first_mother_idx>=0)
   {
     buf = m_event->getParticle(first_mother_idx);
     px2 += buf->getPx();
     py2 += buf->getPy();
     pz2 += buf->getPz();
     e2  += buf->getE();
   }
   
   if(second_mother_idx>=0)
   {
     buf = m_event->getParticle(second_mother_idx);
     px2 += buf->getPx();
     py2 += buf->getPy();
     pz2 += buf->getPz();
     e2  += buf->getE();
   }
   // 3-momentum  // test HepMC style
   double dp = sqrt( (px-px2)*(px-px2) + (py-py2)*(py-py2) + (pz-pz2)*(pz-pz2) );
   // virtuality test as well.
   double m1 = sqrt( fabs( e*e   - px*px   - py*py   - pz*pz   ) );
   double m2 = sqrt( fabs( e2*e2 - px2*px2 - py2*py2 - pz2*pz2 ) );
 
   if( fabs(m1-m2) > 0.0001 || dp > 0.0001*(e+e2))
   {
     Log::RedirectOutput( Log::Warning()<<"Momentum not conserved in vertex: " );
     if(first_mother_idx >=0) m_event->getParticle(first_mother_idx) ->print();
     if(second_mother_idx>=0) m_event->getParticle(second_mother_idx)->print();
     for(int i=m_daughter_start;i<=m_daughter_end;i++) m_event->getParticle(i)->print();
     Log::RevertOutput();
     return false;
   }
   
   return true;
 }

void PhotosHEPEVTParticle::createHistoryEntry ( )

virtual

Creating history entries not implemented in HEPEVT

Implements Photospp::PhotosParticle.

Definition at line 294 of file PhotosHEPEVTParticle.cxx.

References Photospp::Log::Warning().

 {
   Log::Warning()<<"PhotosParticle::createHistoryEntry() not implemented for HEPEVT."<<endl;
 }

void Photospp::PhotosHEPEVTParticle::createHistoryEntry ( )

virtual

Creating history entries not implemented in HEPEVT

Implements Photospp::PhotosParticle.

PhotosHEPEVTParticle * PhotosHEPEVTParticle::createNewParticle	(	int	pdg_id,
		int	status,
		double	mass,
		double	px,
		double	py,
		double	pz,
		double	e
	)

virtual

Creates a new particle of type PhotosHEPEVTParticle, with the given properties. The new particle bares no relations to this particle, but `this particle' provides only a way of creating an instance of this derived class. eg. createNewParticle() is used inside filhep_() so that a PhotosHEPEVTParticle can be created without the method having explicit knowledge of the PhotosHEPEVTParticle class

Implements Photospp::PhotosParticle.

Definition at line 283 of file PhotosHEPEVTParticle.cxx.

References cache, and PhotosHEPEVTParticle().

                                                                   {
 
   // New particles created using this method are added to cache
   // They will be deleted when this particle will be deleted
 
   cache.push_back(new PhotosHEPEVTParticle(pdg_id,status,px,py,pz,e,mass,-1,-1,-1,-1));
   return cache.back();
 }

PhotosHEPEVTParticle* Photospp::PhotosHEPEVTParticle::createNewParticle	(	int	pdg_id,
		int	status,
		double	mass,
		double	px,
		double	py,
		double	pz,
		double	e
	)

virtual

Creates a new particle of type PhotosHEPEVTParticle, with the given properties. The new particle bares no relations to this particle, but `this particle' provides only a way of creating an instance of this derived class. eg. createNewParticle() is used inside filhep_() so that a PhotosHEPEVTParticle can be created without the method having explicit knowledge of the PhotosHEPEVTParticle class

Implements Photospp::PhotosParticle.

void PhotosHEPEVTParticle::createSelfDecayVertex ( PhotosParticle * out )

virtual

Create a self-decay vertex for this particle with 'out' being the outgoing particle in new vertex

Implements Photospp::PhotosParticle.

Definition at line 299 of file PhotosHEPEVTParticle.cxx.

References Photospp::Log::Warning().

 {
   Log::Warning()<<"PhotosHEPEVTParticle::createSelfDecayVertex() not implemented for HEPEVT."<<endl;
 }

void Photospp::PhotosHEPEVTParticle::createSelfDecayVertex ( PhotosParticle * out )

virtual

Create a self-decay vertex for this particle with 'out' being the outgoing particle in new vertex

Implements Photospp::PhotosParticle.

PhotosParticle * PhotosParticle::findLastSelf ( )

inherited

Traverse the event structure and find the final version of this particle which does not have a particle of it's own type as it's daughter. eg. Generally the final stable copy

Definition at line 16 of file PhotosParticle.cxx.

References Photospp::PhotosParticle::getDaughters(), and Photospp::PhotosParticle::getPdgID().

 {
         vector<PhotosParticle*> daughters = getDaughters();
         vector<PhotosParticle*>::iterator pcl_itr = daughters.begin();
 
         //get all daughters and look for stable with same pgd id
         for(;pcl_itr != daughters.end();pcl_itr++)
         {
                 if((*pcl_itr)->getPdgID()==this->getPdgID())
                 return (*pcl_itr)->findLastSelf();
         }
 
         return this;
 }

PhotosParticle* Photospp::PhotosParticle::findLastSelf ( )

inherited

Traverse the event structure and find the final version of this particle which does not have a particle of it's own type as it's daughter. eg. Generally the final stable copy

vector< PhotosParticle * > PhotosParticle::findProductionMothers ( )

inherited

Traverse the event structure and find the first set of mothers which are not of the same type as this particle.

Definition at line 31 of file PhotosParticle.cxx.

References Photospp::PhotosParticle::getMothers(), and Photospp::PhotosParticle::getPdgID().

Referenced by Photospp::PhotosBranch::PhotosBranch().

 {
         vector<PhotosParticle*> mothers = getMothers();
         vector<PhotosParticle*>::iterator pcl_itr = mothers.begin();
 
         //get all mothers and check none have pdg id of this one
         for(;pcl_itr != mothers.end();pcl_itr++)
         {
                 if((*pcl_itr)->getPdgID()==this->getPdgID())
                 return (*pcl_itr)->findProductionMothers();
         }
         return mothers;
 }

std::vector<PhotosParticle *> Photospp::PhotosParticle::findProductionMothers ( )

inherited

Traverse the event structure and find the first set of mothers which are not of the same type as this particle.

std::vector<PhotosParticle*> Photospp::PhotosHEPEVTParticle::getAllDecayProducts ( )

virtual

Returns all particles in the decay tree of this particle via a vector of PhotosParticle

Implements Photospp::PhotosParticle.

std::vector< PhotosParticle * > PhotosHEPEVTParticle::getAllDecayProducts ( )

virtual

Returns all particles in the decay tree of this particle via a vector of PhotosParticle

Implements Photospp::PhotosParticle.

Definition at line 185 of file PhotosHEPEVTParticle.cxx.

References getBarcode(), getDaughters(), Photospp::PhotosParticle::hasDaughters(), and i.

 {
   std::vector<PhotosParticle*> list;
 
   if(!hasDaughters()) // if this particle has no daughters
     return list;
 
   std::vector<PhotosParticle*> daughters = getDaughters();
   
   // copy daughters to list of all decay products
   list.insert(list.end(),daughters.begin(),daughters.end());
   
   // Now, get all daughters recursively, without duplicates.
   // That is, for each daughter:
   // 1)  get list of her daughters
   // 2)  for each particle on this list:
   //  a) check if it is already on the list
   //  b) if it's not, add her to the end of the list
   for(unsigned int i=0;i<list.size();i++)
   {
     std::vector<PhotosParticle*> daughters2 = list[i]->getDaughters();
 
     if(!list[i]->hasDaughters()) continue;
     for(unsigned int j=0;j<daughters2.size();j++)
     {
       bool add=true;
       for(unsigned int k=0;k<list.size();k++)
         if( daughters2[j]->getBarcode() == list[k]->getBarcode() )
         {
           add=false;
           break;
         }
 
       if(add) list.push_back(daughters2[j]);
     }
   }
 
   return list;
 }

int Photospp::PhotosHEPEVTParticle::getBarcode ( )

virtual

Get the barcode (position in list) of this particle

Implements Photospp::PhotosParticle.

int PhotosHEPEVTParticle::getBarcode ( )

virtual

Get the barcode (position in list) of this particle

Implements Photospp::PhotosParticle.

Definition at line 388 of file PhotosHEPEVTParticle.cxx.

References m_barcode.

Referenced by getAllDecayProducts(), isDaughterOf(), isMotherOf(), and setMothers().

                                     {
   return m_barcode;
 }

int PhotosHEPEVTParticle::getDaughterRangeEnd ( )

Get index of last daughter

Definition at line 412 of file PhotosHEPEVTParticle.cxx.

References m_daughter_end.

Referenced by addDaughter().

                                              {
   return m_daughter_end;
 }

int Photospp::PhotosHEPEVTParticle::getDaughterRangeEnd ( )

Get index of last daughter

int Photospp::PhotosHEPEVTParticle::getDaughterRangeStart ( )

Get index of first daughter

int PhotosHEPEVTParticle::getDaughterRangeStart ( )

Get index of first daughter

Definition at line 408 of file PhotosHEPEVTParticle.cxx.

References m_daughter_start.

Referenced by addDaughter().

                                                {
   return m_daughter_start;
 }

std::vector< PhotosParticle * > PhotosHEPEVTParticle::getDaughters ( )

virtual

Returns the daughters of this particle via a vector of PhotosParticle

Implements Photospp::PhotosParticle.

Definition at line 138 of file PhotosHEPEVTParticle.cxx.

References Photospp::PhotosHEPEVTEvent::getParticle(), Photospp::PhotosHEPEVTEvent::getParticleCount(), i, isDaughterOf(), m_daughter_end, m_daughter_start, m_event, m_status, p, and Photospp::Log::Warning().

Referenced by getAllDecayProducts().

                                                              {
 
   std::vector<PhotosParticle*> daughters;
 
   if(!m_event) return daughters;
 
   // Check if m_daughter_start and m_daughter_end are set
   // If not - try to get list of daughters from event
   if(m_daughter_end<0)
   {
     int min_d=65535, max_d=-1;
     for(int i=0;i<m_event->getParticleCount();i++)
     {
       if(m_event->getParticle(i)->isDaughterOf(this))
       {
         if(i<min_d) min_d = i;
         if(i>max_d) max_d = i;
       }
     }
     if(max_d>=0)
     {
       m_daughter_start = min_d;
       m_daughter_end   = max_d;
       m_status         = 2;
     }
   }
 
   // If m_daughter_end is still not set - there are no daughters
   // Otherwsie - get daughters
   if(m_daughter_end>=0)
   {
     for(int i=m_daughter_start;i<=m_daughter_end;i++)
     {
       PhotosParticle *p = m_event->getParticle(i);
       if(p==NULL)
       {
         Log::Warning()<<"PhotosHEPEVTParticle::getDaughters(): No particle with index "<<i<<endl;
         return daughters;
       }
 
       daughters.push_back(p);
     }
   }
 
   return daughters;
 }

std::vector<PhotosParticle*> Photospp::PhotosHEPEVTParticle::getDaughters ( )

virtual

Returns the daughters of this particle via a vector of PhotosParticle

Implements Photospp::PhotosParticle.

vector< PhotosParticle * > PhotosParticle::getDecayTree ( )

inherited

Return whole decay tree starting from this particle

Definition at line 45 of file PhotosParticle.cxx.

References Photospp::PhotosParticle::getBarcode(), Photospp::PhotosParticle::getDaughters(), Photospp::PhotosParticle::getDecayTree(), Photospp::PhotosParticle::getMothers(), i, p, and tree.

Referenced by Photospp::PhotosBranch::createBranches(), Photospp::PhotosParticle::getDecayTree(), and Photospp::Photos::processBranch().

 {
         vector<PhotosParticle *> particles;
         particles.push_back(this);
         vector<PhotosParticle *> daughters = getDaughters();
         for(int i=0;i<(int)daughters.size();i++)
         {
                 // Check if we are the first mother of each daughters
                 // If not - skip this daughter
                 PhotosParticle *p = daughters.at(i);
                 vector<PhotosParticle *> mothers = p->getMothers();
                 if(mothers.size()>1 && mothers.at(0)->getBarcode()!=getBarcode()) continue;
                 vector<PhotosParticle *> tree = p->getDecayTree();
                 particles.insert(particles.end(),tree.begin(),tree.end());
         }
         return particles;
 }

std::vector<PhotosParticle *> Photospp::PhotosParticle::getDecayTree ( )

inherited

Return whole decay tree starting from this particle

double PhotosHEPEVTParticle::getE ( )

virtual

Returns the energy component of the four vector

Implements Photospp::PhotosParticle.

Definition at line 367 of file PhotosHEPEVTParticle.cxx.

References m_e.

Referenced by checkMomentumConservation(), and checkMomentumConservationInEvent().

                                  {
   return m_e;
 }

double Photospp::PhotosHEPEVTParticle::getE ( )

virtual

Returns the energy component of the four vector

Implements Photospp::PhotosParticle.

int Photospp::PhotosHEPEVTParticle::getFirstMotherIndex ( )

Get index of first mother

int PhotosHEPEVTParticle::getFirstMotherIndex ( )

Get index of first mother

Definition at line 400 of file PhotosHEPEVTParticle.cxx.

References m_first_mother.

Referenced by checkMomentumConservation().

                                              {
   return m_first_mother;
 }

double PhotosHEPEVTParticle::getMass ( )

virtual

Get the mass stored (i.e. not calculated from four vector) at generation step

Implements Photospp::PhotosParticle.

Definition at line 351 of file PhotosHEPEVTParticle.cxx.

References m_generated_mass.

                                     {
   return m_generated_mass;
 }

double Photospp::PhotosHEPEVTParticle::getMass ( )

virtual

Get the mass stored (i.e. not calculated from four vector) at generation step

Implements Photospp::PhotosParticle.

std::vector< PhotosParticle * > PhotosHEPEVTParticle::getMothers ( )

virtual

Returns the mothers of this particle via a vector of PhotosParticle

Implements Photospp::PhotosParticle.

Definition at line 120 of file PhotosHEPEVTParticle.cxx.

References Photospp::PhotosHEPEVTEvent::getParticle(), m_event, m_first_mother, m_second_mother, p1, and p2.

                                                            {
 
   std::vector<PhotosParticle*> mothers;
 
   PhotosParticle *p1 = NULL;
   PhotosParticle *p2 = NULL;
   
   if(m_first_mother>=0)  p1 = m_event->getParticle(m_first_mother);
   if(m_second_mother>=0) p2 = m_event->getParticle(m_second_mother);
 
   if(p1) mothers.push_back(p1);
   if(p2) mothers.push_back(p2);
 
   return mothers;
 }

std::vector<PhotosParticle*> Photospp::PhotosHEPEVTParticle::getMothers ( )

virtual

Returns the mothers of this particle via a vector of PhotosParticle

Implements Photospp::PhotosParticle.

double PhotosParticle::getP ( )

inherited

Get scalar momentum

Definition at line 202 of file PhotosParticle.cxx.

References Photospp::PhotosParticle::getPx(), Photospp::PhotosParticle::getPy(), Photospp::PhotosParticle::getPz(), and sqrt().

Referenced by Photospp::PhotosParticle::boostFromRestFrame(), Photospp::PhotosParticle::boostToRestFrame(), Photospp::PhotosParticle::getRotationAngle(), Photospp::PhotosParticle::getVirtuality(), and Photospp::PhotosParticle::rotate().

 {
         return sqrt(getPx()*getPx()+getPy()*getPy()+getPz()*getPz());
 }

double Photospp::PhotosParticle::getP ( )

inherited

Get scalar momentum

double PhotosParticle::getP ( int axis )

inherited

Get momentum component in the direction of "axis" (x,y,z)

Definition at line 207 of file PhotosParticle.cxx.

References Photospp::PhotosParticle::getPx(), Photospp::PhotosParticle::getPy(), Photospp::PhotosParticle::getPz(), Photospp::PhotosParticle::X_AXIS, Photospp::PhotosParticle::Y_AXIS, and Photospp::PhotosParticle::Z_AXIS.

 {
         if(axis==X_AXIS) return getPx();
         if(axis==Y_AXIS) return getPy();
         if(axis==Z_AXIS) return getPz();
         return 0;
 }

double Photospp::PhotosParticle::getP ( int axis )

inherited

Get momentum component in the direction of "axis" (x,y,z)

int PhotosHEPEVTParticle::getPdgID ( )

virtual

Get the PDG ID code of this particle

Implements Photospp::PhotosParticle.

Definition at line 343 of file PhotosHEPEVTParticle.cxx.

References m_pdgid.

                                   {
   return m_pdgid;
 }

int Photospp::PhotosHEPEVTParticle::getPdgID ( )

virtual

Get the PDG ID code of this particle

Implements Photospp::PhotosParticle.

double PhotosHEPEVTParticle::getPx ( )

inlinevirtual

Returns the px component of the four vector

Implements Photospp::PhotosParticle.

Definition at line 355 of file PhotosHEPEVTParticle.cxx.

References m_px.

Referenced by checkMomentumConservation(), and checkMomentumConservationInEvent().

                                          {
   return m_px;
 }

double Photospp::PhotosHEPEVTParticle::getPx ( )

virtual

Returns the px component of the four vector

Implements Photospp::PhotosParticle.

double PhotosHEPEVTParticle::getPy ( )

inlinevirtual

Returns the py component of the four vector

Implements Photospp::PhotosParticle.

Definition at line 359 of file PhotosHEPEVTParticle.cxx.

References m_py.

Referenced by checkMomentumConservation(), and checkMomentumConservationInEvent().

                                          {
   return m_py;
 }

double Photospp::PhotosHEPEVTParticle::getPy ( )

virtual

Returns the py component of the four vector

Implements Photospp::PhotosParticle.

double PhotosHEPEVTParticle::getPz ( )

virtual

Returns the pz component of the four vector

Implements Photospp::PhotosParticle.

Definition at line 363 of file PhotosHEPEVTParticle.cxx.

References m_pz.

Referenced by checkMomentumConservation(), and checkMomentumConservationInEvent().

                                   {
   return m_pz;
 }

double Photospp::PhotosHEPEVTParticle::getPz ( )

virtual

Returns the pz component of the four vector

Implements Photospp::PhotosParticle.

double Photospp::PhotosParticle::getRotationAngle	(	int	axis,
		int	second_axis = `Z_AXIS`
	)

inherited

Returns the angle around the axis "axis" needed to rotate the four momenum is such a way that the non-Z component disappears and Z>0. This is used to in rotating the coordinate system into a frame with only a Z component before calling boostAlongZ().

double PhotosParticle::getRotationAngle	(	int	axis,
		int	second_axis = `Z_AXIS`
	)

inherited

Returns the angle around the axis "axis" needed to rotate the four momenum is such a way that the non-Z component disappears and Z>0. This is used to in rotating the coordinate system into a frame with only a Z component before calling boostAlongZ().

Get the angle needed to rotate the 4 momentum vector so that the x (y) component disapears. (and the Z component is > 0)

if(getP(axis)==0){ if(getPz()>0) return 0; //no rotaion required else return M_PI; }

Definition at line 135 of file PhotosParticle.cxx.

References Photospp::PhotosParticle::getP().

Referenced by Photospp::PhotosParticle::boostFromRestFrame(), and Photospp::PhotosParticle::boostToRestFrame().

 {
         if(getP(second_axis)==0)
         {
                 if(getP(axis)>0) return -M_PI/2.0;
                 else             return  M_PI/2.0;
         }
         if(getP(second_axis)>0) return     -atan(getP(axis)/getP(second_axis));
         else                    return M_PI-atan(getP(axis)/getP(second_axis));
 
 }

int PhotosHEPEVTParticle::getSecondMotherIndex ( )

Get index of second mother

Definition at line 404 of file PhotosHEPEVTParticle.cxx.

References m_second_mother.

Referenced by checkMomentumConservation().

                                               {
   return m_second_mother;
 }

int Photospp::PhotosHEPEVTParticle::getSecondMotherIndex ( )

Get index of second mother

int PhotosHEPEVTParticle::getStatus ( )

virtual

Get the status of this particle

Implements Photospp::PhotosParticle.

Definition at line 347 of file PhotosHEPEVTParticle.cxx.

References m_status.

Referenced by checkMomentumConservationInEvent().

                                    {
   return m_status;
 }

int Photospp::PhotosHEPEVTParticle::getStatus ( )

virtual

Get the status of this particle

Implements Photospp::PhotosParticle.

double PhotosParticle::getVirtuality ( )

virtualinherited

Get sqrt(e^2-p^2)

Definition at line 193 of file PhotosParticle.cxx.

References Photospp::PhotosParticle::getE(), Photospp::PhotosParticle::getP(), and sqrt().

Referenced by Photospp::PH_HEPEVT_Interface::add_particle().

 {
         double e_sq=getE()*getE();
         double p_sq=getP()*getP();
 
         if(e_sq>p_sq) return    sqrt(e_sq-p_sq);
         else          return -1*sqrt(p_sq-e_sq); //if it's negative
 }

virtual double Photospp::PhotosParticle::getVirtuality ( )

virtualinherited

Get sqrt(e^2-p^2)

bool PhotosParticle::hasDaughters ( )

inherited

Return whether the particle has any chidren

Definition at line 10 of file PhotosParticle.cxx.

References Photospp::PhotosParticle::getDaughters().

Referenced by Photospp::PhotosParticle::boostDaughtersFromRestFrame(), Photospp::PhotosParticle::boostDaughtersToRestFrame(), getAllDecayProducts(), Photospp::PhotosHepMCParticle::getAllDecayProducts(), and Photospp::PhotosParticle::rotateDaughters().

 {
         if(getDaughters().size()==0) return false;
         else                         return true;
 }

bool Photospp::PhotosParticle::hasDaughters ( )

inherited

Return whether the particle has any chidren

bool Photospp::PhotosHEPEVTParticle::isDaughterOf ( PhotosHEPEVTParticle * p )

Check if particle 'p' is daughter of this particle

bool PhotosHEPEVTParticle::isDaughterOf ( PhotosHEPEVTParticle * p )

Check if particle 'p' is daughter of this particle

Definition at line 304 of file PhotosHEPEVTParticle.cxx.

References getBarcode(), m_first_mother, and m_second_mother.

Referenced by getDaughters().

 {
   int bc = p->getBarcode();
   if(bc==m_first_mother || bc==m_second_mother) return true;
 
   return false;
 }

bool Photospp::PhotosHEPEVTParticle::isMotherOf ( PhotosHEPEVTParticle * p )

Check if particle 'p' is mother of this particle

bool PhotosHEPEVTParticle::isMotherOf ( PhotosHEPEVTParticle * p )

Check if particle 'p' is mother of this particle

Definition at line 312 of file PhotosHEPEVTParticle.cxx.

References getBarcode(), m_daughter_end, and m_daughter_start.

 {
   int bc = p->getBarcode();
   if(bc>=m_daughter_start && bc<=m_daughter_end) return true;
 
   return false;
 }

void Photospp::PhotosHEPEVTParticle::print ( )

virtual

Print information on this particle into standard output

Implements Photospp::PhotosParticle.

void PhotosHEPEVTParticle::print ( )

virtual

Print information on this particle into standard output

Implements Photospp::PhotosParticle.

Definition at line 320 of file PhotosHEPEVTParticle.cxx.

References m_barcode, m_daughter_end, m_daughter_start, m_e, m_first_mother, m_generated_mass, m_pdgid, m_px, m_py, m_pz, m_second_mother, and m_status.

Referenced by checkMomentumConservation().

                                 {
   char buf[256];
   sprintf(buf,"P: (%2i) %6i %2i | %11.4e %11.4e %11.4e %11.4e | %11.4e | M: %2i %2i | D: %2i %2i\n",
           m_barcode, m_pdgid, m_status, m_px, m_py, m_pz, m_e, m_generated_mass,
           m_first_mother, m_second_mother,   m_daughter_start, m_daughter_end);
 
   cout<<buf;
 }

void PhotosParticle::rotate	(	int	axis,
		double	theta,
		int	second_axis = `Z_AXIS`
	)

inherited

rotate this particles 4-momentum by an angle phi from the axisis "axis" towards the axis "second_axis".

Rotation around an axis X or Y

Definition at line 168 of file PhotosParticle.cxx.

References cos(), Photospp::PhotosParticle::getP(), Photospp::PhotosParticle::setP(), and sin().

Referenced by Photospp::PhotosParticle::boostFromRestFrame(), and Photospp::PhotosParticle::boostToRestFrame().

 {
         double temp_px=getP(axis);
         double temp_pz=getP(second_axis);
         setP(axis,cos(theta)*temp_px + sin(theta)*temp_pz);
         setP(second_axis,-sin(theta)*temp_px + cos(theta)*temp_pz);
 }

void Photospp::PhotosParticle::rotate	(	int	axis,
		double	phi,
		int	second_axis = `Z_AXIS`
	)

inherited

rotate this particles 4-momentum by an angle phi from the axisis "axis" towards the axis "second_axis".

void Photospp::PhotosParticle::rotateDaughters	(	int	axis,
		double	phi,
		int	second_axis = `Z_AXIS`
	)

inherited

rotate 4-momentum of daughters of this particle by an angle phi from the axisis "axis" towards the axis "second_axis".

void PhotosParticle::rotateDaughters	(	int	axis,
		double	phi,
		int	second_axis = `Z_AXIS`
	)

inherited

rotate 4-momentum of daughters of this particle by an angle phi from the axisis "axis" towards the axis "second_axis".

Definition at line 176 of file PhotosParticle.cxx.

References Photospp::PhotosParticle::getDaughters(), and Photospp::PhotosParticle::hasDaughters().

 {
         if(!hasDaughters()) //if there are no daughters
         return;
 
         vector<PhotosParticle*> daughters = getDaughters();
         vector<PhotosParticle*>::iterator pcl_itr = daughters.begin();
 
         //get all daughters then rotate and boost them.
         for(;pcl_itr != daughters.end();pcl_itr++)
         {
                 (*pcl_itr)->rotate(axis,theta,second_axis);
                 (*pcl_itr)->rotateDaughters(axis,theta,second_axis);
         }
         //checkMomentumConservation();
 }

void Photospp::PhotosHEPEVTParticle::setBarcode ( int barcode )

Set barcode (position in list) of this particle

void PhotosHEPEVTParticle::setBarcode ( int barcode )

Set barcode (position in list) of this particle

Definition at line 392 of file PhotosHEPEVTParticle.cxx.

References m_barcode.

Referenced by addDaughter(), and Photospp::PhotosHEPEVTEvent::addParticle().

                                                 {
   m_barcode = barcode;
 }

void Photospp::PhotosHEPEVTParticle::setDaughterRangeEnd ( int i )

inlineprivate

Set index of last daughter

Definition at line 159 of file 3.52/src/eventRecordInterfaces/PhotosHEPEVTParticle.h.

References i, and m_daughter_end.

159 { m_daughter_end =i; }

i

Int_t i

Definition: run_full.C:25

Photospp::PhotosHEPEVTParticle::m_daughter_end

int m_daughter_end

Definition: PhotosHEPEVTParticle.h:171

void Photospp::PhotosHEPEVTParticle::setDaughterRangeEnd ( int i )

inlineprivate

Set index of last daughter

Definition at line 159 of file PhotosHEPEVTParticle.h.

References i, and m_daughter_end.

Referenced by addDaughter().

159 { m_daughter_end =i; }

i

Int_t i

Definition: run_full.C:25

Photospp::PhotosHEPEVTParticle::m_daughter_end

int m_daughter_end

Definition: PhotosHEPEVTParticle.h:171

void Photospp::PhotosHEPEVTParticle::setDaughterRangeStart ( int i )

inlineprivate

Set index of first daughter

Definition at line 156 of file PhotosHEPEVTParticle.h.

References i, and m_daughter_start.

Referenced by addDaughter().

156 { m_daughter_start=i; }

i

Int_t i

Definition: run_full.C:25

Photospp::PhotosHEPEVTParticle::m_daughter_start

int m_daughter_start

Definition: PhotosHEPEVTParticle.h:171

void Photospp::PhotosHEPEVTParticle::setDaughterRangeStart ( int i )

inlineprivate

Set index of first daughter

Definition at line 156 of file 3.52/src/eventRecordInterfaces/PhotosHEPEVTParticle.h.

References i, and m_daughter_start.

156 { m_daughter_start=i; }

i

Int_t i

Definition: run_full.C:25

Photospp::PhotosHEPEVTParticle::m_daughter_start

int m_daughter_start

Definition: PhotosHEPEVTParticle.h:171

void PhotosHEPEVTParticle::setDaughters ( std::vector< PhotosParticle * > daughters )

virtual

Set the daughters of this particle via a vector of PhotosParticle

Implements Photospp::PhotosParticle.

Definition at line 100 of file PhotosHEPEVTParticle.cxx.

References Photospp::Log::Fatal(), i, m_daughter_end, m_daughter_start, and m_event.

                                                                         {
 
   // This particle must be inside some event record to be able to add daughters
   if(m_event==NULL) Log::Fatal("PhotosHEPEVTParticle::setDaughters: particle not inside event record.");
 
   int beg = 65535, end = -1;
 
   for(unsigned int i=0;i<daughters.size();i++)
   {
     int bc = daughters[i]->getBarcode();
     if(bc<0) Log::Fatal("PhotosHEPEVTParticle::setDaughters: all daughters has to be in event record first");
     if(bc<beg) beg = bc;
     if(bc>end) end = bc;
   }
   if(end == -1) beg = -1;
 
   m_daughter_start = beg;
   m_daughter_end   = end;
 }

void Photospp::PhotosHEPEVTParticle::setDaughters ( std::vector< PhotosParticle * > daughters )

virtual

Set the daughters of this particle via a vector of PhotosParticle

Implements Photospp::PhotosParticle.

void PhotosHEPEVTParticle::setE ( double e )

virtual

Set the energy component of the four vector

Implements Photospp::PhotosParticle.

Definition at line 384 of file PhotosHEPEVTParticle.cxx.

References m_e.

                                        {
   m_e = e;
 }

void Photospp::PhotosHEPEVTParticle::setE ( double e )

virtual

Set the energy component of the four vector

Implements Photospp::PhotosParticle.

void PhotosHEPEVTParticle::setEvent ( PhotosHEPEVTEvent * event )

Set event of this particle

Definition at line 396 of file PhotosHEPEVTParticle.cxx.

References m_event.

Referenced by Photospp::PhotosHEPEVTEvent::addParticle().

                                                            {
   m_event = event;
 }

void Photospp::PhotosHEPEVTParticle::setEvent ( PhotosHEPEVTEvent * event )

Set event of this particle

void Photospp::PhotosHEPEVTParticle::setMass ( double mass )

virtual

Set the mass of this particle

Implements Photospp::PhotosParticle.

void PhotosHEPEVTParticle::setMass ( double mass )

virtual

Set the mass of this particle

Implements Photospp::PhotosParticle.

Definition at line 339 of file PhotosHEPEVTParticle.cxx.

References m_generated_mass.

                                              {
   m_generated_mass = mass;
 }

void Photospp::PhotosHEPEVTParticle::setMothers ( std::vector< PhotosParticle * > mothers )

virtual

Set the mothers of this particle via a vector of PhotosParticle

Implements Photospp::PhotosParticle.

void PhotosHEPEVTParticle::setMothers ( std::vector< PhotosParticle * > mothers )

virtual

Set the mothers of this particle via a vector of PhotosParticle

Implements Photospp::PhotosParticle.

Definition at line 84 of file PhotosHEPEVTParticle.cxx.

References Photospp::PhotosHEPEVTEvent::addParticle(), evt, Photospp::Log::Fatal(), getBarcode(), m_barcode, m_first_mother, and m_second_mother.

                                                                     {
 
   // If this particle has not yet been added to the event record
   // then add it to the mothers' event record
   if(m_barcode<0 && mothers.size()>0)
   {
     PhotosHEPEVTEvent *evt = ((PhotosHEPEVTParticle*)mothers[0])->m_event;
     evt->addParticle(this);
   }
 
   if(mothers.size()>2) Log::Fatal("PhotosHEPEVTParticle::setMothers: HEPEVT does not allow more than two mothers!");
   
   if(mothers.size()>0) m_first_mother  = mothers[0]->getBarcode();
   if(mothers.size()>1) m_second_mother = mothers[1]->getBarcode();
 }

void PhotosParticle::setP	(	int	axis,
		double	p_component
	)

inherited

Set momentum component in the direction of "axis" (x,y,z)

Definition at line 215 of file PhotosParticle.cxx.

References Photospp::PhotosParticle::setPx(), Photospp::PhotosParticle::setPy(), Photospp::PhotosParticle::setPz(), Photospp::PhotosParticle::X_AXIS, Photospp::PhotosParticle::Y_AXIS, and Photospp::PhotosParticle::Z_AXIS.

Referenced by Photospp::PhotosParticle::rotate().

 {
         if(axis==X_AXIS) setPx(p_component);
         if(axis==Y_AXIS) setPy(p_component);
         if(axis==Z_AXIS) setPz(p_component);
 }

void Photospp::PhotosParticle::setP	(	int	axis,
		double	p_component
	)

inherited

Set momentum component in the direction of "axis" (x,y,z)

void PhotosHEPEVTParticle::setPdgID ( int pdg_id )

virtual

Set the PDG ID code of this particle

Implements Photospp::PhotosParticle.

Definition at line 331 of file PhotosHEPEVTParticle.cxx.

References m_pdgid.

                                              {
   m_pdgid = pdg_id;
 }

void Photospp::PhotosHEPEVTParticle::setPdgID ( int pdg_id )

virtual

Set the PDG ID code of this particle

Implements Photospp::PhotosParticle.

void PhotosHEPEVTParticle::setPx ( double px )

virtual

Set the px component of the four vector

Implements Photospp::PhotosParticle.

Definition at line 371 of file PhotosHEPEVTParticle.cxx.

References m_px.

                                          {
   m_px = px;
 }

void Photospp::PhotosHEPEVTParticle::setPx ( double px )

virtual

Set the px component of the four vector

Implements Photospp::PhotosParticle.

void PhotosHEPEVTParticle::setPy ( double py )

virtual

Set the px component of the four vector

Implements Photospp::PhotosParticle.

Definition at line 375 of file PhotosHEPEVTParticle.cxx.

References m_py.

                                          {
   m_py = py;
 }

void Photospp::PhotosHEPEVTParticle::setPy ( double py )

virtual

Set the px component of the four vector

Implements Photospp::PhotosParticle.

void Photospp::PhotosHEPEVTParticle::setPz ( double pz )

virtual

Set the pz component of the four vector

Implements Photospp::PhotosParticle.

void PhotosHEPEVTParticle::setPz ( double pz )

virtual

Set the pz component of the four vector

Implements Photospp::PhotosParticle.

Definition at line 380 of file PhotosHEPEVTParticle.cxx.

References m_pz, and pz.

                                          {
   m_pz = pz;
 }

void Photospp::PhotosHEPEVTParticle::setStatus ( int statu )

virtual

Set the status of this particle

Implements Photospp::PhotosParticle.

void PhotosHEPEVTParticle::setStatus ( int statu )

virtual

Set the status of this particle

Implements Photospp::PhotosParticle.

Definition at line 335 of file PhotosHEPEVTParticle.cxx.

References m_status, and status.

                                               {
   m_status = status;
 }

Member Data Documentation

static const int PhotosParticle::ANTIDOWN = -1

staticinherited

anti-down quark

Definition at line 119 of file PhotosParticle.h.

static const int PhotosParticle::ANTIUP = -2

staticinherited

anti-up quark

Definition at line 113 of file PhotosParticle.h.

vector< PhotosHEPEVTParticle * > PhotosHEPEVTParticle::cache

private

List of created particles - if they are not in the event, they will be deleted when no longer needed

Definition at line 187 of file PhotosHEPEVTParticle.h.

Referenced by createNewParticle(), and ~PhotosHEPEVTParticle().

static const int PhotosParticle::DECAYED =2

staticinherited

Decayed particle status

Definition at line 32 of file PhotosParticle.h.

Referenced by Photospp::PhotosHepMCParticle::setMothers().

static const int PhotosParticle::DOWN = 1

staticinherited

down quark

Definition at line 116 of file PhotosParticle.h.

static const int PhotosParticle::ELECTRON = 11

staticinherited

e- particle

Definition at line 101 of file PhotosParticle.h.

static const int PhotosParticle::ELECTRON_ANTINEUTRINO = -12

staticinherited

e antineutrino particle

Definition at line 107 of file PhotosParticle.h.

static const int PhotosParticle::ELECTRON_NEUTRINO = 12

staticinherited

e neutrino particle

Definition at line 104 of file PhotosParticle.h.

static const int PhotosParticle::GAMMA = 22

staticinherited

photon

Definition at line 71 of file PhotosParticle.h.

static const int PhotosParticle::HIGGS = 25

staticinherited

H particle

Definition at line 50 of file PhotosParticle.h.

static const int PhotosParticle::HIGGS_A = 36

staticinherited

A0 particle

Definition at line 56 of file PhotosParticle.h.

static const int PhotosParticle::HIGGS_H = 35

staticinherited

H0 particle

Definition at line 53 of file PhotosParticle.h.

static const int PhotosParticle::HIGGS_MINUS = -37

staticinherited

H- particle

Definition at line 62 of file PhotosParticle.h.

static const int PhotosParticle::HIGGS_PLUS = 37

staticinherited

H+ particle

Definition at line 59 of file PhotosParticle.h.

static const int PhotosParticle::HISTORY =3

staticinherited

History particle status

Definition at line 35 of file PhotosParticle.h.

int PhotosHEPEVTParticle::m_barcode

private

Position in the event record

Definition at line 165 of file PhotosHEPEVTParticle.h.

Referenced by getBarcode(), PhotosHEPEVTParticle(), print(), setBarcode(), setMothers(), and ~PhotosHEPEVTParticle().

int PhotosHEPEVTParticle::m_daughter_end

private

Definition at line 171 of file PhotosHEPEVTParticle.h.

Referenced by addDaughter(), checkMomentumConservation(), getDaughterRangeEnd(), getDaughters(), isMotherOf(), PhotosHEPEVTParticle(), print(), setDaughterRangeEnd(), and setDaughters().

int PhotosHEPEVTParticle::m_daughter_start

private

Range of indexes of daughters (-1 if do not have daughters)

Definition at line 171 of file PhotosHEPEVTParticle.h.

Referenced by addDaughter(), checkMomentumConservation(), getDaughterRangeStart(), getDaughters(), isMotherOf(), PhotosHEPEVTParticle(), print(), setDaughterRangeStart(), and setDaughters().

double PhotosHEPEVTParticle::m_e

private

Definition at line 180 of file PhotosHEPEVTParticle.h.

Referenced by getE(), PhotosHEPEVTParticle(), print(), and setE().

PhotosHEPEVTEvent * PhotosHEPEVTParticle::m_event

private

Event from which this particle is taken

Definition at line 162 of file PhotosHEPEVTParticle.h.

Referenced by addDaughter(), checkMomentumConservation(), getDaughters(), getMothers(), PhotosHEPEVTParticle(), setDaughters(), and setEvent().

int PhotosHEPEVTParticle::m_first_mother

private

Indexes of mothers (-1 if do not have mothers)

Definition at line 168 of file PhotosHEPEVTParticle.h.

Referenced by getFirstMotherIndex(), getMothers(), isDaughterOf(), PhotosHEPEVTParticle(), print(), and setMothers().

double PhotosHEPEVTParticle::m_generated_mass

private

Mass saved at generation step

Definition at line 183 of file PhotosHEPEVTParticle.h.

Referenced by getMass(), PhotosHEPEVTParticle(), print(), and setMass().

int PhotosHEPEVTParticle::m_pdgid

private

PDG ID

Definition at line 174 of file PhotosHEPEVTParticle.h.

Referenced by getPdgID(), PhotosHEPEVTParticle(), print(), and setPdgID().

double PhotosHEPEVTParticle::m_px

private

Momentum

Definition at line 180 of file PhotosHEPEVTParticle.h.

Referenced by getPx(), PhotosHEPEVTParticle(), print(), and setPx().

double PhotosHEPEVTParticle::m_py

private

Definition at line 180 of file PhotosHEPEVTParticle.h.

Referenced by getPy(), PhotosHEPEVTParticle(), print(), and setPy().

double PhotosHEPEVTParticle::m_pz

private

Definition at line 180 of file PhotosHEPEVTParticle.h.

Referenced by getPz(), PhotosHEPEVTParticle(), print(), and setPz().

int PhotosHEPEVTParticle::m_second_mother

private

Definition at line 168 of file PhotosHEPEVTParticle.h.

Referenced by getMothers(), getSecondMotherIndex(), isDaughterOf(), PhotosHEPEVTParticle(), print(), and setMothers().

int PhotosHEPEVTParticle::m_status

private

Status (stable, decayed)

Definition at line 177 of file PhotosHEPEVTParticle.h.

Referenced by getDaughters(), getStatus(), PhotosHEPEVTParticle(), print(), and setStatus().

static const int PhotosParticle::MUON_ANTINEUTRINO = -14

staticinherited

muon antineutrino particle

Definition at line 95 of file PhotosParticle.h.

static const int PhotosParticle::MUON_MINUS = 13

staticinherited

muon- particle

Definition at line 89 of file PhotosParticle.h.

static const int PhotosParticle::MUON_NEUTRINO = 14

staticinherited

muon neutrino particle

Definition at line 92 of file PhotosParticle.h.

static const int PhotosParticle::MUON_PLUS = -13

staticinherited

muon+ particle

Definition at line 86 of file PhotosParticle.h.

static const int PhotosParticle::OTHER = 0

staticinherited

All other particle types

Definition at line 122 of file PhotosParticle.h.

static const int PhotosParticle::POSITRON = -11

staticinherited

e+ particle

Definition at line 98 of file PhotosParticle.h.

static const int PhotosParticle::STABLE =1

staticinherited

Stable particle status

Definition at line 29 of file PhotosParticle.h.

Referenced by Photospp::PhotosEvent::filterParticles(), and Photospp::PhotosHepMCParticle::setMothers().

static const int PhotosParticle::TAU_ANTINEUTRINO = -16

staticinherited

tau antineutrino particle

Definition at line 83 of file PhotosParticle.h.

static const int PhotosParticle::TAU_MINUS = 15

staticinherited

tau- particle

Definition at line 77 of file PhotosParticle.h.

static const int PhotosParticle::TAU_NEUTRINO = 16

staticinherited

tau neutrino particle

Definition at line 80 of file PhotosParticle.h.

static const int PhotosParticle::TAU_PLUS = -15

staticinherited

tau+ particle

Definition at line 74 of file PhotosParticle.h.

static const int PhotosParticle::UP = 2

staticinherited

up quark

Definition at line 110 of file PhotosParticle.h.

static const int PhotosParticle::W_MINUS = -24

staticinherited

W- particle

Definition at line 68 of file PhotosParticle.h.

static const int PhotosParticle::W_PLUS = 24

staticinherited

W+ particle

Definition at line 65 of file PhotosParticle.h.

static const int PhotosParticle::X_AXIS =1

staticinherited

X Axis

Definition at line 38 of file PhotosParticle.h.

Referenced by Photospp::PhotosParticle::boostFromRestFrame(), Photospp::PhotosParticle::boostToRestFrame(), Photospp::PhotosParticle::getP(), and Photospp::PhotosParticle::setP().

static const int PhotosParticle::Y_AXIS =2

staticinherited

Y Axis

Definition at line 41 of file PhotosParticle.h.

Referenced by Photospp::PhotosParticle::boostFromRestFrame(), Photospp::PhotosParticle::boostToRestFrame(), Photospp::PhotosParticle::getP(), and Photospp::PhotosParticle::setP().

static const int PhotosParticle::Z0 = 23

staticinherited

Z0 particle

Definition at line 47 of file PhotosParticle.h.

static const int PhotosParticle::Z_AXIS =3

staticinherited

Z Axis

Definition at line 44 of file PhotosParticle.h.

Referenced by Photospp::PhotosParticle::getP(), and Photospp::PhotosParticle::setP().

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

Public Member Functions

Static Public Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation