XrayFluoPrimaryGeneratorAction Class Reference

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

Detailed Description

Definition at line 56 of file XrayFluoPrimaryGeneratorAction.hh.

Constructor & Destructor Documentation

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

(

)

Definition at line 111 of file XrayFluoPrimaryGeneratorAction.cc.

{
  delete particleGun;
  delete gunMessenger;
}

Member Function Documentation

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:

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

Here is the call graph for this function:

void XrayFluoPrimaryGeneratorAction::SetIsoVert ( G4String  val )

inline

Definition at line 77 of file XrayFluoPrimaryGeneratorAction.hh.

{ isoVert = val  ;}

Here is the caller graph for this function:

void XrayFluoPrimaryGeneratorAction::SetRayleighFlag

(

G4bool

value

)

Definition at line 103 of file XrayFluoPrimaryGeneratorAction.cc.

{
  rayleighFlag = value; 
}

Here is the caller graph for this function:

void XrayFluoPrimaryGeneratorAction::SetRndmFlag ( G4String  val )

inline

Definition at line 68 of file XrayFluoPrimaryGeneratorAction.hh.

{ rndmFlag = val;}

Here is the caller graph for this function:

void XrayFluoPrimaryGeneratorAction::SetRndmVert ( G4String  val )

inline

Definition at line 71 of file XrayFluoPrimaryGeneratorAction.hh.

{ beam = val;}

Here is the caller graph for this function:

void XrayFluoPrimaryGeneratorAction::SetSpectrum ( G4String  val )

inline

Definition at line 74 of file XrayFluoPrimaryGeneratorAction.hh.

{ spectrum= val  ;}

Here is the caller graph for this function:

---

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

• source/geant4.10.03.p03/examples/advanced/xray_fluorescence/include/XrayFluoPrimaryGeneratorAction.hh
• source/geant4.10.03.p03/examples/advanced/xray_fluorescence/src/XrayFluoPrimaryGeneratorAction.cc