#include <B5PrimaryGeneratorAction.hh>

Inheritance diagram for B5PrimaryGeneratorAction:

Collaboration diagram for B5PrimaryGeneratorAction:

Public Member Functions
	B5PrimaryGeneratorAction ()

virtual	~B5PrimaryGeneratorAction ()

virtual void	GeneratePrimaries (G4Event *)

void	SetMomentum (G4double val)

G4double	GetMomentum () const

void	SetSigmaMomentum (G4double val)

G4double	GetSigmaMomentum () const

void	SetSigmaAngle (G4double val)

G4double	GetSigmaAngle () const

void	SetRandomize (G4bool val)

G4bool	GetRandomize () const

Public Member Functions inherited from G4VUserPrimaryGeneratorAction
	G4VUserPrimaryGeneratorAction ()

virtual	~G4VUserPrimaryGeneratorAction ()

Detailed Description

Primary generator

A single particle is generated. User can select

the initial momentum and angle
the momentum and angle spreads
random selection of a particle type from proton, kaon+, pi+, muon+, e+

Definition at line 51 of file B5PrimaryGeneratorAction.hh.

Constructor & Destructor Documentation

B5PrimaryGeneratorAction::B5PrimaryGeneratorAction ( )

Definition at line 43 of file B5PrimaryGeneratorAction.cc.

 : G4VUserPrimaryGeneratorAction(),     
   fParticleGun(nullptr), fMessenger(nullptr), 
   fPositron(nullptr), fMuon(nullptr), fPion(nullptr), 
   fKaon(nullptr), fProton(nullptr),
   fMomentum(1000.*MeV),
   fSigmaMomentum(50.*MeV),
   fSigmaAngle(2.*deg),
   fRandomizePrimary(true)
 {
   G4int nofParticles = 1;
   fParticleGun  = new G4ParticleGun(nofParticles);
   
   auto particleTable = G4ParticleTable::GetParticleTable();
   fPositron = particleTable->FindParticle("e+");
   fMuon = particleTable->FindParticle("mu+");
   fPion = particleTable->FindParticle("pi+");
   fKaon = particleTable->FindParticle("kaon+");
   fProton = particleTable->FindParticle("proton");
   
   // default particle kinematics
   fParticleGun->SetParticlePosition(G4ThreeVector(0.,0.,-8.*m));
   fParticleGun->SetParticleDefinition(fPositron);
   
   // define commands for this class
   DefineCommands();
 }

Here is the call graph for this function:

B5PrimaryGeneratorAction::~B5PrimaryGeneratorAction ( )

virtual

Definition at line 73 of file B5PrimaryGeneratorAction.cc.

 {
   delete fParticleGun;
   delete fMessenger;
 }

Member Function Documentation

void B5PrimaryGeneratorAction::GeneratePrimaries ( G4Event * event )

virtual

Implements G4VUserPrimaryGeneratorAction.

Definition at line 81 of file B5PrimaryGeneratorAction.cc.

 {
   G4ParticleDefinition* particle;  
   if (fRandomizePrimary) {
     G4int i = (int)(5.*G4UniformRand());
     switch(i) {
       case 0:
           particle = fPositron;
           break;
       case 1:
           particle = fMuon;
           break;
       case 2:
           particle = fPion;
           break;
       case 3:
           particle = fKaon;
           break;
       default:
           particle = fProton;
           break;
     }
     fParticleGun->SetParticleDefinition(particle);
   }
   else {
     particle = fParticleGun->GetParticleDefinition();
   }
   
   auto pp = fMomentum + (G4UniformRand()-0.5)*fSigmaMomentum;
   auto mass = particle->GetPDGMass();
   auto ekin = std::sqrt(pp*pp+mass*mass)-mass;
   fParticleGun->SetParticleEnergy(ekin);
   
   auto angle = (G4UniformRand()-0.5)*fSigmaAngle;
   fParticleGun->SetParticleMomentumDirection(
                   G4ThreeVector(std::sin(angle),0.,std::cos(angle)));
   
   fParticleGun->GeneratePrimaryVertex(event);
 }