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 // G4IsotopeProperty class implementation 27 // 28 // Author: H.Kurashige, 5 October 1999 29 // -------------------------------------------------------------------- 30 31 #include "G4IsotopeProperty.hh" 32 33 #include "G4DecayTable.hh" 34 #include "G4PhysicalConstants.hh" 35 #include "G4SystemOfUnits.hh" 36 #include "G4ios.hh" 37 38 #include <iomanip> 39 40 G4IsotopeProperty::~G4IsotopeProperty() 41 { 42 delete fDecayTable; 43 fDecayTable = nullptr; 44 } 45 46 G4IsotopeProperty::G4IsotopeProperty(const G4IsotopeProperty& right) 47 : fAtomicNumber(right.fAtomicNumber), 48 fAtomicMass(right.fAtomicMass), 49 fISpin(right.fISpin), 50 fEnergy(right.fEnergy), 51 fLifeTime(right.fLifeTime), 52 53 fMagneticMoment(right.fMagneticMoment), 54 fIsomerLevel(right.fIsomerLevel), 55 fFloatLevelBase(right.fFloatLevelBase) 56 { 57 // decay table is not copied because G4DecayTable has no copy constructor 58 } 59 60 G4IsotopeProperty& G4IsotopeProperty::operator=(G4IsotopeProperty& right) 61 { 62 if (this != &right) { 63 fAtomicNumber = right.fAtomicNumber; 64 fAtomicMass = right.fAtomicMass; 65 fISpin = right.fISpin; 66 fMagneticMoment = right.fMagneticMoment; 67 fEnergy = right.fEnergy; 68 fLifeTime = right.fLifeTime; 69 fIsomerLevel = right.fIsomerLevel; 70 fFloatLevelBase = right.fFloatLevelBase; 71 // decay table is not copied because G4DecayTable has no copy constructor 72 fDecayTable = nullptr; 73 } 74 return *this; 75 } 76 77 G4bool G4IsotopeProperty::operator==(const G4IsotopeProperty& right) const 78 { 79 G4bool value = true; 80 value = value && (fAtomicNumber == right.fAtomicNumber); 81 value = value && (fAtomicMass == right.fAtomicMass); 82 value = value && (fISpin == right.fISpin); 83 value = value && (fMagneticMoment == right.fMagneticMoment); 84 value = value && (fEnergy == right.fEnergy); 85 value = value && (fLifeTime == right.fLifeTime); 86 value = value && (fIsomerLevel == right.fIsomerLevel); 87 value = value && (fFloatLevelBase == right.fFloatLevelBase); 88 return value; 89 } 90 91 G4bool G4IsotopeProperty::operator!=(const G4IsotopeProperty& right) const 92 { 93 return !(*this == right); 94 } 95 96 void G4IsotopeProperty::DumpInfo() const 97 { 98 #ifdef G4VERBOSE 99 G4cout << "AtomicNumber: " << fAtomicNumber << ", " 100 << "AtomicMass: " << fAtomicMass << G4endl; 101 if ((fISpin % 2) != 0) { 102 G4cout << "Spin: " << fISpin << "/2"; 103 } 104 else { 105 G4cout << "Spin: " << fISpin / 2; 106 } 107 G4cout << ", " 108 << "MagneticMoment: " << fMagneticMoment / MeV * tesla << "[MeV/T]" << G4endl; 109 G4cout << "Isomer Level: " << fIsomerLevel << ", Excited Energy: " << std::setprecision(1) 110 << fEnergy / keV; 111 if (fFloatLevelBase != G4Ions::G4FloatLevelBase::no_Float) { 112 G4cout << " +" << G4Ions::FloatLevelBaseChar(fFloatLevelBase); 113 } 114 G4cout << " [keV]" 115 << ", " << std::setprecision(6) << "Life Time: " << fLifeTime / ns << "[ns]" << G4endl; 116 if (fDecayTable != nullptr) { 117 fDecayTable->DumpInfo(); 118 } 119 else { 120 // G4cout << "Decay Table is not defined !" << G4endl; 121 } 122 #endif 123 } 124