PrimaryGeneratorAction1 Class Reference

#include <PrimaryGeneratorAction1.hh>

Public Member Functions
	PrimaryGeneratorAction1 (G4ParticleGun *)
	~PrimaryGeneratorAction1 ()
void	GeneratePrimaries (G4Event *)

Detailed Description

Definition at line 46 of file PrimaryGeneratorAction1.hh.

Constructor & Destructor Documentation

PrimaryGeneratorAction1::PrimaryGeneratorAction1

(

G4ParticleGun *

gun

)

Definition at line 48 of file PrimaryGeneratorAction1.cc.

: fParticleGun(gun)
{ }

PrimaryGeneratorAction1::~PrimaryGeneratorAction1

(

)

Definition at line 54 of file PrimaryGeneratorAction1.cc.

{ }

Member Function Documentation

void PrimaryGeneratorAction1::GeneratePrimaries

(

G4Event *

anEvent

)

Definition at line 59 of file PrimaryGeneratorAction1.cc.

References python.hepunit::deg, G4UniformRand, G4ParticleGun::GeneratePrimaryVertex(), python.hepunit::GeV, python.hepunit::keV, python.hepunit::MeV, python.hepunit::mm, python.hepunit::pi, G4ParticleGun::SetParticleEnergy(), G4ParticleGun::SetParticleMomentumDirection(), G4VPrimaryGenerator::SetParticlePosition(), python.hepunit::twopi, and z.

{
  const G4double r = 2*mm;
  const G4double zmax = 8*mm;   
  
  //vertex 1 uniform on cylinder
  //
  G4double alpha = twopi*G4UniformRand();  //alpha uniform in (0, 2*pi)
  G4double ux = std::cos(alpha);
  G4double uy = std::sin(alpha);
  G4double z = zmax*(2*G4UniformRand() - 1);  //z uniform in (-zmax, +zmax)
        
  fParticleGun->SetParticlePosition(G4ThreeVector(r*ux,r*uy,z));
  fParticleGun->SetParticleMomentumDirection(G4ThreeVector(ux,uy,0));    
  fParticleGun->SetParticleEnergy(1*MeV);
  fParticleGun->GeneratePrimaryVertex(anEvent);
  
  //vertex 2 at opposite
  //
  alpha += pi;
  ux = std::cos(alpha);
  uy = std::sin(alpha);        
  fParticleGun->SetParticlePosition(G4ThreeVector(r*ux,r*uy,z));
  
  //particle 2 at vertex 2
  //
  const G4double dalpha = 10*deg;
  ux = std::cos(alpha + dalpha);
  uy = std::sin(alpha + dalpha);        
  
  fParticleGun->SetParticleMomentumDirection(G4ThreeVector(ux,uy,0));    
  fParticleGun->SetParticleEnergy(1*keV);
  fParticleGun->GeneratePrimaryVertex(anEvent);
  
  //particle 3 at vertex 2
  //
  ux = std::cos(alpha - dalpha);
  uy = std::sin(alpha - dalpha);        
  
  fParticleGun->SetParticleMomentumDirection(G4ThreeVector(ux,uy,0));    
  fParticleGun->SetParticleEnergy(1*GeV);
  fParticleGun->GeneratePrimaryVertex(anEvent);    
}

---

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

• PrimaryGeneratorAction1.hh
• PrimaryGeneratorAction1.cc