Ex02MCStack Class Reference

Implementation of the TVirtualMCStack interface. More...

#include <Ex02MCStack.h>

Inheritance diagram for Ex02MCStack:

Public Member Functions
	Ex02MCStack (Int_t size)
	Ex02MCStack ()
virtual	~Ex02MCStack ()
virtual void	PushTrack (Int_t toBeDone, Int_t parent, Int_t pdg, Double_t px, Double_t py, Double_t pz, Double_t e, Double_t vx, Double_t vy, Double_t vz, Double_t tof, Double_t polx, Double_t poly, Double_t polz, TMCProcess mech, Int_t &ntr, Double_t weight, Int_t is)
virtual TParticle *	PopNextTrack (Int_t &track)
virtual TParticle *	PopPrimaryForTracking (Int_t i)
virtual void	Print (Option_t *option="") const
void	Reset ()
virtual void	SetCurrentTrack (Int_t track)
void	SetObjectNumber ()
virtual Int_t	GetNtrack () const
virtual Int_t	GetNprimary () const
virtual TParticle *	GetCurrentTrack () const
virtual Int_t	GetCurrentTrackNumber () const
virtual Int_t	GetCurrentParentTrackNumber () const
Ex02Particle *	GetParticle (Int_t id) const

Private Attributes
std::stack< Ex02Particle * >	fStack
	The stack of particles (transient)
TObjArray *	fParticles
	The array of particle (persistent)
Int_t	fCurrentTrack
	The current track number.
Int_t	fNPrimary
	The number of primaries.
Int_t	fObjectNumber
	The Root object number counter.

Detailed Description

Implementation of the TVirtualMCStack interface.

Date

05/04/2002

Author

I. Hrivnacova; IPN, Orsay

Definition at line 32 of file Ex02MCStack.h.

Constructor & Destructor Documentation

Ex02MCStack() [1/2]

Ex02MCStack::Ex02MCStack

(

Int_t

size

)

Standard constructor

Parameters

size	The stack size

Definition at line 33 of file Ex02MCStack.cxx.

  : fParticles(0), fCurrentTrack(-1), fNPrimary(0), fObjectNumber(0)
{
  /// Standard constructor
  /// \param size  The stack size
 
  fParticles = new TObjArray(size);
}

Ex02MCStack() [2/2]

Ex02MCStack::Ex02MCStack

(

)

Default constructor

Definition at line 43 of file Ex02MCStack.cxx.

  : fParticles(0), fCurrentTrack(-1), fNPrimary(0), fObjectNumber(0)
{
  /// Default constructor
}

~Ex02MCStack()

Ex02MCStack::~Ex02MCStack ( )

virtual

Destructor

Definition at line 50 of file Ex02MCStack.cxx.

{
  /// Destructor
 
  if (fParticles) fParticles->Delete();
  delete fParticles;
}

Member Function Documentation

PushTrack()

void Ex02MCStack::PushTrack ( Int_t toBeDone, Int_t parent, Int_t pdg, Double_t px, Double_t py, Double_t pz, Double_t e, Double_t vx, Double_t vy, Double_t vz, Double_t tof, Double_t polx, Double_t poly, Double_t polz, TMCProcess mech, Int_t & ntr, Double_t weight, Int_t is )

virtual

Create a new particle and push into stack; adds it to the particles array (fParticles) and if not done to the stack (fStack).

Parameters

toBeDone	1 if particles should go to tracking, 0 otherwise
parent	number of the parent track, -1 if track is primary
pdg	PDG encoding
px	particle momentum - x component [GeV/c]
py	particle momentum - y component [GeV/c]
pz	particle momentum - z component [GeV/c]
e	total energy [GeV]
vx	position - x component [cm]
vy	position - y component [cm]
vz	position - z component [cm]
tof	time of flight [s]
polx	polarization - x component
poly	polarization - y component
polz	polarization - z component
mech	creator process VMC code
ntr	track number (is filled by the stack
weight	particle weight
is	generation status code

Definition at line 63 of file Ex02MCStack.cxx.

{
  /// Create a new particle and push into stack;
  /// adds it to the particles array (fParticles) and if not done to the
  /// stack (fStack).
  /// \param toBeDone  1 if particles should go to tracking, 0 otherwise
  /// \param parent    number of the parent track, -1 if track is primary
  /// \param pdg       PDG encoding
  /// \param px        particle momentum - x component [GeV/c]
  /// \param py        particle momentum - y component [GeV/c]
  /// \param pz        particle momentum - z component [GeV/c]
  /// \param e         total energy [GeV]
  /// \param vx        position - x component [cm]
  /// \param vy        position - y component  [cm]
  /// \param vz        position - z component  [cm]
  /// \param tof       time of flight [s]
  /// \param polx      polarization - x component
  /// \param poly      polarization - y component
  /// \param polz      polarization - z component
  /// \param mech      creator process VMC code
  /// \param ntr       track number (is filled by the stack
  /// \param weight    particle weight
  /// \param is        generation status code
 
  const Int_t kFirstDaughter = -1;
  const Int_t kLastDaughter = -1;
 
  TParticle* particleDef = new TParticle(pdg, is, parent, -1, kFirstDaughter,
    kLastDaughter, px, py, pz, e, vx, vy, vz, tof);
 
  particleDef->SetPolarisation(polx, poly, polz);
  particleDef->SetWeight(weight);
  particleDef->SetUniqueID(mech);
 
  Ex02Particle* mother = 0;
  if (parent >= 0)
    mother = GetParticle(parent);
  else
    fNPrimary++;
 
  Ex02Particle* particle = new Ex02Particle(GetNtrack(), particleDef, mother);
  if (mother) mother->AddDaughter(particle);
  fParticles->Add(particle);
 
  if (toBeDone) fStack.push(particle);
 
  ntr = GetNtrack() - 1;
}

PopNextTrack()

TParticle * Ex02MCStack::PopNextTrack ( Int_t & track )

virtual

Get next particle for tracking from the stack.

Returns

The popped particle object

Parameters

track

The index of the popped track

Definition at line 116 of file Ex02MCStack.cxx.

{
  /// Get next particle for tracking from the stack.
  /// \return       The popped particle object
  /// \param track  The index of the popped track
 
  itrack = -1;
  if (fStack.empty()) return 0;
 
  Ex02Particle* particle = fStack.top();
  fStack.pop();
 
  if (!particle) return 0;
 
  itrack = particle->GetID();
  fCurrentTrack = itrack;
 
  return particle->GetParticle();
}

PopPrimaryForTracking()

TParticle * Ex02MCStack::PopPrimaryForTracking ( Int_t i )

virtual

Return i -th particle in fParticles.

Returns

The popped primary particle object

Parameters

i	The index of primary particle to be popped

Definition at line 137 of file Ex02MCStack.cxx.

{
  /// Return \em i -th particle in fParticles.
  /// \return   The popped primary particle object
  /// \param i  The index of primary particle to be popped
 
  if (i < 0 || i >= fNPrimary)
    Fatal("GetPrimaryForTracking", "Index out of range");
 
  return ((Ex02Particle*)fParticles->At(i))->GetParticle();
}

Print()

void Ex02MCStack::Print ( Option_t * option = "" ) const

virtual

Print info for all particles.

Definition at line 150 of file Ex02MCStack.cxx.

{
  /// Print info for all particles.
 
  cout << "Ex02MCStack Info  " << endl;
  cout << "Total number of particles:   " << GetNtrack() << endl;
  cout << "Number of primary particles: " << GetNprimary() << endl;
 
  for (Int_t i = 0; i < GetNtrack(); i++) {
    GetParticle(i)->Print();
    // GetParticle(i)->PrintDaughters();
  }
}

Reset()

void Ex02MCStack::Reset

(

)

Delete contained particles, reset particles array and stack.

Definition at line 165 of file Ex02MCStack.cxx.

{
  /// Delete contained particles, reset particles array and stack.
 
  // reset fStack
  fCurrentTrack = -1;
  fNPrimary = 0;
  // fParticles->Delete();
  fParticles->Clear();
 
  // Restore Object count
  // To save space in the table keeping track of all referenced objects
  // we assume that our events do not address each other. We reset the
  // object count to what it was at the beginning of the event
  TProcessID::SetObjectCount(fObjectNumber);
}

SetCurrentTrack()

void Ex02MCStack::SetCurrentTrack ( Int_t track )

virtual

Set the current track number to a given value.

Parameters

track

The current track number

Definition at line 183 of file Ex02MCStack.cxx.

{
  /// Set the current track number to a given value.
  /// \param  track The current track number
 
  fCurrentTrack = track;
}

SetObjectNumber()

void Ex02MCStack::SetObjectNumber

(

)

Set the current value of Root object counter into fObjectNumber. Tis value will be restored in Reset.

Definition at line 192 of file Ex02MCStack.cxx.

{
  /// Set the current value of Root object counter into fObjectNumber.
  /// Tis value will be restored in Reset.
 
  fObjectNumber = TProcessID::GetObjectCount();
}

GetNtrack()

Int_t Ex02MCStack::GetNtrack ( ) const

virtual

Returns

The total number of all tracks.

Definition at line 201 of file Ex02MCStack.cxx.

{
  /// \return  The total number of all tracks.
 
  return fParticles->GetEntriesFast();
}

GetNprimary()

Int_t Ex02MCStack::GetNprimary ( ) const

virtual

Returns

The total number of primary tracks.

Definition at line 209 of file Ex02MCStack.cxx.

{
  /// \return  The total number of primary tracks.
 
  return fNPrimary;
}

GetCurrentTrack()

TParticle * Ex02MCStack::GetCurrentTrack ( ) const

virtual

Returns

The current track particle

Definition at line 217 of file Ex02MCStack.cxx.

{
  /// \return  The current track particle
 
  Ex02Particle* current = GetParticle(fCurrentTrack);
 
  if (current)
    return current->GetParticle();
  else
    return 0;
}

GetCurrentTrackNumber()

Int_t Ex02MCStack::GetCurrentTrackNumber ( ) const

virtual

Returns

The current track number

Definition at line 230 of file Ex02MCStack.cxx.

{
  /// \return  The current track number
 
  return fCurrentTrack;
}

GetCurrentParentTrackNumber()

Int_t Ex02MCStack::GetCurrentParentTrackNumber ( ) const

virtual

Returns

The current track parent ID.

Definition at line 238 of file Ex02MCStack.cxx.

{
  /// \return  The current track parent ID.
 
  Ex02Particle* current = GetParticle(fCurrentTrack);
 
  if (!current) return -1;
 
  Ex02Particle* mother = current->GetMother();
 
  if (!mother) return -1;
 
  return mother->GetID();
}

GetParticle()

Ex02Particle * Ex02MCStack::GetParticle

(

Int_t

id

)

const

Returns

The id -th particle in fParticles

Parameters

id	The index of the particle to be returned

Definition at line 254 of file Ex02MCStack.cxx.

{
  /// \return   The \em id -th particle in fParticles
  /// \param id The index of the particle to be returned
 
  if (id < 0 || id >= fParticles->GetEntriesFast())
    Fatal("GetParticle", "Index out of range");
 
  return (Ex02Particle*)fParticles->At(id);
}

Member Data Documentation

fStack

std::stack<Ex02Particle*> Ex02MCStack::fStack

private

The stack of particles (transient)

Definition at line 63 of file Ex02MCStack.h.

fParticles

TObjArray* Ex02MCStack::fParticles

private

The array of particle (persistent)

Definition at line 64 of file Ex02MCStack.h.

fCurrentTrack

Int_t Ex02MCStack::fCurrentTrack

private

The current track number.

Definition at line 65 of file Ex02MCStack.h.

fNPrimary

Int_t Ex02MCStack::fNPrimary

private

The number of primaries.

Definition at line 66 of file Ex02MCStack.h.

fObjectNumber

Int_t Ex02MCStack::fObjectNumber

private

The Root object number counter.

Definition at line 67 of file Ex02MCStack.h.

---

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

• examples/E02/include/Ex02MCStack.h
• examples/E02/src/Ex02MCStack.cxx