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 /// \file medical/dna/range/src/PrimaryGenerat 27 /// \brief Implementation of the PrimaryGenera 28 // 29 // $Id: PrimaryGeneratorAction.cc 73024 2013-0 30 // 31 //....oooOO0OOooo........oooOO0OOooo........oo 32 //....oooOO0OOooo........oooOO0OOooo........oo 33 34 #include "PrimaryGeneratorAction.hh" 35 36 #include "DetectorConstruction.hh" 37 38 #include "G4Electron.hh" 39 #include "G4Event.hh" 40 #include "G4Gamma.hh" 41 #include "G4ParticleDefinition.hh" 42 #include "G4ParticleGun.hh" 43 #include "G4ParticleTable.hh" 44 #include "G4PhysicalConstants.hh" 45 #include "G4RandomDirection.hh" 46 #include "G4RunManager.hh" 47 #include "G4StateManager.hh" 48 #include "G4SystemOfUnits.hh" 49 #include "Randomize.hh" 50 51 //....oooOO0OOooo........oooOO0OOooo........oo 52 53 PrimaryGeneratorAction::PrimaryGeneratorAction 54 : G4VUserPrimaryGeneratorAction(), G4VStateD 55 { 56 fDetector = dynamic_cast<const DetectorConst 57 G4RunManager::GetRunManager()->GetUserDete 58 G4int n_particle = 1; 59 fParticleGun = new G4ParticleGun(n_particle) 60 } 61 62 //....oooOO0OOooo........oooOO0OOooo........oo 63 64 PrimaryGeneratorAction::~PrimaryGeneratorActio 65 { 66 G4StateManager::GetStateManager()->Deregiste 67 if (fParticleGun) delete fParticleGun; 68 } 69 70 //....oooOO0OOooo........oooOO0OOooo........oo 71 72 void PrimaryGeneratorAction::GeneratePrimaries 73 { 74 G4double R = fDetector->GetNPRadius(); 75 G4double Ryz = 9999; 76 77 // G4ParticleDefinition *particle = G4Electr 78 // fParticleGun->SetParticleDefinition(parti 79 80 G4double scale = 1.; 81 82 G4double xpos = 0, ypos = 0, zpos = 0; 83 while (!(Ryz < R)) { 84 ypos = (2. * G4UniformRand() - 1.) * R; 85 zpos = (2. * G4UniformRand() - 1.) * R; 86 Ryz = std::sqrt(ypos * ypos + zpos * zpos) 87 } 88 xpos = std::sqrt(R * R - (ypos * ypos + zpos 89 fParticleGun->SetParticlePosition(G4ThreeVec 90 fParticleGun->SetParticleMomentumDirection(G 91 92 fParticleGun->GeneratePrimaryVertex(anEvent) 93 } 94 95 //....oooOO0OOooo........oooOO0OOooo........oo 96 97 G4bool PrimaryGeneratorAction::Notify(G4Applic 98 { 99 if (requestedState == G4State_Idle) { 100 if (fParticleGun != 0) return true; 101 G4double R = fDetector->GetNPRadius(); 102 G4double Ryz = 9999; 103 ; 104 105 // G4ParticleDefinition *particle = G4Elec 106 // fParticleGun->SetParticleDefinition(par 107 108 fParticleGun->SetParticleMomentumDirection 109 110 while (!(Ryz < R)) { 111 G4double ypos = (2. * G4UniformRand() - 112 G4double zpos = (2. * G4UniformRand() - 113 Ryz = std::sqrt(ypos * ypos + zpos * zpo 114 G4double xpos = std::sqrt(R * R - (ypos 115 fParticleGun->SetParticlePosition(G4Thre 116 } 117 // fParticleGun->SetParticlePosition(G4Thr 118 // fParticleGun->GeneratePrimaryVertex(anE 119 } 120 121 return true; 122 } 123