The primary generator action class with particle gun. More...

#include <Doxymodules_visualization.h>

Inheritance diagram for VtkVis::PrimaryGeneratorAction:

Public Member Functions
	PrimaryGeneratorAction ()

	~PrimaryGeneratorAction () override

void	GeneratePrimaries (G4Event *) override

const G4ParticleGun *	GetParticleGun () const

Private Attributes
G4ParticleGun *	fParticleGun = nullptr

G4Box *	fEnvelopeBox = nullptr

Detailed Description

The primary generator action class with particle gun.

The default kinematic is a 6 MeV gamma, randomly distribued in front of the phantom across 80% of the (X,Y) phantom size.

Definition at line 52 of file Doxymodules_visualization.h.

Constructor & Destructor Documentation

◆ PrimaryGeneratorAction()

VtkVis::PrimaryGeneratorAction::PrimaryGeneratorAction ( )

Definition at line 47 of file PrimaryGeneratorAction.cc.

{
  G4int n_particle = 1;
  fParticleGun  = new G4ParticleGun(n_particle);
 
  // default particle kinematic
  G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
  G4String particleName;
  G4ParticleDefinition* particle
    = particleTable->FindParticle(particleName="e-");
  fParticleGun->SetParticleDefinition(particle);
  fParticleGun->SetParticleMomentumDirection(G4ThreeVector(0.,0.,1.));
  fParticleGun->SetParticleEnergy(16.5*GeV);
}

◆ ~PrimaryGeneratorAction()

VtkVis::PrimaryGeneratorAction::~PrimaryGeneratorAction ( )

override

Definition at line 64 of file PrimaryGeneratorAction.cc.

{
  delete fParticleGun;
}

Member Function Documentation

◆ GeneratePrimaries()

void VtkVis::PrimaryGeneratorAction::GeneratePrimaries ( G4Event * anEvent )

override

Definition at line 71 of file PrimaryGeneratorAction.cc.

{
  //this function is called at the begining of ecah event
  //
 
  // In order to avoid dependence of PrimaryGeneratorAction
  // on DetectorConstruction class we get Envelope volume
  // from G4LogicalVolumeStore.
 
  G4double envSizeXY = 0;
  G4double envSizeZ = 0;
 
  if (!fEnvelopeBox)
  {
    G4LogicalVolume* envLV
      = G4LogicalVolumeStore::GetInstance()->GetVolume("Envelope");
    if ( envLV ) fEnvelopeBox = dynamic_cast<G4Box*>(envLV->GetSolid());
  }
 
  if ( fEnvelopeBox ) {
    envSizeXY = fEnvelopeBox->GetXHalfLength()*2.;
    envSizeZ = fEnvelopeBox->GetZHalfLength()*2.;
  }
  else  {
    G4ExceptionDescription msg;
    msg << "Envelope volume of box shape not found.\n";
    msg << "Perhaps you have changed geometry.\n";
    msg << "The gun will be place at the center.";
    G4Exception("PrimaryGeneratorAction::GeneratePrimaries()",
     "MyCode0002",JustWarning,msg);
  }
 
  G4double size = 0.8;
  G4double x0 = size * envSizeXY * (G4UniformRand()-0.5);
  G4double y0 = size * envSizeXY * (G4UniformRand()-0.5);
  G4double z0 = 0.5*envSizeZ;
 
  fParticleGun->SetParticlePosition(G4ThreeVector(0,0,-z0));
 
  fParticleGun->GeneratePrimaryVertex(anEvent);
}

◆ GetParticleGun()

const G4ParticleGun * VtkVis::PrimaryGeneratorAction::GetParticleGun ( ) const

inline

Definition at line 59 of file PrimaryGeneratorAction.hh.

59{ return fParticleGun; }

Member Data Documentation

◆ fParticleGun

G4ParticleGun* VtkVis::PrimaryGeneratorAction::fParticleGun = nullptr

private

Definition at line 62 of file PrimaryGeneratorAction.hh.

◆ fEnvelopeBox

G4Box* VtkVis::PrimaryGeneratorAction::fEnvelopeBox = nullptr