Geant4 Cross Reference |
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 // =============== Begin Documentation Comments =============== 28 //! 29 //! \file FFPrimaryGeneratorAction.cc 30 //! \author B. Wendt (brycen.linn.wendt@cern.ch) 31 //! \date June 06, 2014 32 //! 33 //! \brief Implementation of the FFPrimaryGeneratorAction class 34 //! 35 //! \details Generates 4.5 MeV neutrons from the solid volume defined in 36 //! FFDetectorConstruction as "NeturonSource" and shoots them off 37 //! in a isotropically sampled direction 38 //! 39 // ================ End Documentation Comments ================ 40 // 41 // Modified: 42 // 43 // 23-06-14 BWendt 44 // Implemented "GetNeutronSourceCenter()" and added code to correctly sample 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........oooOO0OOooo........oooOO0OOooo...... 64 FFPrimaryGeneratorAction::FFPrimaryGeneratorAction() 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(G4Neutron::Definition()); 76 fParticleGun->SetParticleEnergy(4.5 * MeV); 77 } 78 79 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 80 void FFPrimaryGeneratorAction::GeneratePrimaries(G4Event* event) 81 { 82 #ifndef NDEBUG 83 G4cout << "Shooting event " << ++fEventNumber << G4endl; 84 #endif // NDEBUG 85 const G4ThreeVector sourceCenter = GetNeutronSourceCenter(); 86 87 // Sample the neutron source location 88 const G4double radius = fNeutronSolid->GetOuterRadius(); 89 const G4double z = fNeutronSolid->GetZHalfLength() * 2; 90 G4ThreeVector randomLocation; 91 randomLocation.setRThetaPhi(radius * std::sqrt(G4UniformRand()), G4UniformRand() * 180 * deg, 0); 92 randomLocation.setZ(z * (G4UniformRand() - 0.5)); 93 G4ThreeVector location(randomLocation.x() + sourceCenter.x(), 94 randomLocation.y() + sourceCenter.y(), 95 randomLocation.z() + sourceCenter.z()); 96 #ifndef NDEBUG 97 G4cout << "Emission Location: r: " << location << G4endl; 98 #endif // NDEBUG 99 100 // Sample the neutron emission direction 101 G4ThreeVector direction; 102 direction.setRThetaPhi(1.0, std::acos(G4UniformRand() * 2 - 1), 103 (G4UniformRand() * 2 - 1) * 180 * deg); 104 #ifndef NDEBUG 105 G4cout << "Emission Direction: r: " << direction << G4endl; 106 #endif // NDEBUG 107 108 // Load the event 109 fParticleGun->SetParticlePosition(location); 110 fParticleGun->SetParticleMomentumDirection(direction); 111 fParticleGun->GeneratePrimaryVertex(event); 112 } 113 114 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 115 G4ThreeVector FFPrimaryGeneratorAction::GetNeutronSourceCenter(void) 116 { 117 // Get the dimensions of the neutron source 118 if (fNeutronSolid == NULL) { 119 G4LogicalVolume* temp = G4LogicalVolumeStore::GetInstance()->GetVolume("NeutronSource"); 120 if (temp != NULL) { 121 fNeutronSolid = dynamic_cast<G4Tubs*>(temp->GetSolid()); 122 } 123 124 if (fNeutronSolid == NULL) { 125 G4Exception( 126 "FFPrimaryGeneratorAction::" 127 "GeneratePrimaries(G4Event*)", 128 "Neutron source solid volume not found", EventMustBeAborted, "This run will be aborted"); 129 } 130 } 131 132 // Get the position of the neutron source within the water 133 if (fNeutronPhysical == NULL) { 134 fNeutronPhysical = G4PhysicalVolumeStore::GetInstance()->GetVolume("NeutronSource"); 135 } 136 137 if (fNeutronPhysical == NULL) { 138 G4Exception("FFPrimaryGeneratorAction::GetNeutronSourceCenter(void)", 139 "Neutron source physical volume not found", EventMustBeAborted, 140 "This run will be aborted"); 141 } 142 143 // Get the position of the water within the tank 144 if (fH2OPhysical == NULL) { 145 fH2OPhysical = G4PhysicalVolumeStore::GetInstance()->GetVolume("Tank_H2O"); 146 147 if (fH2OPhysical == NULL) { 148 G4Exception( 149 "FFPrimaryGeneratorAction::" 150 "GetNeutronSourceCenter(void)", 151 "Tank H2O physical volume not found", EventMustBeAborted, "This run will be aborted"); 152 } 153 } 154 155 // Get the position of the tank within the world 156 if (fTankPhysical == NULL) { 157 fTankPhysical = G4PhysicalVolumeStore::GetInstance()->GetVolume("Tank_Wall"); 158 159 if (fTankPhysical == NULL) { 160 G4Exception( 161 "FFPrimaryGeneratorAction::" 162 "GetNeutronSourceCenter(void)", 163 "Tank physical volume not found", EventMustBeAborted, "This run will be aborted"); 164 } 165 } 166 167 return fNeutronPhysical->GetTranslation() + fH2OPhysical->GetTranslation() 168 + fTankPhysical->GetTranslation(); 169 } 170 171 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 172 FFPrimaryGeneratorAction::~FFPrimaryGeneratorAction() 173 { 174 delete fParticleGun; 175 } 176