G4RichTrajectory Class Reference

#include <G4RichTrajectory.hh>

Inheritance diagram for G4RichTrajectory:

Public Member Functions
	G4RichTrajectory ()
	G4RichTrajectory (const G4Track *aTrack)
	G4RichTrajectory (G4RichTrajectory &)
virtual	~G4RichTrajectory ()
void *	operator new (size_t)
void	operator delete (void *)
int	operator== (const G4RichTrajectory &right) const
void	ShowTrajectory (std::ostream &os=G4cout) const
void	DrawTrajectory () const
void	AppendStep (const G4Step *aStep)
void	MergeTrajectory (G4VTrajectory *secondTrajectory)
int	GetPointEntries () const
G4VTrajectoryPoint *	GetPoint (G4int i) const
virtual const std::map < G4String, G4AttDef > *	GetAttDefs () const
virtual std::vector< G4AttValue > *	CreateAttValues () const
Public Member Functions inherited from G4Trajectory
	G4Trajectory ()
	G4Trajectory (const G4Track *aTrack)
	G4Trajectory (G4Trajectory &)
virtual	~G4Trajectory ()
void *	operator new (size_t)
void	operator delete (void *)
int	operator== (const G4Trajectory &right) const
G4int	GetTrackID () const
G4int	GetParentID () const
G4String	GetParticleName () const
G4double	GetCharge () const
G4int	GetPDGEncoding () const
G4double	GetInitialKineticEnergy () const
G4ThreeVector	GetInitialMomentum () const
G4ParticleDefinition *	GetParticleDefinition ()
Public Member Functions inherited from G4VTrajectory
	G4VTrajectory ()
virtual	~G4VTrajectory ()
G4bool	operator== (const G4VTrajectory &right) const

Detailed Description

Definition at line 65 of file G4RichTrajectory.hh.

Constructor & Destructor Documentation

G4RichTrajectory::G4RichTrajectory

(

)

Definition at line 62 of file G4RichTrajectory.cc.

                                  :
  fpRichPointsContainer(0),
  fpCreatorProcess(0),
  fpEndingProcess(0),
  fFinalKineticEnergy(0.)
{
}

G4RichTrajectory::G4RichTrajectory

(

const G4Track *

aTrack

)

Definition at line 70 of file G4RichTrajectory.cc.

References G4Track::GetCreatorProcess(), G4Track::GetKineticEnergy(), G4Track::GetNextTouchableHandle(), and G4Track::GetTouchableHandle().

                                                       :
  G4Trajectory(aTrack)  // Note: this initialises the base class data
            // members and, unfortunately but never mind,
            // creates a G4TrajectoryPoint in
            // TrajectoryPointContainer that we cannot
            // access because it's private.  We store the
            // same information (plus more) in a
            // G4RichTrajectoryPoint in the
            // RichTrajectoryPointsContainer
{
  fpInitialVolume = aTrack->GetTouchableHandle();
  fpInitialNextVolume = aTrack->GetNextTouchableHandle();
  fpCreatorProcess = aTrack->GetCreatorProcess();
  // On construction, set final values to initial values.
  // Final values are updated at the addition of every step - see AppendStep.
  fpFinalVolume = aTrack->GetTouchableHandle();
  fpFinalNextVolume = aTrack->GetNextTouchableHandle();
  fpEndingProcess = aTrack->GetCreatorProcess();
  fFinalKineticEnergy = aTrack->GetKineticEnergy();
  // Insert the first rich trajectory point (see note above)...
  fpRichPointsContainer = new RichTrajectoryPointsContainer;
  fpRichPointsContainer->push_back(new G4RichTrajectoryPoint(aTrack));
}

G4RichTrajectory::G4RichTrajectory

(

G4RichTrajectory &

right

)

Definition at line 94 of file G4RichTrajectory.cc.

                                                          :
  G4Trajectory(right)
{
  fpInitialVolume = right.fpInitialVolume;
  fpInitialNextVolume = right.fpInitialNextVolume;
  fpCreatorProcess = right.fpCreatorProcess;
  fpFinalVolume = right.fpFinalVolume;
  fpFinalNextVolume = right.fpFinalNextVolume;
  fpEndingProcess = right.fpEndingProcess;
  fFinalKineticEnergy = right.fFinalKineticEnergy;
  fpRichPointsContainer = new RichTrajectoryPointsContainer;
  for(size_t i=0;i<right.fpRichPointsContainer->size();i++)
  {
    G4RichTrajectoryPoint* rightPoint =
      (G4RichTrajectoryPoint*)((*(right.fpRichPointsContainer))[i]);
    fpRichPointsContainer->push_back(new G4RichTrajectoryPoint(*rightPoint));
  }
}

G4RichTrajectory::~G4RichTrajectory ( )

virtual

Definition at line 113 of file G4RichTrajectory.cc.

{
  if (fpRichPointsContainer) {
    //  fpRichPointsContainer->clearAndDestroy();
    size_t i;
    for(i=0;i<fpRichPointsContainer->size();i++){
      delete  (*fpRichPointsContainer)[i];
    }
    fpRichPointsContainer->clear();
    delete fpRichPointsContainer;
  }
}

Member Function Documentation

void G4RichTrajectory::AppendStep ( const G4Step *  aStep )

virtual

Reimplemented from G4Trajectory.

Definition at line 126 of file G4RichTrajectory.cc.

References G4Track::GetCurrentStepNumber(), G4StepPoint::GetKineticEnergy(), G4Track::GetNextTouchableHandle(), G4Step::GetPostStepPoint(), G4Step::GetPreStepPoint(), G4StepPoint::GetProcessDefinedStep(), G4Step::GetTotalEnergyDeposit(), G4Track::GetTouchableHandle(), and G4Step::GetTrack().

{
  fpRichPointsContainer->push_back(new G4RichTrajectoryPoint(aStep));
  // Except for first step, which is a sort of virtual step to start
  // the track, compute the final values...
  const G4Track* track = aStep->GetTrack();
  const G4StepPoint* postStepPoint = aStep->GetPostStepPoint();
  if (track->GetCurrentStepNumber() > 0) {
    fpFinalVolume = track->GetTouchableHandle();
    fpFinalNextVolume = track->GetNextTouchableHandle();
    fpEndingProcess = postStepPoint->GetProcessDefinedStep();
    fFinalKineticEnergy =
      aStep->GetPreStepPoint()->GetKineticEnergy() -
      aStep->GetTotalEnergyDeposit();
  }
}

std::vector< G4AttValue > * G4RichTrajectory::CreateAttValues ( ) const

virtual

Reimplemented from G4Trajectory.

Definition at line 237 of file G4RichTrajectory.cc.

References G4Trajectory::CreateAttValues(), G4BestUnit, G4cout, GetAttDefs(), G4VProcess::GetProcessName(), G4VProcess::GetProcessType(), G4VProcess::GetProcessTypeName(), and G4VTouchable::GetVolume().

{
  // Create base class att values...
  std::vector<G4AttValue>* values = G4Trajectory::CreateAttValues();

  if (fpInitialVolume && fpInitialVolume->GetVolume()) {
    values->push_back(G4AttValue("IVPath",Path(fpInitialVolume),""));
  } else {
    values->push_back(G4AttValue("IVPath","None",""));
  }

  if (fpInitialNextVolume && fpInitialNextVolume->GetVolume()) {
    values->push_back(G4AttValue("INVPath",Path(fpInitialNextVolume),""));
  } else {
    values->push_back(G4AttValue("INVPath","None",""));
  }

  if (fpCreatorProcess) {
    values->push_back(G4AttValue("CPN",fpCreatorProcess->GetProcessName(),""));
    G4ProcessType type = fpCreatorProcess->GetProcessType();
    values->push_back(G4AttValue("CPTN",G4VProcess::GetProcessTypeName(type),""));
  } else {
    values->push_back(G4AttValue("CPN","None",""));
    values->push_back(G4AttValue("CPTN","None",""));
  }

  if (fpFinalVolume && fpFinalVolume->GetVolume()) {
    values->push_back(G4AttValue("FVPath",Path(fpFinalVolume),""));
  } else {
    values->push_back(G4AttValue("FVPath","None",""));
  }

  if (fpFinalNextVolume && fpFinalNextVolume->GetVolume()) {
    values->push_back(G4AttValue("FNVPath",Path(fpFinalNextVolume),""));
  } else {
    values->push_back(G4AttValue("FNVPath","None",""));
  }

  if (fpEndingProcess) {
    values->push_back(G4AttValue("EPN",fpEndingProcess->GetProcessName(),""));
    G4ProcessType type = fpEndingProcess->GetProcessType();
    values->push_back(G4AttValue("EPTN",G4VProcess::GetProcessTypeName(type),""));
  } else {
    values->push_back(G4AttValue("EPN","None",""));
    values->push_back(G4AttValue("EPTN","None",""));
  }

  values->push_back
    (G4AttValue("FKE",G4BestUnit(fFinalKineticEnergy,"Energy"),""));

#ifdef G4ATTDEBUG
  G4cout << G4AttCheck(values,GetAttDefs());
#endif

  return values;
}

void G4RichTrajectory::DrawTrajectory ( ) const

virtual

Reimplemented from G4Trajectory.

Definition at line 165 of file G4RichTrajectory.cc.

References G4VTrajectory::DrawTrajectory().

{
  // Invoke the default implementation in G4VTrajectory...
  G4VTrajectory::DrawTrajectory();
  // ... or override with your own code here.
}

const std::map< G4String, G4AttDef > * G4RichTrajectory::GetAttDefs ( ) const

virtual

Reimplemented from G4Trajectory.

Definition at line 172 of file G4RichTrajectory.cc.

References G4Trajectory::GetAttDefs(), and G4AttDefStore::GetInstance().

Referenced by CreateAttValues(), G4VisCommandList::SetNewValue(), and G4VisCommandSceneAddTrajectories::SetNewValue().

{
  G4bool isNew;
  std::map<G4String,G4AttDef>* store
    = G4AttDefStore::GetInstance("G4RichTrajectory",isNew);
  if (isNew) {

    // Get att defs from base class...
    *store = *(G4Trajectory::GetAttDefs());

    G4String ID;

    ID = "IVPath";
    (*store)[ID] = G4AttDef(ID,"Initial Volume Path",
                "Physics","","G4String");

    ID = "INVPath";
    (*store)[ID] = G4AttDef(ID,"Initial Next Volume Path",
                "Physics","","G4String");

    ID = "CPN";
    (*store)[ID] = G4AttDef(ID,"Creator Process Name",
                "Physics","","G4String");

    ID = "CPTN";
    (*store)[ID] = G4AttDef(ID,"Creator Process Type Name",
                "Physics","","G4String");

    ID = "FVPath";
    (*store)[ID] = G4AttDef(ID,"Final Volume Path",
                "Physics","","G4String");

    ID = "FNVPath";
    (*store)[ID] = G4AttDef(ID,"Final Next Volume Path",
                "Physics","","G4String");

    ID = "EPN";
    (*store)[ID] = G4AttDef(ID,"Ending Process Name",
                "Physics","","G4String");

    ID = "EPTN";
    (*store)[ID] = G4AttDef(ID,"Ending Process Type Name",
                "Physics","","G4String");

    ID = "FKE";
    (*store)[ID] = G4AttDef(ID,"Final kinetic energy",
                "Physics","G4BestUnit","G4double");

  }

  return store;
}

G4VTrajectoryPoint* G4RichTrajectory::GetPoint ( G4int  i ) const

inlinevirtual

Reimplemented from G4Trajectory.

Definition at line 91 of file G4RichTrajectory.hh.

   { return (*fpRichPointsContainer)[i]; }

int G4RichTrajectory::GetPointEntries ( ) const

inlinevirtual

Reimplemented from G4Trajectory.

Definition at line 90 of file G4RichTrajectory.hh.

Referenced by MergeTrajectory().

{ return fpRichPointsContainer->size(); }

void G4RichTrajectory::MergeTrajectory ( G4VTrajectory *  secondTrajectory )

virtual

Reimplemented from G4Trajectory.

Definition at line 143 of file G4RichTrajectory.cc.

References GetPointEntries().

{
  if(!secondTrajectory) return;

  G4RichTrajectory* seco = (G4RichTrajectory*)secondTrajectory;
  G4int ent = seco->GetPointEntries();
  for(G4int i=1;i<ent;i++) {
    // initial point of the second trajectory should not be merged
    fpRichPointsContainer->push_back((*(seco->fpRichPointsContainer))[i]);
    //    fpRichPointsContainer->push_back(seco->fpRichPointsContainer->removeAt(1));
  }
  delete (*seco->fpRichPointsContainer)[0];
  seco->fpRichPointsContainer->clear();
}

void G4RichTrajectory::operator delete ( void *  aRichTrajectory )

inline

Definition at line 122 of file G4RichTrajectory.hh.

References aRichTrajectoryAllocator.

{
  aRichTrajectoryAllocator->FreeSingle((G4RichTrajectory*)aRichTrajectory);
}

void * G4RichTrajectory::operator new ( size_t  )

inline

Definition at line 115 of file G4RichTrajectory.hh.

References aRichTrajectoryAllocator.

{
  if (!aRichTrajectoryAllocator)
  { aRichTrajectoryAllocator = new G4Allocator<G4RichTrajectory>; }
  return (void*)aRichTrajectoryAllocator->MallocSingle();
}

int G4RichTrajectory::operator== ( const G4RichTrajectory &  right ) const

inline

Definition at line 82 of file G4RichTrajectory.hh.

  {return (this==&right);} 

void G4RichTrajectory::ShowTrajectory ( std::ostream &  os = G4cout ) const

virtual

Reimplemented from G4Trajectory.

Definition at line 158 of file G4RichTrajectory.cc.

References G4VTrajectory::ShowTrajectory().

{
  // Invoke the default implementation in G4VTrajectory...
  G4VTrajectory::ShowTrajectory(os);
  // ... or override with your own code here.
}

---

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

• G4RichTrajectory.hh
• G4RichTrajectory.cc