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 // G4MaterialCutsCouple 26 // G4MaterialCutsCouple 27 // 27 // 28 // Class description: 28 // Class description: 29 // 29 // 30 // The same material can be used in regions wi 30 // The same material can be used in regions with different cut values, 31 // and physics processes must prepare several 31 // and physics processes must prepare several different cross-sections 32 // for that material. 32 // for that material. 33 // The G4ProductionCutsTable has G4MaterialCut 33 // The G4ProductionCutsTable has G4MaterialCutsCouple objects, each of 34 // which consists of a material paired with a 34 // which consists of a material paired with a cut value. 35 // G4MaterialCutsCouples objects are numbered 35 // G4MaterialCutsCouples objects are numbered with an index which is the 36 // same as the index of a G4PhysicsVector for 36 // same as the index of a G4PhysicsVector for the corresponding 37 // G4MaterialCutsCouple in the G4PhysicsTable. 37 // G4MaterialCutsCouple in the G4PhysicsTable. 38 // The list of G4MaterialCutsCouple objects us 38 // The list of G4MaterialCutsCouple objects used in the geometry setup 39 // is updated before starting the event loop i 39 // is updated before starting the event loop in each run. 40 40 41 // Author: H.Kurashige, 17 September 2002 - Fi 41 // Author: H.Kurashige, 17 September 2002 - First implementation 42 // ------------------------------------------- 42 // -------------------------------------------------------------------- 43 #ifndef G4MaterialCutsCouple_hh 43 #ifndef G4MaterialCutsCouple_hh 44 #define G4MaterialCutsCouple_hh 1 44 #define G4MaterialCutsCouple_hh 1 45 45 46 #include "globals.hh" 46 #include "globals.hh" 47 #include "G4ios.hh" 47 #include "G4ios.hh" 48 #include "G4ProductionCuts.hh" 48 #include "G4ProductionCuts.hh" 49 49 50 class G4Material; 50 class G4Material; 51 51 52 class G4MaterialCutsCouple 52 class G4MaterialCutsCouple 53 { 53 { 54 public: 54 public: 55 55 56 G4MaterialCutsCouple(); 56 G4MaterialCutsCouple(); 57 G4MaterialCutsCouple(const G4Material*, G4 57 G4MaterialCutsCouple(const G4Material*, G4ProductionCuts* cut = nullptr); 58 // Constructors 58 // Constructors 59 59 60 virtual ~G4MaterialCutsCouple(); 60 virtual ~G4MaterialCutsCouple(); 61 // Destructor 61 // Destructor 62 62 63 G4MaterialCutsCouple(const G4MaterialCutsC 63 G4MaterialCutsCouple(const G4MaterialCutsCouple& right); 64 G4MaterialCutsCouple& operator=(const G4Ma 64 G4MaterialCutsCouple& operator=(const G4MaterialCutsCouple& right); 65 // Copy constructor & assignment operato 65 // Copy constructor & assignment operator 66 66 67 G4bool operator==(const G4MaterialCutsCoup 67 G4bool operator==(const G4MaterialCutsCouple& right) const; 68 G4bool operator!=(const G4MaterialCutsCoup 68 G4bool operator!=(const G4MaterialCutsCouple& right) const; 69 // Equality operators 69 // Equality operators 70 70 71 void SetMaterial(const G4Material*); 71 void SetMaterial(const G4Material*); 72 // Set pointer to material 72 // Set pointer to material 73 73 74 const G4Material* GetMaterial() const; 74 const G4Material* GetMaterial() const; 75 // Get pointer to material 75 // Get pointer to material 76 76 77 void SetProductionCuts(G4ProductionCuts*); 77 void SetProductionCuts(G4ProductionCuts*); 78 // Set pointer to production cuts 78 // Set pointer to production cuts 79 79 80 G4ProductionCuts* GetProductionCuts() cons 80 G4ProductionCuts* GetProductionCuts() const; 81 // Get pointer to production cuts 81 // Get pointer to production cuts 82 82 83 G4bool IsRecalcNeeded() const; 83 G4bool IsRecalcNeeded() const; 84 // Return true if cut and/or material ha 84 // Return true if cut and/or material has been modified 85 // after last calculation of PhysicsTabl 85 // after last calculation of PhysicsTable 86 86 87 void PhysicsTableUpdated(); 87 void PhysicsTableUpdated(); 88 // Inform end of calculation of Physics 88 // Inform end of calculation of Physics Table 89 89 90 void SetIndex(G4int idx); 90 void SetIndex(G4int idx); 91 G4int GetIndex() const; 91 G4int GetIndex() const; 92 // Set/Get the index number in G4Product 92 // Set/Get the index number in G4ProductionCutsTable 93 93 94 void SetUseFlag(G4bool flg = true); 94 void SetUseFlag(G4bool flg = true); 95 G4bool IsUsed() const; 95 G4bool IsUsed() const; 96 96 97 private: 97 private: 98 98 99 G4bool isMaterialModified = fal 99 G4bool isMaterialModified = false; 100 const G4Material* fMaterial = nullptr; 100 const G4Material* fMaterial = nullptr; 101 G4ProductionCuts* fCuts = nullptr; 101 G4ProductionCuts* fCuts = nullptr; 102 G4int indexNumber = -1; 102 G4int indexNumber = -1; 103 G4bool isUsedInGeometry = false 103 G4bool isUsedInGeometry = false; 104 }; 104 }; 105 105 106 // ------------------ 106 // ------------------ 107 // Inline methods 107 // Inline methods 108 // ------------------ 108 // ------------------ 109 109 110 inline 110 inline 111 void G4MaterialCutsCouple::SetIndex(G4int idx) 111 void G4MaterialCutsCouple::SetIndex(G4int idx) 112 { 112 { 113 indexNumber = idx; 113 indexNumber = idx; 114 } 114 } 115 115 116 inline 116 inline 117 G4int G4MaterialCutsCouple::GetIndex() const 117 G4int G4MaterialCutsCouple::GetIndex() const 118 { 118 { 119 return indexNumber; 119 return indexNumber; 120 } 120 } 121 121 122 inline 122 inline 123 void G4MaterialCutsCouple::SetUseFlag(G4bool f 123 void G4MaterialCutsCouple::SetUseFlag(G4bool flg) 124 { 124 { 125 isUsedInGeometry = flg; 125 isUsedInGeometry = flg; 126 } 126 } 127 127 128 inline 128 inline 129 G4bool G4MaterialCutsCouple::IsUsed() const 129 G4bool G4MaterialCutsCouple::IsUsed() const 130 { 130 { 131 return isUsedInGeometry; 131 return isUsedInGeometry; 132 } 132 } 133 133 134 inline 134 inline 135 void G4MaterialCutsCouple::SetProductionCuts(G 135 void G4MaterialCutsCouple::SetProductionCuts(G4ProductionCuts* aCut) 136 { 136 { 137 fCuts = aCut; 137 fCuts = aCut; 138 } 138 } 139 139 140 inline 140 inline 141 G4ProductionCuts* G4MaterialCutsCouple::GetPro 141 G4ProductionCuts* G4MaterialCutsCouple::GetProductionCuts() const 142 { 142 { 143 return fCuts; 143 return fCuts; 144 } 144 } 145 145 146 inline 146 inline 147 G4bool G4MaterialCutsCouple::operator==(const 147 G4bool G4MaterialCutsCouple::operator==(const G4MaterialCutsCouple& right) const 148 { 148 { 149 return (this == &right); 149 return (this == &right); 150 } 150 } 151 151 152 inline 152 inline 153 G4bool G4MaterialCutsCouple::operator!=(const 153 G4bool G4MaterialCutsCouple::operator!=(const G4MaterialCutsCouple& right) const 154 { 154 { 155 return (this != &right); 155 return (this != &right); 156 } 156 } 157 157 158 inline 158 inline 159 void G4MaterialCutsCouple::SetMaterial(const 159 void G4MaterialCutsCouple::SetMaterial(const G4Material* material) 160 { 160 { 161 fMaterial = material; 161 fMaterial = material; 162 isMaterialModified = true; 162 isMaterialModified = true; 163 } 163 } 164 164 165 inline 165 inline 166 const G4Material* G4MaterialCutsCouple::GetMat 166 const G4Material* G4MaterialCutsCouple::GetMaterial() const 167 { 167 { 168 return fMaterial; 168 return fMaterial; 169 } 169 } 170 170 171 inline 171 inline 172 G4bool G4MaterialCutsCouple::IsRecalcNeeded() 172 G4bool G4MaterialCutsCouple::IsRecalcNeeded() const 173 { 173 { 174 G4bool isCutModified = false; 174 G4bool isCutModified = false; 175 if (fCuts != nullptr ) isCutModified = fCuts 175 if (fCuts != nullptr ) isCutModified = fCuts->IsModified(); 176 return (isMaterialModified || isCutModified) 176 return (isMaterialModified || isCutModified); 177 } 177 } 178 178 179 inline 179 inline 180 void G4MaterialCutsCouple::PhysicsTableUpdated 180 void G4MaterialCutsCouple::PhysicsTableUpdated() 181 { 181 { 182 if (fCuts != nullptr ) fCuts->PhysicsTableUp 182 if (fCuts != nullptr ) fCuts->PhysicsTableUpdated(); 183 isMaterialModified = false; 183 isMaterialModified = false; 184 } 184 } 185 185 186 #endif 186 #endif 187 187