#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 ()

Private Member Functions
void	SetDefaultPrimaryParticle ()

Private Attributes
G4double	meanKineticEnergy

G4double	sigmaEnergy

G4double	X0

G4double	Y0

G4double	Z0

G4double	sigmaY

G4double	sigmaZ

G4double	Theta

G4ParticleGun *	particleGun

IORTPrimaryGeneratorMessenger *	gunMessenger

Detailed Description

Definition at line 50 of file IORTPrimaryGeneratorAction.hh.

Constructor & Destructor Documentation

◆ IORTPrimaryGeneratorAction()

IORTPrimaryGeneratorAction::IORTPrimaryGeneratorAction ( )

Definition at line 55 of file IORTPrimaryGeneratorAction.cc.

 {
   // Define the messenger
   gunMessenger = new IORTPrimaryGeneratorMessenger(this);
 
   particleGun  = new G4ParticleGun();
 
   SetDefaultPrimaryParticle();  
 }  

Here is the call graph for this function:

◆ ~IORTPrimaryGeneratorAction()

IORTPrimaryGeneratorAction::~IORTPrimaryGeneratorAction ( )

Definition at line 65 of file IORTPrimaryGeneratorAction.cc.

 {
   delete particleGun;
 
   delete gunMessenger;
 }

Member Function Documentation

◆ GeneratePrimaries()

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:

◆ GetmeanKineticEnergy()

G4double IORTPrimaryGeneratorAction::GetmeanKineticEnergy ( void )

Definition at line 273 of file IORTPrimaryGeneratorAction.cc.

274 { return meanKineticEnergy;}

IORTPrimaryGeneratorAction::meanKineticEnergy

G4double meanKineticEnergy

Definition: IORTPrimaryGeneratorAction.hh:75

◆ GetParticleGun()

G4ParticleGun* IORTPrimaryGeneratorAction::GetParticleGun ( void )

inline

Definition at line 71 of file IORTPrimaryGeneratorAction.hh.

71 {return particleGun;}

IORTPrimaryGeneratorAction::particleGun

G4ParticleGun * particleGun

Definition: IORTPrimaryGeneratorAction.hh:87

Here is the call graph for this function:

◆ SetDefaultPrimaryParticle()

void IORTPrimaryGeneratorAction::SetDefaultPrimaryParticle ( )

private

Definition at line 72 of file IORTPrimaryGeneratorAction.cc.

 {    
   // ****************************
   // Default primary particle
   // ****************************
   
   // Define primary particles: electrons // protons
   G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
   G4ParticleDefinition* particle = particleTable -> FindParticle("e-");   // ("proton")
   particleGun -> SetParticleDefinition(particle); 
 
   // Define the energy of primary particles:
   // gaussian distribution with mean energy = 10.0 *MeV   
   // and sigma = 400.0 *keV
   G4double defaultMeanKineticEnergy = 10.0 *CLHEP::MeV;   
   meanKineticEnergy = defaultMeanKineticEnergy;
 
   G4double defaultsigmaEnergy = 100.0 *CLHEP::keV;   
   sigmaEnergy = defaultsigmaEnergy;
   
   // Write these values into the analysis if needed. Have to be written separately on change.
   IORTAnalysisManager::GetInstance()->
     setBeamMetaData(meanKineticEnergy, sigmaEnergy);
 
 
   // Define the parameters of the initial position: 
   // the y, z coordinates have a gaussian distribution
   
   
   G4double defaultX0 = -862.817 *CLHEP::mm;                 
   X0 = defaultX0;
 
   G4double defaultY0 = 0.0 *CLHEP::mm;  
   Y0 = defaultY0;
 
   G4double defaultZ0 = 0.0 *CLHEP::mm;  
   Z0 = defaultZ0;
 
   G4double defaultsigmaY = 1. *CLHEP::mm;  
   sigmaY = defaultsigmaY;
 
   G4double defaultsigmaZ = 1. *CLHEP::mm;  
   sigmaZ = defaultsigmaZ;
 
   // Define the parameters of the momentum of primary particles: 
   // The momentum along the y and z axis has a gaussian distribution
  
 /* 
   G4double defaultsigmaMomentumY = 0.0;  
   sigmaMomentumY = defaultsigmaMomentumY;
 
   G4double defaultsigmaMomentumZ = 0.0;  
   sigmaMomentumZ = defaultsigmaMomentumZ;
 */
 
   G4double defaultTheta = 6.0 *CLHEP::deg;  
   Theta = defaultTheta;
 
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ SetmeanKineticEnergy()

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:

◆ SetsigmaEnergy()

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