| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 // -------------------------------------------
27 // =============== Begin Documentation Commen
28 //!
29 //! \file FFPrimaryGeneratorAction.cc
30 //! \author B. Wendt (brycen.linn.wendt@ce
31 //! \date June 06, 2014
32 //!
33 //! \brief Implementation of the FFPrimar
34 //!
35 //! \details Generates 4.5 MeV neutrons fro
36 //! FFDetectorConstruction as "Net
37 //! in a isotropically sampled dir
38 //!
39 // ================ End Documentation Comment
40 //
41 // Modified:
42 //
43 // 23-06-14
44 // Implemented "GetNeutronSourceCenter()" and
45 // the neutron starting location
46 //
47 // -------------------------------------------
48
49 #include "FFPrimaryGeneratorAction.hh"
50
51 #include "G4Event.hh"
52 #include "G4LogicalVolume.hh"
53 #include "G4LogicalVolumeStore.hh"
54 #include "G4Neutron.hh"
55 #include "G4ParticleGun.hh"
56 #include "G4PhysicalVolumeStore.hh"
57 #include "G4SystemOfUnits.hh"
58 #include "G4Tubs.hh"
59 #include "G4VPhysicalVolume.hh"
60 #include "Randomize.hh"
61 #include "globals.hh"
62
63 //....oooOO0OOooo........oooOO0OOooo........oo
64 FFPrimaryGeneratorAction::FFPrimaryGeneratorAc
65 : G4VUserPrimaryGeneratorAction(),
66 #ifndef NDEBUG
67 fEventNumber(0),
68 #endif // NDEBUG
69 fH2OPhysical(NULL),
70 fNeutronPhysical(NULL),
71 fNeutronSolid(NULL),
72 fParticleGun(new G4ParticleGun(1)),
73 fTankPhysical(NULL)
74 {
75 fParticleGun->SetParticleDefinition(G4Neutro
76 fParticleGun->SetParticleEnergy(4.5 * MeV);
77 }
78
79 //....oooOO0OOooo........oooOO0OOooo........oo
80 void FFPrimaryGeneratorAction::GeneratePrimari
81 {
82 #ifndef NDEBUG
83 G4cout << "Shooting event " << ++fEventNumbe
84 #endif // NDEBUG
85 const G4ThreeVector sourceCenter = GetNeutro
86
87 // Sample the neutron source location
88 const G4double radius = fNeutronSolid->GetOu
89 const G4double z = fNeutronSolid->GetZHalfLe
90 G4ThreeVector randomLocation;
91 randomLocation.setRThetaPhi(radius * std::sq
92 randomLocation.setZ(z * (G4UniformRand() - 0
93 G4ThreeVector location(randomLocation.x() +
94 randomLocation.y() +
95 randomLocation.z() +
96 #ifndef NDEBUG
97 G4cout << "Emission Location: r: " << locati
98 #endif // NDEBUG
99
100 // Sample the neutron emission direction
101 G4ThreeVector direction;
102 direction.setRThetaPhi(1.0, std::acos(G4Unif
103 (G4UniformRand() * 2
104 #ifndef NDEBUG
105 G4cout << "Emission Direction: r: " << direc
106 #endif // NDEBUG
107
108 // Load the event
109 fParticleGun->SetParticlePosition(location);
110 fParticleGun->SetParticleMomentumDirection(d
111 fParticleGun->GeneratePrimaryVertex(event);
112 }
113
114 //....oooOO0OOooo........oooOO0OOooo........oo
115 G4ThreeVector FFPrimaryGeneratorAction::GetNeu
116 {
117 // Get the dimensions of the neutron source
118 if (fNeutronSolid == NULL) {
119 G4LogicalVolume* temp = G4LogicalVolumeSto
120 if (temp != NULL) {
121 fNeutronSolid = dynamic_cast<G4Tubs*>(te
122 }
123
124 if (fNeutronSolid == NULL) {
125 G4Exception(
126 "FFPrimaryGeneratorAction::"
127 "GeneratePrimaries(G4Event*)",
128 "Neutron source solid volume not found
129 }
130 }
131
132 // Get the position of the neutron source wi
133 if (fNeutronPhysical == NULL) {
134 fNeutronPhysical = G4PhysicalVolumeStore::
135 }
136
137 if (fNeutronPhysical == NULL) {
138 G4Exception("FFPrimaryGeneratorAction::Get
139 "Neutron source physical volum
140 "This run will be aborted");
141 }
142
143 // Get the position of the water within the
144 if (fH2OPhysical == NULL) {
145 fH2OPhysical = G4PhysicalVolumeStore::GetI
146
147 if (fH2OPhysical == NULL) {
148 G4Exception(
149 "FFPrimaryGeneratorAction::"
150 "GetNeutronSourceCenter(void)",
151 "Tank H2O physical volume not found",
152 }
153 }
154
155 // Get the position of the tank within the w
156 if (fTankPhysical == NULL) {
157 fTankPhysical = G4PhysicalVolumeStore::Get
158
159 if (fTankPhysical == NULL) {
160 G4Exception(
161 "FFPrimaryGeneratorAction::"
162 "GetNeutronSourceCenter(void)",
163 "Tank physical volume not found", Even
164 }
165 }
166
167 return fNeutronPhysical->GetTranslation() +
168 + fTankPhysical->GetTranslation();
169 }
170
171 //....oooOO0OOooo........oooOO0OOooo........oo
172 FFPrimaryGeneratorAction::~FFPrimaryGeneratorA
173 {
174 delete fParticleGun;
175 }
176