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 // G4ParticleDefinition inline methods impleme 26 // G4ParticleDefinition inline methods implementation 27 // 27 // 28 // Authors: G.Cosmo, 2 December 1995 - Design, 28 // Authors: G.Cosmo, 2 December 1995 - Design, based on object model 29 // M.Asai, 29 January 1996 - First im 29 // M.Asai, 29 January 1996 - First implementation 30 // ------------------------------------------- 30 // -------------------------------------------------------------------- 31 31 32 inline G4int G4ParticleDefinition::GetInstance 32 inline G4int G4ParticleDefinition::GetInstanceID() const 33 { 33 { 34 return g4particleDefinitionInstanceID; 34 return g4particleDefinitionInstanceID; 35 } 35 } 36 36 37 inline G4bool G4ParticleDefinition::GetPDGStab 37 inline G4bool G4ParticleDefinition::GetPDGStable() const 38 { 38 { 39 if (IsGeneralIon()) { 39 if (IsGeneralIon()) { 40 return (GetIonLifeTime() < 0.); 40 return (GetIonLifeTime() < 0.); 41 } 41 } 42 42 43 return thePDGStable; 43 return thePDGStable; 44 } 44 } 45 45 46 inline G4double G4ParticleDefinition::GetPDGLi 46 inline G4double G4ParticleDefinition::GetPDGLifeTime() const 47 { 47 { 48 return thePDGLifeTime; 48 return thePDGLifeTime; 49 } 49 } 50 50 51 inline G4double G4ParticleDefinition::GetIonLi 51 inline G4double G4ParticleDefinition::GetIonLifeTime() const 52 { 52 { 53 //-- No longer needed to access to G4IonTabl 53 //-- No longer needed to access to G4IonTable. 54 //-- Method GetIonLifeTime() itself is kept 54 //-- Method GetIonLifeTime() itself is kept for compatibility 55 55 56 return thePDGLifeTime; 56 return thePDGLifeTime; 57 } 57 } 58 58 59 inline G4ParticleTable* G4ParticleDefinition:: 59 inline G4ParticleTable* G4ParticleDefinition::GetParticleTable() const 60 { 60 { 61 return theParticleTable; 61 return theParticleTable; 62 } 62 } 63 63 64 inline G4DecayTable* G4ParticleDefinition::Get 64 inline G4DecayTable* G4ParticleDefinition::GetDecayTable() const 65 { 65 { 66 return theDecayTable; 66 return theDecayTable; 67 } 67 } 68 68 69 inline void G4ParticleDefinition::SetDecayTabl 69 inline void G4ParticleDefinition::SetDecayTable(G4DecayTable* aDecayTable) 70 { 70 { 71 theDecayTable = aDecayTable; 71 theDecayTable = aDecayTable; 72 } 72 } 73 73 74 inline void G4ParticleDefinition::SetVerboseLe 74 inline void G4ParticleDefinition::SetVerboseLevel(G4int value) 75 { 75 { 76 verboseLevel = value; 76 verboseLevel = value; 77 } 77 } 78 78 79 inline G4int G4ParticleDefinition::GetVerboseL 79 inline G4int G4ParticleDefinition::GetVerboseLevel() const 80 { 80 { 81 return verboseLevel; 81 return verboseLevel; 82 } 82 } 83 83 84 inline G4ProcessManager* G4ParticleDefinition: 84 inline G4ProcessManager* G4ParticleDefinition::GetMasterProcessManager() const 85 { 85 { 86 return theProcessManagerShadow; 86 return theProcessManagerShadow; 87 } 87 } 88 88 89 inline void G4ParticleDefinition::SetMasterPro 89 inline void G4ParticleDefinition::SetMasterProcessManager(G4ProcessManager* aNewPM) 90 { 90 { 91 theProcessManagerShadow = aNewPM; 91 theProcessManagerShadow = aNewPM; 92 } 92 } 93 93 94 inline G4int G4ParticleDefinition::GetQuarkCon 94 inline G4int G4ParticleDefinition::GetQuarkContent(G4int flavor) const 95 { 95 { 96 G4int content = 0; 96 G4int content = 0; 97 if ((flavor > 0) && (flavor <= NumberOfQuark 97 if ((flavor > 0) && (flavor <= NumberOfQuarkFlavor)) { 98 content = theQuarkContent[flavor - 1]; 98 content = theQuarkContent[flavor - 1]; 99 } 99 } 100 else { 100 else { 101 #ifdef G4VERBOSE 101 #ifdef G4VERBOSE 102 if (verboseLevel > 0) { 102 if (verboseLevel > 0) { 103 std::ostringstream message; 103 std::ostringstream message; 104 message << "Invalid Quark Flavor: " << f 104 message << "Invalid Quark Flavor: " << flavor; 105 G4Exception("G4ParticleDefinition::GetQu 105 G4Exception("G4ParticleDefinition::GetQuarkContent()", "InvalidFlavor", JustWarning, message); 106 } 106 } 107 #endif 107 #endif 108 } 108 } 109 return content; 109 return content; 110 } 110 } 111 111 112 inline G4int G4ParticleDefinition::GetAntiQuar 112 inline G4int G4ParticleDefinition::GetAntiQuarkContent(G4int flavor) const 113 { 113 { 114 G4int content = 0; 114 G4int content = 0; 115 if ((flavor > 0) && (flavor <= NumberOfQuark 115 if ((flavor > 0) && (flavor <= NumberOfQuarkFlavor)) { 116 content = theAntiQuarkContent[flavor - 1]; 116 content = theAntiQuarkContent[flavor - 1]; 117 } 117 } 118 else { 118 else { 119 #ifdef G4VERBOSE 119 #ifdef G4VERBOSE 120 if (verboseLevel > 0) { 120 if (verboseLevel > 0) { 121 std::ostringstream message; 121 std::ostringstream message; 122 message << "Invalid Quark Flavor: " << f 122 message << "Invalid Quark Flavor: " << flavor; 123 G4Exception("G4ParticleDefinition::GetAn 123 G4Exception("G4ParticleDefinition::GetAntiQuarkContent()", "InvalidFlavor", JustWarning, 124 message); 124 message); 125 } 125 } 126 #endif 126 #endif 127 } 127 } 128 return content; 128 return content; 129 } 129 } 130 130 131 inline void G4ParticleDefinition::SetParticleS 131 inline void G4ParticleDefinition::SetParticleSubType(const G4String& subtype) 132 { 132 { 133 theParticleSubType = subtype; 133 theParticleSubType = subtype; 134 } 134 } 135 135 136 inline void G4ParticleDefinition::SetAntiPDGEn 136 inline void G4ParticleDefinition::SetAntiPDGEncoding(G4int aEncoding) 137 { 137 { 138 theAntiPDGEncoding = aEncoding; 138 theAntiPDGEncoding = aEncoding; 139 } 139 } 140 140 141 inline G4bool G4ParticleDefinition::GetApplyCu 141 inline G4bool G4ParticleDefinition::GetApplyCutsFlag() const 142 { 142 { 143 return fApplyCutsFlag; 143 return fApplyCutsFlag; 144 } 144 } 145 145 146 inline void G4ParticleDefinition::SetAtomicNum 146 inline void G4ParticleDefinition::SetAtomicNumber(G4int i) 147 { 147 { 148 theAtomicNumber = i; 148 theAtomicNumber = i; 149 } 149 } 150 150 151 inline G4int G4ParticleDefinition::GetAtomicNu 151 inline G4int G4ParticleDefinition::GetAtomicNumber() const 152 { 152 { 153 return theAtomicNumber; 153 return theAtomicNumber; 154 } 154 } 155 155 156 inline void G4ParticleDefinition::SetAtomicMas 156 inline void G4ParticleDefinition::SetAtomicMass(G4int i) 157 { 157 { 158 theAtomicMass = i; 158 theAtomicMass = i; 159 } 159 } 160 160 161 inline G4int G4ParticleDefinition::GetAtomicMa 161 inline G4int G4ParticleDefinition::GetAtomicMass() const 162 { 162 { 163 return theAtomicMass; 163 return theAtomicMass; 164 } 164 } 165 165 166 inline void G4ParticleDefinition::SetPDGMagnet 166 inline void G4ParticleDefinition::SetPDGMagneticMoment(G4double magneticMoment) 167 { 167 { 168 thePDGMagneticMoment = magneticMoment; 168 thePDGMagneticMoment = magneticMoment; 169 } 169 } 170 170 171 inline G4bool G4ParticleDefinition::IsGeneralI 171 inline G4bool G4ParticleDefinition::IsGeneralIon() const 172 { 172 { 173 return isGeneralIon; 173 return isGeneralIon; 174 } 174 } 175 175 176 inline G4bool G4ParticleDefinition::IsMuonicAt 176 inline G4bool G4ParticleDefinition::IsMuonicAtom() const 177 { 177 { 178 return isMuonicAtom; 178 return isMuonicAtom; 179 } 179 } 180 180 181 inline G4int G4ParticleDefinition::GetParticle 181 inline G4int G4ParticleDefinition::GetParticleDefinitionID() const 182 { 182 { 183 return g4particleDefinitionInstanceID; 183 return g4particleDefinitionInstanceID; 184 } 184 } 185 185 186 inline G4bool G4ParticleDefinition::IsHypernuc 186 inline G4bool G4ParticleDefinition::IsHypernucleus() const 187 { 187 { 188 return GetNumberOfLambdasInHypernucleus() > 188 return GetNumberOfLambdasInHypernucleus() > 0; 189 } 189 } 190 190 191 inline G4int G4ParticleDefinition::GetNumberOf 191 inline G4int G4ParticleDefinition::GetNumberOfLambdasInHypernucleus() const 192 { 192 { 193 G4int numberOfLambdas = 0; 193 G4int numberOfLambdas = 0; 194 // PDG code of hypernuclei: 10LZZZAAAI 194 // PDG code of hypernuclei: 10LZZZAAAI 195 if (thePDGEncoding > 0 && thePDGEncoding / 1 195 if (thePDGEncoding > 0 && thePDGEncoding / 1000000000 != 0) { 196 numberOfLambdas = (thePDGEncoding / 100000 196 numberOfLambdas = (thePDGEncoding / 10000000) % 100; 197 } 197 } 198 return numberOfLambdas; 198 return numberOfLambdas; 199 } 199 } 200 200 201 inline G4bool G4ParticleDefinition::IsAntiHype 201 inline G4bool G4ParticleDefinition::IsAntiHypernucleus() const 202 { 202 { 203 return GetNumberOfAntiLambdasInAntiHypernucl 203 return GetNumberOfAntiLambdasInAntiHypernucleus() > 0; 204 } 204 } 205 205 206 inline G4int G4ParticleDefinition::GetNumberOf 206 inline G4int G4ParticleDefinition::GetNumberOfAntiLambdasInAntiHypernucleus() const 207 { 207 { 208 G4int numberOfAntiLambdas = 0; 208 G4int numberOfAntiLambdas = 0; 209 // PDG code of anti-hypernuclei: -10LZZZAAAI 209 // PDG code of anti-hypernuclei: -10LZZZAAAI 210 if (thePDGEncoding < 0 && thePDGEncoding / 1 210 if (thePDGEncoding < 0 && thePDGEncoding / 1000000000 != 0) { 211 numberOfAntiLambdas = (std::abs(thePDGEnco 211 numberOfAntiLambdas = (std::abs(thePDGEncoding) / 10000000) % 100; 212 } 212 } 213 return numberOfAntiLambdas; 213 return numberOfAntiLambdas; 214 } 214 } 215 215