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

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 54 of file RE01Trajectory.hh.

Constructor & Destructor Documentation

RE01Trajectory()

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::~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

AppendStep()

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

CreateAttValues()

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:

DrawTrajectory()

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:

GetAttDefs()

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:

GetCharge()

virtual G4double RE01Trajectory::GetCharge ( ) const

inlinevirtual

Implements G4VTrajectory.

Definition at line 75 of file RE01Trajectory.hh.

{ return fPDGCharge; }

GetInitialMomentum()

virtual G4ThreeVector RE01Trajectory::GetInitialMomentum ( ) const

inlinevirtual

Implements G4VTrajectory.

Definition at line 77 of file RE01Trajectory.hh.

{ return fMomentum; }

GetParentID()

virtual G4int RE01Trajectory::GetParentID ( ) const

inlinevirtual

Implements G4VTrajectory.

Definition at line 73 of file RE01Trajectory.hh.

{ return fParentID; }

GetParticleName()

virtual G4String RE01Trajectory::GetParticleName ( ) const

inlinevirtual

Implements G4VTrajectory.

Definition at line 74 of file RE01Trajectory.hh.

{ return fParticleName; }

GetPDGEncoding()

virtual G4int RE01Trajectory::GetPDGEncoding ( ) const

inlinevirtual

Implements G4VTrajectory.

Definition at line 76 of file RE01Trajectory.hh.

{ return fPDGEncoding; }

GetPoint()

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

inlinevirtual

Implements G4VTrajectory.

Definition at line 79 of file RE01Trajectory.hh.

  { return (*fPositionRecord)[i]; }

GetPointEntries()

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:

GetTrackID()

virtual G4int RE01Trajectory::GetTrackID ( ) const

inlinevirtual

Implements G4VTrajectory.

Definition at line 72 of file RE01Trajectory.hh.

{ return fTrackID; }

MergeTrajectory()

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:

operator delete()

void RE01Trajectory::operator delete ( void *  aTrajectory )

inline

Definition at line 106 of file RE01Trajectory.hh.

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

operator new()

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

operator==()

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

inline

Definition at line 69 of file RE01Trajectory.hh.

  {return (this==&right);} 

ShowTrajectory()

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:

Member Data Documentation

fGlobalTime

G4double RE01Trajectory::fGlobalTime

private

Definition at line 93 of file RE01Trajectory.hh.

fMomentum

G4ThreeVector RE01Trajectory::fMomentum

private

Definition at line 91 of file RE01Trajectory.hh.

fParentID

G4int RE01Trajectory::fParentID

private

Definition at line 85 of file RE01Trajectory.hh.

fParticleDefinition

G4ParticleDefinition* RE01Trajectory::fParticleDefinition

private

Definition at line 87 of file RE01Trajectory.hh.

fParticleName

G4String RE01Trajectory::fParticleName

private

Definition at line 88 of file RE01Trajectory.hh.

fPDGCharge

G4double RE01Trajectory::fPDGCharge

private

Definition at line 89 of file RE01Trajectory.hh.

fPDGEncoding

G4int RE01Trajectory::fPDGEncoding

private

Definition at line 90 of file RE01Trajectory.hh.

fPositionRecord

RE01TrajectoryPointContainer* RE01Trajectory::fPositionRecord

private

Definition at line 83 of file RE01Trajectory.hh.

fTrackID

G4int RE01Trajectory::fTrackID

private

Definition at line 84 of file RE01Trajectory.hh.

fTrackStatus

G4int RE01Trajectory::fTrackStatus

private

Definition at line 86 of file RE01Trajectory.hh.

fVertexPosition

G4ThreeVector RE01Trajectory::fVertexPosition

private

Definition at line 92 of file RE01Trajectory.hh.

---

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

• RE01Trajectory.hh
• RE01Trajectory.cc