UltraPrimaryGeneratorAction Class Reference

#include <UltraPrimaryGeneratorAction.hh>

Inheritance diagram for UltraPrimaryGeneratorAction:

Public Member Functions
	UltraPrimaryGeneratorAction ()
	~UltraPrimaryGeneratorAction ()
void	GeneratePrimaries (G4Event *anEvent)
G4GeneralParticleSource *	GetParticleGun ()
Public Member Functions inherited from G4VUserPrimaryGeneratorAction
	G4VUserPrimaryGeneratorAction ()
virtual	~G4VUserPrimaryGeneratorAction ()

Detailed Description

Definition at line 53 of file UltraPrimaryGeneratorAction.hh.

Constructor & Destructor Documentation

UltraPrimaryGeneratorAction::UltraPrimaryGeneratorAction

(

)

Definition at line 61 of file UltraPrimaryGeneratorAction.cc.

References python.hepunit::cm, python.hepunit::eV, G4ParticleTable::FindParticle(), G4SingleParticleSource::GetAngDist(), G4GeneralParticleSource::GetCurrentSource(), G4SingleParticleSource::GetEneDist(), G4ParticleTable::GetParticleTable(), G4SingleParticleSource::GetPosDist(), G4SPSPosDistribution::SetCentreCoords(), G4SPSEneDistribution::SetEnergyDisType(), G4SPSEneDistribution::SetMonoEnergy(), G4GeneralParticleSource::SetParticleDefinition(), G4SPSAngDistribution::SetParticleMomentumDirection(), G4GeneralParticleSource::SetParticlePolarization(), G4SPSPosDistribution::SetPosDisShape(), G4SPSPosDistribution::SetPosDisType(), G4SPSPosDistribution::SetPosRot1(), G4SPSPosDistribution::SetPosRot2(), and G4SPSPosDistribution::SetRadius().

{

  particleGun = new G4GeneralParticleSource();


  // Define here the user default properties for the General Particle Source (GPS)
  // Can be modified through the GPS Messenger (/gps/... commands)

  G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable();
  G4String particleName;

  G4ParticleDefinition* opticalphoton = particleTable->FindParticle(particleName="opticalphoton");

  //....PARTICLE DEFINITIONS
  particleGun->SetParticleDefinition(opticalphoton);

  G4ThreeVector Polarization = G4ThreeVector(1.,1.,0.) ;
  particleGun->SetParticlePolarization(Polarization);
  
  // DEFINE A MONO-ENERGETIC SOURCE
  G4SPSEneDistribution *eneDist = particleGun->GetCurrentSource()->GetEneDist() ;
  eneDist->SetEnergyDisType("Mono");
  eneDist->SetMonoEnergy(3.0*eV);

  // SET POSITION DISTRIBUTION 
  G4SPSPosDistribution *posDist = particleGun->GetCurrentSource()->GetPosDist() ;
  posDist->SetPosDisType("Plane");
  posDist->SetPosDisShape("Circle");
  posDist->SetRadius(20.0*cm);

#ifdef ULTRA_MIRROR_USE
#define ULTRA_REFLECTION_USE
#endif

#ifdef ULTRA_GROUND_USE
#define ULTRA_REFLECTION_USE
#endif

  G4SPSAngDistribution *angDist = particleGun->GetCurrentSource()->GetAngDist() ;

#ifdef ULTRA_REFLECTION_USE
  angDist->SetParticleMomentumDirection(G4ThreeVector(0.0,-1.0,0.0)) ;
  posDist->SetPosRot1(G4ThreeVector(1.,0.,0.));
  posDist->SetPosRot2(G4ThreeVector(0.,0.,-1.));
  posDist->SetCentreCoords(G4ThreeVector(0.0*cm,90.0*cm,150.0*cm));

#else
  angDist->SetParticleMomentumDirection(G4ThreeVector(0.0,0.0,-1.0)) ;
  posDist->SetCentreCoords(G4ThreeVector(0.0*cm,0.0*cm,150.0*cm));

#endif

}

UltraPrimaryGeneratorAction::~UltraPrimaryGeneratorAction

(

)

Definition at line 118 of file UltraPrimaryGeneratorAction.cc.

{
  delete particleGun;
}

Member Function Documentation

void UltraPrimaryGeneratorAction::GeneratePrimaries ( G4Event *  anEvent )

virtual

Implements G4VUserPrimaryGeneratorAction.

Definition at line 125 of file UltraPrimaryGeneratorAction.cc.

References python.hepunit::c_light, energy(), python.hepunit::eV, FatalException, G4cerr, G4cout, G4endl, G4Exception(), G4GeneralParticleSource::GeneratePrimaryVertex(), G4GeneralParticleSource::GetCurrentSource(), G4SingleParticleSource::GetEneDist(), G4SPSEneDistribution::GetEnergyDisType(), G4Event::GetEventID(), UltraDetectorConstruction::GetLambdaMax(), UltraDetectorConstruction::GetLambdaMin(), G4GeneralParticleSource::GetParticleDefinition(), G4GeneralParticleSource::GetParticleEnergy(), G4ParticleDefinition::GetParticleName(), G4SingleParticleSource::GetPosDist(), G4SPSPosDistribution::GetPosDisType(), G4RunManager::GetRunManager(), and python.hepunit::h_Planck.

{

  G4int iEvent = anEvent->GetEventID() ;
  if ( iEvent == 0 ){

    G4cout << particleGun->GetParticleDefinition()->GetParticleName()           << " " ;
    G4cout << particleGun->GetCurrentSource()->GetEneDist()->GetEnergyDisType() << " " ;
    G4cout << particleGun->GetCurrentSource()->GetPosDist()->GetPosDisType()    << G4endl ;


// Check if optical photon wavelength is within limits set for material optical properties tables. 

   


}
  particleGun->GeneratePrimaryVertex(anEvent);

    if (particleGun->GetParticleDefinition()->GetParticleName() == "opticalphoton"){
     
        const UltraDetectorConstruction * detector =  
        dynamic_cast<const UltraDetectorConstruction *>((G4RunManager::GetRunManager())->GetUserDetectorConstruction()) ;

    G4double lambda_min = detector->GetLambdaMin() ;
    G4double lambda_max = detector->GetLambdaMax() ;

        G4double energy = particleGun->GetParticleEnergy() ;

    if (h_Planck*c_light/energy > lambda_max || h_Planck*c_light/energy < lambda_min){
           G4cerr << "Error ! Optical photon energy (" << energy/eV << " eV) out of limits set by material optical properties tables. \n" 
              << "Please check that photon wavelength is within the following interval: [" 
              << lambda_min/nm << "," 
              << lambda_max/nm << "] nm" 
              << ", i.e., ["
              << h_Planck*c_light/lambda_max/eV << ","
              << h_Planck*c_light/lambda_min/eV << "] eV"
              << G4endl ;
    
           G4Exception("UltraPrimaryGeneratorAction::GeneratePrimaries()","AirSh005",
               FatalException,"Wavelength outside the valid range") ;
    }
 }

}

G4GeneralParticleSource* UltraPrimaryGeneratorAction::GetParticleGun ( )

inline

Definition at line 64 of file UltraPrimaryGeneratorAction.hh.

{return particleGun;}    

---

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

• UltraPrimaryGeneratorAction.hh
• UltraPrimaryGeneratorAction.cc