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 // G4ParticlePropertyData class implementation << 27 // 26 // 28 // Author: H.Kurashige, 9 June 2003 << 27 // class G4ParticlePropertyData 29 // ------------------------------------------- << 28 // 30 << 29 // Implementation 31 #include "G4ParticlePropertyData.hh" << 30 // >> 31 // History: >> 32 // first implementation by H Kurashige 9 June 2003 >> 33 // Add magnetic moment by H Kurashige Mar 2007 32 34 33 #include "G4PhysicalConstants.hh" << 34 #include "G4SystemOfUnits.hh" << 35 #include "G4ios.hh" 35 #include "G4ios.hh" 36 #include "globals.hh" 36 #include "globals.hh" >> 37 #include "G4PhysicalConstants.hh" >> 38 #include "G4SystemOfUnits.hh" >> 39 #include "G4ParticlePropertyData.hh" >> 40 37 41 38 G4ParticlePropertyData::G4ParticlePropertyData << 42 ///////////////////////////////////////////////////////////// 39 : theParticleName(particleName) << 43 G4ParticlePropertyData::~G4ParticlePropertyData() 40 { 44 { 41 for (std::size_t flv = 0; flv < NumberOfQuar << 42 theQuarkContent[flv] = 0; << 43 theAntiQuarkContent[flv] = 0; << 44 } << 45 } 45 } 46 46 47 G4ParticlePropertyData::G4ParticlePropertyData << 47 ///////////////////////////////////////////////////////////// 48 { << 48 G4ParticlePropertyData::G4ParticlePropertyData(const G4String& particleName): 49 verboseLevel = right.verboseLevel; << 49 theParticleName(particleName), 50 theParticleName = right.theParticleName; << 50 thePDGMass(0.0), 51 thePDGMass = right.thePDGMass; << 51 thePDGWidth(0.0), 52 thePDGWidth = right.thePDGWidth; << 52 thePDGCharge(0.0), 53 thePDGCharge = right.thePDGCharge; << 53 thePDGiSpin(0), 54 thePDGiSpin = right.thePDGiSpin; << 54 thePDGiParity(0), 55 thePDGiParity = right.thePDGiParity; << 55 thePDGiConjugation(0), 56 thePDGiConjugation = right.thePDGiConjugatio << 56 thePDGiGParity(0), 57 thePDGiGParity = right.thePDGiGParity; << 57 thePDGiIsospin(0), 58 thePDGiIsospin = right.thePDGiIsospin; << 58 thePDGiIsospin3(0), 59 thePDGiIsospin3 = right.thePDGiIsospin3; << 59 thePDGMagneticMoment(0.0), 60 thePDGMagneticMoment = right.thePDGMagneticM << 60 theLeptonNumber(0), 61 theLeptonNumber = right.theLeptonNumber; << 61 theBaryonNumber(0), 62 theBaryonNumber = right.theBaryonNumber; << 62 thePDGEncoding(0), 63 thePDGEncoding = right.thePDGEncoding; << 63 theAntiPDGEncoding(0), 64 theAntiPDGEncoding = right.theAntiPDGEncodin << 64 thePDGLifeTime(-1.0), 65 for (std::size_t flv = 0; flv < NumberOfQuar << 65 fPDGMassModified(false), 66 theQuarkContent[flv] = right.theQuarkConte << 66 fPDGWidthModified(false), 67 theAntiQuarkContent[flv] = right.theAntiQu << 67 fPDGChargeModified(false), >> 68 fPDGiSpinModified(false), >> 69 fPDGiParityModified(false), >> 70 fPDGiConjugationModified(false), >> 71 fPDGiGParityModified(false), >> 72 fPDGiIsospinModified(false), >> 73 fPDGiIsospin3Modified(false), >> 74 fPDGIsospinModified(false), >> 75 fPDGIsospin3Modified(false), >> 76 fPDGMagneticMomentModified(false), >> 77 fLeptonNumberModified(false), >> 78 fBaryonNumberModified(false), >> 79 fPDGEncodingModified(false), >> 80 fAntiPDGEncodingModified(false), >> 81 fQuarkContentModified(false), >> 82 fAntiQuarkContentModified(false), >> 83 fPDGLifeTimeModified(false), >> 84 verboseLevel(1) >> 85 { >> 86 for (size_t flv=0; flv<NumberOfQuarkFlavor; ++flv) { >> 87 theQuarkContent[flv] =0; >> 88 theAntiQuarkContent[flv]=0; 68 } 89 } 69 thePDGLifeTime = right.thePDGLifeTime; << 70 } 90 } 71 91 72 G4ParticlePropertyData& G4ParticlePropertyData << 92 >> 93 //////////////////////// >> 94 G4ParticlePropertyData::G4ParticlePropertyData(const G4ParticlePropertyData &right): >> 95 fPDGMassModified(false), >> 96 fPDGWidthModified(false), >> 97 fPDGChargeModified(false), >> 98 fPDGiSpinModified(false), >> 99 fPDGiParityModified(false), >> 100 fPDGiConjugationModified(false), >> 101 fPDGiGParityModified(false), >> 102 fPDGiIsospinModified(false), >> 103 fPDGiIsospin3Modified(false), >> 104 fPDGIsospinModified(false), >> 105 fPDGIsospin3Modified(false), >> 106 fPDGMagneticMomentModified(false), >> 107 fLeptonNumberModified(false), >> 108 fBaryonNumberModified(false), >> 109 fPDGEncodingModified(false), >> 110 fAntiPDGEncodingModified(false), >> 111 fQuarkContentModified(false), >> 112 fAntiQuarkContentModified(false), >> 113 fPDGLifeTimeModified(false) >> 114 { >> 115 verboseLevel = right.verboseLevel; >> 116 theParticleName = right.theParticleName; >> 117 thePDGMass = right.thePDGMass; >> 118 thePDGWidth = right. thePDGWidth; >> 119 thePDGCharge = right.thePDGCharge; >> 120 thePDGiSpin = right.thePDGiSpin; >> 121 thePDGiParity = right.thePDGiParity; >> 122 thePDGiConjugation = right.thePDGiConjugation; >> 123 thePDGiGParity = right.thePDGiGParity; >> 124 thePDGiIsospin = right.thePDGiIsospin; >> 125 thePDGiIsospin3 = right.thePDGiIsospin3; >> 126 thePDGMagneticMoment = right.thePDGMagneticMoment; >> 127 theLeptonNumber = right.theLeptonNumber; >> 128 theBaryonNumber = right.theBaryonNumber; >> 129 thePDGEncoding = right.thePDGEncoding; >> 130 theAntiPDGEncoding = right.theAntiPDGEncoding; >> 131 for (size_t flv=0; flv<NumberOfQuarkFlavor; ++flv) { >> 132 theQuarkContent[flv] = right.theQuarkContent[flv]; >> 133 theAntiQuarkContent[flv]= right.theAntiQuarkContent[flv]; >> 134 } >> 135 thePDGLifeTime = right.thePDGLifeTime; >> 136 } >> 137 >> 138 //////////////////////// >> 139 G4ParticlePropertyData & G4ParticlePropertyData::operator=(const G4ParticlePropertyData &right) 73 { 140 { 74 if (this != &right) { 141 if (this != &right) { 75 verboseLevel = right.verboseLevel; << 142 verboseLevel = right.verboseLevel; 76 theParticleName = right.theParticleName; << 143 theParticleName = right.theParticleName; 77 thePDGMass = right.thePDGMass; << 144 thePDGMass = right.thePDGMass; 78 thePDGWidth = right.thePDGWidth; << 145 thePDGWidth = right. thePDGWidth; 79 thePDGCharge = right.thePDGCharge; << 146 thePDGCharge = right.thePDGCharge; 80 thePDGiSpin = right.thePDGiSpin; << 147 thePDGiSpin = right.thePDGiSpin; 81 thePDGiParity = right.thePDGiParity; << 148 thePDGiParity = right.thePDGiParity; 82 thePDGiConjugation = right.thePDGiConjugat << 149 thePDGiConjugation = right.thePDGiConjugation; 83 thePDGiGParity = right.thePDGiGParity; << 150 thePDGiGParity = right.thePDGiGParity; 84 thePDGiIsospin = right.thePDGiIsospin; << 151 thePDGiIsospin = right.thePDGiIsospin; 85 thePDGiIsospin3 = right.thePDGiIsospin3; << 152 thePDGiIsospin3 = right.thePDGiIsospin3; 86 thePDGMagneticMoment = right.thePDGMagneti << 153 thePDGMagneticMoment = right.thePDGMagneticMoment; 87 theLeptonNumber = right.theLeptonNumber; << 154 theLeptonNumber = right.theLeptonNumber; 88 theBaryonNumber = right.theBaryonNumber; << 155 theBaryonNumber = right.theBaryonNumber; 89 thePDGEncoding = right.thePDGEncoding; << 156 thePDGEncoding = right.thePDGEncoding; 90 theAntiPDGEncoding = right.theAntiPDGEncod << 157 theAntiPDGEncoding = right.theAntiPDGEncoding; 91 for (std::size_t flv = 0; flv < NumberOfQu << 158 for (size_t flv=0; flv<NumberOfQuarkFlavor; ++flv) { 92 theQuarkContent[flv] = right.theQuarkCon << 159 theQuarkContent[flv] = right.theQuarkContent[flv]; 93 theAntiQuarkContent[flv] = right.theAnti << 160 theAntiQuarkContent[flv]= right.theAntiQuarkContent[flv]; 94 } 161 } 95 thePDGLifeTime = right.thePDGLifeTime; << 162 thePDGLifeTime = right.thePDGLifeTime; 96 fPDGMassModified = true; << 163 fPDGMassModified = true; 97 fPDGWidthModified = true; << 164 fPDGWidthModified = true; 98 fPDGChargeModified = true; << 165 fPDGChargeModified = true; 99 fPDGiSpinModified = true; << 166 fPDGiSpinModified = true; 100 fPDGiParityModified = true; << 167 fPDGiParityModified = true; 101 fPDGiConjugationModified = true; << 168 fPDGiConjugationModified = true; 102 fPDGiGParityModified = true; << 169 fPDGiGParityModified = true; 103 fPDGiIsospinModified = true; << 170 fPDGiIsospinModified = true; 104 fPDGiIsospin3Modified = true; << 171 fPDGiIsospin3Modified = true; 105 fPDGIsospinModified = true; << 172 fPDGIsospinModified = true; 106 fPDGIsospin3Modified = true; << 173 fPDGIsospin3Modified = true; 107 fPDGMagneticMomentModified = true; 174 fPDGMagneticMomentModified = true; 108 fLeptonNumberModified = true; << 175 fLeptonNumberModified = true; 109 fBaryonNumberModified = true; << 176 fBaryonNumberModified = true; 110 fPDGEncodingModified = true; << 177 fPDGEncodingModified = true; 111 fAntiPDGEncodingModified = true; << 178 fAntiPDGEncodingModified = true; 112 fQuarkContentModified = true; << 179 fQuarkContentModified = true; 113 fAntiQuarkContentModified = true; << 180 fAntiQuarkContentModified = true; 114 fPDGLifeTimeModified = true; << 181 fPDGLifeTimeModified = true; 115 } 182 } 116 return *this; 183 return *this; 117 } 184 } 118 << 185 119 G4bool G4ParticlePropertyData::operator==(cons << 186 //////////////////////// >> 187 G4bool G4ParticlePropertyData::operator==(const G4ParticlePropertyData &right) const 120 { 188 { 121 return (this == &right); 189 return (this == &right); 122 } 190 } 123 191 124 G4bool G4ParticlePropertyData::operator!=(cons << 192 //////////////////////// >> 193 G4bool G4ParticlePropertyData::operator!=(const G4ParticlePropertyData &right) const 125 { 194 { 126 return (this != &right); 195 return (this != &right); 127 } 196 } 128 197 >> 198 //////////////////////// 129 void G4ParticlePropertyData::Print() const 199 void G4ParticlePropertyData::Print() const 130 { 200 { 131 #ifdef G4VERBOSE 201 #ifdef G4VERBOSE 132 G4cout << " Particle Name : " << theParticle 202 G4cout << " Particle Name : " << theParticleName << G4endl; 133 G4cout << " PDG particle code : " << thePDGE 203 G4cout << " PDG particle code : " << thePDGEncoding; 134 G4cout << " [PDG anti-particle code: " << th << 204 G4cout << " [PDG anti-particle code: " << this->GetAntiPDGEncoding() << "]"<< G4endl; 135 G4cout << " Mass [GeV/c2] : " << thePDGMass << 205 G4cout << " Mass [GeV/c2] : " << thePDGMass/GeV ; 136 G4cout << " Width : " << thePDGWidth / G << 206 G4cout << " Width : " << thePDGWidth/GeV << G4endl; 137 G4cout << " Lifetime [nsec] : " << thePDGLif << 207 G4cout << " Lifetime [nsec] : " << thePDGLifeTime/ns << G4endl; 138 G4cout << " Charge [e]: " << thePDGCharge / << 208 G4cout << " Charge [e]: " << thePDGCharge/eplus << G4endl; 139 G4cout << " Spin : " << thePDGiSpin << "/2" 209 G4cout << " Spin : " << thePDGiSpin << "/2" << G4endl; 140 G4cout << " Parity : " << thePDGiParity << G 210 G4cout << " Parity : " << thePDGiParity << G4endl; 141 G4cout << " Charge conjugation : " << thePDG 211 G4cout << " Charge conjugation : " << thePDGiConjugation << G4endl; 142 G4cout << " Isospin : (I,Iz): (" << thePDGiI << 212 G4cout << " Isospin : (I,Iz): (" << thePDGiIsospin <<"/2"; 143 G4cout << " , " << thePDGiIsospin3 << "/2 ) 213 G4cout << " , " << thePDGiIsospin3 << "/2 ) " << G4endl; 144 G4cout << " GParity : " << thePDGiGParity << 214 G4cout << " GParity : " << thePDGiGParity << G4endl; 145 G4cout << " MagneticMoment [MeV/T]: "; 215 G4cout << " MagneticMoment [MeV/T]: "; 146 if (thePDGMagneticMoment != 0.0) { 216 if (thePDGMagneticMoment != 0.0) { 147 G4cout << thePDGMagneticMoment / MeV * tes << 217 G4cout << thePDGMagneticMoment/MeV*tesla << G4endl; 148 } << 218 }else { 149 else { << 150 G4cout << "not defined " << G4endl; 219 G4cout << "not defined " << G4endl; 151 } 220 } 152 G4cout << " Lepton number : " << theLeptonNu 221 G4cout << " Lepton number : " << theLeptonNumber; 153 G4cout << " Baryon number : " << theBaryonNu << 222 G4cout << " Baryon number : " << theBaryonNumber << G4endl; 154 G4cout << " Quark contents (d,u,s,c,b,t) 223 G4cout << " Quark contents (d,u,s,c,b,t) : " << theQuarkContent[0]; 155 G4cout << ", " << theQuarkContent[1]; 224 G4cout << ", " << theQuarkContent[1]; 156 G4cout << ", " << theQuarkContent[2]; 225 G4cout << ", " << theQuarkContent[2]; 157 G4cout << ", " << theQuarkContent[3]; 226 G4cout << ", " << theQuarkContent[3]; 158 G4cout << ", " << theQuarkContent[4]; 227 G4cout << ", " << theQuarkContent[4]; 159 G4cout << ", " << theQuarkContent[5] << G4en 228 G4cout << ", " << theQuarkContent[5] << G4endl; 160 G4cout << " AntiQuark contents 229 G4cout << " AntiQuark contents : " << theAntiQuarkContent[0]; 161 G4cout << ", " << theAntiQuarkContent[1]; 230 G4cout << ", " << theAntiQuarkContent[1]; 162 G4cout << ", " << theAntiQuarkContent[2]; 231 G4cout << ", " << theAntiQuarkContent[2]; 163 G4cout << ", " << theAntiQuarkContent[3]; 232 G4cout << ", " << theAntiQuarkContent[3]; 164 G4cout << ", " << theAntiQuarkContent[4]; 233 G4cout << ", " << theAntiQuarkContent[4]; 165 G4cout << ", " << theAntiQuarkContent[5] << 234 G4cout << ", " << theAntiQuarkContent[5] << G4endl; 166 #endif 235 #endif 167 } 236 } >> 237 >> 238 >> 239 >> 240 >> 241 >> 242 >> 243 >> 244 >> 245 168 246