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 // GEANT4 class header file 28 // 29 // Class: G4IonStoppingData 30 // 31 // Base class: G4VIonDEDXTable 32 // 33 // Author: Anton Lechner (Anton.Lechner@cern.ch) 34 // 35 // First implementation: 03. 11. 2009 36 // 37 // Modifications: 38 // 39 // 40 // Class description: Class which can read ion stopping power data from 41 // $G4LEDATA/ion_stopping_data 42 // 43 // Comments: 44 // 45 // =========================================================================== 46 47 #ifndef G4IONSTOPPINGDATA_HH 48 #define G4IONSTOPPINGDATA_HH 49 50 #include "G4VIonDEDXTable.hh" 51 #include "globals.hh" 52 53 #include <map> 54 #include <utility> 55 #include <vector> 56 57 class G4IonStoppingData : public G4VIonDEDXTable 58 { 59 public: 60 G4IonStoppingData(const G4String& dir, G4bool icru); 61 ~G4IonStoppingData() override; 62 G4IonStoppingData(G4IonStoppingData&) = delete; 63 const G4IonStoppingData& operator=(const G4IonStoppingData&) = delete; 64 65 // Function for checking the availability of stopping power tables 66 // for a given ion-material couple, where the material consists of 67 // a single element only. 68 G4bool IsApplicable(G4int atomicNumberIon, // Atomic number of ion 69 G4int atomicNumberElem // Atomic number of elemental material 70 ) override; 71 72 // Function for checking the availability of stopping power tables 73 // for given ion-material couples. 74 G4bool IsApplicable(G4int atomicNumberIon, // Atomic number of ion 75 const G4String& matIdentifier // Name or chemical formula of material 76 ) override; 77 78 // Function which invokes the read/build process of physics vectors from 79 // files in G4LEDATA 80 G4bool BuildPhysicsVector(G4int ionZ, const G4String& matName) override; 81 82 // Function which invokes the read/build process of physics vectors from 83 // files in G4LEDATA 84 G4bool BuildPhysicsVector(G4int ionZ, G4int matZ) override; 85 86 // Function returning the stopping power vector for given ion-material 87 // couples, where the material consists of a single element only. 88 G4PhysicsVector* GetPhysicsVector(G4int atomicNumberIon, // Atomic number of ion 89 G4int atomicNumberElem // Atomic number of elemental material 90 ) override; 91 92 // Function returning the stopping power vector for given ion-material 93 // couples. 94 G4PhysicsVector* GetPhysicsVector(G4int atomicNumberIon, // Atomic number of ion 95 const G4String& matIdenfier // Name or chemical formula of material 96 ) override; 97 98 // Function returning the stopping power value for given ion-material 99 // couples, where the material consists of a single element only, and 100 // given energy. 101 G4double GetDEDX(G4double kinEnergyPerNucleon, // Kinetic energy per nucleon 102 G4int atomicNumberIon, // Atomic number of ion 103 G4int atomicNumberElem // Atomic number of elemental material 104 ); 105 106 // Function returning the stopping power value for given ion-material 107 // couples and given energy. 108 G4double GetDEDX(G4double kinEnergyPerNucleon, // Kinetic energy per nucleon 109 G4int atomicNumberIon, // Atomic number of ion 110 const G4String& matIdentifier // Name or chemical formula of material 111 ); 112 113 // Function for adding dE/dx vector for an elemental materials. The last 114 // argument only applies to elemental materials. 115 G4bool AddPhysicsVector(G4PhysicsVector* physicsVector, // Physics vector 116 G4int atomicNumberIon, // Atomic number of ion 117 const G4String& matIdentifier // Name or chemical formula of material 118 ); 119 120 // Function for adding dE/dx vector for an elemental materials. The last 121 // argument only applies to elemental materials. 122 G4bool AddPhysicsVector(G4PhysicsVector* physicsVector, // Physics vector 123 G4int atomicNumberIon, // Atomic number of ion 124 G4int atomicNumberElem // Atomic number of elemental material 125 ); 126 127 // Function for removing dE/dx vector for a compound materials 128 G4bool RemovePhysicsVector(G4int atomicNumberIon, // Atomic number of ion 129 const G4String& matIdentifier // Name or chemical formula of material 130 ); 131 // Function for removing dE/dx vector for a compound materials 132 G4bool RemovePhysicsVector(G4int atomicNumberIon, // Atomic number of ion 133 G4int atomicNumberElem // Atomic number of elemental material 134 ); 135 // Function deleting all physics vectors and clearing the maps 136 void ClearTable(); 137 138 // Function printing the ion-material pairs of available vectors to stdout 139 void DumpMap(); 140 141 private: 142 // Subdirectory of G4LEDATA 143 G4String subDir; 144 145 using G4IonDEDXKeyElem = std::pair<G4int, G4int>; 146 using G4IonDEDXKeyMat = std::pair<G4int, G4String>; 147 148 using G4IonDEDXMapElem = std::map<G4IonDEDXKeyElem, G4PhysicsVector*>; 149 using G4IonDEDXMapMat = std::map<G4IonDEDXKeyMat, G4PhysicsVector*>; 150 151 G4IonDEDXMapElem dedxMapElements; 152 G4IonDEDXMapMat dedxMapMaterials; 153 154 G4bool fICRU90; 155 }; 156 157 #endif // G4IONSTOPPINGDATA_HH 158