FissionFragment/src/FFPrimaryGeneratorAction.cc

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 // -------------------------------------------------------------
0027 //  =============== Begin Documentation Comments ===============
0028 //!
0029 //! \file       FFPrimaryGeneratorAction.cc
0030 //! \author     B. Wendt (brycen.linn.wendt@cern.ch)
0031 //! \date       June 06, 2014
0032 //!
0033 //! \brief      Implementation of the FFPrimaryGeneratorAction class
0034 //!
0035 //! \details    Generates 4.5 MeV neutrons from the solid volume defined in
0036 //!             FFDetectorConstruction as "NeturonSource" and shoots them off
0037 //!             in a isotropically sampled direction
0038 //!
0039 //  ================ End Documentation Comments ================
0040 //
0041 //  Modified:
0042 //
0043 //  23-06-14                                              BWendt
0044 //  Implemented "GetNeutronSourceCenter()" and added code to correctly sample
0045 //  the neutron starting location
0046 //
0047 // -------------------------------------------------------------
0048
0049 #include "FFPrimaryGeneratorAction.hh"
0050
0051 #include "G4Event.hh"
0052 #include "G4LogicalVolume.hh"
0053 #include "G4LogicalVolumeStore.hh"
0054 #include "G4Neutron.hh"
0055 #include "G4ParticleGun.hh"
0056 #include "G4PhysicalVolumeStore.hh"
0057 #include "G4SystemOfUnits.hh"
0058 #include "G4Tubs.hh"
0059 #include "G4VPhysicalVolume.hh"
0060 #include "Randomize.hh"
0061 #include "globals.hh"
0062
0063 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0064 FFPrimaryGeneratorAction::FFPrimaryGeneratorAction()
0065   : G4VUserPrimaryGeneratorAction(),
0066 #ifndef NDEBUG
0067     fEventNumber(0),
0068 #endif  // NDEBUG
0069     fH2OPhysical(NULL),
0070     fNeutronPhysical(NULL),
0071     fNeutronSolid(NULL),
0072     fParticleGun(new G4ParticleGun(1)),
0073     fTankPhysical(NULL)
0074 {
0075   fParticleGun->SetParticleDefinition(G4Neutron::Definition());
0076   fParticleGun->SetParticleEnergy(4.5 * MeV);
0077 }
0078
0079 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0080 void FFPrimaryGeneratorAction::GeneratePrimaries(G4Event* event)
0081 {
0082 #ifndef NDEBUG
0083   G4cout << "Shooting event " << ++fEventNumber << G4endl;
0084 #endif  // NDEBUG
0085   const G4ThreeVector sourceCenter = GetNeutronSourceCenter();
0086
0087   // Sample the neutron source location
0088   const G4double radius = fNeutronSolid->GetOuterRadius();
0089   const G4double z = fNeutronSolid->GetZHalfLength() * 2;
0090   G4ThreeVector randomLocation;
0091   randomLocation.setRThetaPhi(radius * std::sqrt(G4UniformRand()), G4UniformRand() * 180 * deg, 0);
0092   randomLocation.setZ(z * (G4UniformRand() - 0.5));
0093   G4ThreeVector location(randomLocation.x() + sourceCenter.x(),
0094                          randomLocation.y() + sourceCenter.y(),
0095                          randomLocation.z() + sourceCenter.z());
0096 #ifndef NDEBUG
0097   G4cout << "Emission Location: r: " << location << G4endl;
0098 #endif  // NDEBUG
0099
0100   // Sample the neutron emission direction
0101   G4ThreeVector direction;
0102   direction.setRThetaPhi(1.0, std::acos(G4UniformRand() * 2 - 1),
0103                          (G4UniformRand() * 2 - 1) * 180 * deg);
0104 #ifndef NDEBUG
0105   G4cout << "Emission Direction: r: " << direction << G4endl;
0106 #endif  // NDEBUG
0107
0108   // Load the event
0109   fParticleGun->SetParticlePosition(location);
0110   fParticleGun->SetParticleMomentumDirection(direction);
0111   fParticleGun->GeneratePrimaryVertex(event);
0112 }
0113
0114 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0115 G4ThreeVector FFPrimaryGeneratorAction::GetNeutronSourceCenter(void)
0116 {
0117   // Get the dimensions of the neutron source
0118   if (fNeutronSolid == NULL) {
0119     G4LogicalVolume* temp = G4LogicalVolumeStore::GetInstance()->GetVolume("NeutronSource");
0120     if (temp != NULL) {
0121       fNeutronSolid = dynamic_cast<G4Tubs*>(temp->GetSolid());
0122     }
0123
0124     if (fNeutronSolid == NULL) {
0125       G4Exception(
0126         "FFPrimaryGeneratorAction::"
0127         "GeneratePrimaries(G4Event*)",
0128         "Neutron source solid volume not found", EventMustBeAborted, "This run will be aborted");
0129     }
0130   }
0131
0132   // Get the position of the neutron source within the water
0133   if (fNeutronPhysical == NULL) {
0134     fNeutronPhysical = G4PhysicalVolumeStore::GetInstance()->GetVolume("NeutronSource");
0135   }
0136
0137   if (fNeutronPhysical == NULL) {
0138     G4Exception("FFPrimaryGeneratorAction::GetNeutronSourceCenter(void)",
0139                 "Neutron source physical volume not found", EventMustBeAborted,
0140                 "This run will be aborted");
0141   }
0142
0143   // Get the position of the water within the tank
0144   if (fH2OPhysical == NULL) {
0145     fH2OPhysical = G4PhysicalVolumeStore::GetInstance()->GetVolume("Tank_H2O");
0146
0147     if (fH2OPhysical == NULL) {
0148       G4Exception(
0149         "FFPrimaryGeneratorAction::"
0150         "GetNeutronSourceCenter(void)",
0151         "Tank H2O physical volume not found", EventMustBeAborted, "This run will be aborted");
0152     }
0153   }
0154
0155   // Get the position of the tank within the world
0156   if (fTankPhysical == NULL) {
0157     fTankPhysical = G4PhysicalVolumeStore::GetInstance()->GetVolume("Tank_Wall");
0158
0159     if (fTankPhysical == NULL) {
0160       G4Exception(
0161         "FFPrimaryGeneratorAction::"
0162         "GetNeutronSourceCenter(void)",
0163         "Tank physical volume not found", EventMustBeAborted, "This run will be aborted");
0164     }
0165   }
0166
0167   return fNeutronPhysical->GetTranslation() + fH2OPhysical->GetTranslation()
0168          + fTankPhysical->GetTranslation();
0169 }
0170
0171 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
0172 FFPrimaryGeneratorAction::~FFPrimaryGeneratorAction()
0173 {
0174   delete fParticleGun;
0175 }