PndStack Class Reference

#include <PndStack.h>

Inheritance diagram for PndStack:

Public Member Functions
	PndStack (Int_t size=100)
virtual	~PndStack ()
virtual void	PushTrack (Int_t toBeDone, Int_t parentID, Int_t pdgCode, Double_t px, Double_t py, Double_t pz, Double_t e, Double_t vx, Double_t vy, Double_t vz, Double_t time, Double_t polx, Double_t poly, Double_t polz, TMCProcess proc, Int_t &ntr, Double_t weight, Int_t is)
virtual void	PushTrack (Int_t toBeDone, Int_t parentID, Int_t pdgCode, Double_t px, Double_t py, Double_t pz, Double_t e, Double_t vx, Double_t vy, Double_t vz, Double_t time, Double_t polx, Double_t poly, Double_t polz, TMCProcess proc, Int_t &ntr, Double_t weight, Int_t is, Int_t secondParentId)
virtual TParticle *	PopNextTrack (Int_t &iTrack)
virtual TParticle *	PopPrimaryForTracking (Int_t iPrim)
virtual void	SetCurrentTrack (Int_t iTrack)
virtual Int_t	GetNtrack () const
virtual Int_t	GetNprimary () const
virtual TParticle *	GetCurrentTrack () const
virtual Int_t	GetCurrentTrackNumber () const
virtual Int_t	GetCurrentParentTrackNumber () const
virtual void	AddParticle (TParticle *part)
virtual void	FillTrackArray ()
virtual void	UpdateTrackIndex (TRefArray *detArray)
virtual void	Reset ()
virtual void	Register ()
virtual void	Print (Int_t iVerbose=0) const
void	StoreSecondaries (Bool_t choice=kTRUE)
void	SetMinPoints (Int_t min)
void	SetEnergyCut (Double_t eMin)
void	StoreMothers (Bool_t choice=kTRUE)
void	SetPersistence (Bool_t choice)
void	AddPoint (DetectorId iDet)
void	AddPoint (DetectorId iDet, Int_t iTrack)
TParticle *	GetParticle (Int_t trackId) const
TClonesArray *	GetListOfParticles ()

Private Member Functions
void	SelectTracks ()
void	SetGeneratorFlags (Int_t myid)
	PndStack (const PndStack &L)
PndStack &	operator= (const PndStack &)

Private Attributes
std::stack< TParticle * >	fStack
TClonesArray *	fParticles
TClonesArray *	fTracks
std::map< Int_t, Bool_t >	fStoreMap
std::map< Int_t, Bool_t >::iterator	fStoreIter
std::map< Int_t, Int_t >	fIndexMap
std::map< Int_t, Int_t >::iterator	fIndexIter
std::map< std::pair< Int_t, Int_t >, Int_t >	fPointsMap
Int_t	fCurrentTrack
Int_t	fNPrimaries
	Index of current track. More...
Int_t	fNParticles
	Number of primary particles. More...
Int_t	fNTracks
	Number of entries in fParticles. More...
Int_t	fIndex
	Number of entries in fTracks. More...
Bool_t	fStoreSecondaries
	Used for merging. More...
Int_t	fMinPoints
Double32_t	fEnergyCut
Bool_t	fStoreMothers
Bool_t	fPersistence

Detailed Description

Definition at line 45 of file PndStack.h.

Constructor & Destructor Documentation

PndStack::PndStack

(

Int_t

size = 100

)

Default constructor param size Estimated track number

Definition at line 27 of file PndStack.cxx.

References fIndexMap, fParticles, fPointsMap, fStoreMap, and fTracks.

                            :
fStack(), fParticles(), fTracks(), fStoreMap(), fStoreIter(), fIndexMap(), fIndexIter(), fPointsMap(), fCurrentTrack(-1), fNPrimaries(0), fNParticles(0), fNTracks(0), fIndex(0), fStoreSecondaries(kTRUE), fMinPoints(1), fEnergyCut(0), fStoreMothers(kTRUE), fPersistence(kTRUE)
{
  fStoreMap.clear();
  fIndexMap.clear();
  fPointsMap.clear();
  fParticles = new TClonesArray("TParticle", size);
  fTracks    = new TClonesArray("PndMCTrack", size);
}

PndStack::~PndStack ( )

virtual

Destructor

Definition at line 42 of file PndStack.cxx.

References fParticles, and fTracks.

                    {
  if (fParticles) {
    fParticles->Delete();
    delete fParticles;
  }
  if (fTracks) {
    fTracks->Delete();
    delete fTracks;
  }
}

PndStack::PndStack ( const PndStack &  L )

private

Member Function Documentation

void PndStack::AddParticle ( TParticle *  part )

virtual

Add a TParticle to the fParticles array

Definition at line 186 of file PndStack.cxx.

References fIndex, and fParticles.

                                             {
  TClonesArray& array = *fParticles;
  TParticle* newPart = new(array[fIndex]) TParticle(*oldPart);
  newPart->SetWeight(oldPart->GetWeight());
  newPart->SetUniqueID(oldPart->GetUniqueID());
  fIndex++;
}

void PndStack::AddPoint

(

DetectorId

iDet

)

Increment number of points for the current track in a given detector

Parameters

iDet	Detector unique identifier

Definition at line 408 of file PndStack.cxx.

References a, fCurrentTrack, and fPointsMap.

Referenced by PndRich::ProcessHits(), PndFts2::ProcessHits(), PndFts::ProcessHits(), PndFtof::ProcessHits(), PndEmc::ProcessHits(), PndGemDetector::ProcessHits(), PndSciT::ProcessHits(), PndStt::ProcessHits(), PndDrc::ProcessHits(), PndSdsDetector::ProcessHits(), PndDsk::ProcessHitsCerenkov(), PndDsk::ProcessHitsCerenkov_FLG(), PndMdt::ProcessHitsFast(), PndDsk::ProcessHitsParticle(), and PndMdt::ProcessHitsRoot().

                                        {
  Int_t iDet = detId;
  pair<Int_t, Int_t> a(fCurrentTrack, iDet);
  if ( fPointsMap.find(a) == fPointsMap.end() ) fPointsMap[a] = 1;
  else fPointsMap[a]++;
}

void PndStack::AddPoint	(	DetectorId	iDet,
		Int_t	iTrack
	)

Increment number of points for an arbitrary track in a given detector

Parameters

iDet	Detector unique identifier
iTrack	Track number

Definition at line 419 of file PndStack.cxx.

References a, and fPointsMap.

                                                      {
  if ( iTrack < 0 ) return;
  Int_t iDet = detId;
  pair<Int_t, Int_t> a(iTrack, iDet);
  if ( fPointsMap.find(a) == fPointsMap.end() ) fPointsMap[a] = 1;
  else fPointsMap[a]++;
}

void PndStack::FillTrackArray ( )

virtual

Fill the MCTrack output array, applying filter criteria

Definition at line 198 of file PndStack.cxx.

References a, Bool_t, DEBUG, fIndexMap, fNParticles, fNTracks, fPointsMap, fStoreIter, fStoreMap, GetParticle(), kHYP, kRICH, Print(), SelectTracks(), SetGeneratorFlags(), PndMCTrack::SetNPoints(), store(), and track.

                              {

  LOG(DEBUG) << "PndStack: Filling MCTrack array..."  ;

  // --> Reset index map and number of output tracks
  fIndexMap.clear();
  fNTracks = 0;

  // --> Check tracks for selection criteria
  SelectTracks();

  // --> Loop over fParticles array and copy selected tracks
  for (Int_t iPart=0; iPart<fNParticles; iPart++) {

    fStoreIter = fStoreMap.find(iPart);
    if (fStoreIter == fStoreMap.end() ) {
      LOG(ERROR) << "PndStack: Particle  "<<iPart<<"  not found in storage map!"  ;
      Fatal("PndStack::FillTrackArray","Particle not found in storage map.");
    }
    Bool_t store = (*fStoreIter).second;

    if (store) {
      PndMCTrack* track =
        new( (*fTracks)[fNTracks]) PndMCTrack(GetParticle(iPart));
      fIndexMap[iPart] = fNTracks;
      // --> Set the number of points in the detectors for this track
      for (Int_t iDet=kRICH; iDet<=kHYP; iDet++) {
        pair<Int_t, Int_t> a(iPart, iDet);
        track->SetNPoints(iDet, fPointsMap[a]);
      }

          SetGeneratorFlags(iPart);

      fNTracks++;

    }else{
      fIndexMap[iPart] = -2;
    }

  }

  // --> Map index for primary mothers
  fIndexMap[-1] = -1;

  // --> Screen output
  Print(0);

}

Int_t PndStack::GetCurrentParentTrackNumber ( ) const

virtual

Get the track number of the parent of the current track Declared in TVirtualMCStack

Definition at line 432 of file PndStack.cxx.

References GetCurrentTrack().

                                                  { 
  TParticle* currentPart = GetCurrentTrack();
  if ( currentPart ) return currentPart->GetFirstMother();
  else               return -1;
}

TParticle * PndStack::GetCurrentTrack ( ) const

virtual

Get the current track's particle Declared in TVirtualMCStack

Definition at line 173 of file PndStack.cxx.

References fCurrentTrack, and GetParticle().

Referenced by GetCurrentParentTrackNumber().

                                           {
  TParticle* currentPart = GetParticle(fCurrentTrack);
  if ( ! currentPart) {
    LOG(WARNING) << "PndStack: Current track not found in stack!"  ;
    Warning("PndStack::GetCurrentTrack", "Track not found in stack");
  }
  return currentPart;
}

virtual Int_t PndStack::GetCurrentTrackNumber ( ) const

inlinevirtual

Get the number of the current track Declared in TVirtualMCStack

Definition at line 134 of file PndStack.h.

References fCurrentTrack.

{ return fCurrentTrack; }

TClonesArray* PndStack::GetListOfParticles ( )

inline

Definition at line 192 of file PndStack.h.

References fParticles.

{ return fParticles; }

virtual Int_t PndStack::GetNprimary ( ) const

inlinevirtual

Get number of primary tracks Declared in TVirtualMCStack

Definition at line 122 of file PndStack.h.

References fNPrimaries.

{ return fNPrimaries; }

virtual Int_t PndStack::GetNtrack ( ) const

inlinevirtual

Get total number of tracks Declared in TVirtualMCStack

Definition at line 116 of file PndStack.h.

References fNParticles.

Referenced by PndFastSim::Exec().

{ return fNParticles; }

TParticle * PndStack::GetParticle

(

Int_t

trackId

)

const

Accessors

Definition at line 442 of file PndStack.cxx.

References fNParticles, and fParticles.

Referenced by PndFastSim::Exec(), FillTrackArray(), GetCurrentTrack(), and SelectTracks().

                                                    {
  if (trackID < 0 || trackID >= fNParticles) {
    LOG(ERROR) << "PndStack: Particle index out of range." << trackID  ;
    Fatal("PndStack::GetParticle", "Index out of range");
  }
  return (TParticle*)fParticles->At(trackID);
}

PndStack& PndStack::operator= ( const PndStack &  )

inlineprivate

Definition at line 248 of file PndStack.h.

{return *this;};

TParticle * PndStack::PopNextTrack ( Int_t &  iTrack )

virtual

Get next particle for tracking from the stack. Declared in TVirtualMCStack

Parameters

iTrack

index of popped track (return)

Returns

Pointer to the TParticle of the track

Definition at line 116 of file PndStack.cxx.

References fCurrentTrack, and fStack.

                                               {

  // If end of stack: Return empty pointer
  if (fStack.empty()) {
    iTrack = -1;
    return NULL;
  }

  // If not, get next particle from stack
  TParticle* thisParticle = fStack.top();
  fStack.pop();

  if ( !thisParticle) {
    iTrack = 0;
    return NULL;
  }

  fCurrentTrack = thisParticle->GetStatusCode();
  iTrack = fCurrentTrack;

  return thisParticle;

}

TParticle * PndStack::PopPrimaryForTracking ( Int_t  iPrim )

virtual

Get primary particle by index for tracking from stack Declared in TVirtualMCStack

Parameters

iPrim

index of primary particle

Returns

Pointer to the TParticle of the track

Definition at line 144 of file PndStack.cxx.

References fNPrimaries, fParticles, and kDoneBit.

                                                      {

  // Get the iPrimth particle from the fStack TClonesArray. This
  // should be a primary (if the index is correct).

  // Test for index
  if (iPrim < 0 || iPrim >= fNPrimaries) {
    LOG(ERROR) << "PndStack: Primary index out of range! " << iPrim  ;
    Fatal("PndStack::PopPrimaryForTracking", "Index out of range");
  }

  // Return the iPrim-th TParticle from the fParticle array. This should be
  // a primary.
  TParticle* part = (TParticle*)fParticles->At(iPrim);

  if ( ! (part->GetMother(0) < 0) ) {
    LOG(ERROR) << "PndStack:: Not a primary track! ,  " <<iPrim  ;
    Fatal("PndStack::PopPrimaryForTracking", "Not a primary track");
  }

  if(!part->TestBit(kDoneBit)) return NULL;
  else return part;

}

void PndStack::Print ( Int_t  iVerbose = 0 ) const

virtual

Output to screen

Parameters

iVerbose

0=events summary, 1=track info

Definition at line 393 of file PndStack.cxx.

References DEBUG, fNParticles, fNPrimaries, fNTracks, and fTracks.

Referenced by FillTrackArray().

                                         {
  LOG(DEBUG) << "  PndStack: Number of primaries  = " <<fNPrimaries
             << "\n  Total number of particles  = " << fNParticles
             << "\n  Number of tracks in output = " << fNTracks  ;

  if (iVerbose) {
    for (Int_t iTrack=0; iTrack<fNTracks; iTrack++)
      ((PndMCTrack*) fTracks->At(iTrack))->Print(iTrack);
  }
}

void PndStack::PushTrack ( Int_t  toBeDone, Int_t  parentID, Int_t  pdgCode, Double_t  px, Double_t  py, Double_t  pz, Double_t  e, Double_t  vx, Double_t  vy, Double_t  vz, Double_t  time, Double_t  polx, Double_t  poly, Double_t  polz, TMCProcess  proc, Int_t &  ntr, Double_t  weight, Int_t  is  )

virtual

Add a TParticle to the stack. Declared in TVirtualMCStack

Parameters

toBeDone	Flag for tracking
parentID	Index of mother particle
pdgCode	Particle type (PDG encoding)
px,py,pz	Momentum components at start vertex [GeV]
e	Total energy at start vertex [GeV]
vx,vy,vz	Coordinates of start vertex [cm]
time	Start time of track [s]
polx,poly,polz	Polarisation vector
proc	Production mechanism (VMC encoding)
ntr	Track number (filled by the stack)
weight	Particle weight
is	Generation status code (whatever that means)

Definition at line 57 of file PndStack.cxx.

Referenced by PndFastSim::Exec().

                                                    {

        PushTrack( toBeDone, parentId, pdgCode,
                                  px,  py,  pz,
                                  e,  vx,  vy,  vz,
                                  time,  polx,  poly,
                                  polz, proc, ntr,
                                  weight, is, -1);
}

void PndStack::PushTrack ( Int_t  toBeDone, Int_t  parentID, Int_t  pdgCode, Double_t  px, Double_t  py, Double_t  pz, Double_t  e, Double_t  vx, Double_t  vy, Double_t  vz, Double_t  time, Double_t  polx, Double_t  poly, Double_t  polz, TMCProcess  proc, Int_t &  ntr, Double_t  weight, Int_t  is, Int_t  secondParentId  )

virtual

Definition at line 73 of file PndStack.cxx.

References fNParticles, fNPrimaries, fParticles, fStack, kDoneBit, and particle.

                                                                         {

  (void)is; // To remove "unused" warnings
  // --> Get TParticle array
  TClonesArray& partArray = *fParticles;

  // --> Create new TParticle and add it to the TParticle array
  Int_t trackId = fNParticles;
  Int_t nPoints = 0;
  Int_t daughter1Id = -1;
  Int_t daughter2Id = -1;
  TParticle* particle =
    new(partArray[fNParticles++]) TParticle(pdgCode, trackId, parentId,
                                            nPoints, daughter1Id,
                                            daughter2Id, px, py, pz, e,
                                            vx, vy, vz, time);
  particle->SetLastMother(secondparentID);
  particle->SetPolarisation(polx, poly, polz);
  particle->SetWeight(weight);
  particle->SetUniqueID(proc);

  // --> Increment counter
  if (parentId < 0) fNPrimaries++;

  // --> Set argument variable
  ntr = trackId;

  // --> Push particle on the stack if toBeDone is set
   if (toBeDone == 1) {
      particle->SetBit(kDoneBit);
      fStack.push(particle);
   }
}

void PndStack::Register ( )

virtual

Register the MCTrack array to the Root Manager

Definition at line 385 of file PndStack.cxx.

References fPersistence, and fTracks.

                        {
  FairRootManager::Instance()->Register("MCTrack", "Stack", fTracks, fPersistence);
}

void PndStack::Reset ( )

virtual

Resets arrays and stack and deletes particles and tracks

Definition at line 371 of file PndStack.cxx.

References fCurrentTrack, fIndex, fNParticles, fNPrimaries, fNTracks, fParticles, fPointsMap, fStack, and fTracks.

                     {
  fIndex = 0;
  fCurrentTrack = -1;
  fNPrimaries = fNParticles = fNTracks = 0;
  while (! fStack.empty() ) fStack.pop();
  fParticles->Clear();
  fTracks->Clear();
  fPointsMap.clear();
}

void PndStack::SelectTracks ( )

private

Mark tracks for output using selection criteria

Definition at line 454 of file PndStack.cxx.

References a, Bool_t, Double_t, energy, fEnergyCut, fMinPoints, fNParticles, fPointsMap, fStoreMap, fStoreMothers, fStoreSecondaries, GetParticle(), i, kHYP, kRICH, p, and store().

Referenced by FillTrackArray().

                            {

  // --> Clear storage map
  fStoreMap.clear();

  // --> Check particles in the fParticle array
  for (Int_t i=0; i<fNParticles; i++) {

    TParticle* thisPart = GetParticle(i);
    Bool_t store = kTRUE;

    // --> Get track parameters
    Int_t iMother   = thisPart->GetMother(0);
    TLorentzVector p;
    thisPart->Momentum(p);
    Double_t energy = p.E();
    Double_t mass   = p.M();
//    Double_t mass   = thisPart->GetMass();// Why?? Mass (given by generator) is inside by Lorentzvector!!! I dont care about PSG mass!
    Double_t eKin = energy - mass;
    if(eKin < 0.0) eKin=0.0; // sometimes due to different PDG masses between ROOT and G4!!!!!!
    // --> Calculate number of points
    Int_t nPoints = 0;
    for (Int_t iDet=kRICH; iDet<=kHYP; iDet++) {
      pair<Int_t, Int_t> a(i, iDet);
      if ( fPointsMap.find(a) != fPointsMap.end() )
        nPoints += fPointsMap[a];
    }

    // --> Check for cuts (store primaries in any case)
    if (iMother < 0)            store = kTRUE;
    else {
      if (!fStoreSecondaries)   store = kFALSE;
      if (nPoints < fMinPoints) store = kFALSE;
      if (eKin < fEnergyCut)    store = kFALSE;
    }

    // --> Set storage flag
    fStoreMap[i] = store;

  }

  // --> If flag is set, flag recursively mothers of selected tracks
  if (fStoreMothers) {
    for (Int_t i=0; i<fNParticles; i++) {
      if (fStoreMap[i]) {
        Int_t iMother = GetParticle(i)->GetMother(0);
        while(iMother >= 0) {
          fStoreMap[iMother] = kTRUE;
          iMother = GetParticle(iMother)->GetMother(0);
        }
      }
    }
  }

}

virtual void PndStack::SetCurrentTrack ( Int_t  iTrack )

inlinevirtual

Set the current track number Declared in TVirtualMCStack

Parameters

iTrack

track number

Definition at line 110 of file PndStack.h.

References fCurrentTrack.

{ fCurrentTrack     = iTrack; }

void PndStack::SetEnergyCut ( Double_t  eMin )

inline

Definition at line 172 of file PndStack.h.

References fEnergyCut.

{ fEnergyCut        = eMin;   }

void PndStack::SetGeneratorFlags ( Int_t  myid )

private

Definition at line 248 of file PndStack.cxx.

References fIndexMap, fParticles, fTracks, PndMCTrack::GetMotherID(), i, m, n, PndMCTrack::SetGeneratorCreated(), and PndMCTrack::SetGeneratorDecayed().

Referenced by FillTrackArray().

{
        if(myid<0) return;

        PndMCTrack* mytrack;
        {
                Int_t myid2=fIndexMap[myid];
                if(myid2<0){
                        LOG(ERROR)<<"=== This should not happen negative index in MAP!!"  ;
                        return;
                }

                mytrack = (PndMCTrack*)fTracks->At(myid2);
        }

        Int_t n;
        Int_t daughters=0, daughtersp=0;
        // fParticles; // TParticle
        n=fParticles->GetEntries();
        for (Int_t i=0; i<n; i++) {
                TParticle* part = (TParticle*)fParticles->At(i);
                Int_t m;
                m=part->GetMother(0);
                if(myid==m){
                        daughters++;
                }else if(m==-1){
                        m=part->GetMother(1);
                        if(myid==m){
                                daughtersp++;
                        }
                }else if(m==-2){
                        // removed should not happen before this is called
                        // and anyway not on the TParticle Level
                        LOG(ERROR)<<"=== Problem!!! part mother -2"  ;
                }
        }

        Int_t mymo1=mytrack->GetMotherID();

        if( ((TParticle*)fParticles->At(myid))->GetMother(0)!=mymo1){
                LOG(ERROR)<<"=== Problem: Mothers != "<< myid  ;
    }
        if(mymo1==-1){
                if(daughters!=0 && daughtersp!=0){
                        LOG(ERROR)<<"=== Problem: particle with index "<<myid<<" has  daughters= "<<daughters<<"  && daughtersp= "<<daughtersp  ;
                }

                mytrack->SetGeneratorCreated();
                if(daughtersp>0) mytrack->SetGeneratorDecayed();
//              cout << myid <<" ("<<mytrack->GetPdgCode()<<"): "<<daughters<<","<<daughtersp<<" ==> " <<mytrack->IsGeneratorCreated()<<" "<<mytrack->IsGeneratorDecayed()<<" "<<mytrack->IsGeneratorLast()<<endl;
        }
}

void PndStack::SetMinPoints ( Int_t  min )

inline

Definition at line 171 of file PndStack.h.

References fMinPoints, and min().

Referenced by Config().

{ fMinPoints        = min;    }

void PndStack::SetPersistence ( Bool_t  choice )

inline

Definition at line 174 of file PndStack.h.

References choice, and fPersistence.

Referenced by Config().

{ fPersistence      = choice; }

void PndStack::StoreMothers ( Bool_t  choice = kTRUE )

inline

Definition at line 173 of file PndStack.h.

References choice, and fStoreMothers.

{ fStoreMothers     = choice; }

void PndStack::StoreSecondaries ( Bool_t  choice = kTRUE )

inline

Modifiers

Definition at line 170 of file PndStack.h.

References choice, and fStoreSecondaries.

Referenced by Config().

{ fStoreSecondaries = choice; }

void PndStack::UpdateTrackIndex ( TRefArray *  detArray )

virtual

Update the track index in the MCTracks and MCPoints

Definition at line 303 of file PndStack.cxx.

References DEBUG, fIndexIter, fIndexMap, fNTracks, fTracks, PndMCTrack::GetMotherID(), PndMCTrack::GetSecondMotherID(), i, point, PndMCTrack::SetMotherID(), PndMCTrack::SetSecondMotherID(), and track.

                                                  {

  LOG(DEBUG) << "PndStack: Updating track indizes..."  ;
  Int_t nColl = 0;

  FairMCEventHeader* header = (FairMCEventHeader*)FairRootManager::Instance()->GetObject("MCEventHeader.");

  // First update mother ID in MCTracks
  for (Int_t i=0; i<fNTracks; i++) {
    PndMCTrack* track = (PndMCTrack*)fTracks->At(i);
    Int_t iMotherOld = track->GetMotherID();
    fIndexIter = fIndexMap.find(iMotherOld);
    if (fIndexIter == fIndexMap.end()) {
      LOG(ERROR) << "PndStack: Particle index  "<<iMotherOld<<" not found in dex map! "  ;
      Fatal("PndStack::UpdateTrackIndex","Particle index not found in map");
    }
    track->SetMotherID( (*fIndexIter).second );
    if(iMotherOld==-1){
        iMotherOld = track->GetSecondMotherID();
        fIndexIter = fIndexMap.find(iMotherOld);
        if (fIndexIter == fIndexMap.end()) {
           LOG(ERROR) << "PndStack: Particle index  "<<iMotherOld<<" not found in dex map! (second mother id)"  ;
           Fatal("PndStack::UpdateTrackIndex","Particle index not found in map");
        }
        track->SetSecondMotherID( (*fIndexIter).second );
    }
  }

  // Now iterate through all active detectors
  TIterator* detIter = detList->MakeIterator();
  detIter->Reset();
  FairDetector* det = NULL;
        while ((det = (FairDetector*) detIter->Next())) {

                // --> Get hit collections from detector
                Int_t iColl = 0;
                TClonesArray* hitArray;
                while ((hitArray = det->GetCollection(iColl++))) {
                        nColl++;
                        Int_t nPoints = hitArray->GetEntriesFast();

                        // --> Update track index for all MCPoints in the collection
                        for (Int_t iPoint = 0; iPoint < nPoints; iPoint++) {
                                FairMCPoint* point = (FairMCPoint*) hitArray->At(iPoint);
                                Int_t iTrack = point->GetTrackID();

                                fIndexIter = fIndexMap.find(iTrack);
                                if (fIndexIter == fIndexMap.end()) {
                                        LOG(ERROR) << "PndStack: Particle index "<<iTrack<<" not found in index map! "  ;
                                        Fatal("PndStack::UpdateTrackIndex",
                                                        "Particle index not found in map");
                                }
                                point->SetTrackID((*fIndexIter).second);
//                              std::cout << "Header->GetEventID() " << header->GetEventID() << std::endl;
                                point->SetLink(FairLink(-1, (header->GetEventID()-1), "MCTrack", (*fIndexIter).second));
                        }

                }   // Collections of this detector
        }     // List of active detectors

  LOG(DEBUG) << "...stack and "<<nColl<<" collections updated."  ;
  delete detIter;
}

Member Data Documentation

Int_t PndStack::fCurrentTrack

private

Some indizes and counters

Definition at line 228 of file PndStack.h.

Referenced by AddPoint(), GetCurrentTrack(), GetCurrentTrackNumber(), PopNextTrack(), Reset(), and SetCurrentTrack().

Double32_t PndStack::fEnergyCut

private

Definition at line 238 of file PndStack.h.

Referenced by SelectTracks(), and SetEnergyCut().

Int_t PndStack::fIndex

private

Number of entries in fTracks.

Definition at line 232 of file PndStack.h.

Referenced by AddParticle(), and Reset().

std::map<Int_t, Int_t>::iterator PndStack::fIndexIter

private

Definition at line 220 of file PndStack.h.

Referenced by UpdateTrackIndex().

std::map<Int_t, Int_t> PndStack::fIndexMap

private

STL map from particle index to track index

Definition at line 219 of file PndStack.h.

Referenced by FillTrackArray(), PndStack(), SetGeneratorFlags(), and UpdateTrackIndex().

Int_t PndStack::fMinPoints

private

Definition at line 237 of file PndStack.h.

Referenced by SelectTracks(), and SetMinPoints().

Int_t PndStack::fNParticles

private

Number of primary particles.

Definition at line 230 of file PndStack.h.

Referenced by FillTrackArray(), GetNtrack(), GetParticle(), Print(), PushTrack(), Reset(), and SelectTracks().

Int_t PndStack::fNPrimaries

private

Index of current track.

Definition at line 229 of file PndStack.h.

Referenced by GetNprimary(), PopPrimaryForTracking(), Print(), PushTrack(), and Reset().

Int_t PndStack::fNTracks

private

Number of entries in fParticles.

Definition at line 231 of file PndStack.h.

Referenced by FillTrackArray(), Print(), Reset(), and UpdateTrackIndex().

TClonesArray* PndStack::fParticles

private

Array of TParticles (contains all TParticles put into or created by the transport

Definition at line 206 of file PndStack.h.

Referenced by AddParticle(), GetListOfParticles(), GetParticle(), PndStack(), PopPrimaryForTracking(), PushTrack(), Reset(), SetGeneratorFlags(), and ~PndStack().

Bool_t PndStack::fPersistence

private

Definition at line 240 of file PndStack.h.

Referenced by Register(), and SetPersistence().

std::map<std::pair<Int_t, Int_t>, Int_t> PndStack::fPointsMap

private

STL map from track index and detector ID to number of MCPoints

Definition at line 224 of file PndStack.h.

Referenced by AddPoint(), FillTrackArray(), PndStack(), Reset(), and SelectTracks().

std::stack<TParticle*> PndStack::fStack

private

STL stack (FILO) used to handle the TParticles for tracking

Definition at line 200 of file PndStack.h.

Referenced by PopNextTrack(), PushTrack(), and Reset().

std::map<Int_t, Bool_t>::iterator PndStack::fStoreIter

private

Definition at line 215 of file PndStack.h.

Referenced by FillTrackArray().

std::map<Int_t, Bool_t> PndStack::fStoreMap

private

STL map from particle index to storage flag

Definition at line 214 of file PndStack.h.

Referenced by FillTrackArray(), PndStack(), and SelectTracks().

Bool_t PndStack::fStoreMothers

private

Definition at line 239 of file PndStack.h.

Referenced by SelectTracks(), and StoreMothers().

Bool_t PndStack::fStoreSecondaries

private

Used for merging.

Variables defining the criteria for output selection

Definition at line 236 of file PndStack.h.

Referenced by SelectTracks(), and StoreSecondaries().

TClonesArray* PndStack::fTracks

private

Array of PndMCTracks containg the tracks written to the output

Definition at line 210 of file PndStack.h.

Referenced by PndStack(), Print(), Register(), Reset(), SetGeneratorFlags(), UpdateTrackIndex(), and ~PndStack().

---

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

• PndStack.h
• PndStack.cxx