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

#include <Doxymodules_analysis.h>

Inheritance diagram for B1PrimaryGeneratorAction:

Public Member Functions
	B1PrimaryGeneratorAction ()

virtual	~B1PrimaryGeneratorAction ()

virtual void	GeneratePrimaries (G4Event *)

const G4ParticleGun *	GetParticleGun () const

Private Attributes
G4ParticleGun *	fParticleGun

G4Box *	fEnvelopeBox

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 50 of file Doxymodules_analysis.h.

Constructor & Destructor Documentation

◆ B1PrimaryGeneratorAction()

B1PrimaryGeneratorAction::B1PrimaryGeneratorAction ( )

Definition at line 44 of file B1PrimaryGeneratorAction.cc.

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

◆ ~B1PrimaryGeneratorAction()

B1PrimaryGeneratorAction::~B1PrimaryGeneratorAction ( )

virtual

Definition at line 64 of file B1PrimaryGeneratorAction.cc.

{
  delete fParticleGun;
}

Member Function Documentation

◆ GeneratePrimaries()

void B1PrimaryGeneratorAction::GeneratePrimaries ( G4Event * anEvent )

virtual

Definition at line 71 of file B1PrimaryGeneratorAction.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("B1PrimaryGeneratorAction::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(x0,y0,z0));
 
  fParticleGun->GeneratePrimaryVertex(anEvent);
}

◆ GetParticleGun()

const G4ParticleGun * B1PrimaryGeneratorAction::GetParticleGun ( ) const

inline

Definition at line 56 of file B1PrimaryGeneratorAction.hh.

56{ return fParticleGun; }

Member Data Documentation

◆ fParticleGun

G4ParticleGun* B1PrimaryGeneratorAction::fParticleGun

private

Definition at line 59 of file B1PrimaryGeneratorAction.hh.

◆ fEnvelopeBox

G4Box* B1PrimaryGeneratorAction::fEnvelopeBox