IORTPrimaryGeneratorAction Class Reference

#include <IORTPrimaryGeneratorAction.hh>

Inheritance diagram for IORTPrimaryGeneratorAction:

Collaboration diagram for IORTPrimaryGeneratorAction:

Public Member Functions
	IORTPrimaryGeneratorAction ()
	~IORTPrimaryGeneratorAction ()
void	SetsigmaEnergy (G4double)
void	SetmeanKineticEnergy (G4double)
void	GeneratePrimaries (G4Event *)
void	SetXposition (G4double)
void	SetYposition (G4double)
void	SetZposition (G4double)
void	SetsigmaY (G4double)
void	SetsigmaZ (G4double)
void	SetTheta (G4double)
G4double	GetmeanKineticEnergy (void)
G4ParticleGun *	GetParticleGun (void)
Public Member Functions inherited from G4VUserPrimaryGeneratorAction
	G4VUserPrimaryGeneratorAction ()
virtual	~G4VUserPrimaryGeneratorAction ()

Detailed Description

Definition at line 50 of file IORTPrimaryGeneratorAction.hh.

Constructor & Destructor Documentation

IORTPrimaryGeneratorAction::IORTPrimaryGeneratorAction

(

)

Definition at line 55 of file IORTPrimaryGeneratorAction.cc.

{
  // Define the messenger
  gunMessenger = new IORTPrimaryGeneratorMessenger(this);

  particleGun  = new G4ParticleGun();

  SetDefaultPrimaryParticle();  
}  

IORTPrimaryGeneratorAction::~IORTPrimaryGeneratorAction

(

)

Definition at line 65 of file IORTPrimaryGeneratorAction.cc.

{
  delete particleGun;

  delete gunMessenger;
}

Member Function Documentation

void IORTPrimaryGeneratorAction::GeneratePrimaries ( G4Event *  anEvent )

virtual

Implements G4VUserPrimaryGeneratorAction.

Definition at line 132 of file IORTPrimaryGeneratorAction.cc.

{
  // Increment the event counter
  IORTAnalysisManager::GetInstance()->startNewEvent();

  // ****************************************
  // Set the beam angular apread 
  // and spot size
  // beam spot size
  // ****************************************

  // Set the position of the primary particles
  G4double x = X0;
  G4double y = Y0;
  G4double z = Z0;

  if ( sigmaY > 0.0 )
    {
      y += G4RandGauss::shoot( Y0, sigmaY );
    }
 
  if ( sigmaZ > 0.0 )
    {
      z += G4RandGauss::shoot( Z0, sigmaZ );
    }

  particleGun -> SetParticlePosition(G4ThreeVector( x , y , z ) );
 
  // ********************************************
  // Set the beam energy and energy spread
  // ********************************************

  G4double kineticEnergy = G4RandGauss::shoot( meanKineticEnergy, sigmaEnergy );
  particleGun -> SetParticleEnergy ( kineticEnergy );

  // Set the direction of the primary particles  
  

  /*                          
  G4double momentumX = 1.0;
  G4double momentumY = 0.0;
  G4double momentumZ = 0.0;

  if ( sigmaMomentumY  > 0.0 )
    {
      momentumY += G4RandGauss::shoot( 0., sigmaMomentumY );
    }
  if ( sigmaMomentumZ  > 0.0 )
    {
      momentumZ += G4RandGauss::shoot( 0., sigmaMomentumZ );
    }
 
  particleGun -> SetParticleMomentumDirection( G4ThreeVector(momentumX,momentumY,momentumZ) );

  */
  
               
  G4double Mx;    
  G4double My;
  G4double Mz;
  G4double condizione;
  
while (true)  {

  //Mx =  CLHEP::RandFlat::shoot(0.9,1);
  //My =  CLHEP::RandFlat::shoot(-0.1,0.1);
  //Mz =  CLHEP::RandFlat::shoot(-0.1,0.1);

  Mx =  CLHEP::RandFlat::shoot(0.7,1);
  My =  CLHEP::RandFlat::shoot(-0.3,0.3); // ranges good for 0<Theta<20
  Mz =  CLHEP::RandFlat::shoot(-0.3,0.3);

  condizione = std::sqrt(Mx*Mx + My*My + Mz*Mz);

 
  if (condizione < 1)  {
    Mx = Mx/condizione;
    My = My/condizione;
    Mz = Mz/condizione;


    if (Mx > std::cos(Theta)) { 
      break;
        }
    }
}
  
 
  particleGun -> SetParticleMomentumDirection( G4ThreeVector(Mx,My,Mz) );
  

  // Generate a primary particle
  particleGun -> GeneratePrimaryVertex( anEvent ); 
} 

Here is the call graph for this function:

G4double IORTPrimaryGeneratorAction::GetmeanKineticEnergy

(

void

)

Definition at line 273 of file IORTPrimaryGeneratorAction.cc.

{ return meanKineticEnergy;}

G4ParticleGun* IORTPrimaryGeneratorAction::GetParticleGun ( void  )

inline

Definition at line 71 of file IORTPrimaryGeneratorAction.hh.

{return particleGun;}

void IORTPrimaryGeneratorAction::SetmeanKineticEnergy

(

G4double

val

)

Definition at line 227 of file IORTPrimaryGeneratorAction.cc.

{
    meanKineticEnergy = val;

  // Update the beam-data in the analysis manager
  IORTAnalysisManager::GetInstance()->
    setBeamMetaData(meanKineticEnergy, sigmaEnergy);


} 

Here is the call graph for this function:

void IORTPrimaryGeneratorAction::SetsigmaEnergy

(

G4double

val

)

Definition at line 238 of file IORTPrimaryGeneratorAction.cc.

{ 
  sigmaEnergy = val;
  // Update the sigmaenergy in the metadata.
  IORTAnalysisManager::GetInstance()->
    setBeamMetaData(meanKineticEnergy, sigmaEnergy);
}

Here is the call graph for this function:

void IORTPrimaryGeneratorAction::SetsigmaY

(

G4double

val

)

Definition at line 255 of file IORTPrimaryGeneratorAction.cc.

{ sigmaY = val;}

void IORTPrimaryGeneratorAction::SetsigmaZ

(

G4double

val

)

Definition at line 258 of file IORTPrimaryGeneratorAction.cc.

{ sigmaZ = val;}

void IORTPrimaryGeneratorAction::SetTheta

(

G4double

val

)

Definition at line 269 of file IORTPrimaryGeneratorAction.cc.

{ Theta = val;}

void IORTPrimaryGeneratorAction::SetXposition

(

G4double

val

)

Definition at line 246 of file IORTPrimaryGeneratorAction.cc.

{ X0 = val;}

void IORTPrimaryGeneratorAction::SetYposition

(

G4double

val

)

Definition at line 249 of file IORTPrimaryGeneratorAction.cc.

{ Y0 = val;}

void IORTPrimaryGeneratorAction::SetZposition

(

G4double

val

)

Definition at line 252 of file IORTPrimaryGeneratorAction.cc.

{ Z0 = val;}

---

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

• source/geant4.10.03.p03/examples/advanced/iort_therapy/include/IORTPrimaryGeneratorAction.hh
• source/geant4.10.03.p03/examples/advanced/iort_therapy/src/IORTPrimaryGeneratorAction.cc