PrimaryGeneratorAction1 Class Reference

#include <PrimaryGeneratorAction1.hh>

Collaboration diagram for PrimaryGeneratorAction1:

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

Private Attributes
G4ParticleGun *	fParticleGun

Detailed Description

Definition at line 46 of file PrimaryGeneratorAction1.hh.

Constructor & Destructor Documentation

PrimaryGeneratorAction1()

PrimaryGeneratorAction1::PrimaryGeneratorAction1

(

G4ParticleGun *

gun

)

Definition at line 48 of file PrimaryGeneratorAction1.cc.

: fParticleGun(gun)
{ }

~PrimaryGeneratorAction1()

PrimaryGeneratorAction1::~PrimaryGeneratorAction1

(

)

Definition at line 54 of file PrimaryGeneratorAction1.cc.

{ }

Member Function Documentation

GeneratePrimaries()

void PrimaryGeneratorAction1::GeneratePrimaries

(

G4Event *

anEvent

)

Definition at line 59 of file PrimaryGeneratorAction1.cc.

{
  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 3 (same as 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);
  
  // randomize time zero of anEvent
  //
  G4double tmin = 0*s, tmax = 10*s;
  G4double t0 = tmin + (tmax - tmin)*G4UniformRand();
  G4PrimaryVertex* aVertex = anEvent->GetPrimaryVertex();
  while (aVertex) {
    aVertex->SetT0(t0);
    aVertex = aVertex->GetNext();
  }
}

Here is the call graph for this function:

Member Data Documentation

fParticleGun

G4ParticleGun* PrimaryGeneratorAction1::fParticleGun

private

Definition at line 56 of file PrimaryGeneratorAction1.hh.

---

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

• PrimaryGeneratorAction1.hh
• PrimaryGeneratorAction1.cc