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 SteppingAction.hh 27 /// \brief Definition of the SteppingAction cl 28 // 29 // 30 31 //....oooOO0OOooo........oooOO0OOooo........oo 32 //....oooOO0OOooo........oooOO0OOooo........oo 33 34 #ifndef SteppingAction_H 35 #define SteppingAction_H 1 36 37 #include "G4ThreeVector.hh" 38 #include "G4UserSteppingAction.hh" 39 #include "globals.hh" 40 41 #include <array> 42 43 class Run; 44 45 //....oooOO0OOooo........oooOO0OOooo........oo 46 47 class SteppingAction : public G4UserSteppingAc 48 { 49 public: 50 SteppingAction(); 51 ~SteppingAction() override = default; 52 53 void UserSteppingAction(const G4Step*) ove 54 // This is the main method where the step 55 // the scoring volumes are collected, and 56 // are filled up in the Run object where t 57 // printed out at the end of the Run). 58 // (For simplicity and brevity, we avoid h 59 // some statistics ourself, which will be 60 61 void Initialize(); 62 // This method is called by RunAction::Beg 63 // initialization of the stepping-action a 64 // This is necessary because different run 65 // types, kinetic energies, and detector c 66 67 void SetRunPointer(Run* inputValue = nullp 68 // This method is called by RunAction::Beg 69 // stepping-action the pointer to the run 70 // This pointer is then used to pass the i 71 // to the run object. 72 73 G4double GetCubicVolumeScoringUpDown() con 74 G4double GetCubicVolumeScoringSide() const 75 // The cubic-volumes of the scoring volume 76 // the sum of step lengths in those scorin 77 // Notice that two of the three scoring vo 78 // have the same cubic-volume, that we cal 79 80 static const G4int fkNumberScoringVolumes 81 static const G4int fkNumberKinematicRegion 82 static const G4int fkNumberParticleTypes = 83 84 static const G4int fkNumberCombinations = 85 fkNumberScoringVolumes * fkNumberKinemat 86 static const std::array<G4String, fkNumber 87 static const std::array<G4String, fkNumber 88 static const std::array<G4String, fkNumber 89 static G4int GetIndex(const G4int iScoring 90 const G4int iParticl 91 92 private: 93 Run* fRunPtr; // Pointer to the Run objec 94 G4int fPrimaryParticleId; 95 G4double fPrimaryParticleEnergy; 96 G4ThreeVector fPrimaryParticleDirection; 97 G4String fAbsorberMaterialName; 98 G4String fActiveMaterialName; 99 G4bool fIsFirstStepOfTheEvent; 100 G4bool fIsFirstStepInAbsorberLayer; 101 G4bool fIsFirstStepInActiveLayer; 102 G4bool fIsFirstStepInScoringUpDown; 103 G4bool fIsFirstStepInScoringSide; 104 G4double fCubicVolumeScoringUpDown; 105 G4double fCubicVolumeScoringSide; 106 107 std::array<G4double, fkNumberCombinations> 108 // Array to collect the sum of step length 109 // according to the various cases (kinemat 110 // Note that the fluence in a scoring volu 111 // in that scoring volume divided by the c 112 }; 113 114 //....oooOO0OOooo........oooOO0OOooo........oo 115 116 #endif 117