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 /// \file VG01PrimaryGeneratorAction.cc 28 /// \brief Implementation of the VG01PrimaryGeneratorAction class 29 // 30 // Authors: J. Apostolakis & S. Wenzel (CERN) 2018-2021 31 // 32 // Started from FullCMS code by Mihaly Novak (CERN) 2017 33 34 #include "VG01PrimaryGeneratorAction.hh" 35 36 #include "G4Event.hh" 37 #include "G4ParticleDefinition.hh" 38 #include "G4ParticleGun.hh" 39 #include "G4ParticleTable.hh" 40 #include "G4RandomDirection.hh" 41 #include "G4SystemOfUnits.hh" 42 43 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 44 45 VG01PrimaryGeneratorAction::VG01PrimaryGeneratorAction(EGeneratorMode mode) 46 : G4VUserPrimaryGeneratorAction(), fParticleGun(0), fMode(mode) 47 { 48 G4int n_particle = 1; 49 fParticleGun = new G4ParticleGun(n_particle); 50 51 G4ParticleTable* particleTable = G4ParticleTable::GetParticleTable(); 52 G4String particleName; 53 fParticleGun->SetParticleDefinition(particleTable->FindParticle(particleName = "geantino")); 54 fParticleGun->SetParticleEnergy(1.0 * GeV); 55 fParticleGun->SetParticlePosition(G4ThreeVector(0, 0, 0)); // -1.2*m, 0.1, 0.1)); 56 } 57 58 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 59 60 VG01PrimaryGeneratorAction::~VG01PrimaryGeneratorAction() 61 { 62 delete fParticleGun; 63 } 64 65 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 66 67 void VG01PrimaryGeneratorAction::GeneratePrimaries(G4Event* anEvent) 68 { 69 G4ThreeVector dir(1.0, 0.0, 0.0); 70 71 // G4ThreeVector FixedMode[3] = 72 // { G4ThreeVector(1,0,0), G4ThreeVector(1,0.1,0), G4ThreeVector( 1, 0, 0.1)}; 73 G4ThreeVector AxisModeVectors[3] = {G4ThreeVector(1, 0, 0), G4ThreeVector(0, 1, 0), 74 G4ThreeVector(0, 0, 1.0)}; 75 76 if (fMode == kFreeMode) { 77 // Do not fix the direction 78 // -- UI commands to particle mode can change it! 79 fParticleGun->GeneratePrimaryVertex(anEvent); 80 return; 81 //---- 82 } 83 else if (fMode == kFixedMode) { 84 G4int i = anEvent->GetEventID() % 3; 85 dir = G4ThreeVector(1.0, 0.0, 0.0); 86 switch (i) { 87 case 0: 88 break; 89 case 1: 90 dir.setY(0.1); 91 break; 92 case 2: 93 dir.setZ(0.1); 94 break; 95 } 96 } 97 else if (fMode == kUniformMode) { 98 dir = G4RandomDirection(); 99 } 100 else if (fMode == kAxisMode) { 101 G4int i = anEvent->GetEventID() % 3; 102 dir = AxisModeVectors[i]; 103 } 104 105 fParticleGun->SetParticleMomentumDirection(dir); 106 fParticleGun->GeneratePrimaryVertex(anEvent); 107 } 108