LBNETrajectory Class Reference

#include <LBNETrajectory.hh>

Inheritance diagram for LBNETrajectory:

Public Member Functions
	LBNETrajectory ()
	LBNETrajectory (const G4Track *aTrack)
	LBNETrajectory (LBNETrajectory &)
	LBNETrajectory (const LBNETrajectory *right)
virtual	~LBNETrajectory ()
void *	operator new (size_t)
void	operator delete (void *)
int	operator== (const LBNETrajectory &right) const
G4int	GetTrackID () const
G4int	GetParentID () const
G4String	GetParticleName () const
G4double	GetCharge () const
G4double	GetMass () const
G4int	GetPDGEncoding () const
const G4ThreeVector &	GetVertexPosition () const
virtual int	GetPointEntries () const
virtual G4VTrajectoryPoint *	GetPoint (G4int i) const
virtual G4ThreeVector	GetMomentum (G4int i) const
virtual G4String	GetPreStepVolumeName (G4int i) const
G4ThreeVector	GetInitialMomentum () const
virtual G4int	GetTgen () const
G4int	GetDecayCode () const
virtual G4double	GetNImpWt () const
virtual G4double	GetStepLength (G4int i) const
G4int	GetMaterialNumber1rst () const
void	SetMaterialNumber1rst (G4double matNum)
G4int	GetMaterialNumberLast () const
void	SetMaterialNumberLast (G4double matNum)
G4String	GetVolName1rst () const
G4double	GetTimeStart () const
G4ThreeVector	GetPolarization () const
G4int	GetPDGNucleus () const
G4String	GetProcessName () const
G4String	GetMaterialName1rst () const
virtual G4String	GetMaterialName (G4int i) const
const G4ThreeVector	GetParentMomentumAtThisProduction () const
virtual void	ShowTrajectory () const
virtual void	ShowTrajectory (std::ostream &o) 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)
G4ParticleDefinition *	GetParticleDefinition ()

Private Attributes
LBNETrajectoryMomentumContainer *	fMomentumRecord
LBNETrajectoryPointContainer *	fPositionRecord
G4int	fDecayCode
G4int	fEventNo
G4int	fTgen
G4double	fNImpWt
G4int	fTrackID
G4int	fParentID
G4ParticleDefinition *	fParticleDefinition
G4String	fParticleName
G4double	fPDGCharge
G4int	fPDGEncoding
G4ThreeVector	fMomentum
G4ThreeVector	fVertexPosition
G4double	fParticleMass
LBNETrajectoryVolumeName *	fPreStepVolume
DVec *	fStepLength
G4int	fMaterialNumber1rst
G4String	fVolName1rst
G4int	fMaterialNumberLast
G4double	fTimeStart
G4ThreeVector	fPolarization
G4int	fPDGNucleus
G4String	fProcessName
G4ThreeVector	fParentMomentumAtThisProduction
LBNETrajectoryMaterialName *	fMaterialName
G4String	fMaterialName1rst

Detailed Description

Definition at line 32 of file LBNETrajectory.hh.

Constructor & Destructor Documentation

LBNETrajectory::LBNETrajectory

(

)

Definition at line 30 of file LBNETrajectory.cc.

{
   fParticleDefinition = 0;
   fParticleName       = "";
   fPDGCharge          = 0;
   fPDGEncoding        = 0;
   fTrackID            = 0;
   fParentID           = 0;
   fPositionRecord     = 0;
   fMomentum           = G4ThreeVector(0.,0.,0.);
   fMomentumRecord     = 0;
   fVertexPosition     = G4ThreeVector(0.,0.,0.);
   fParticleMass       = 0.;
   fDecayCode          = 0;
   fTgen               = 0;
   fNImpWt             = 1.;
   fPreStepVolume      = 0;
   fStepLength         = 0;
   fMaterialNumber1rst = 0;
   fVolName1rst        = std::string("");
   fMaterialNumberLast = 0;
   fTimeStart = 0.;
   fPolarization       = G4ThreeVector(0.,0.,0.);
   fParentMomentumAtThisProduction = G4ThreeVector(0.,0.,0.); //for ppfx....A. Bashyal
   fPDGNucleus         = 0;
   fProcessName        = std::string("");
   fMaterialName1rst   = std::string("");

}

LBNETrajectory::LBNETrajectory

(

const G4Track *

aTrack

)

Definition at line 60 of file LBNETrajectory.cc.

{
   fParticleDefinition = aTrack->GetDefinition();
   fParticleName       = fParticleDefinition->GetParticleName();
   fPDGCharge          = fParticleDefinition->GetPDGCharge();
   fPDGEncoding        = fParticleDefinition->GetPDGEncoding();
   fTrackID            = aTrack->GetTrackID();
   fParentID           = aTrack->GetParentID();
   fPositionRecord     = new LBNETrajectoryPointContainer();
   fPositionRecord -> push_back(new G4TrajectoryPoint(aTrack->GetPosition()));
   fMomentumRecord     = new LBNETrajectoryMomentumContainer();
   fMomentumRecord -> push_back(aTrack->GetMomentum());
   fPreStepVolume      = new LBNETrajectoryVolumeName();  
   fPreStepVolume  -> push_back(aTrack->GetVolume()->GetName());
   fMaterialName       = new LBNETrajectoryMaterialName();
   fMaterialName   -> push_back(aTrack->GetVolume()->GetLogicalVolume()->GetMaterial()->GetName());
   fStepLength         = new DVec();  
   fStepLength     -> push_back(aTrack->GetStepLength());
   fMomentum           = aTrack->GetMomentum();
   fVertexPosition     = aTrack->GetPosition();
   fParticleMass       = aTrack->GetDefinition()->GetPDGMass();
   if (!aTrack->GetLogicalVolumeAtVertex()) {
//     std::cerr << " LBNETrajectory::LBNETrajectory No logical volume at vertex... ???    " << std::endl;
     fMaterialNumber1rst = 0; 
     fMaterialName1rst = std::string("Unknown"); 
     fVolName1rst = std::string("Unknown"); 
   } else {
     fMaterialNumber1rst = aTrack->GetLogicalVolumeAtVertex()->GetMaterial()->GetIndex(); 
     fMaterialName1rst = aTrack->GetLogicalVolumeAtVertex()->GetMaterial()->GetName(); 
     fVolName1rst = aTrack->GetLogicalVolumeAtVertex()->GetName(); 
   }
   fTimeStart = aTrack->GetGlobalTime();
   fPolarization = aTrack->GetPolarization();
   
   LBNETrackInformation* info = (LBNETrackInformation*)(aTrack->GetUserInformation());
   if (info!=0) 
   {
     fDecayCode = info->GetDecayCode();
     fTgen      = info->GetTgen();
     fNImpWt    = info->GetNImpWt();
     fParentMomentumAtThisProduction = info->GetParentMomentumAtThisProduction();     
//     std::cerr << " Transfering ImpWeight = " <<  fNImpWt << " from trInfo, track id  " <<  aTrack->GetTrackID() << std::endl;
   }
   else 
   { 
   fDecayCode = 0;
     fTgen      = 0; 
     fNImpWt    = 1.;
     fParentMomentumAtThisProduction = G4ThreeVector(0.,0.,0.);     
   }
   // Now the tricky bit: Dk2Nu package wants the process name and the nucleus 
   if (aTrack->GetCreatorProcess() != 0)
     fProcessName = aTrack->GetCreatorProcess()->GetProcessName();
   else fProcessName = std::string ("BeamParticle");
   // Well, that one was indeed easy. The PDG of the target nucleus now..
   fPDGNucleus = 0;
//   if ((aTrack->GetCreatorProcess() != 0) && (aTrack->GetCreatorProcess()->GetProcessType() ==  fHadronic)) {
//      std::cerr << " Hadronic interaction, what do I do to get the nucleus that got hit?!!!" << std::endl;
//     exit(2);
// The goal is dig through the creator process ... but the info is lost.. 
// 
//       const G4HadronicProcess *hProc = static_cast<const G4HadronicProcess* >(aTrack->GetCreatorProcess());
//       std::cerr << " Current interaction length " << hProc->GetCurrentInteractionLength() 
//                 << " name " << hProc->GetProcessName() << " Process sub type " << hProc->GetProcessSubType() << std::endl;
//       G4IsoParticleChange *ic =  G4HadronicProcess::GetIsotopeProductionInfo();
//       if (ic !=0) {
//         std::cerr << " .... No Isotope production info.. " << std::endl;
//       } else { 
//           std::cerr << " .... Isotope is " << ic->GetIsotope() << std::endl; 
// This access actually causes a crash.. So, the trail is cold regarding getting info about the isotope in questions..
//  Paul Lebrun, April 2014. 
//       }
//   }
   
}

LBNETrajectory::LBNETrajectory

(

LBNETrajectory &

right

)

Definition at line 137 of file LBNETrajectory.cc.

  : G4VTrajectory()
{
  fParticleName       = right.fParticleName;
  fParticleDefinition = right.fParticleDefinition;
  fPDGCharge          = right.fPDGCharge;
  fPDGEncoding        = right.fPDGEncoding;
  fTrackID            = right.fTrackID;
  fParentID           = right.fParentID;
  fMomentum           = right.fMomentum;
  fVertexPosition     = right.fVertexPosition;
  fParticleMass       = right.fParticleMass;
  fDecayCode          = right.fDecayCode;
  fTgen               = right.fTgen;
  fNImpWt             = right.fNImpWt;
  fPositionRecord     = new LBNETrajectoryPointContainer();
  fMomentumRecord     = new LBNETrajectoryMomentumContainer();
  fPreStepVolume      = new LBNETrajectoryVolumeName();  
  fStepLength         = new DVec();
  fParentMomentumAtThisProduction = right.fParentMomentumAtThisProduction;
  fMaterialName       = new LBNETrajectoryMaterialName();  

  for(size_t i=0;i<right.fPositionRecord->size();i++)
  {
     G4TrajectoryPoint* rightPoint = (G4TrajectoryPoint*)((*(right.fPositionRecord))[i]);
     fPositionRecord->push_back(new G4TrajectoryPoint(*rightPoint));
  }
  for(size_t i=0;i<right.fMomentumRecord->size();i++)
  {
     G4ThreeVector rightMomentum = (G4ThreeVector)((*(right.fMomentumRecord))[i]);
     fMomentumRecord->push_back(rightMomentum);
  }
  for(size_t i=0;i<right.fPreStepVolume->size();i++)
  {
     G4String rightPreStepVolume=(G4String)((*(right.fPreStepVolume))[i]);
     fPreStepVolume->push_back(rightPreStepVolume);
  }
  for(size_t i=0;i<right.fStepLength->size();i++)
  {
     G4double rightsteplength =(G4double)((*(right.fStepLength))[i]);
     fStepLength->push_back(rightsteplength);
  }
    for(size_t i=0;i<right.fMaterialName->size();i++)
  {
    G4String rightMaterialName = (G4String)((*(right.fMaterialName))[i]);
      fMaterialName->push_back(rightMaterialName);
  }
  fMaterialNumber1rst = right.fMaterialNumber1rst;
  fVolName1rst        = right.fVolName1rst;
  fMaterialNumberLast = right.fMaterialNumberLast;
  fTimeStart = right.fTimeStart;
  fPolarization       = right.fPolarization;
  fPDGNucleus = right.fPDGNucleus;
  fProcessName = right.fProcessName;
  fMaterialName1rst = right.fMaterialName1rst;
  fMaterialName     = right.fMaterialName;
  

}

LBNETrajectory::LBNETrajectory

(

const LBNETrajectory *

right

)

Definition at line 198 of file LBNETrajectory.cc.

  : G4VTrajectory()
{
  fParticleName       = right -> fParticleName;
  fParticleDefinition = right -> fParticleDefinition;
  fPDGCharge          = right -> fPDGCharge;
  fPDGEncoding        = right -> fPDGEncoding;
  fTrackID            = right -> fTrackID;
  fParentID           = right -> fParentID;
  fMomentum           = right -> fMomentum;
  fVertexPosition     = right -> fVertexPosition;
  fParticleMass       = right -> fParticleMass;
  fDecayCode          = right -> fDecayCode;
  fTgen               = right -> fTgen;
  fNImpWt             = right -> fNImpWt;
  fParentMomentumAtThisProduction = right->fParentMomentumAtThisProduction;
  fMaterialName       = new LBNETrajectoryMaterialName();
  

  fPositionRecord     = new LBNETrajectoryPointContainer();
  for(size_t i=0; i < (right->fPositionRecord)->size(); ++i)
  {
     G4TrajectoryPoint* rightPoint = (G4TrajectoryPoint*)((*(right->fPositionRecord))[i]);
     fPositionRecord->push_back(new G4TrajectoryPoint(*rightPoint));
  }

  fMomentumRecord = new LBNETrajectoryMomentumContainer();
  for(size_t i=0; i < (right->fMomentumRecord)->size(); ++i)
  {
     G4ThreeVector rightMomentum = (G4ThreeVector)((*(right->fMomentumRecord))[i]);
     fMomentumRecord->push_back(rightMomentum);
  }

  fPreStepVolume = new LBNETrajectoryVolumeName();  
  for(size_t i=0; i < (right->fPreStepVolume)->size(); ++i)
  {
     G4String rightPreStepVolume =(G4String)((*(right->fPreStepVolume))[i]);
     fPreStepVolume->push_back(rightPreStepVolume);
  }

  fStepLength = new DVec();  
  for(size_t i=0; i < (right->fStepLength)->size(); ++i)
  {
     G4double rightsteplength =(G4double)((*(right->fStepLength))[i]);
     fStepLength->push_back(rightsteplength);
  }
  for(size_t i=0;i<right->fMaterialName->size();i++)
  {
      G4String rightMaterialName = (G4String)((*(right->fMaterialName))[i]);
      fMaterialName->push_back(rightMaterialName);
  }
  
  fMaterialNumber1rst = right->fMaterialNumber1rst;
  fVolName1rst        = right->fVolName1rst;
  fMaterialNumberLast = right->fMaterialNumberLast;
  fTimeStart = right->fTimeStart;
  fPolarization       = right->fPolarization;
  fPDGNucleus = right->fPDGNucleus;
  fProcessName = right->fProcessName;
  fMaterialName1rst = right->fMaterialName1rst;
  fMaterialName  = right->fMaterialName;
  
  
}

LBNETrajectory::~LBNETrajectory ( )

virtual

Definition at line 264 of file LBNETrajectory.cc.

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


  fMomentumRecord->clear();  
  delete fMomentumRecord;

  fPreStepVolume->clear();
  delete fPreStepVolume;

  fStepLength->clear();
  delete fStepLength;
  
  fMaterialName->clear();
  delete fMaterialName;
 
}

Member Function Documentation

void LBNETrajectory::AppendStep ( const G4Step *  aStep )

virtual

Definition at line 484 of file LBNETrajectory.cc.

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

   G4Track* aTrack=aStep->GetTrack();
   LBNETrackInformation* info=(LBNETrackInformation*)(aTrack->GetUserInformation());
   if (info!=0) {
     fDecayCode=info->GetDecayCode();
     fParentMomentumAtThisProduction = info->GetParentMomentumAtThisProduction();     
     fTgen=info->GetTgen();
   }
   else {fDecayCode=-1;
        fParentMomentumAtThisProduction = G4ThreeVector(0.,0.,0.);}  
   G4StepPoint * steppoint=aStep->GetPreStepPoint(); 
   G4String PreVolumeName=steppoint->GetPhysicalVolume()->GetName(); 
   G4String MaterialName =
   steppoint->GetPhysicalVolume()->GetLogicalVolume()->GetMaterial()->GetName();
   fPreStepVolume->push_back(PreVolumeName); 
   fStepLength->push_back(aStep->GetStepLength());
   fMaterialName->push_back(MaterialName);
}

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

virtual

Definition at line 452 of file LBNETrajectory.cc.

{
   std::vector<G4AttValue>* values = new std::vector<G4AttValue>;
   
   values->push_back
      (G4AttValue("TrkID",G4UIcommand::ConvertToString(fTrackID),""));
   
   values->push_back
      (G4AttValue("ParTrkID",G4UIcommand::ConvertToString(fParentID),""));
   
   values->push_back(G4AttValue("Name",fParticleName,""));
   
   values->push_back
      (G4AttValue("Q",G4UIcommand::ConvertToString(fPDGCharge),""));
   
   values->push_back
      (G4AttValue("PDG",G4UIcommand::ConvertToString(fPDGEncoding),""));
   
   values->push_back
      (G4AttValue("IMom",G4BestUnit(GetInitialMomentum(),"Energy"),""));
   
   values->push_back
      (G4AttValue("IMomMag",G4BestUnit((GetInitialMomentum()).mag(),"Energy"),""));
   
   values->push_back
      (G4AttValue("NPts",G4UIcommand::ConvertToString(GetPointEntries()),""));
   
   return values;
}

void LBNETrajectory::DrawTrajectory ( ) const

virtual

Definition at line 314 of file LBNETrajectory.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 color(1.0,1.0,1.0); //white

   if(fPDGCharge ==  0.)   color = G4Colour::White();
   if(fPDGCharge ==  1.0)  color = G4Colour(1.0, 0.5, 0.0); //orange
   if(fPDGCharge == -1.0)  color = G4Colour(0.5, 0.0, 0.5); //purple


   if (fParticleDefinition==G4Proton::ProtonDefinition()) color = G4Color::Green();
   if (fParticleDefinition==G4Neutron::NeutronDefinition()) color = G4Color::Grey();
   if (fParticleDefinition==G4PionMinus::PionMinusDefinition()) color = G4Color::Cyan();
   if (fParticleDefinition==G4PionPlus::PionPlusDefinition()) color = G4Color::Blue();
   if (fParticleDefinition==G4KaonMinus::KaonMinusDefinition()) color = G4Color::Yellow();
   if (fParticleDefinition==G4KaonPlus::KaonPlusDefinition()) color = G4Color::Yellow();
   if (fParticleDefinition==G4MuonMinus::MuonMinusDefinition()) color = G4Color::Red();
   if (fParticleDefinition==G4MuonPlus::MuonPlusDefinition()) color = G4Color::Red();

   if ((fParticleDefinition == G4NeutrinoE::NeutrinoEDefinition())  ||
       (fParticleDefinition == G4NeutrinoMu::NeutrinoMuDefinition()) ||
       (fParticleDefinition == G4NeutrinoTau::NeutrinoTauDefinition()) ||
       (fParticleDefinition == G4AntiNeutrinoE::AntiNeutrinoEDefinition()) ||
       (fParticleDefinition == G4AntiNeutrinoMu::AntiNeutrinoMuDefinition()) ||
       (fParticleDefinition == G4AntiNeutrinoTau::AntiNeutrinoTauDefinition()))
      color = G4Color::Magenta();







G4VisAttributes attribs(color);
   pPolyline.SetVisAttributes(attribs);
   if(pVVisManager) pVVisManager->Draw(pPolyline);
      
/*
   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);
   }

   //

   //draw only protons,pi+ and pi-
   G4VisAttributes attribs;
   if (fParticleDefinition==G4Proton::ProtonDefinition()) {
     colour=G4Colour(0.,0.,1.);
     attribs=G4VisAttributes(colour);
     pPolyline.SetVisAttributes(attribs);
     if(pVVisManager) pVVisManager->Draw(pPolyline);
   }
   if (fParticleDefinition==G4PionMinus::PionMinusDefinition()) {
     colour=G4Colour(1.,0.,0.);
     attribs=G4VisAttributes(colour);
     pPolyline.SetVisAttributes(attribs);
     if(pVVisManager) pVVisManager->Draw(pPolyline);
   }
   if (fParticleDefinition==G4PionPlus::PionPlusDefinition()) {
     colour=G4Colour(0.,1.,0.);
     attribs=G4VisAttributes(colour);
     pPolyline.SetVisAttributes(attribs);
     if(pVVisManager) pVVisManager->Draw(pPolyline);
   }

*/

}

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

virtual

Definition at line 413 of file LBNETrajectory.cc.

{
    G4bool isNew;
    std::map<G4String,G4AttDef>* store
       = G4AttDefStore::GetInstance("LBNETrajectory",isNew);
    if (isNew) 
    {
       
       G4String ID("TrkID");
       (*store)[ID] = G4AttDef(ID,"Track ID","Physics","","G4int");
       
       G4String PID("ParTrkID");
       (*store)[PID] = G4AttDef(PID,"Parent ID","Physics","","G4int");
 
       G4String PN("Name");
       (*store)[PN] = G4AttDef(PN,"Particle Name","Physics","","G4String");
       
       G4String Ch("Q");
       (*store)[Ch] = G4AttDef(Ch,"Charge","Physics","e+","G4double");
       
       G4String PDG("PDG");
       (*store)[PDG] = G4AttDef(PDG,"PDG Encoding","Physics","","G4int");
       
       G4String IMom("IMom");
       (*store)[IMom] = G4AttDef(IMom, "Momentum at start of trajectory",
                                 "Physics","G4BestUnit","G4ThreeVector");
       
       G4String IMag("IMomMag");
       (*store)[IMag] = 
          G4AttDef(IMag, "Magnitude of momentum at start of trajectory",
                   "Physics","G4BestUnit","G4double");
       
       G4String NTP("NPts");
       (*store)[NTP] = G4AttDef(NTP,"No. of points","Physics","","G4int");
       
    }
    return store;
}

G4double LBNETrajectory::GetCharge ( ) const

inline

Definition at line 52 of file LBNETrajectory.hh.

   { return fPDGCharge; }

G4int LBNETrajectory::GetDecayCode ( ) const

inline

Definition at line 72 of file LBNETrajectory.hh.

   { return fDecayCode;}

G4ThreeVector LBNETrajectory::GetInitialMomentum ( ) const

inline

Definition at line 68 of file LBNETrajectory.hh.

   { return fMomentum; }                            

G4double LBNETrajectory::GetMass ( ) const

inline

Definition at line 54 of file LBNETrajectory.hh.

   { return fParticleMass; }

virtual G4String LBNETrajectory::GetMaterialName ( G4int  i ) const

inlinevirtual

Definition at line 106 of file LBNETrajectory.hh.

   {return (*fMaterialName)[i];}

G4String LBNETrajectory::GetMaterialName1rst ( ) const

inline

Definition at line 104 of file LBNETrajectory.hh.

   { return fMaterialName1rst; }

G4int LBNETrajectory::GetMaterialNumber1rst ( ) const

inline

Definition at line 79 of file LBNETrajectory.hh.

   {return fMaterialNumber1rst;}

G4int LBNETrajectory::GetMaterialNumberLast ( ) const

inline

Definition at line 84 of file LBNETrajectory.hh.

   {return fMaterialNumberLast;}

virtual G4ThreeVector LBNETrajectory::GetMomentum ( G4int  i ) const

inlinevirtual

Definition at line 64 of file LBNETrajectory.hh.

   { return (*fMomentumRecord)[i]; }

virtual G4double LBNETrajectory::GetNImpWt ( ) const

inlinevirtual

Definition at line 74 of file LBNETrajectory.hh.

   { return fNImpWt;}

G4int LBNETrajectory::GetParentID ( ) const

inline

Definition at line 48 of file LBNETrajectory.hh.

   { return fParentID; }

const G4ThreeVector LBNETrajectory::GetParentMomentumAtThisProduction ( ) const

inline

Definition at line 108 of file LBNETrajectory.hh.

                                                                        {
        return fParentMomentumAtThisProduction;
    }  

G4ParticleDefinition * LBNETrajectory::GetParticleDefinition

(

)

Definition at line 509 of file LBNETrajectory.cc.

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

G4String LBNETrajectory::GetParticleName ( ) const

inline

Definition at line 50 of file LBNETrajectory.hh.

   { return fParticleName; }

G4int LBNETrajectory::GetPDGEncoding ( ) const

inline

Definition at line 56 of file LBNETrajectory.hh.

   { return fPDGEncoding; }

G4int LBNETrajectory::GetPDGNucleus ( ) const

inline

Definition at line 98 of file LBNETrajectory.hh.

                                      {
    return fPDGNucleus; }

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

inlinevirtual

Definition at line 62 of file LBNETrajectory.hh.

   { return (*fPositionRecord)[i]; }

virtual int LBNETrajectory::GetPointEntries ( ) const

inlinevirtual

Definition at line 60 of file LBNETrajectory.hh.

   { return fPositionRecord->size(); }

G4ThreeVector LBNETrajectory::GetPolarization ( ) const

inline

Definition at line 95 of file LBNETrajectory.hh.

   { return fPolarization;}

virtual G4String LBNETrajectory::GetPreStepVolumeName ( G4int  i ) const

inlinevirtual

Definition at line 66 of file LBNETrajectory.hh.

   { return (*fPreStepVolume)[i]; }

G4String LBNETrajectory::GetProcessName ( ) const

inline

Definition at line 101 of file LBNETrajectory.hh.

   { return fProcessName; } 

virtual G4double LBNETrajectory::GetStepLength ( G4int  i ) const

inlinevirtual

Definition at line 76 of file LBNETrajectory.hh.

      {return (*fStepLength)[i];}

virtual G4int LBNETrajectory::GetTgen ( ) const

inlinevirtual

Definition at line 70 of file LBNETrajectory.hh.

   { return fTgen;}

G4double LBNETrajectory::GetTimeStart ( ) const

inline

Definition at line 92 of file LBNETrajectory.hh.

   {return fTimeStart;}

G4int LBNETrajectory::GetTrackID ( ) const

inline

Definition at line 46 of file LBNETrajectory.hh.

   { return fTrackID; }

const G4ThreeVector& LBNETrajectory::GetVertexPosition ( ) const

inline

Definition at line 58 of file LBNETrajectory.hh.

   { return fVertexPosition; }

G4String LBNETrajectory::GetVolName1rst ( ) const

inline

Definition at line 89 of file LBNETrajectory.hh.

   {return fVolName1rst;}

void LBNETrajectory::MergeTrajectory ( G4VTrajectory *  secondTrajectory )

virtual

Definition at line 514 of file LBNETrajectory.cc.

{
  if(!secondTrajectory) return;

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

}

void LBNETrajectory::operator delete ( void *  aTrajectory )

inline

Definition at line 165 of file LBNETrajectory.hh.

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

void * LBNETrajectory::operator new ( size_t  )

inline

Definition at line 158 of file LBNETrajectory.hh.

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

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

inline

Definition at line 43 of file LBNETrajectory.hh.

   {return (this==&right);} 

void LBNETrajectory::SetMaterialNumber1rst ( G4double  matNum )

inline

Definition at line 81 of file LBNETrajectory.hh.

   {fMaterialNumber1rst = matNum;}

void LBNETrajectory::SetMaterialNumberLast ( G4double  matNum )

inline

Definition at line 86 of file LBNETrajectory.hh.

    {fMaterialNumberLast = matNum;}

void LBNETrajectory::ShowTrajectory ( ) const

virtual

Definition at line 290 of file LBNETrajectory.cc.

{
   G4cout << G4endl << "TrackID =" << fTrackID 
        << " : ParentID=" << fParentID << G4endl;
   G4cout << "Particle name : " << fParticleName 
        << "  Charge : " << fPDGCharge << G4endl;
   G4cout << "Original momentum : " <<
G4BestUnit(fMomentum,"Energy") << G4endl;
   G4cout << "Vertex : " << G4BestUnit(fVertexPosition,"Length") << G4endl;
   G4cout << "  Current trajectory has " << fPositionRecord->size() 
        << " points." << G4endl;

   for( size_t i=0 ; i < fPositionRecord->size() ; i++){
       G4TrajectoryPoint* aTrajectoryPoint = (G4TrajectoryPoint*)((*fPositionRecord)[i]);
       G4cout << "Point[" << i << "]" 
            << " Position= " << aTrajectoryPoint->GetPosition() << G4endl;
   }
}

void LBNETrajectory::ShowTrajectory ( std::ostream &  o ) const

virtual

Definition at line 309 of file LBNETrajectory.cc.

{
    G4VTrajectory::ShowTrajectory(o);
}

Member Data Documentation

G4int LBNETrajectory::fDecayCode

private

Definition at line 126 of file LBNETrajectory.hh.

G4int LBNETrajectory::fEventNo

private

Definition at line 127 of file LBNETrajectory.hh.

LBNETrajectoryMaterialName* LBNETrajectory::fMaterialName

private

Definition at line 151 of file LBNETrajectory.hh.

G4String LBNETrajectory::fMaterialName1rst

private

Definition at line 152 of file LBNETrajectory.hh.

G4int LBNETrajectory::fMaterialNumber1rst

private

Definition at line 142 of file LBNETrajectory.hh.

G4int LBNETrajectory::fMaterialNumberLast

private

Definition at line 144 of file LBNETrajectory.hh.

G4ThreeVector LBNETrajectory::fMomentum

private

Definition at line 136 of file LBNETrajectory.hh.

LBNETrajectoryMomentumContainer* LBNETrajectory::fMomentumRecord

private

Definition at line 124 of file LBNETrajectory.hh.

G4double LBNETrajectory::fNImpWt

private

Definition at line 129 of file LBNETrajectory.hh.

G4int LBNETrajectory::fParentID

private

Definition at line 131 of file LBNETrajectory.hh.

G4ThreeVector LBNETrajectory::fParentMomentumAtThisProduction

private

Definition at line 150 of file LBNETrajectory.hh.

G4ParticleDefinition* LBNETrajectory::fParticleDefinition

private

Definition at line 132 of file LBNETrajectory.hh.

G4double LBNETrajectory::fParticleMass

private

Definition at line 138 of file LBNETrajectory.hh.

G4String LBNETrajectory::fParticleName

private

Definition at line 133 of file LBNETrajectory.hh.

G4double LBNETrajectory::fPDGCharge

private

Definition at line 134 of file LBNETrajectory.hh.

G4int LBNETrajectory::fPDGEncoding

private

Definition at line 135 of file LBNETrajectory.hh.

G4int LBNETrajectory::fPDGNucleus

private

Definition at line 147 of file LBNETrajectory.hh.

G4ThreeVector LBNETrajectory::fPolarization

private

Definition at line 146 of file LBNETrajectory.hh.

LBNETrajectoryPointContainer* LBNETrajectory::fPositionRecord

private

Definition at line 125 of file LBNETrajectory.hh.

LBNETrajectoryVolumeName* LBNETrajectory::fPreStepVolume

private

Definition at line 139 of file LBNETrajectory.hh.

G4String LBNETrajectory::fProcessName

private

Definition at line 149 of file LBNETrajectory.hh.

DVec* LBNETrajectory::fStepLength

private

Definition at line 140 of file LBNETrajectory.hh.

G4int LBNETrajectory::fTgen

private

Definition at line 128 of file LBNETrajectory.hh.

G4double LBNETrajectory::fTimeStart

private

Definition at line 145 of file LBNETrajectory.hh.

G4int LBNETrajectory::fTrackID

private

Definition at line 130 of file LBNETrajectory.hh.

G4ThreeVector LBNETrajectory::fVertexPosition

private

Definition at line 137 of file LBNETrajectory.hh.

G4String LBNETrajectory::fVolName1rst

private

Definition at line 143 of file LBNETrajectory.hh.

---

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

• g4lbne/include/LBNETrajectory.hh
• g4lbne/src/LBNETrajectory.cc