RE01Trajectory Class Reference

#include <RE01Trajectory.hh>

Inheritance diagram for RE01Trajectory:

Collaboration 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
Public Member Functions inherited from G4VTrajectory
	G4VTrajectory ()
virtual	~G4VTrajectory ()
G4bool	operator== (const G4VTrajectory &right) const

Detailed Description

Definition at line 54 of file RE01Trajectory.hh.

Constructor & Destructor Documentation

RE01Trajectory::RE01Trajectory

(

const G4Track *

aTrack

)

Definition at line 53 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();
}

Here is the call graph for this function:

RE01Trajectory::~RE01Trajectory ( )

virtual

Definition at line 93 of file RE01Trajectory.cc.

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

  delete fPositionRecord;
}

Member Function Documentation

void RE01Trajectory::AppendStep ( const G4Step *  aStep )

virtual

Implements G4VTrajectory.

Definition at line 254 of file RE01Trajectory.cc.

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

Here is the call graph for this function:

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

virtual

Reimplemented from G4VTrajectory.

Definition at line 217 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;
}

Here is the call graph for this function:

void RE01Trajectory::DrawTrajectory ( ) const

virtual

Reimplemented from G4VTrajectory.

Definition at line 127 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);
}

Here is the call graph for this function:

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

virtual

Reimplemented from G4VTrajectory.

Definition at line 171 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;
}

Here is the call graph for this function:

virtual G4double RE01Trajectory::GetCharge ( ) const

inlinevirtual

Implements G4VTrajectory.

Definition at line 75 of file RE01Trajectory.hh.

{ return fPDGCharge; }

virtual G4ThreeVector RE01Trajectory::GetInitialMomentum ( ) const

inlinevirtual

Implements G4VTrajectory.

Definition at line 77 of file RE01Trajectory.hh.

{ return fMomentum; }

virtual G4int RE01Trajectory::GetParentID ( ) const

inlinevirtual

Implements G4VTrajectory.

Definition at line 73 of file RE01Trajectory.hh.

{ return fParentID; }

virtual G4String RE01Trajectory::GetParticleName ( ) const

inlinevirtual

Implements G4VTrajectory.

Definition at line 74 of file RE01Trajectory.hh.

{ return fParticleName; }

virtual G4int RE01Trajectory::GetPDGEncoding ( ) const

inlinevirtual

Implements G4VTrajectory.

Definition at line 76 of file RE01Trajectory.hh.

{ return fPDGEncoding; }

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

inlinevirtual

Implements G4VTrajectory.

Definition at line 79 of file RE01Trajectory.hh.

  { return (*fPositionRecord)[i]; }

virtual int RE01Trajectory::GetPointEntries ( ) const

inlinevirtual

Implements G4VTrajectory.

Definition at line 78 of file RE01Trajectory.hh.

{ return fPositionRecord->size(); }

Here is the caller graph for this function:

virtual G4int RE01Trajectory::GetTrackID ( ) const

inlinevirtual

Implements G4VTrajectory.

Definition at line 72 of file RE01Trajectory.hh.

{ return fTrackID; }

void RE01Trajectory::MergeTrajectory ( G4VTrajectory *  secondTrajectory )

virtual

Implements G4VTrajectory.

Definition at line 261 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();

}

Here is the call graph for this function:

void RE01Trajectory::operator delete ( void *  aTrajectory )

inline

Definition at line 106 of file RE01Trajectory.hh.

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

void * RE01Trajectory::operator new ( size_t  )

inline

Definition at line 99 of file RE01Trajectory.hh.

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

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

inline

Definition at line 69 of file RE01Trajectory.hh.

  {return (this==&right);} 

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

virtual

Reimplemented from G4VTrajectory.

Definition at line 105 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;
   }
}

Here is the call graph for this function:

Here is the caller graph for this function:

---

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

• source/geant4.10.03.p02/examples/extended/runAndEvent/RE01/include/RE01Trajectory.hh
• source/geant4.10.03.p02/examples/extended/runAndEvent/RE01/src/RE01Trajectory.cc