#include <XrayFluoPrimaryGeneratorAction.hh>

Inheritance diagram for XrayFluoPrimaryGeneratorAction:

Collaboration diagram for XrayFluoPrimaryGeneratorAction:

Public Member Functions
	XrayFluoPrimaryGeneratorAction (const XrayFluoDetectorConstruction *)

	~XrayFluoPrimaryGeneratorAction ()

void	GeneratePrimaries (G4Event *)

void	SetRndmFlag (G4String val)

void	SetRndmVert (G4String val)

void	SetSpectrum (G4String val)

void	SetIsoVert (G4String val)

void	ActivatePhaseSpace (G4String)

void	SetRayleighFlag (G4bool)

Public Member Functions inherited from G4VUserPrimaryGeneratorAction
	G4VUserPrimaryGeneratorAction ()

virtual	~G4VUserPrimaryGeneratorAction ()

Private Attributes
G4ParticleGun *	particleGun

const XrayFluoDetectorConstruction *	XrayFluoDetector

const XrayFluoRunAction *	runAction

XrayFluoPrimaryGeneratorMessenger *	gunMessenger

G4String	rndmFlag

G4String	beam

G4String	spectrum

G4String	isoVert

G4bool	phaseSpaceGunFlag

G4bool	rayleighFlag

G4ThreeVector	detectorPosition

Detailed Description

Definition at line 56 of file XrayFluoPrimaryGeneratorAction.hh.

Constructor & Destructor Documentation

◆ XrayFluoPrimaryGeneratorAction()

XrayFluoPrimaryGeneratorAction::XrayFluoPrimaryGeneratorAction ( const XrayFluoDetectorConstruction * XrayFluoDC )

Definition at line 56 of file XrayFluoPrimaryGeneratorAction.cc.

   :rndmFlag("off"),beam("off"),spectrum("off"),isoVert("off"),phaseSpaceGunFlag(false), 
    rayleighFlag(true), detectorPosition(0) 
 {
   runAction = 0;
   XrayFluoDetector = XrayFluoDC;
 
   G4int n_particle = 1;
   particleGun  = new G4ParticleGun(n_particle);
   
   //create a messenger for this class
   gunMessenger = new XrayFluoPrimaryGeneratorMessenger(this); 
 
   // default particle kinematic
   G4ParticleDefinition* particle
     = G4Gamma::Definition();
   particleGun->SetParticleDefinition(particle);
   particleGun->SetParticleMomentumDirection(G4ThreeVector(0.,0.,1.));
   particleGun->SetParticleEnergy(10. * keV);
 
   G4double position = -0.5*(XrayFluoDetector->GetWorldSizeZ());
   particleGun->SetParticlePosition(G4ThreeVector(0.*cm,0.*cm,position));
 
   G4cout << "XrayFluoPrimaryGeneratorAction created" << G4endl;
   
 }

Here is the call graph for this function:

◆ ~XrayFluoPrimaryGeneratorAction()

XrayFluoPrimaryGeneratorAction::~XrayFluoPrimaryGeneratorAction ( )

Definition at line 111 of file XrayFluoPrimaryGeneratorAction.cc.

 {
   delete particleGun;
   delete gunMessenger;
 }

Member Function Documentation

◆ ActivatePhaseSpace()

void XrayFluoPrimaryGeneratorAction::ActivatePhaseSpace ( G4String fileName )

Definition at line 86 of file XrayFluoPrimaryGeneratorAction.cc.

                                                                          {
 
   // load phase-space
   phaseSpaceGunFlag = true;
 
   // reads the data stored on disk form previous runs 
   // and get these data to data members
 
   XrayFluoAnalysisManager* analysis =  XrayFluoAnalysisManager::getInstance();
   analysis->LoadGunData(fileName, rayleighFlag);
   detectorPosition = XrayFluoDetector->GetDetectorPosition();
   detectorPosition.setR(detectorPosition.r()-(5.*cm)); // 5 cm before the detector, so in front of it.
 
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ GeneratePrimaries()

void XrayFluoPrimaryGeneratorAction::GeneratePrimaries ( G4Event * anEvent )

virtual

Implements G4VUserPrimaryGeneratorAction.

Definition at line 119 of file XrayFluoPrimaryGeneratorAction.cc.

 {
   //retrieve runAction, if not done
   if (!runAction)
     {
       //Sequential runaction
       if (G4RunManager::GetRunManager()->GetRunManagerType() == 
       G4RunManager::sequentialRM)
     runAction = static_cast<const XrayFluoRunAction*>
       (G4RunManager::GetRunManager()->GetUserRunAction());  
       else //MT master runaction
     runAction = static_cast<const XrayFluoRunAction*>
       (G4MTRunManager::GetMasterRunManager()->GetUserRunAction());  
       if (!runAction)
     G4cout << "Something wrong here!" << G4endl;
     }
  
   //this function is called at the begining of event
   // 
   G4double z0 = -0.5*(XrayFluoDetector->GetWorldSizeZ());
   G4double y0 = 0.*cm, x0 = 0.*cm;
   if (rndmFlag == "on")
     {
       y0 = (XrayFluoDetector->GetDia3SizeXY())/std::sqrt(2.)*(G4UniformRand()-0.5); // it was GetSampleSizeXY(), 
       x0 = (XrayFluoDetector->GetDia3SizeXY())/std::sqrt(2.)*(G4UniformRand()-0.5); // not divided by std::sqrt(2.)
     } 
   particleGun->SetParticlePosition(G4ThreeVector(x0,y0,z0));
   
   //randomize starting point
   if (beam == "on")
     {
       G4double radius = 0.5 * mm;
       G4double rho = radius*std::sqrt(G4UniformRand());
       G4double theta = 2*pi*G4UniformRand()*rad;
       G4double position = -0.5*(XrayFluoDetector->GetWorldSizeZ());
       
       G4double y = rho * std::sin(theta);
       G4double x = rho * std::cos(theta);
       
       particleGun->SetParticlePosition(G4ThreeVector(x,y,position));
     }
   //shoot particles according to a certain spectrum
   if (spectrum =="on")
     {
       G4String particle =  particleGun->GetParticleDefinition()
     ->GetParticleName();
       if(particle == "proton"|| particle == "alpha")
     {
       G4DataVector* energies =  runAction->GetEnergies();
       G4DataVector* data =  runAction->GetData();
      
       G4double sum = runAction->GetDataSum();
       G4double partSum = 0;
       G4int j = 0;
       G4double random= sum*G4UniformRand();
       while (partSum<random)
         {
           partSum += (*data)[j];
           j++;
         }
      
       particleGun->SetParticleEnergy((*energies)[j]);
     
     }
       else if (particle == "gamma")
     {
       const XrayFluoDataSet* dataSet = runAction->GetGammaSet();
       
       G4int i = 0;
       G4int id = 0;
       G4double minEnergy = 0. * keV;
       G4double particleEnergy= 0.;
       G4double maxEnergy = 10. * keV;
       G4double energyRange = maxEnergy - minEnergy;
 
        while ( i == 0)
         {
           G4double random = G4UniformRand();
           
           G4double randomNum = G4UniformRand(); //*5.0E6;
           
           particleEnergy = (random*energyRange) + minEnergy;
           
           if ((dataSet->FindValue(particleEnergy,id)) > randomNum)
         {
           i = 1;
           
         }
         }
        particleGun->SetParticleEnergy(particleEnergy);
     }
     }
 
   // Randomize starting point and direction
   
   if (isoVert == "on")
     {
       G4double rho = 1. *m;
       //theta in [0;pi/2]
       G4double theta = (pi/2)*G4UniformRand();
       //phi in [-pi;pi]
       G4double phi = (G4UniformRand()*2*pi)- pi;
       G4double x = rho*std::sin(theta)*std::sin(phi);
       G4double y = rho*std::sin(theta)*std::cos(phi);
       G4double z = -(rho*std::cos(theta));
       particleGun->SetParticlePosition(G4ThreeVector(x,y,z));
       
       G4double Xdim = XrayFluoDetector->GetSampleSizeXY();
       G4double Ydim = XrayFluoDetector->GetSampleSizeXY();
       
       G4double Dx = Xdim*(G4UniformRand()-0.5);
       
       G4double Dy = Ydim*(G4UniformRand()-0.5);
       
       particleGun->SetParticleMomentumDirection(G4ThreeVector(-x+Dx,-y+Dy,-z));
       
     }
 
   // using prevoiously genereated emissions from sample.....
 
   if (phaseSpaceGunFlag){
 
     particleGun->SetParticlePosition(detectorPosition); 
     particleGun->SetParticleMomentumDirection(detectorPosition);
    
     G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
    
     const std::pair<G4double,G4String> kine = 
       XrayFluoAnalysisManager::getInstance()->GetEmittedParticleEnergyAndType();
 
     G4double energy = kine.first;
     G4ParticleDefinition* particle = particleTable->FindParticle(kine.second);
 
     particleGun->SetParticleEnergy(energy);
     particleGun->SetParticleDefinition(particle);
 
 
   }
 
   G4double partEnergy = particleGun->GetParticleEnergy();
   XrayFluoAnalysisManager* analysis =  XrayFluoAnalysisManager::getInstance();
   analysis->analysePrimaryGenerator(partEnergy/keV);
 
  
   particleGun->GeneratePrimaryVertex(anEvent);
 }