#include <XrayFluoPrimaryGeneratorAction.hh>

Inheritance diagram for XrayFluoPrimaryGeneratorAction:

Public Member Functions
	XrayFluoPrimaryGeneratorAction (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 ( XrayFluoDetectorConstruction * XrayFluoDC )

Definition at line 56 of file XrayFluoPrimaryGeneratorAction.cc.

References python.hepunit::cm, G4ParticleTable::FindParticle(), G4cout, G4endl, G4ParticleTable::GetParticleTable(), XrayFluoDetectorConstruction::GetWorldSizeZ(), python.hepunit::keV, G4ParticleGun::SetParticleDefinition(), G4ParticleGun::SetParticleEnergy(), G4ParticleGun::SetParticleMomentumDirection(), and G4VPrimaryGenerator::SetParticlePosition().

   :rndmFlag("off"),beam("off"),spectrum("off"),isoVert("off"),phaseSpaceGunFlag(false), 
    rayleighFlag(true), particleEnergies(0),  particleTypes(0), detectorPosition(0) 
 {
 
   XrayFluoDetector = XrayFluoDC;
 
   G4int n_particle = 1;
   particleGun  = new G4ParticleGun(n_particle);
   
   //create a messenger for this class
   gunMessenger = new XrayFluoPrimaryGeneratorMessenger(this);
   runManager = new XrayFluoRunAction();
   
   // default particle kinematic
   
   G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
   G4String particleName;
   G4ParticleDefinition* particle
     = particleTable->FindParticle(particleName="gamma");
   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;
   
 }

XrayFluoPrimaryGeneratorAction::~XrayFluoPrimaryGeneratorAction ( )

Definition at line 123 of file XrayFluoPrimaryGeneratorAction.cc.

References G4cout, and G4endl.

 {
   delete particleGun;
   delete gunMessenger;
   delete runManager;
 
   G4cout << "XrayFluoPrimaryGeneratorAction deleted" << G4endl;
 
 }

Member Function Documentation

void XrayFluoPrimaryGeneratorAction::ActivatePhaseSpace ( G4String )

Definition at line 106 of file XrayFluoPrimaryGeneratorAction.cc.

Referenced by XrayFluoPrimaryGeneratorMessenger::SetNewValue().

void XrayFluoPrimaryGeneratorAction::GeneratePrimaries ( G4Event * anEvent )

virtual

Implements G4VUserPrimaryGeneratorAction.

Definition at line 135 of file XrayFluoPrimaryGeneratorAction.cc.

References python.hepunit::cm, energy(), G4ParticleTable::FindParticle(), XrayFluoDataSet::FindValue(), G4UniformRand, G4ParticleGun::GeneratePrimaryVertex(), XrayFluoRunAction::GetData(), XrayFluoRunAction::GetDataSum(), XrayFluoDetectorConstruction::GetDia3SizeXY(), XrayFluoRunAction::GetEnergies(), G4Event::GetEventID(), XrayFluoRunAction::GetGammaSet(), G4ParticleGun::GetParticleDefinition(), G4ParticleGun::GetParticleEnergy(), G4ParticleDefinition::GetParticleName(), G4ParticleTable::GetParticleTable(), XrayFluoDetectorConstruction::GetSampleSizeXY(), XrayFluoDetectorConstruction::GetWorldSizeZ(), python.hepunit::keV, python.hepunit::m, python.hepunit::mm, python.hepunit::pi, python.hepunit::rad, G4ParticleGun::SetParticleDefinition(), G4ParticleGun::SetParticleEnergy(), G4ParticleGun::SetParticleMomentumDirection(), G4VPrimaryGenerator::SetParticlePosition(), test::x, z, and G4InuclParticleNames::z0.

 {
   //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 =  runManager->GetEnergies();
       G4DataVector* data =  runManager->GetData();
      
       G4double sum = runManager->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 = runManager->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){
 
     size_t i = anEvent->GetEventID();
 
     particleGun->SetParticlePosition(detectorPosition); 
     particleGun->SetParticleMomentumDirection(detectorPosition);
     G4double energy;
     G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
     G4ParticleDefinition* particle;
     if (i < particleEnergies->size()) {      
       energy = (*particleEnergies)[i];
       particle = particleTable->FindParticle((*particleTypes)[i]);
     }
 
     else {
       energy = 0.;
       particle = particleTable->FindParticle("gamma");
     }
 
     particleGun->SetParticleEnergy(energy);
     particleGun->SetParticleDefinition(particle);
 
 
   }
 
 #ifdef G4ANALYSIS_USE 
 
   G4double partEnergy = particleGun->GetParticleEnergy();
   XrayFluoAnalysisManager* analysis =  XrayFluoAnalysisManager::getInstance();
   analysis->analysePrimaryGenerator(partEnergy/keV);
 
 #endif
  
   particleGun->GeneratePrimaryVertex(anEvent);
 }

void XrayFluoPrimaryGeneratorAction::SetIsoVert ( G4String val )

inline

Definition at line 77 of file XrayFluoPrimaryGeneratorAction.hh.

Referenced by XrayFluoPrimaryGeneratorMessenger::SetNewValue().

77 { isoVert = val ;}

void XrayFluoPrimaryGeneratorAction::SetRayleighFlag ( G4bool value )

Definition at line 112 of file XrayFluoPrimaryGeneratorAction.cc.

References G4cout, and G4endl.

Referenced by XrayFluoPrimaryGeneratorMessenger::SetNewValue().

 {
 
   G4cout <<  "value: " << value << G4endl; 
   rayleighFlag = value;
   
 }

void XrayFluoPrimaryGeneratorAction::SetRndmFlag ( G4String val )

inline

Definition at line 68 of file XrayFluoPrimaryGeneratorAction.hh.

Referenced by XrayFluoPrimaryGeneratorMessenger::SetNewValue().

68 { rndmFlag = val;}

void XrayFluoPrimaryGeneratorAction::SetRndmVert ( G4String val )

inline

Definition at line 71 of file XrayFluoPrimaryGeneratorAction.hh.

Referenced by XrayFluoPrimaryGeneratorMessenger::SetNewValue().

71 { beam = val;}

void XrayFluoPrimaryGeneratorAction::SetSpectrum ( G4String val )

inline

Definition at line 74 of file XrayFluoPrimaryGeneratorAction.hh.

Referenced by XrayFluoPrimaryGeneratorMessenger::SetNewValue().

74 { spectrum= val ;}

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

Public Member Functions

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation