RE01Trajectory Class Reference

#include <RE01Trajectory.hh>

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
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.

References DBL_MAX, G4PrimaryParticle::GetCharge(), G4Track::GetDefinition(), G4Track::GetDynamicParticle(), G4PrimaryParticle::GetG4code(), G4Track::GetGlobalTime(), G4Track::GetMomentum(), G4Track::GetParentID(), G4ParticleDefinition::GetParticleName(), G4ParticleDefinition::GetPDGCharge(), G4PrimaryParticle::GetPDGcode(), G4ParticleDefinition::GetPDGEncoding(), G4Track::GetPosition(), G4DynamicParticle::GetPrimaryParticle(), RE01TrackInformation::GetSourceTrackID(), G4Track::GetTrackID(), RE01TrackInformation::GetTrackingStatus(), G4Track::GetUserInformation(), and G4InuclParticleNames::pp.

: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 ( )

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.

References G4StepPoint::GetPosition(), and G4Step::GetPostStepPoint().

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

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

virtual

Reimplemented from G4VTrajectory.

Definition at line 217 of file RE01Trajectory.cc.

References G4UIcommand::ConvertToString(), G4BestUnit, GetPointEntries(), and CLHEP::Hep3Vector::mag().

{
  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;
}

void RE01Trajectory::DrawTrajectory ( ) const

virtual

Reimplemented from G4VTrajectory.

Definition at line 127 of file RE01Trajectory.cc.

References G4VVisManager::Draw(), G4Electron::ElectronDefinition(), G4Gamma::GammaDefinition(), G4VVisManager::GetConcreteInstance(), G4ParticleDefinition::GetParticleType(), G4TrajectoryPoint::GetPosition(), G4MuonMinus::MuonMinusDefinition(), G4MuonPlus::MuonPlusDefinition(), G4Positron::PositronDefinition(), and G4Visible::SetVisAttributes().

{

   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);
}

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

virtual

Reimplemented from G4VTrajectory.

Definition at line 171 of file RE01Trajectory.cc.

References G4AttDefStore::GetInstance(), G4InuclParticleNames::pn, and status.

{
  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;
}

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.

Referenced by CreateAttValues(), and MergeTrajectory().

{ return fPositionRecord->size(); }

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.

References GetPointEntries().

{
  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();

}

void RE01Trajectory::operator delete ( void *  aTrajectory )

inline

Definition at line 106 of file RE01Trajectory.hh.

References myTrajectoryAllocator.

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

void * RE01Trajectory::operator new ( size_t  )

inline

Definition at line 99 of file RE01Trajectory.hh.

References myTrajectoryAllocator.

{
  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.

References G4BestUnit, G4endl, and G4TrajectoryPoint::GetPosition().

Referenced by RE01EventAction::EndOfEventAction().

{
   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;
   }
}

---

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

• RE01Trajectory.hh
• RE01Trajectory.cc