#include <B5PrimaryGeneratorAction.hh>

Inheritance 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.

References G4ParticleTable::FindParticle(), G4ParticleTable::GetParticleTable(), python.hepunit::m, G4ParticleGun::SetParticleDefinition(), and G4VPrimaryGenerator::SetParticlePosition().

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

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.

References G4UniformRand, G4ParticleGun::GeneratePrimaryVertex(), G4ParticleGun::GetParticleDefinition(), G4ParticleDefinition::GetPDGMass(), int(), G4InuclParticleNames::pp, G4ParticleGun::SetParticleDefinition(), G4ParticleGun::SetParticleEnergy(), and G4ParticleGun::SetParticleMomentumDirection().

 {
     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();
     }
     
     G4double pp = fMomentum + (G4UniformRand()-0.5)*fSigmaMomentum;
     G4double mass = particle->GetPDGMass();
     G4double Ekin = std::sqrt(pp*pp+mass*mass)-mass;
     fParticleGun->SetParticleEnergy(Ekin);
     
     G4double angle = (G4UniformRand()-0.5)*fSigmaAngle;
     fParticleGun->SetParticleMomentumDirection(G4ThreeVector(std::sin(angle),0.,
                                                              std::cos(angle)));
     
     fParticleGun->GeneratePrimaryVertex(event);
 }