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