XrayFluoMercuryPrimaryGeneratorAction Class Reference

#include <XrayFluoMercuryPrimaryGeneratorAction.hh>

Inheritance diagram for XrayFluoMercuryPrimaryGeneratorAction:

Public Member Functions
	XrayFluoMercuryPrimaryGeneratorAction (XrayFluoMercuryDetectorConstruction *)
	~XrayFluoMercuryPrimaryGeneratorAction ()
void	GeneratePrimaries (G4Event *)
void	SetGlobalFlag (G4bool val)
void	SetSpectrum (G4String val)
Public Member Functions inherited from G4VUserPrimaryGeneratorAction
	G4VUserPrimaryGeneratorAction ()
virtual	~G4VUserPrimaryGeneratorAction ()

Detailed Description

Definition at line 53 of file XrayFluoMercuryPrimaryGeneratorAction.hh.

Constructor & Destructor Documentation

XrayFluoMercuryPrimaryGeneratorAction::XrayFluoMercuryPrimaryGeneratorAction

(

XrayFluoMercuryDetectorConstruction *

XrayFluoDC

)

Definition at line 55 of file XrayFluoMercuryPrimaryGeneratorAction.cc.

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

  :globalFlag(false),spectrum("off")
{

  XrayFluoDetector = XrayFluoDC;

  G4int n_particle = 1;
  particleGun  = new G4ParticleGun(n_particle);
  
  //create a messenger for this class
  gunMessenger = new XrayFluoMercuryPrimaryGeneratorMessenger(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 << "XrayFluoMercuryPrimaryGeneratorAction created" << G4endl;
  
}

XrayFluoMercuryPrimaryGeneratorAction::~XrayFluoMercuryPrimaryGeneratorAction

(

)

Definition at line 89 of file XrayFluoMercuryPrimaryGeneratorAction.cc.

References G4cout, and G4endl.

{
  delete particleGun;
  delete gunMessenger;
  delete runManager;

  G4cout << "XrayFluoMercuryPrimaryGeneratorAction deleted" << G4endl;

}

Member Function Documentation

void XrayFluoMercuryPrimaryGeneratorAction::GeneratePrimaries ( G4Event *  anEvent )

virtual

Implements G4VUserPrimaryGeneratorAction.

Definition at line 101 of file XrayFluoMercuryPrimaryGeneratorAction.cc.

References test::a, python.hepunit::deg, XrayFluoDataSet::FindValue(), G4UniformRand, G4ParticleGun::GeneratePrimaryVertex(), XrayFluoRunAction::GetData(), XrayFluoRunAction::GetDataSum(), XrayFluoRunAction::GetEnergies(), XrayFluoRunAction::GetGammaSet(), XrayFluoMercuryDetectorConstruction::GetMercuryDia(), XrayFluoMercuryDetectorConstruction::GetOpticAperture(), XrayFluoMercuryDetectorConstruction::GetOrbitInclination(), G4ParticleGun::GetParticleDefinition(), G4ParticleGun::GetParticleEnergy(), G4ParticleDefinition::GetParticleName(), XrayFluoMercuryDetectorConstruction::GetSunDia(), XrayFluoMercuryDetectorConstruction::GetWorldSizeZ(), python.hepunit::keV, python.hepunit::m, python.hepunit::pi, python.hepunit::rad, G4ParticleGun::SetParticleEnergy(), G4ParticleGun::SetParticleMomentumDirection(), G4VPrimaryGenerator::SetParticlePosition(), and G4InuclParticleNames::z0.

{
  //this function is called at the begining of event
  // 

  // Conidering the sunas a Poin-like source.

  G4double z0 = -0.5*(XrayFluoDetector->GetWorldSizeZ());
  G4double y0 = 0.*m, x0 = 0.*m;


  // Let's try to illuminate only the prtion of Mercury surface that can be seen by the detector.

  G4double spacecraftLatitude = XrayFluoDetector->GetOrbitInclination();
  G4double mercuryDia = XrayFluoDetector->GetMercuryDia();
  G4double sunDia = XrayFluoDetector->GetSunDia();
  G4double opticField = XrayFluoDetector->GetOpticAperture();
  
 
  G4double a = 2*std::tan(opticField/2);
  
  //  if (!pointLikeFlag) {

    // let's decide from wich point of the sun surface the particle is coming:

  G4double theta = std::acos(2.*G4UniformRand() - 1.0);
    G4double phi = 2. * pi * G4UniformRand();
    G4double rho = sunDia/2;                         

    G4double sunPosX = x0 + rho * std::sin(theta) * std::cos(phi);
    G4double sunPosY = y0 + rho * std::sin(theta) * std::sin(phi);
    G4double sunPosZ = z0 + rho * std::cos(theta);           

    particleGun->SetParticlePosition(G4ThreeVector(sunPosX,sunPosY,sunPosZ));

    // the angle at the center of Mercury subtending the area seen by the optics:
    G4double alpha = 2 * a/mercuryDia; 
    
    if(!globalFlag){
      theta = alpha * G4UniformRand() + (180.*deg - spacecraftLatitude)-alpha/2.;
      phi = alpha * G4UniformRand() + 90. * deg - alpha/2.;     
    }    
    
    else if(globalFlag){
      theta = pi/2. * rad * G4UniformRand() + 90.*deg ; //was 900., probably an error
      phi = 2*pi*rad * G4UniformRand() ;
    }
    
    rho = mercuryDia/2.;                        
    
    G4double mercuryPosX = rho * std::sin(theta) * std::cos(phi);
    G4double mercuryPosY = rho * std::sin(theta) * std::sin(phi);
    G4double mercuryPosZ = rho * std::cos(theta);           
    
    particleGun->SetParticleMomentumDirection(
                G4ThreeVector(mercuryPosX-sunPosX ,mercuryPosY-sunPosY,mercuryPosZ-sunPosZ));

    //  }
//   if (pointLikeFlag) {

//   // theta is the angle that the mean direction of the incident light (on the desired 
//   // point of the surface of Mercury) makes with  the Z-axis
//   G4double theta = std::asin( mercuryDia/2. * std::sin(spacecraftLatitude) / 
//           std::sqrt(std::pow(z0,2)+std::pow(mercuryDia/2.,2)-2*mercuryDia/2.*z0*std::cos(spacecraftLatitude)) );  

//   // on the y axis, the light emitted from the Sun must be in [theta-phi;theta+phi]
//   G4double phi = std::asin( mercuryDia/2.*std::sin(spacecraftLatitude) + a*std::cos(spacecraftLatitude) /
//             std::sqrt( std::pow(mercuryDia/2.*std::sin(spacecraftLatitude) + a*std::cos(spacecraftLatitude) , 2) +
//               std::pow(z0 - mercuryDia/2.*std::cos(spacecraftLatitude) - a*std::sin(spacecraftLatitude) , 2)) ) 
//     - theta;  
  
//   // on the x axis, the light emitted from the Sun must be in [-zeta;zeta]
//   G4double zeta = std::atan( a/std::sqrt(std::pow(z0,2)+std::pow(mercuryDia,2)-2*mercuryDia*z0*std::cos(spacecraftLatitude)) );
  
  
  
//   //alpha in [-zeta;zeta]
//   G4double alpha = (2*zeta)*G4UniformRand() - zeta;
//   //beta in [theta-phi;theta+phi]
//   G4double beta = (G4UniformRand()*2*phi) - phi + theta;
  
//   G4double dirY = std::sin(beta);
//   G4double dirX = std::sin(alpha);
  
//   particleGun->SetParticleMomentumDirection(G4ThreeVector(dirX.,dirY,1.));
  
//   particleGun->SetParticlePosition(G4ThreeVector(x0,y0,z0));

//   }


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

#ifdef G4ANALYSIS_USE 

  G4double partEnergy = particleGun->GetParticleEnergy();
  XrayFluoAnalysisManager* analysis =  XrayFluoAnalysisManager::getInstance();
  analysis->analysePrimaryGenerator(partEnergy/keV);

#endif

  particleGun->GeneratePrimaryVertex(anEvent);
}

void XrayFluoMercuryPrimaryGeneratorAction::SetGlobalFlag ( G4bool  val )

inline

Definition at line 65 of file XrayFluoMercuryPrimaryGeneratorAction.hh.

Referenced by XrayFluoMercuryPrimaryGeneratorMessenger::SetNewValue().

{globalFlag = val;};

void XrayFluoMercuryPrimaryGeneratorAction::SetSpectrum ( G4String  val )

inline

Definition at line 75 of file XrayFluoMercuryPrimaryGeneratorAction.hh.

Referenced by XrayFluoMercuryPrimaryGeneratorMessenger::SetNewValue().

{ spectrum = val;};

---

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

• XrayFluoMercuryPrimaryGeneratorAction.hh
• XrayFluoMercuryPrimaryGeneratorAction.cc