#include <UltraPrimaryGeneratorAction.hh>

Inheritance diagram for UltraPrimaryGeneratorAction:

Collaboration diagram for UltraPrimaryGeneratorAction:

Public Member Functions
	UltraPrimaryGeneratorAction ()

	~UltraPrimaryGeneratorAction ()

void	GeneratePrimaries (G4Event *anEvent)

G4GeneralParticleSource *	GetParticleGun ()

Public Member Functions inherited from G4VUserPrimaryGeneratorAction
	G4VUserPrimaryGeneratorAction ()

virtual	~G4VUserPrimaryGeneratorAction ()

Private Attributes
G4GeneralParticleSource *	particleGun

Detailed Description

Definition at line 53 of file UltraPrimaryGeneratorAction.hh.

Constructor & Destructor Documentation

◆ UltraPrimaryGeneratorAction()

UltraPrimaryGeneratorAction::UltraPrimaryGeneratorAction ( )

Definition at line 61 of file UltraPrimaryGeneratorAction.cc.

 {
 
   particleGun = new G4GeneralParticleSource();
 
 }

◆ ~UltraPrimaryGeneratorAction()

UltraPrimaryGeneratorAction::~UltraPrimaryGeneratorAction ( )

Definition at line 72 of file UltraPrimaryGeneratorAction.cc.

 {
   delete particleGun;
 }

Member Function Documentation

◆ GeneratePrimaries()

void UltraPrimaryGeneratorAction::GeneratePrimaries ( G4Event * anEvent )

virtual

Implements G4VUserPrimaryGeneratorAction.

Definition at line 79 of file UltraPrimaryGeneratorAction.cc.

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