QweakSimTrajectory Class Reference

Stores the information about the various tracks. More...

#include <QweakSimTrajectory.hh>

Inherits G4VTrajectory.

Public Member Functions
	QweakSimTrajectory ()
	QweakSimTrajectory (const G4Track *aTrack)
	QweakSimTrajectory (QweakSimTrajectory &)
virtual	~QweakSimTrajectory ()
void *	operator new (size_t)
void	operator delete (void *)
int	operator== (const QweakSimTrajectory &right) const
G4int	GetTrackID () const
G4int	GetParentID () const
G4String	GetParticleName () const
G4double	GetCharge () const
G4int	GetPDGEncoding () const
G4ThreeVector	GetInitialMomentum () const
virtual int	GetPointEntries () const
virtual G4VTrajectoryPoint *	GetPoint (G4int i) const
virtual void	ShowTrajectory (std::ostream &os=G4cout) const
virtual void	DrawTrajectory (G4int i_mode=0) const
virtual void	AppendStep (const G4Step *aStep)
virtual void	MergeTrajectory (G4VTrajectory *secondTrajectory)
void	SetDrawTrajectory (G4bool b)
void	SetForceDrawTrajectory (G4bool b)
void	SetForceNoDrawTrajectory (G4bool b)
G4ParticleDefinition *	GetParticleDefinition ()

Private Attributes
QweakSimTrajectoryPointContainer *	positionRecord
G4int	fTrackID
G4int	fParentID
G4ParticleDefinition *	fpParticleDefinition
G4String	ParticleName
G4double	PDGCharge
G4int	PDGEncoding
G4ThreeVector	InitialMomentum
G4ThreeVector	momentum
G4ThreeVector	vertexPosition
G4double	globalTime
G4bool	drawit
G4bool	forceNoDraw
G4bool	forceDraw

Detailed Description

Stores the information about the various tracks.

It can be viewed as a perfect track reconstruction in a hypothetical perfect Qweak setup

Placeholder for a long explaination

Definition at line 71 of file QweakSimTrajectory.hh.

Constructor & Destructor Documentation

QweakSimTrajectory::QweakSimTrajectory

(

)

Definition at line 60 of file QweakSimTrajectory.cc.

References drawit, forceDraw, forceNoDraw, fParentID, fpParticleDefinition, fTrackID, globalTime, momentum, ParticleName, PDGCharge, PDGEncoding, positionRecord, and vertexPosition.

{
  G4cout << G4endl << "###### Calling QweakSimTrajectory::QweakSimTrajectory()" << G4endl << G4endl;

    drawit       = FALSE;
    forceNoDraw  = FALSE;
    forceDraw    = FALSE;


   fpParticleDefinition = 0;
   ParticleName         = "";
   PDGCharge            = 0;
   PDGEncoding          = 0;
   fTrackID             = 0;
   fParentID            = 0;
   positionRecord       = 0;

   momentum             = G4ThreeVector(0.,0.,0.);
   vertexPosition       = G4ThreeVector(0.,0.,0.);
   globalTime           = 0.;

   G4cout << G4endl << "###### Leaving QweakSimTrajectory::QweakSimTrajectory()" << G4endl << G4endl;
}

QweakSimTrajectory::QweakSimTrajectory

(

const G4Track *

aTrack

)

Definition at line 85 of file QweakSimTrajectory.cc.

References fParentID, fpParticleDefinition, fTrackID, globalTime, momentum, ParticleName, PDGCharge, PDGEncoding, positionRecord, and vertexPosition.

{
 //   G4cout << G4endl << "###### Calling QweakSimTrajectory::QweakSimTrajectory(aTrack)" << G4endl << G4endl;

   fpParticleDefinition = aTrack->GetDefinition();
   ParticleName         = fpParticleDefinition->GetParticleName();
   PDGCharge            = fpParticleDefinition->GetPDGCharge();
   PDGEncoding          = fpParticleDefinition->GetPDGEncoding();
   fTrackID             = aTrack->GetTrackID();
   fParentID            = aTrack->GetParentID();

   positionRecord =  new QweakSimTrajectoryPointContainer();
   positionRecord -> push_back(new G4TrajectoryPoint(aTrack->GetPosition()));

   momentum       = aTrack->GetMomentum();
   vertexPosition = aTrack->GetPosition();
   globalTime     = aTrack->GetGlobalTime();


   // G4cout << G4endl << "###### Leaving QweakSimTrajectory::QweakSimTrajectory(aTrack)" << G4endl << G4endl;
}

QweakSimTrajectory::QweakSimTrajectory

(

QweakSimTrajectory &

right

)

Definition at line 108 of file QweakSimTrajectory.cc.

References fParentID, fpParticleDefinition, fTrackID, globalTime, momentum, ParticleName, PDGCharge, PDGEncoding, positionRecord, and vertexPosition.

    : G4VTrajectory()
{
  ParticleName          = right.ParticleName;
  fpParticleDefinition  = right.fpParticleDefinition;
  PDGCharge             = right.PDGCharge;
  PDGEncoding           = right.PDGEncoding;
  fTrackID              = right.fTrackID;
  fParentID             = right.fParentID;

  positionRecord = new QweakSimTrajectoryPointContainer();

  for(int i=0;i<(int)right.positionRecord->size();i++)
  {
    G4TrajectoryPoint* rightPoint = (G4TrajectoryPoint*)((*(right.positionRecord))[i]);

    positionRecord->push_back(new G4TrajectoryPoint(*rightPoint));
  }

  momentum       = right.momentum;
  vertexPosition = right.vertexPosition;
  globalTime     = right.globalTime;
}

QweakSimTrajectory::~QweakSimTrajectory ( )

virtual

Definition at line 133 of file QweakSimTrajectory.cc.

References positionRecord.

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

  delete positionRecord;
}

Member Function Documentation

void QweakSimTrajectory::AppendStep ( const G4Step *  aStep )

virtual

Definition at line 288 of file QweakSimTrajectory.cc.

References positionRecord.

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

void QweakSimTrajectory::DrawTrajectory ( G4int  i_mode = 0 ) const

virtual

Definition at line 172 of file QweakSimTrajectory.cc.

Referenced by QweakSimEventAction::EndOfEventAction().

{
 //  G4cout << G4endl << "###### Calling QweakSimTrajectory::DrawTrajectory()" << G4endl << G4endl;

    //if(!forceDraw && (!drawit || forceNoDraw))

//    if( (forceNoDraw == true))
//         return;

  // If i_mode>=0, draws a trajectory as a polyline (blue for
  // positive, red for negative, green for neutral) and, if i_mode!=0,
  // adds markers - yellow circles for step points and magenta squares
  // for auxiliary points, if any - whose screen size in pixels is
  // given by abs(i_mode)/1000.  E.g: i_mode = 5000 gives easily
  // visible markers.

   G4VVisManager* pVVisManager = G4VVisManager::GetConcreteInstance();

   G4ThreeVector pos;
   G4Polyline    pPolyline;

   for (int i = 0; i < (int) positionRecord->size() ; i++) {

       G4TrajectoryPoint* aTrajectoryPoint = (G4TrajectoryPoint*)((*positionRecord)[i]);
       pos = aTrajectoryPoint->GetPosition();
       pPolyline.push_back( pos );
   }


   G4Colour colour(0.75,0.75,0.75);    // LightGray

   G4Colour red      ( 255/255.,   0/255.,   0/255.);
   G4Colour blue     (   0/255.,   0/255., 255/255.);   
   G4Colour green    (   0/255., 255/255.,   0/255.);
   G4Colour yellow   ( 255/255., 255/255.,   0/255.);

   G4Colour white    ( 255/255., 255/255., 255/255.);
   
   G4Colour orange   ( 255/255., 127/255.,   0/255.);
   G4Colour magenta  ( 237/255., 173/255., 255/255.);   
   G4Colour magenta1 ( 104/255.,  49/255.,  94/255.);



   // Default colors for particlea
        if ( PDGCharge < 0.) {colour = red;}
   else if ( PDGCharge > 0.) {colour = blue;}
   else                      
   {   
       //G4cout << "%%%%% We have something neutral here !!!" << G4endl;
       colour = green;
   }
   
        if (fpParticleDefinition == G4Gamma     ::GammaDefinition()     )  {
      //G4cout << "%%%%% We have a green Gamma here !!!" << G4endl;
      colour = green;} 
        if (fpParticleDefinition == G4Electron  ::ElectronDefinition()  )  
      {
      //G4cout << "%%%%% We have a red Electron here !!!" << G4endl;
      colour = red;  
      //colour = white;  
      } 
        if (fpParticleDefinition == G4Positron  ::PositronDefinition()  )  
      {
      //G4cout << "%%%%% We have a blue Positron here !!!" << G4endl;
      colour = blue; 
      } 

//    else if (fpParticleDefinition == G4MuonMinus ::MuonMinusDefinition()
//           ||fpParticleDefinition == G4MuonPlus  ::MuonPlusDefinition()) ) colour = G4Colour(1.,0.,1.);      // Magenta
//
//    else if(fpParticleDefinition->GetParticleType()=="meson")
//    {
//       if(PDGCharge!=0.)
//          colour = yellow;
//       else
//          colour = G4Colour(0.5,0.,0.);  // HalfRed
//    }
//    else if(fpParticleDefinition->GetParticleType()=="baryon")
//    {
//       if(PDGCharge!=0.)
//          colour = orange; // Orange
//       else
//          colour = G4Colour(0.5,0.39,0.);// HalfOrange
//    }

      if( fpParticleDefinition == G4Neutron       ::NeutronDefinition())       
      { 
    //G4cout << "%%%%% We have a white Neutron here !!!" << G4endl;
    colour = white;   
      }
      if( fpParticleDefinition == G4Proton        ::ProtonDefinition())        
      {
    //G4cout << "%%%%% We have an orange Proton here !!!" << G4endl;
    colour = orange;  
      }
      if( fpParticleDefinition == G4OpticalPhoton ::OpticalPhotonDefinition()) 
      { 
    //G4cout << "%%%%% We have a magenta OpticalPhoton here !!!" << G4endl;  
    colour = magenta; 
      }

   G4VisAttributes attribs(colour);
   pPolyline.SetVisAttributes(attribs);

  // if(pVVisManager) 
   pVVisManager->Draw(pPolyline);


  // G4cout << G4endl << "###### Leaving QweakSimTrajectory::DrawTrajectory()" << G4endl << G4endl;
}

Here is the caller graph for this function:

G4double QweakSimTrajectory::GetCharge ( ) const

inline

Definition at line 99 of file QweakSimTrajectory.hh.

References PDGCharge.

Referenced by QweakSimCerenkov_PMTSD::GetParentTrajectory(), and QweakSimPMTOnly_PMTSD::GetParentTrajectory().

   { return PDGCharge; }

Here is the caller graph for this function:

G4ThreeVector QweakSimTrajectory::GetInitialMomentum ( ) const

inline

Definition at line 103 of file QweakSimTrajectory.hh.

References InitialMomentum.

   { return InitialMomentum; }

G4int QweakSimTrajectory::GetParentID ( ) const

inline

Definition at line 95 of file QweakSimTrajectory.hh.

References fParentID.

   { return fParentID; }

G4ParticleDefinition * QweakSimTrajectory::GetParticleDefinition

(

)

Definition at line 294 of file QweakSimTrajectory.cc.

References ParticleName.

{
   return (G4ParticleTable::GetParticleTable()->FindParticle(ParticleName));
}

G4String QweakSimTrajectory::GetParticleName ( ) const

inline

Definition at line 97 of file QweakSimTrajectory.hh.

References ParticleName.

   { return ParticleName; }

G4int QweakSimTrajectory::GetPDGEncoding ( ) const

inline

Definition at line 101 of file QweakSimTrajectory.hh.

References PDGEncoding.

   { return PDGEncoding; }

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

inlinevirtual

Definition at line 108 of file QweakSimTrajectory.hh.

References positionRecord.

       { return (*positionRecord)[i]; }

virtual int QweakSimTrajectory::GetPointEntries ( ) const

inlinevirtual

Definition at line 106 of file QweakSimTrajectory.hh.

References positionRecord.

Referenced by MergeTrajectory().

       { return positionRecord->size(); }

Here is the caller graph for this function:

G4int QweakSimTrajectory::GetTrackID ( ) const

inline

Definition at line 93 of file QweakSimTrajectory.hh.

References fTrackID.

Referenced by QweakSimCerenkov_PMTSD::GetParentTrajectory(), and QweakSimPMTOnly_PMTSD::GetParentTrajectory().

   { return fTrackID; }

Here is the caller graph for this function:

void QweakSimTrajectory::MergeTrajectory ( G4VTrajectory *  secondTrajectory )

virtual

Definition at line 300 of file QweakSimTrajectory.cc.

References GetPointEntries(), and positionRecord.

{
  if(!secondTrajectory) return;

  QweakSimTrajectory* seco = (QweakSimTrajectory*)secondTrajectory;
  G4int ent = seco->GetPointEntries();

  for(int i=1;i<ent;i++) // initial point of the second trajectory should not be merged
  {
    positionRecord->push_back((*(seco->positionRecord))[i]);
  }

  delete (*seco->positionRecord)[0];

  seco->positionRecord->clear();

}

Here is the call graph for this function:

void QweakSimTrajectory::operator delete ( void *  aTrajectory )

inline

Definition at line 166 of file QweakSimTrajectory.hh.

References myTrajectoryAllocator.

{
  myTrajectoryAllocator.FreeSingle((QweakSimTrajectory*)aTrajectory);
}

void * QweakSimTrajectory::operator new ( size_t  )

inline

Definition at line 158 of file QweakSimTrajectory.hh.

References myTrajectoryAllocator.

{
  void* aTrajectory;
  aTrajectory = (void*)myTrajectoryAllocator.MallocSingle();
  return aTrajectory;
}

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

inline

Definition at line 89 of file QweakSimTrajectory.hh.

   {return (this==&right);}

void QweakSimTrajectory::SetDrawTrajectory ( G4bool  b )

inline

Definition at line 123 of file QweakSimTrajectory.hh.

References drawit.

{drawit=b;}

void QweakSimTrajectory::SetForceDrawTrajectory ( G4bool  b )

inline

Definition at line 124 of file QweakSimTrajectory.hh.

References forceDraw.

{forceDraw=b;}

void QweakSimTrajectory::SetForceNoDrawTrajectory ( G4bool  b )

inline

Definition at line 125 of file QweakSimTrajectory.hh.

References forceNoDraw.

{forceNoDraw=b;}

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

virtual

Definition at line 145 of file QweakSimTrajectory.cc.

References fParentID, fTrackID, globalTime, momentum, ParticleName, PDGCharge, PDGEncoding, positionRecord, and vertexPosition.

{

  os << G4endl << "TrackID =" << fTrackID 
     << " : ParentID=" << fParentID << G4endl;

  os << "Particle name : " << ParticleName << "  PDG code : " << PDGEncoding
     << "  Charge : " << PDGCharge << G4endl;

  os << "Original momentum : " <<
      G4BestUnit(momentum,"Energy") << G4endl;

  os << "Vertex : "        << G4BestUnit(vertexPosition,"Length")
     << "  Global time : " << G4BestUnit(globalTime,"Time") << G4endl;

  os << "  Current trajectory has " << positionRecord->size() 
     << " points." << G4endl;

   for( size_t i=0 ; i < positionRecord->size() ; i++){
       G4TrajectoryPoint* aTrajectoryPoint = (G4TrajectoryPoint*)((*positionRecord)[i]);

       G4cout << "Point[" << i << "]" << " Position= " << aTrajectoryPoint->GetPosition() << G4endl;
   }
}

Field Documentation

G4bool QweakSimTrajectory::drawit

private

Definition at line 148 of file QweakSimTrajectory.hh.

Referenced by QweakSimTrajectory(), and SetDrawTrajectory().

G4bool QweakSimTrajectory::forceDraw

private

Definition at line 150 of file QweakSimTrajectory.hh.

Referenced by QweakSimTrajectory(), and SetForceDrawTrajectory().

G4bool QweakSimTrajectory::forceNoDraw

private

Definition at line 149 of file QweakSimTrajectory.hh.

Referenced by QweakSimTrajectory(), and SetForceNoDrawTrajectory().

G4int QweakSimTrajectory::fParentID

private

Definition at line 135 of file QweakSimTrajectory.hh.

Referenced by GetParentID(), QweakSimTrajectory(), and ShowTrajectory().

G4ParticleDefinition* QweakSimTrajectory::fpParticleDefinition

private

Definition at line 136 of file QweakSimTrajectory.hh.

Referenced by QweakSimTrajectory().

G4int QweakSimTrajectory::fTrackID

private

Definition at line 134 of file QweakSimTrajectory.hh.

Referenced by GetTrackID(), QweakSimTrajectory(), and ShowTrajectory().

G4double QweakSimTrajectory::globalTime

private

Definition at line 145 of file QweakSimTrajectory.hh.

Referenced by QweakSimTrajectory(), and ShowTrajectory().

G4ThreeVector QweakSimTrajectory::InitialMomentum

private

Definition at line 141 of file QweakSimTrajectory.hh.

Referenced by GetInitialMomentum().

G4ThreeVector QweakSimTrajectory::momentum

private

Definition at line 143 of file QweakSimTrajectory.hh.

Referenced by QweakSimTrajectory(), and ShowTrajectory().

G4String QweakSimTrajectory::ParticleName

private

Definition at line 137 of file QweakSimTrajectory.hh.

Referenced by GetParticleDefinition(), GetParticleName(), QweakSimTrajectory(), and ShowTrajectory().

G4double QweakSimTrajectory::PDGCharge

private

Definition at line 138 of file QweakSimTrajectory.hh.

Referenced by GetCharge(), QweakSimTrajectory(), and ShowTrajectory().

G4int QweakSimTrajectory::PDGEncoding

private

Definition at line 139 of file QweakSimTrajectory.hh.

Referenced by GetPDGEncoding(), QweakSimTrajectory(), and ShowTrajectory().

QweakSimTrajectoryPointContainer* QweakSimTrajectory::positionRecord

private

Definition at line 133 of file QweakSimTrajectory.hh.

Referenced by AppendStep(), GetPoint(), GetPointEntries(), MergeTrajectory(), QweakSimTrajectory(), ShowTrajectory(), and ~QweakSimTrajectory().

G4ThreeVector QweakSimTrajectory::vertexPosition

private

Definition at line 144 of file QweakSimTrajectory.hh.

Referenced by QweakSimTrajectory(), and ShowTrajectory().

---

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

• QweakSimTrajectory.hh
• QweakSimTrajectory.cc