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 // original by H.P. Wellisch 27 // modified by J.L. Chuma, TRIUMF, 19-Nov-1996 28 // last modified: 27-Mar-1997 29 // Chr. Volcker, 10-Nov-1997: new methods and class variables. 30 // M.G. Pia, 2 Oct 1998: modified GetFermiMomentum (original design was 31 // the source of memory leaks) 32 // G.Folger, spring 2010: add integer A/Z interface 33 // A. Ribon, autumn 2021: extended to hypernuclei 34 35 #ifndef G4Nucleus_h 36 #define G4Nucleus_h 1 37 // Class Description 38 // This class knows how to describe a nucleus; 39 // to be used in your physics implementation (not physics list) in case you need this physics. 40 // Class Description - End 41 42 43 #include "globals.hh" 44 #include "G4ThreeVector.hh" 45 #include "G4ParticleTypes.hh" 46 #include "G4ReactionProduct.hh" 47 #include "G4DynamicParticle.hh" 48 #include "G4ReactionProductVector.hh" 49 #include "Randomize.hh" 50 51 class G4Nucleus 52 { 53 public: 54 55 G4Nucleus(); 56 G4Nucleus(const G4double A, const G4double Z, const G4int numberOfLambdas = 0); 57 G4Nucleus(const G4int A, const G4int Z, const G4int numberOfLambdas = 0); 58 G4Nucleus(const G4Material* aMaterial); 59 60 ~G4Nucleus(); 61 62 inline G4Nucleus( const G4Nucleus &right ) 63 { *this = right; } 64 65 inline G4Nucleus& operator = (const G4Nucleus& right) 66 { 67 if (this != &right) { 68 theA=right.theA; 69 theZ=right.theZ; 70 theL=right.theL; 71 aEff=right.aEff; 72 zEff=right.zEff; 73 fIsotope = right.fIsotope; 74 pnBlackTrackEnergy=right.pnBlackTrackEnergy; 75 dtaBlackTrackEnergy=right.dtaBlackTrackEnergy; 76 pnBlackTrackEnergyfromAnnihilation = 77 right.pnBlackTrackEnergyfromAnnihilation; 78 dtaBlackTrackEnergyfromAnnihilation = 79 right.dtaBlackTrackEnergyfromAnnihilation; 80 theTemp = right.theTemp; 81 excitationEnergy = right.excitationEnergy; 82 momentum = right.momentum; 83 fermiMomentum = right.fermiMomentum; 84 } 85 return *this; 86 } 87 88 inline G4bool operator==( const G4Nucleus &right ) const 89 { return ( this == (G4Nucleus *) &right ); } 90 91 inline G4bool operator!=( const G4Nucleus &right ) const 92 { return ( this != (G4Nucleus *) &right ); } 93 94 void ChooseParameters( const G4Material *aMaterial ); 95 96 void SetParameters( const G4double A, const G4double Z, const G4int numberOfLambdas = 0 ); 97 void SetParameters( const G4int A, const G4int Z, const G4int numberOfLambdas = 0 ); 98 99 inline G4int GetA_asInt() const 100 { return theA; } 101 102 inline G4int GetN_asInt() const 103 { return theA-theZ-theL; } 104 105 inline G4int GetZ_asInt() const 106 { return theZ; } 107 108 inline G4int GetL() const // Number of Lambdas (in the case of a hypernucleus) 109 { return theL; } 110 111 inline const G4Isotope* GetIsotope() 112 { return fIsotope; } 113 114 inline void SetIsotope(const G4Isotope* iso) 115 { 116 fIsotope = iso; 117 if(iso) { 118 theZ = iso->GetZ(); 119 theA = iso->GetN(); 120 theL = 0; 121 aEff = theA; 122 zEff = theZ; 123 } 124 } 125 126 G4DynamicParticle *ReturnTargetParticle() const; 127 128 G4double AtomicMass( const G4double A, const G4double Z, const G4int numberOfLambdas = 0 ) const; 129 G4double AtomicMass( const G4int A, const G4int Z, const G4int numberOfLambdas = 0 ) const; 130 131 G4double GetThermalPz( const G4double mass, const G4double temp ) const; 132 133 G4ReactionProduct GetThermalNucleus(G4double aMass, G4double temp=-1) const; 134 135 G4ReactionProduct GetBiasedThermalNucleus(G4double aMass, G4ThreeVector aVelocity, G4double temp=-1) const; 136 137 void DoKinematicsOfThermalNucleus(const G4double mu, const G4double vT_norm, const G4ThreeVector& aVelocity, 138 G4ReactionProduct& result) const; 139 140 G4double Cinema( G4double kineticEnergy ); 141 142 G4double EvaporationEffects( G4double kineticEnergy ); 143 144 G4double AnnihilationEvaporationEffects(G4double kineticEnergy, G4double ekOrg); 145 146 inline G4double GetPNBlackTrackEnergy() const 147 { return pnBlackTrackEnergy; } 148 149 inline G4double GetDTABlackTrackEnergy() const 150 { return dtaBlackTrackEnergy; } 151 152 inline G4double GetAnnihilationPNBlackTrackEnergy() const 153 { return pnBlackTrackEnergyfromAnnihilation; } 154 155 inline G4double GetAnnihilationDTABlackTrackEnergy() const 156 { return dtaBlackTrackEnergyfromAnnihilation; } 157 158 // ****************** methods introduced by ChV *********************** 159 // return fermi momentum 160 G4ThreeVector GetFermiMomentum(); 161 162 /* 163 // return particle to be absorbed. 164 G4DynamicParticle* ReturnAbsorbingParticle(G4double weight); 165 */ 166 167 // final nucleus fragmentation. Return List of particles 168 // which should be used for further tracking. 169 G4ReactionProductVector* Fragmentate(); 170 171 172 // excitation Energy... 173 void AddExcitationEnergy(G4double anEnergy); 174 175 176 // momentum of absorbed Particles .. 177 void AddMomentum(const G4ThreeVector aMomentum); 178 179 // return excitation Energy 180 G4double GetEnergyDeposit() {return excitationEnergy; } 181 182 183 184 // ****************************** end ChV ****************************** 185 186 187 private: 188 189 G4int theA; 190 G4int theZ; 191 G4int theL; // Number of Lambdas (in the case of hypernucleus) 192 G4double aEff; // effective atomic weight 193 G4double zEff; // effective atomic number 194 195 const G4Isotope* fIsotope; 196 197 G4double pnBlackTrackEnergy; // the kinetic energy available for 198 // proton/neutron black track particles 199 G4double dtaBlackTrackEnergy; // the kinetic energy available for 200 // deuteron/triton/alpha particles 201 G4double pnBlackTrackEnergyfromAnnihilation; 202 // kinetic energy available for proton/neutron black 203 // track particles based on baryon annihilation 204 G4double dtaBlackTrackEnergyfromAnnihilation; 205 // kinetic energy available for deuteron/triton/alpha 206 // black track particles based on baryon annihilation 207 208 209 // ************************** member variables by ChV ******************* 210 // Excitation Energy leading to evaporation or deexcitation. 211 G4double excitationEnergy; 212 213 // Momentum, accumulated by absorbing Particles 214 G4ThreeVector momentum; 215 216 // Fermi Gas model: at present, we assume constant nucleon density for all 217 // nuclei. The radius of a nucleon is taken to be 1 fm. 218 // see for example S.Fl"ugge, Encyclopedia of Physics, Vol XXXIX, 219 // Structure of Atomic Nuclei (Berlin-Gottingen-Heidelberg, 1957) page 426. 220 221 // maximum momentum possible from fermi gas model: 222 G4double fermiMomentum; 223 G4double theTemp; // temperature 224 // ****************************** end ChV ****************************** 225 226 }; 227 228 #endif 229 230