GarfieldG4FastSimulationModel Class Reference

#include <GarfieldG4FastSimulationModel.h>

Inheritance diagram for GarfieldG4FastSimulationModel:

Public Member Functions
	GarfieldG4FastSimulationModel (G4String, G4Region *)
	GarfieldG4FastSimulationModel (G4String)
	~GarfieldG4FastSimulationModel ()
void	SetPhysics (GarfieldPhysics *fGarfieldPhysics)
virtual G4bool	IsApplicable (const G4ParticleDefinition &)
virtual G4bool	ModelTrigger (const G4FastTrack &)
virtual void	DoIt (const G4FastTrack &, G4FastStep &)

Private Attributes
GarfieldPhysics *	fGarfieldPhysics

Detailed Description

Definition at line 36 of file GarfieldG4FastSimulationModel.h.

Constructor & Destructor Documentation

GarfieldG4FastSimulationModel() [1/2]

GarfieldG4FastSimulationModel::GarfieldG4FastSimulationModel	(	G4String	modelName,
		G4Region *	envelope )

Definition at line 38 of file GarfieldG4FastSimulationModel.cxx.

  : G4VFastSimulationModel(modelName, envelope)
{
  fGarfieldPhysics = GarfieldPhysics::GetInstance();
  fGarfieldPhysics->InitializePhysics();
}

GarfieldG4FastSimulationModel() [2/2]

GarfieldG4FastSimulationModel::GarfieldG4FastSimulationModel

(

G4String

modelName

)

Definition at line 46 of file GarfieldG4FastSimulationModel.cxx.

  : G4VFastSimulationModel(modelName)
{
  // I.H.
  fGarfieldPhysics = GarfieldPhysics::GetInstance();
  fGarfieldPhysics->InitializePhysics();
}

~GarfieldG4FastSimulationModel()

GarfieldG4FastSimulationModel::~GarfieldG4FastSimulationModel

(

)

Definition at line 54 of file GarfieldG4FastSimulationModel.cxx.

{}

Member Function Documentation

SetPhysics()

void GarfieldG4FastSimulationModel::SetPhysics

(

GarfieldPhysics *

fGarfieldPhysics

)

IsApplicable()

G4bool GarfieldG4FastSimulationModel::IsApplicable ( const G4ParticleDefinition & particleType )

virtual

Definition at line 66 of file GarfieldG4FastSimulationModel.cxx.

{
 
  G4String particleName = particleType.GetParticleName();
 
  if (fGarfieldPhysics->FindParticleName(particleName)) {
    return true;
  }
  return false;
}

ModelTrigger()

G4bool GarfieldG4FastSimulationModel::ModelTrigger ( const G4FastTrack & fastTrack )

virtual

Definition at line 78 of file GarfieldG4FastSimulationModel.cxx.

{
  double ekin_MeV = fastTrack.GetPrimaryTrack()->GetKineticEnergy() / MeV;
  G4String particleName =
    fastTrack.GetPrimaryTrack()->GetParticleDefinition()->GetParticleName();
  if (fGarfieldPhysics->FindParticleNameEnergy(particleName, ekin_MeV)) {
    return true;
  }
  return false;
}

DoIt()

void GarfieldG4FastSimulationModel::DoIt ( const G4FastTrack & fastTrack, G4FastStep & fastStep )

virtual

Definition at line 89 of file GarfieldG4FastSimulationModel.cxx.

{
 
  G4TouchableHandle theTouchable =
    fastTrack.GetPrimaryTrack()->GetTouchableHandle();
  G4String name = theTouchable->GetVolume()->GetName();
 
  G4ThreeVector pdirection = fastTrack.GetPrimaryTrack()->GetMomentum().unit();
  G4ThreeVector localdir = fastTrack.GetPrimaryTrackLocalDirection();
 
  G4ThreeVector worldPosition = fastTrack.GetPrimaryTrack()->GetPosition();
  G4ThreeVector localPosition = fastTrack.GetPrimaryTrackLocalPosition();
 
  double ekin_MeV = fastTrack.GetPrimaryTrack()->GetKineticEnergy() / MeV;
  G4double globalTime = fastTrack.GetPrimaryTrack()->GetGlobalTime();
 
  G4String particleName =
    fastTrack.GetPrimaryTrack()->GetParticleDefinition()->GetParticleName();
 
  fastStep.KillPrimaryTrack();
  fastStep.SetPrimaryTrackPathLength(0.0);
 
  if (particleName == "kaon+") {
    particleName = "K+";
  }
  else if (particleName == "kaon-") {
    particleName = "K-";
  }
  else if (particleName == "anti_proton") {
    particleName = "anti-proton";
  }
 
  fGarfieldPhysics->DoIt(particleName, ekin_MeV, globalTime,
    localPosition.x() / CLHEP::cm, localPosition.y() / CLHEP::cm,
    localPosition.z() / CLHEP::cm, localdir.x(), localdir.y(), localdir.z());
 
  fastStep.SetTotalEnergyDeposited(fGarfieldPhysics->GetEnergyDeposit_MeV());
 
  if (fGarfieldPhysics->GetCreateSecondariesInGeant4()) {
    std::vector<GarfieldParticle*>* secondaryParticles =
      fGarfieldPhysics->GetSecondaryParticles();
 
    if (secondaryParticles->size() > 0) {
      fastStep.SetNumberOfSecondaryTracks(secondaryParticles->size());
 
      G4double totalEnergySecondaries_MeV = 0;
 
      for (std::vector<GarfieldParticle*>::iterator it =
             secondaryParticles->begin();
           it != secondaryParticles->end(); ++it) {
        G4double x = (*it)->getX_mm();
        G4double y = (*it)->getY_mm();
        G4double z = (*it)->getZ_mm();
        G4double eKin_MeV = (*it)->getEkin_MeV();
        G4double dx = (*it)->getDX();
        G4double dy = (*it)->getDY();
        G4double dz = (*it)->getDZ();
        G4double time = (*it)->getTime();
        G4ThreeVector momentumDirection(dx, dy, dz);
        G4ThreeVector position(x, y, z);
        if ((*it)->getParticleName() == "e-") {
          G4DynamicParticle particle(
            G4Electron::ElectronDefinition(), momentumDirection, eKin_MeV);
          fastStep.CreateSecondaryTrack(particle, position, time, true);
          totalEnergySecondaries_MeV += eKin_MeV;
        }
        else if ((*it)->getParticleName() == "gamma") {
          G4DynamicParticle particle(
            G4Gamma::GammaDefinition(), momentumDirection, eKin_MeV);
          fastStep.CreateSecondaryTrack(particle, position, time, true);
          totalEnergySecondaries_MeV += eKin_MeV;
        }
      }
    }
  }
}

Member Data Documentation

fGarfieldPhysics

GarfieldPhysics* GarfieldG4FastSimulationModel::fGarfieldPhysics

private

Definition at line 55 of file GarfieldG4FastSimulationModel.h.

---

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

• examples/ExGarfield/geant4/include/GarfieldG4FastSimulationModel.h
• examples/ExGarfield/geant4/src/GarfieldG4FastSimulationModel.cxx