Geant4  10.00.p01
UltraPrimaryGeneratorAction.cc
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 //
27 // --------------------------------------------------------------
28 // GEANT 4 - ULTRA experiment example
29 // --------------------------------------------------------------
30 //
31 // Code developed by:
32 // B. Tome, M.C. Espirito-Santo, A. Trindade, P. Rodrigues
33 //
34 // ****************************************************
35 // * UltraPrimaryGeneratorAction.cc
36 // ****************************************************
37 //
38 // Class used in the definition of the optical photons source
39 // A plane, circular source is used. Depending on the source position, optical
40 // photons may reach the UVscope directly or after reflection. By default direct
41 // incidence is used. The source parameters can be set directly in this class
42 // or through the GeneralParticleSource messenger class.
43 //
46 
47 #include "G4PhysicalConstants.hh"
48 #include "G4SystemOfUnits.hh"
49 #include "G4RunManager.hh"
50 #include "G4Event.hh"
52 #include "G4SPSAngDistribution.hh"
53 #include "G4SPSEneDistribution.hh"
54 #include "G4SPSPosDistribution.hh"
55 #include "G4ParticleTable.hh"
56 #include "G4ParticleDefinition.hh"
57 #include "G4ThreeVector.hh"
58 
59 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
60 
62 {
63 
65 
66 
67  // Define here the user default properties for the General Particle Source (GPS)
68  // Can be modified through the GPS Messenger (/gps/... commands)
69 
71  G4String particleName;
72 
73  G4ParticleDefinition* opticalphoton = particleTable->FindParticle(particleName="opticalphoton");
74 
75  //....PARTICLE DEFINITIONS
76  particleGun->SetParticleDefinition(opticalphoton);
77 
78  G4ThreeVector Polarization = G4ThreeVector(1.,1.,0.) ;
79  particleGun->SetParticlePolarization(Polarization);
80 
81  // DEFINE A MONO-ENERGETIC SOURCE
83  eneDist->SetEnergyDisType("Mono");
84  eneDist->SetMonoEnergy(3.0*eV);
85 
86  // SET POSITION DISTRIBUTION
88  posDist->SetPosDisType("Plane");
89  posDist->SetPosDisShape("Circle");
90  posDist->SetRadius(20.0*cm);
91 
92 #ifdef ULTRA_MIRROR_USE
93 #define ULTRA_REFLECTION_USE
94 #endif
95 
96 #ifdef ULTRA_GROUND_USE
97 #define ULTRA_REFLECTION_USE
98 #endif
99 
101 
102 #ifdef ULTRA_REFLECTION_USE
103  angDist->SetParticleMomentumDirection(G4ThreeVector(0.0,-1.0,0.0)) ;
104  posDist->SetPosRot1(G4ThreeVector(1.,0.,0.));
105  posDist->SetPosRot2(G4ThreeVector(0.,0.,-1.));
106  posDist->SetCentreCoords(G4ThreeVector(0.0*cm,90.0*cm,150.0*cm));
107 
108 #else
109  angDist->SetParticleMomentumDirection(G4ThreeVector(0.0,0.0,-1.0)) ;
110  posDist->SetCentreCoords(G4ThreeVector(0.0*cm,0.0*cm,150.0*cm));
111 
112 #endif
113 
114 }
115 
116 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
117 
119 {
120  delete particleGun;
121 }
122 
123 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
124 
126 {
127 
128  G4int iEvent = anEvent->GetEventID() ;
129  if ( iEvent == 0 ){
130 
134 
135 
136 // Check if optical photon wavelength is within limits set for material optical properties tables.
137 
138 
139 
140 
141 }
143 
144  if (particleGun->GetParticleDefinition()->GetParticleName() == "opticalphoton"){
145 
146  const UltraDetectorConstruction * detector =
147  dynamic_cast<const UltraDetectorConstruction *>((G4RunManager::GetRunManager())->GetUserDetectorConstruction()) ;
148 
149  G4double lambda_min = detector->GetLambdaMin() ;
150  G4double lambda_max = detector->GetLambdaMax() ;
151 
153 
154  if (h_Planck*c_light/energy > lambda_max || h_Planck*c_light/energy < lambda_min){
155  G4cerr << "Error ! Optical photon energy (" << energy/eV << " eV) out of limits set by material optical properties tables. \n"
156  << "Please check that photon wavelength is within the following interval: ["
157  << lambda_min/nm << ","
158  << lambda_max/nm << "] nm"
159  << ", i.e., ["
160  << h_Planck*c_light/lambda_max/eV << ","
161  << h_Planck*c_light/lambda_min/eV << "] eV"
162  << G4endl ;
163 
164  G4Exception("UltraPrimaryGeneratorAction::GeneratePrimaries()","AirSh005",
165  FatalException,"Wavelength outside the valid range") ;
166  }
167  }
168 
169 }
170 
171 
static const double cm
Definition: G4SIunits.hh:106
G4ParticleDefinition * GetParticleDefinition()
G4SPSEneDistribution * GetEneDist()
G4ParticleDefinition * FindParticle(G4int PDGEncoding)
CLHEP::Hep3Vector G4ThreeVector
G4SPSPosDistribution * GetPosDist()
G4GeneralParticleSource * particleGun
int G4int
Definition: G4Types.hh:78
const G4String & GetParticleName() const
G4int GetEventID() const
Definition: G4Event.hh:140
void SetParticleDefinition(G4ParticleDefinition *aParticleDefinition)
G4GLOB_DLL std::ostream G4cout
void SetCentreCoords(G4ThreeVector)
void SetEnergyDisType(G4String)
void SetParticlePolarization(G4ThreeVector aVal)
void SetPosRot1(G4ThreeVector)
G4SingleParticleSource * GetCurrentSource()
G4SPSAngDistribution * GetAngDist()
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
Definition: G4Exception.cc:41
static G4RunManager * GetRunManager()
Definition: G4RunManager.cc:74
static const double eV
Definition: G4SIunits.hh:194
static G4ParticleTable * GetParticleTable()
G4double energy(const ThreeVector &p, const G4double m)
void SetPosRot2(G4ThreeVector)
void SetParticleMomentumDirection(G4ParticleMomentum aMomentumDirection)
#define G4endl
Definition: G4ios.hh:61
double G4double
Definition: G4Types.hh:76
G4GLOB_DLL std::ostream G4cerr