EDepSim::TrajectoryMap Class Reference

#include <EDepSimTrajectoryMap.hh>

Public Member Functions
	~TrajectoryMap ()

Static Public Member Functions
static G4VTrajectory *	Get (int trackId)
	Provide a map between the track id and the trajectory object. More...
static void	Add (G4VTrajectory *traj)
	Add a trajectory to the map. More...
static void	Clear ()
static int	FindPrimaryId (int trackId)

Private Member Functions
	TrajectoryMap ()

Static Private Attributes
static std::map< int, G4VTrajectory * >	fMap

Detailed Description

Definition at line 15 of file EDepSimTrajectoryMap.hh.

Constructor & Destructor Documentation

EDepSim::TrajectoryMap::~TrajectoryMap ( )

inline

Definition at line 17 of file EDepSimTrajectoryMap.hh.

{}

EDepSim::TrajectoryMap::TrajectoryMap ( )

inlineprivate

The constructor is private so that it can only be created using the static get method.

Definition at line 42 of file EDepSimTrajectoryMap.hh.

{}

Member Function Documentation

void EDepSim::TrajectoryMap::Add ( G4VTrajectory *  traj )

static

Add a trajectory to the map.

Definition at line 17 of file EDepSimTrajectoryMap.cc.

                                                  {
        int trackId = traj->GetTrackID();
        fMap[trackId] = traj;
}

void EDepSim::TrajectoryMap::Clear ( )

static

Clear the trajectory map. This must be done in the EDepSim::UserEventAction::BeginOfEventAction() method.

Definition at line 13 of file EDepSimTrajectoryMap.cc.

                                 {
    fMap.clear();
}

int EDepSim::TrajectoryMap::FindPrimaryId ( int  trackId )

static

Find the primary track ID for the current track. This is the primary that is the ultimate parent of the current track.

Definition at line 22 of file EDepSimTrajectoryMap.cc.

                                                   {
    int currentId = trackId;
    int parentId = trackId;
    int loopCount=0;
    for (loopCount=0;loopCount<10000;++loopCount) {
        G4VTrajectory* t = Get(currentId);
        if (!t) break;
        parentId = t->GetParentID();
        // Check to see the search loop should terminate.
        G4VTrajectory* p = Get(parentId);
        // There is no parent so break so this is a primary trajectory.
        if (!p) break;
        // Decay products are primary trajectories since they should be
        // independently reconstructed
        EDepSim::Trajectory * edepTraj = dynamic_cast<EDepSim::Trajectory*>(t);
        if (!edepTraj) EDepSimThrow("Invalid Trajectory");
        if (edepTraj->GetProcessName() == "Decay") break;
        // A parent ID of zero means that this particle is a primary particle,
        // so that makes it a primary trajectory too.
        if (parentId == 0) break;
        currentId = parentId;
    }
    if (loopCount>9999) {
        EDepSimLog("Infinite Loop in EDepSim::TrajectoryMap::FindPrimaryId(): "
                 << "Track Id: " << trackId);
    }
    
    return currentId;
}

G4VTrajectory * EDepSim::TrajectoryMap::Get ( int  trackId )

static

Provide a map between the track id and the trajectory object.

Definition at line 52 of file EDepSimTrajectoryMap.cc.

                                                    {
    std::map<int,G4VTrajectory*>::iterator t = fMap.find(trackId);
    if (t == fMap.end()) {
        return NULL;
    }
    return t->second;
}

Member Data Documentation

std::map< int, G4VTrajectory * > EDepSim::TrajectoryMap::fMap

staticprivate

A map to the trajectories information indexed the the track id. Be careful since the trajectory information is owned by the event, so if you try to use this after a trajectory has been deleted... bad things will happen.

Definition at line 38 of file EDepSimTrajectoryMap.hh.

---

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

• edep-sim/src/EDepSimTrajectoryMap.hh
• edep-sim/src/EDepSimTrajectoryMap.cc