PrimaryGeneratorAction2 Class Reference

#include <PrimaryGeneratorAction2.hh>

Public Member Functions
	PrimaryGeneratorAction2 (G4ParticleGun *)
	~PrimaryGeneratorAction2 ()
void	GeneratePrimaries (G4Event *)
G4double	RejectAccept ()
G4double	InverseCumul ()

Detailed Description

Definition at line 47 of file PrimaryGeneratorAction2.hh.

Constructor & Destructor Documentation

PrimaryGeneratorAction2::PrimaryGeneratorAction2

(

G4ParticleGun *

gun

)

Definition at line 48 of file PrimaryGeneratorAction2.cc.

: particleGun(gun)
{    
  // energy distribution
  //
  InitFunction();
}

PrimaryGeneratorAction2::~PrimaryGeneratorAction2

(

)

Definition at line 58 of file PrimaryGeneratorAction2.cc.

{ }

Member Function Documentation

void PrimaryGeneratorAction2::GeneratePrimaries

(

G4Event *

anEvent

)

Definition at line 63 of file PrimaryGeneratorAction2.cc.

References energy(), G4UniformRand, G4ParticleGun::GeneratePrimaryVertex(), InverseCumul(), G4ParticleGun::SetParticleEnergy(), G4ParticleGun::SetParticleMomentumDirection(), and python.hepunit::twopi.

{
  //cosAlpha uniform in [cos(0), cos(pi)]
  G4double cosAlpha = 1. - 2*G4UniformRand();
  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);

  particleGun->SetParticleMomentumDirection(dir);
  
  //set energy from a tabulated distribution
  //
  //G4double energy = RejectAccept();
  G4double energy = InverseCumul();  
  particleGun->SetParticleEnergy(energy);    

  //create vertex
  //   
  particleGun->GeneratePrimaryVertex(anEvent);
}

G4double PrimaryGeneratorAction2::InverseCumul

(

)

Definition at line 148 of file PrimaryGeneratorAction2.cc.

References test::b, test::c, G4UniformRand, and G4INCL::Math::sign().

Referenced by GeneratePrimaries().

{
  // tabulated function
  // f is assumed positive, linear per segment, continuous 
  // --> cumulative function is second order polynomial
  
  //choose y randomly
  G4double y_rndm = G4UniformRand()*Fc[nPoints-1];
  //find bin
  G4int j = nPoints-2;
  while ((Fc[j] > y_rndm) && (j > 0)) j--;
  //y_rndm --> x_rndm :  Fc(x) is second order polynomial
  G4double x_rndm = x[j];
  G4double aa = a[j];
  if (aa != 0.) {
    G4double b = f[j]/aa, c = 2*(y_rndm - Fc[j])/aa;
    G4double delta = b*b + c;
    G4int sign = 1; if (aa < 0.) sign = -1;
    x_rndm += sign*std::sqrt(delta) - b;    
  } else if (f[j] > 0.) {
    x_rndm += (y_rndm - Fc[j])/f[j];
  };
  return x_rndm;
}

G4double PrimaryGeneratorAction2::RejectAccept

(

)

Definition at line 126 of file PrimaryGeneratorAction2.cc.

References G4UniformRand.

{
  // tabulated function 
  // f is assumed positive, linear per segment, continuous
  //  
  G4double x_rndm = 0., y_rndm = 0., f_inter = -1.;
  
  while (y_rndm > f_inter) {
    //choose a point randomly
    x_rndm = x[0] + G4UniformRand()*(x[nPoints-1] - x[0]);
    y_rndm = G4UniformRand()*fMax;
    //find bin
    G4int j = nPoints-2;
    while ((x[j] > x_rndm) && (j > 0)) j--;
    //compute f(x_rndm) by linear interpolation
    f_inter = f[j] + a[j]*(x_rndm - x[j]);
  };
  return x_rndm;
}

---

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

• PrimaryGeneratorAction2.hh
• PrimaryGeneratorAction2.cc