RE01Trajectory Class Reference

#include <Doxymodules_runAndEvent.h>

Inheritance diagram for RE01Trajectory:

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

Private Attributes
RE01TrajectoryPointContainer *	fPositionRecord
G4int	fTrackID
G4int	fParentID
G4int	fTrackStatus
G4ParticleDefinition *	fParticleDefinition
G4String	fParticleName
G4double	fPDGCharge
G4int	fPDGEncoding
G4ThreeVector	fMomentum
G4ThreeVector	fVertexPosition
G4double	fGlobalTime

Detailed Description

Definition at line 38 of file Doxymodules_runAndEvent.h.

Constructor & Destructor Documentation

RE01Trajectory()

RE01Trajectory::RE01Trajectory

(

const G4Track *

aTrack

)

Definition at line 52 of file RE01Trajectory.cc.

:G4VTrajectory(),
 fPositionRecord(0),fParticleDefinition(0)
{
   fParticleDefinition = aTrack->GetDefinition();
   fParticleName = fParticleDefinition->GetParticleName();
   fPDGCharge = fParticleDefinition->GetPDGCharge();
   fPDGEncoding = fParticleDefinition->GetPDGEncoding();
   if(fParticleName=="unknown")
   {
     G4PrimaryParticle*pp = aTrack->GetDynamicParticle()->GetPrimaryParticle();
     if(pp)
     {
       if(pp->GetCharge()<DBL_MAX) fPDGCharge = pp->GetCharge();
       fPDGEncoding = pp->GetPDGcode();
       if(pp->GetG4code()!=0)
       {
         fParticleName += " : ";
         fParticleName += pp->GetG4code()->GetParticleName();
       }
     }
   }
   fTrackID = aTrack->GetTrackID();
   RE01TrackInformation* trackInfo
    = (RE01TrackInformation*)(aTrack->GetUserInformation());
   fTrackStatus = trackInfo->GetTrackingStatus();
   if(fTrackStatus == 1)
   { fParentID = aTrack->GetParentID(); }
   else if(fTrackStatus == 2)
   { fParentID = trackInfo->GetSourceTrackID(); }
   else
   { fParentID = -1; }
   fPositionRecord = new RE01TrajectoryPointContainer();
   fPositionRecord->push_back(new G4TrajectoryPoint(aTrack->GetPosition()));
   fMomentum = aTrack->GetMomentum();
   fVertexPosition = aTrack->GetPosition();
   fGlobalTime = aTrack->GetGlobalTime();
}

~RE01Trajectory()

RE01Trajectory::~RE01Trajectory ( )

virtual

Definition at line 92 of file RE01Trajectory.cc.

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

Member Function Documentation

ShowTrajectory()

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

virtual

Definition at line 104 of file RE01Trajectory.cc.

{
   os << G4endl << "TrackID =" << fTrackID 
        << " : ParentID=" << fParentID << " : TrackStatus=" << fTrackStatus << G4endl;
   os << "Particle name : " << fParticleName << "  PDG code : " << fPDGEncoding
        << "  Charge : " << fPDGCharge << G4endl;
   os << "Original momentum : " <<
        G4BestUnit(fMomentum,"Energy") << G4endl;
   os << "Vertex : " << G4BestUnit(fVertexPosition,"Length")
        << "  Global time : " << G4BestUnit(fGlobalTime,"Time") << G4endl;
   os << "  Current trajectory has " << fPositionRecord->size() 
        << " points." << G4endl;
 
   for( size_t i=0 ; i < fPositionRecord->size() ; i++){
       G4TrajectoryPoint* aTrajectoryPoint = 
         (G4TrajectoryPoint*)((*fPositionRecord)[i]);
       os << "Point[" << i << "]" 
            << " Position= " << aTrajectoryPoint->GetPosition() << G4endl;
   }
}

DrawTrajectory()

void RE01Trajectory::DrawTrajectory ( ) const

virtual

Definition at line 126 of file RE01Trajectory.cc.

{
 
   G4VVisManager* pVVisManager = G4VVisManager::GetConcreteInstance();
   G4ThreeVector pos;
 
   G4Polyline pPolyline;
   for (size_t i = 0; i < fPositionRecord->size() ; i++) {
     G4TrajectoryPoint* aTrajectoryPoint = 
       (G4TrajectoryPoint*)((*fPositionRecord)[i]);
     pos = aTrajectoryPoint->GetPosition();
     pPolyline.push_back( pos );
   }
 
   G4Colour colour(0.2,0.2,0.2);
   if(fParticleDefinition==G4Gamma::GammaDefinition())
      colour = G4Colour(0.,0.,1.);
   else if(fParticleDefinition==G4Electron::ElectronDefinition()
         ||fParticleDefinition==G4Positron::PositronDefinition())
      colour = G4Colour(1.,1.,0.);
   else if(fParticleDefinition==G4MuonMinus::MuonMinusDefinition()
         ||fParticleDefinition==G4MuonPlus::MuonPlusDefinition())
      colour = G4Colour(0.,1.,0.);
   else if(fParticleDefinition->GetParticleType()=="meson")
   {
      if(fPDGCharge!=0.)
         colour = G4Colour(1.,0.,0.);
      else
         colour = G4Colour(0.5,0.,0.);
   }
   else if(fParticleDefinition->GetParticleType()=="baryon")
   {
      if(fPDGCharge!=0.)
         colour = G4Colour(0.,1.,1.);
      else
         colour = G4Colour(0.,0.5,0.5);
   }
 
   G4VisAttributes attribs(colour);
   pPolyline.SetVisAttributes(attribs);
   if(pVVisManager) pVVisManager->Draw(pPolyline);
}

GetAttDefs()

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

virtual

Definition at line 170 of file RE01Trajectory.cc.

{
  G4bool isNew;
  std::map<G4String,G4AttDef>* store
    = G4AttDefStore::GetInstance("RE01Trajectory",isNew);
  if (isNew) {
 
    G4String id("ID");
    (*store)[id] = G4AttDef(id,"Track ID","Bookkeeping","","G4int");
 
    G4String pid("PID");
    (*store)[pid] = G4AttDef(pid,"Parent ID","Bookkeeping","","G4int");
 
    G4String status("Status");
    (*store)[status] = G4AttDef(status,"Track Status","Bookkeeping","","G4int");
 
    G4String pn("PN");
    (*store)[pn] = G4AttDef(pn,"Particle Name","Bookkeeping","","G4String");
 
    G4String ch("Ch");
    (*store)[ch] = G4AttDef(ch,"Charge","Physics","e+","G4double");
 
    G4String pdg("PDG");
    (*store)[pdg] = G4AttDef(pdg,"PDG Encoding","Bookkeeping","","G4int");
 
    G4String imom("IMom");
    (*store)[imom] = G4AttDef(imom, "Momentum of track at start of trajectory",
                              "Physics","G4BestUnit","G4ThreeVector");
 
    G4String imag("IMag");
    (*store)[imag] = 
      G4AttDef(imag, "Magnitude of momentum of track at start of trajectory",
               "Physics","G4BestUnit","G4double");
 
    G4String vtxPos("VtxPos");
    (*store)[vtxPos] = G4AttDef(vtxPos, "Vertex position",
                              "Physics","G4BestUnit","G4ThreeVector");
 
    G4String ntp("NTP");
    (*store)[ntp] = G4AttDef(ntp,"No. of points","Bookkeeping","","G4int");
 
  }
  return store;
}

CreateAttValues()

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

virtual

Definition at line 216 of file RE01Trajectory.cc.

{
  std::vector<G4AttValue>* values = new std::vector<G4AttValue>;
 
  values->push_back
    (G4AttValue("ID",G4UIcommand::ConvertToString(fTrackID),""));
 
  values->push_back
    (G4AttValue("PID",G4UIcommand::ConvertToString(fParentID),""));
 
  values->push_back
    (G4AttValue("Status",G4UIcommand::ConvertToString(fTrackStatus),""));
 
  values->push_back(G4AttValue("PN",fParticleName,""));
 
  values->push_back
    (G4AttValue("Ch",G4UIcommand::ConvertToString(fPDGCharge),""));
 
  values->push_back
    (G4AttValue("PDG",G4UIcommand::ConvertToString(fPDGEncoding),""));
 
  values->push_back
    (G4AttValue("IMom",G4BestUnit(fMomentum,"Energy"),""));
 
  values->push_back
    (G4AttValue("IMag",G4BestUnit(fMomentum.mag(),"Energy"),""));
 
  values->push_back
    (G4AttValue("VtxPos",G4BestUnit(fVertexPosition,"Length"),""));
 
  values->push_back
    (G4AttValue("NTP",G4UIcommand::ConvertToString(GetPointEntries()),""));
 
  return values;
}

AppendStep()

void RE01Trajectory::AppendStep ( const G4Step *  aStep )

virtual

Definition at line 253 of file RE01Trajectory.cc.

{
   fPositionRecord->push_back( 
       new G4TrajectoryPoint(aStep->GetPostStepPoint()->GetPosition() ));
}

MergeTrajectory()

void RE01Trajectory::MergeTrajectory ( G4VTrajectory *  secondTrajectory )

virtual

Definition at line 260 of file RE01Trajectory.cc.

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

operator new()

void * RE01Trajectory::operator new ( size_t  )

inline

Definition at line 98 of file RE01Trajectory.hh.

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

operator delete()

void RE01Trajectory::operator delete ( void *  aTrajectory )

inline

Definition at line 105 of file RE01Trajectory.hh.

{
  myTrajectoryAllocator->FreeSingle((RE01Trajectory*)aTrajectory);
}

operator==()

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

inline

Definition at line 68 of file RE01Trajectory.hh.

  {return (this==&right);} 

GetTrackID()

virtual G4int RE01Trajectory::GetTrackID ( ) const

inlinevirtual

Definition at line 71 of file RE01Trajectory.hh.

{ return fTrackID; }

GetParentID()

virtual G4int RE01Trajectory::GetParentID ( ) const

inlinevirtual

Definition at line 72 of file RE01Trajectory.hh.

{ return fParentID; }

GetParticleName()

virtual G4String RE01Trajectory::GetParticleName ( ) const

inlinevirtual

Definition at line 73 of file RE01Trajectory.hh.

{ return fParticleName; }

GetCharge()

virtual G4double RE01Trajectory::GetCharge ( ) const

inlinevirtual

Definition at line 74 of file RE01Trajectory.hh.

{ return fPDGCharge; }

GetPDGEncoding()

virtual G4int RE01Trajectory::GetPDGEncoding ( ) const

inlinevirtual

Definition at line 75 of file RE01Trajectory.hh.

{ return fPDGEncoding; }

GetInitialMomentum()

virtual G4ThreeVector RE01Trajectory::GetInitialMomentum ( ) const

inlinevirtual

Definition at line 76 of file RE01Trajectory.hh.

{ return fMomentum; }

GetPointEntries()

virtual int RE01Trajectory::GetPointEntries ( ) const

inlinevirtual

Definition at line 77 of file RE01Trajectory.hh.

{ return fPositionRecord->size(); }

GetPoint()

virtual G4VTrajectoryPoint * RE01Trajectory::GetPoint ( G4int  i ) const

inlinevirtual

Definition at line 78 of file RE01Trajectory.hh.

  { return (*fPositionRecord)[i]; }

Member Data Documentation

fPositionRecord

RE01TrajectoryPointContainer* RE01Trajectory::fPositionRecord

private

Definition at line 82 of file RE01Trajectory.hh.

fTrackID

G4int RE01Trajectory::fTrackID

private

Definition at line 83 of file RE01Trajectory.hh.

fParentID

G4int RE01Trajectory::fParentID

private

Definition at line 84 of file RE01Trajectory.hh.

fTrackStatus

G4int RE01Trajectory::fTrackStatus

private

Definition at line 85 of file RE01Trajectory.hh.

fParticleDefinition

G4ParticleDefinition* RE01Trajectory::fParticleDefinition

private

Definition at line 86 of file RE01Trajectory.hh.

fParticleName

G4String RE01Trajectory::fParticleName

private

Definition at line 87 of file RE01Trajectory.hh.

fPDGCharge

G4double RE01Trajectory::fPDGCharge

private

Definition at line 88 of file RE01Trajectory.hh.

fPDGEncoding

G4int RE01Trajectory::fPDGEncoding

private

Definition at line 89 of file RE01Trajectory.hh.

fMomentum

G4ThreeVector RE01Trajectory::fMomentum

private

Definition at line 90 of file RE01Trajectory.hh.

fVertexPosition

G4ThreeVector RE01Trajectory::fVertexPosition

private

Definition at line 91 of file RE01Trajectory.hh.

fGlobalTime

G4double RE01Trajectory::fGlobalTime

private

Definition at line 92 of file RE01Trajectory.hh.

---

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

• .doxygen/Doxymodules_runAndEvent.h
• extended/runAndEvent/RE01/include/RE01Trajectory.hh
• extended/runAndEvent/RE01/src/RE01Trajectory.cc