PrimaryGeneratorAction4 Class Reference

#include <PrimaryGeneratorAction4.hh>

Collaboration diagram for PrimaryGeneratorAction4:

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

Private Attributes
G4double	fRmin3
G4double	fRmax3
G4double	fCosAlphaMin
G4double	fCosAlphaMax
G4ParticleGun *	fParticleGun

Detailed Description

Definition at line 47 of file PrimaryGeneratorAction4.hh.

Constructor & Destructor Documentation

PrimaryGeneratorAction4()

PrimaryGeneratorAction4::PrimaryGeneratorAction4

(

G4ParticleGun *

gun

)

Definition at line 48 of file PrimaryGeneratorAction4.cc.

: 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 = 60.*deg;
  fCosAlphaMin = std::cos(alphaMin);
  fCosAlphaMax = std::cos(alphaMax);  
}

~PrimaryGeneratorAction4()

PrimaryGeneratorAction4::~PrimaryGeneratorAction4

(

)

Definition at line 68 of file PrimaryGeneratorAction4.cc.

{ }

Member Function Documentation

GeneratePrimaries()

void PrimaryGeneratorAction4::GeneratePrimaries

(

G4Event *

anEvent

)

Definition at line 73 of file PrimaryGeneratorAction4.cc.

{  
  //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);
}

Here is the call graph for this function:

GetCosAlphaMax()

G4double PrimaryGeneratorAction4::GetCosAlphaMax ( )

inline

Definition at line 59 of file PrimaryGeneratorAction4.hh.

{return fCosAlphaMax;};

GetCosAlphaMin()

G4double PrimaryGeneratorAction4::GetCosAlphaMin ( )

inline

Definition at line 58 of file PrimaryGeneratorAction4.hh.

{return fCosAlphaMin;};

GetRmax3()

G4double PrimaryGeneratorAction4::GetRmax3 ( )

inline

Definition at line 57 of file PrimaryGeneratorAction4.hh.

{return fRmax3;};

GetRmin3()

G4double PrimaryGeneratorAction4::GetRmin3 ( )

inline

Definition at line 56 of file PrimaryGeneratorAction4.hh.

{return fRmin3;};

Member Data Documentation

fCosAlphaMax

G4double PrimaryGeneratorAction4::fCosAlphaMax

private

Definition at line 65 of file PrimaryGeneratorAction4.hh.

fCosAlphaMin

G4double PrimaryGeneratorAction4::fCosAlphaMin

private

Definition at line 64 of file PrimaryGeneratorAction4.hh.

fParticleGun

G4ParticleGun* PrimaryGeneratorAction4::fParticleGun

private

Definition at line 67 of file PrimaryGeneratorAction4.hh.

fRmax3

G4double PrimaryGeneratorAction4::fRmax3

private

Definition at line 63 of file PrimaryGeneratorAction4.hh.

fRmin3

G4double PrimaryGeneratorAction4::fRmin3

private

Definition at line 59 of file PrimaryGeneratorAction4.hh.

---

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

• PrimaryGeneratorAction4.hh
• PrimaryGeneratorAction4.cc