PrimaryGeneratorAction4 Class Reference

#include <PrimaryGeneratorAction4.hh>

Public Member Functions
	PrimaryGeneratorAction4 (G4ParticleGun *)
	~PrimaryGeneratorAction4 ()
void	GeneratePrimaries (G4Event *)
G4double	GetRmin3 ()
G4double	GetRmax3 ()
G4double	GetCosAlphaMin ()
G4double	GetCosAlphaMax ()

Detailed Description

Definition at line 47 of file PrimaryGeneratorAction4.hh.

Constructor & Destructor Documentation

PrimaryGeneratorAction4::PrimaryGeneratorAction4

(

G4ParticleGun *

gun

)

Definition at line 48 of file PrimaryGeneratorAction4.cc.

References python.hepunit::deg, and python.hepunit::mm.

: fParticleGun(gun)
{
  // vertex volume
  //  
  G4double Rmin = 2.*mm; 
  G4double Rmax = 8.*mm;
  fRmin3 = Rmin*Rmin*Rmin;
  fRmax3 = Rmax*Rmax*Rmax;
  
  //opening angle
  //
  G4double alphaMin =  0.*deg;
  G4double alphaMax = 10.*deg;
  fCosAlphaMin = std::cos(alphaMin);
  fCosAlphaMax = std::cos(alphaMax);  
}

PrimaryGeneratorAction4::~PrimaryGeneratorAction4

(

)

Definition at line 68 of file PrimaryGeneratorAction4.cc.

{ }

Member Function Documentation

void PrimaryGeneratorAction4::GeneratePrimaries

(

G4Event *

anEvent

)

Definition at line 73 of file PrimaryGeneratorAction4.cc.

References G4UniformRand, G4ParticleGun::GeneratePrimaryVertex(), python.hepunit::MeV, CLHEP::Hep3Vector::rotateUz(), G4ParticleGun::SetParticleEnergy(), G4ParticleGun::SetParticleMomentumDirection(), G4VPrimaryGenerator::SetParticlePosition(), and python.hepunit::twopi.

{  
  //vertex position uniform in spherical shell
  //
  G4double cosTheta = 2*G4UniformRand() - 1;  //cosTheta uniform in [0, pi]
  G4double sinTheta = std::sqrt(1. - cosTheta*cosTheta);
  G4double phi      = twopi*G4UniformRand();  //phi uniform in [0, 2*pi]
  G4ThreeVector ur(sinTheta*std::cos(phi),sinTheta*std::sin(phi),cosTheta);
  
  G4double R3 = fRmin3 + G4UniformRand()*(fRmax3 - fRmin3);
  G4double R  = std::pow(R3, 1./3);  
        
  fParticleGun->SetParticlePosition(R*ur);

  //particle direction uniform around ur 
  //    
  //1- in World frame
  //cosAlpha uniform in [cos(alphaMin), cos(alphaMax)]
  G4double cosAlpha = fCosAlphaMin - G4UniformRand()*(fCosAlphaMin - fCosAlphaMax); 
  G4double sinAlpha = std::sqrt(1. - cosAlpha*cosAlpha);
  G4double psi      = twopi*G4UniformRand();  //psi uniform in (0,2*pi)  
  G4ThreeVector dir(sinAlpha*std::cos(psi),sinAlpha*std::sin(psi),cosAlpha);
  
  //2- rotate dir   (rotateUz transforms uz to ur)
  dir.rotateUz(ur);           

  fParticleGun->SetParticleMomentumDirection(dir);
  
  //energy
  //  
  fParticleGun->SetParticleEnergy(1*MeV);
  
  //create vertex
  //   
  fParticleGun->GeneratePrimaryVertex(anEvent);
}

G4double PrimaryGeneratorAction4::GetCosAlphaMax ( )

inline

Definition at line 59 of file PrimaryGeneratorAction4.hh.

{return fCosAlphaMax;};

G4double PrimaryGeneratorAction4::GetCosAlphaMin ( )

inline

Definition at line 58 of file PrimaryGeneratorAction4.hh.

{return fCosAlphaMin;};

G4double PrimaryGeneratorAction4::GetRmax3 ( )

inline

Definition at line 57 of file PrimaryGeneratorAction4.hh.

{return fRmax3;};

G4double PrimaryGeneratorAction4::GetRmin3 ( )

inline

Definition at line 56 of file PrimaryGeneratorAction4.hh.

{return fRmin3;};

---

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

• PrimaryGeneratorAction4.hh
• PrimaryGeneratorAction4.cc