Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer 3 // * License and Disclaimer * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. 9 // * include a list of copyright holders. * 10 // * 10 // * * 11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. * 17 // * 17 // * * 18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************* 24 // ******************************************************************** 25 // 25 // 26 /// \file TrackingAction.hh 26 /// \file TrackingAction.hh 27 /// \brief Definition of the TrackingAction cl 27 /// \brief Definition of the TrackingAction class 28 // 28 // 29 // << 29 // 30 30 31 //....oooOO0OOooo........oooOO0OOooo........oo 31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 32 //....oooOO0OOooo........oooOO0OOooo........oo 32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 33 33 34 #ifndef TrackingAction_h << 34 #ifndef TrackingAction_h 35 #define TrackingAction_h 1 35 #define TrackingAction_h 1 36 36 37 #include "G4UserTrackingAction.hh" << 38 #include "globals.hh" 37 #include "globals.hh" 39 << 38 #include "G4UserTrackingAction.hh" 40 #include <array> 39 #include <array> 41 40 42 class Run; 41 class Run; 43 42 44 //....oooOO0OOooo........oooOO0OOooo........oo 43 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 45 44 46 class TrackingAction : public G4UserTrackingAc << 45 class TrackingAction : public G4UserTrackingAction { 47 { << 46 // We are using this class to monitor the average multiplicity, the average 48 // We are using this class to monitor the << 47 // kinetic energy, and the average total energy flow (i.e. sum of the 49 // kinetic energy, and the average total e << 48 // kinetic energies) of different particle types as they are produced 50 // kinetic energies) of different particle << 49 // inside the calorimeter. 51 // inside the calorimeter. << 50 // The aim is then to try to correlate some changes in these (more primitive) 52 // The aim is then to try to correlate som << 51 // quantities with the observed changes in the (more indirect and complex) 53 // quantities with the observed changes in << 52 // particle fluences. 54 // particle fluences. << 55 public: 53 public: 56 TrackingAction(); 54 TrackingAction(); 57 ~TrackingAction() override = default; 55 ~TrackingAction() override = default; 58 void PreUserTrackingAction(const G4Track*) << 56 void PreUserTrackingAction( const G4Track* ) override; 59 void PostUserTrackingAction(const G4Track* << 57 void PostUserTrackingAction( const G4Track* ) override; 60 58 61 void Initialize(); 59 void Initialize(); 62 // This method is called by RunAction::Beg 60 // This method is called by RunAction::BeginOfRunAction for the 63 // initialization of the tracking-action a 61 // initialization of the tracking-action at the beginning of each Run. 64 62 65 void SetRunPointer(Run* inputValue = nullp << 63 void SetRunPointer( Run* inputValue = nullptr ) { fRunPtr = inputValue; } 66 // This method is called by RunAction::Beg 64 // This method is called by RunAction::BeginOfRunAction for providing to the 67 // tracking-action the pointer to the run 65 // tracking-action the pointer to the run object at the beginning of each Run. 68 // This pointer is then used to pass the i 66 // This pointer is then used to pass the information collected by the tracking-action 69 // to the run object. 67 // to the run object. 70 << 68 71 static const G4int fkNumberScoringVolumes << 69 static const G4int fkNumberScoringVolumes = 1; // calorimeter 72 static const G4int fkNumberKinematicRegion 70 static const G4int fkNumberKinematicRegions = 3; // all, below 20 MeV, above 20 MeV 73 static const G4int fkNumberParticleTypes = << 71 static const G4int fkNumberParticleTypes = 11; // all, e, gamma, mu, nu, pi, n, p, ions, 74 72 // other-mesons, other-baryons 75 static const G4int fkNumberCombinations = 73 static const G4int fkNumberCombinations = 76 fkNumberScoringVolumes * fkNumberKinemat << 74 fkNumberScoringVolumes*fkNumberKinematicRegions*fkNumberParticleTypes; 77 static const std::array<G4String, fkNumber << 75 static const std::array< G4String, fkNumberScoringVolumes > fkArrayScoringVolumeNames; 78 static const std::array<G4String, fkNumber << 76 static const std::array< G4String, fkNumberKinematicRegions > fkArrayKinematicRegionNames; 79 static const std::array<G4String, fkNumber << 77 static const std::array< G4String, fkNumberParticleTypes > fkArrayParticleTypeNames; 80 static G4int GetIndex(const G4int iScoring << 78 static G4int GetIndex( const G4int iScoringVolume, const G4int iKinematicRegion, 81 const G4int iParticl << 79 const G4int iParticleType ); 82 80 83 private: 81 private: 84 Run* fRunPtr; // Pointer to the Run objec 82 Run* fRunPtr; // Pointer to the Run object 85 << 83 86 std::array<G4long, fkNumberCombinations> f << 84 std::array< G4long, fkNumberCombinations > fArrayMultiplicities; 87 std::array<G4double, fkNumberCombinations> << 85 std::array< G4double, fkNumberCombinations > fArraySumKineticEnergies; 88 // Keep record of the fkNumber of particle 86 // Keep record of the fkNumber of particles and their kinetic energy at production, 89 // according to the particle type and thei 87 // according to the particle type and their kinetic energy range (below/above 20 MeV). 90 }; 88 }; 91 89 92 //....oooOO0OOooo........oooOO0OOooo........oo 90 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 93 91 94 #endif 92 #endif 95 93