Geant4 Cross Reference |
>> 1 // This code implementation is the intellectual property of >> 2 // the GEANT4 collaboration. 1 // 3 // 2 // ******************************************* << 4 // By copying, distributing or modifying the Program (or any work 3 // * License and Disclaimer << 5 // based on the Program) you indicate your acceptance of this statement, 4 // * << 6 // and all its terms. 5 // * The Geant4 software is copyright of th << 6 // * the Geant4 Collaboration. It is provided << 7 // * conditions of the Geant4 Software License << 8 // * LICENSE and available at http://cern.ch/ << 9 // * include a list of copyright holders. << 10 // * << 11 // * Neither the authors of this software syst << 12 // * institutes,nor the agencies providing fin << 13 // * work make any representation or warran << 14 // * regarding this software system or assum << 15 // * use. Please see the license in the file << 16 // * for the full disclaimer and the limitatio << 17 // * << 18 // * This code implementation is the result << 19 // * technical work of the GEANT4 collaboratio << 20 // * By using, copying, modifying or distri << 21 // * any work based on the software) you ag << 22 // * use in resulting scientific publicati << 23 // * acceptance of all terms of the Geant4 Sof << 24 // ******************************************* << 25 // 7 // 26 // G4ElectronOccupancy class implementation << 8 // $Id: G4ElectronOccupancy.cc,v 1.4.6.1.2.1 1999/12/08 17:34:08 gunter Exp $ >> 9 // GEANT4 tag $Name: geant4-01-01 $ 27 // 10 // 28 // Author: Hisaya Kurashige, 17 Aug 1999 << 11 // >> 12 // ---------------------------------------------------------------------- >> 13 // GEANT 4 class implementation file >> 14 // >> 15 // For information related to this code contact: >> 16 // CERN, CN Division, ASD group >> 17 // History: first implementation, based on object model of >> 18 // Hisaya Kurashige, 17 Aug 1999 >> 19 // ---------------------------------------------------------------- >> 20 // This class has information of occupation of electrons >> 21 // in atomic orbits 29 // ------------------------------------------- 22 // --------------------------------------------------------------- 30 23 31 #include "G4ElectronOccupancy.hh" << 32 24 33 #include <sstream> << 25 #include "G4ElectronOccupancy.hh" 34 26 35 G4Allocator<G4ElectronOccupancy>*& aElectronOc << 27 G4Allocator<G4ElectronOccupancy> aElectronOccupancyAllocator; 36 { << 37 G4ThreadLocalStatic G4Allocator<G4ElectronOc << 38 return _instance; << 39 } << 40 28 41 G4ElectronOccupancy::G4ElectronOccupancy(G4int << 29 G4ElectronOccupancy::G4ElectronOccupancy(G4int sizeOrbit ): >> 30 theSizeOfOrbit(sizeOrbit) 42 { 31 { 43 // check size 32 // check size 44 if ((theSizeOfOrbit < 1) || (theSizeOfOrbit << 33 if ( (theSizeOfOrbit <1 ) || (theSizeOfOrbit > MaxSizeOfOrbit) ) { 45 theSizeOfOrbit = MaxSizeOfOrbit; 34 theSizeOfOrbit = MaxSizeOfOrbit; 46 } 35 } 47 36 48 // allocate and clear the array of theOccupa << 37 // allocate and clear the array of theOccupancies 49 theOccupancies = new G4int[theSizeOfOrbit]; 38 theOccupancies = new G4int[theSizeOfOrbit]; 50 for (G4int index = 0; index < theSizeOfOrbit << 39 G4int index =0; 51 theOccupancies[index] = 0; << 40 for (index = 0; index < theSizeOfOrbit; index++) { >> 41 theOccupancies[index] =0; 52 } 42 } 53 43 54 theTotalOccupancy = 0; << 44 theTotalOccupancy =0; 55 } 45 } 56 46 57 G4ElectronOccupancy::~G4ElectronOccupancy() 47 G4ElectronOccupancy::~G4ElectronOccupancy() 58 { 48 { 59 theSizeOfOrbit = -1; << 49 theSizeOfOrbit = -1; 60 50 61 delete[] theOccupancies; << 51 delete [] theOccupancies; 62 theOccupancies = nullptr; << 52 theOccupancies =0; 63 theTotalOccupancy = 0; << 53 theTotalOccupancy =0; 64 } 54 } 65 55 >> 56 66 G4ElectronOccupancy::G4ElectronOccupancy(const 57 G4ElectronOccupancy::G4ElectronOccupancy(const G4ElectronOccupancy& right) 67 { 58 { 68 theSizeOfOrbit = right.theSizeOfOrbit; 59 theSizeOfOrbit = right.theSizeOfOrbit; 69 60 70 // allocate and clear the array of theOccupa << 61 // allocate and clear the array of theOccupancies 71 theOccupancies = new G4int[theSizeOfOrbit]; 62 theOccupancies = new G4int[theSizeOfOrbit]; 72 for (G4int index = 0; index < theSizeOfOrbit << 63 G4int index =0; >> 64 for (index = 0; index < theSizeOfOrbit; index++) { 73 theOccupancies[index] = right.theOccupanci 65 theOccupancies[index] = right.theOccupancies[index]; 74 } 66 } 75 67 76 theTotalOccupancy = right.theTotalOccupancy; 68 theTotalOccupancy = right.theTotalOccupancy; 77 } 69 } 78 70 79 G4ElectronOccupancy& G4ElectronOccupancy::oper 71 G4ElectronOccupancy& G4ElectronOccupancy::operator=(const G4ElectronOccupancy& right) 80 { 72 { 81 if (this != &right) { << 73 if ( this != &right) { 82 theSizeOfOrbit = right.theSizeOfOrbit; 74 theSizeOfOrbit = right.theSizeOfOrbit; 83 << 75 84 // allocate and clear the array of theOccu << 76 // allocate and clear the array of theOccupancies 85 delete[] theOccupancies; << 77 if ( theOccupancies != 0 ) delete [] theOccupancies; 86 theOccupancies = new G4int[theSizeOfOrbit] 78 theOccupancies = new G4int[theSizeOfOrbit]; 87 for (G4int index = 0; index < theSizeOfOrb << 79 G4int index =0; >> 80 for (index = 0; index < theSizeOfOrbit; index++) { 88 theOccupancies[index] = right.theOccupan 81 theOccupancies[index] = right.theOccupancies[index]; 89 } 82 } 90 << 83 91 theTotalOccupancy = right.theTotalOccupanc 84 theTotalOccupancy = right.theTotalOccupancy; 92 } 85 } 93 return *this; 86 return *this; 94 } 87 } 95 88 96 G4bool G4ElectronOccupancy::operator==(const G << 89 G4int G4ElectronOccupancy::operator==(const G4ElectronOccupancy& right) const 97 { 90 { >> 91 G4int index; 98 G4bool value = true; 92 G4bool value = true; 99 for (G4int index = 0; index < MaxSizeOfOrbit << 93 for (index = 0; index < MaxSizeOfOrbit; index++) { 100 if ((index < theSizeOfOrbit) && (index < r << 94 if ( (index < theSizeOfOrbit ) && ( index < right.theSizeOfOrbit) ) { 101 value = value && (theOccupancies[index] << 95 value = value && 102 } << 96 (theOccupancies[index] == right.theOccupancies[index]) ; 103 else if ((index < theSizeOfOrbit) && (inde << 97 } else if ((index < theSizeOfOrbit ) && ( index >= right.theSizeOfOrbit)) { 104 value = value && (!static_cast<G4bool>(t << 98 value = value && (theOccupancies[index] == 0); 105 } << 99 } else if ((index >= theSizeOfOrbit ) && ( index <right.theSizeOfOrbit)) { 106 else if ((index >= theSizeOfOrbit) && (ind << 100 value = value && (right.theOccupancies[index] == 0); 107 value = value && (!static_cast<G4bool>(r << 108 } 101 } 109 } 102 } 110 return value; 103 return value; 111 } 104 } 112 105 113 G4bool G4ElectronOccupancy::operator!=(const G << 106 G4int G4ElectronOccupancy::operator!=(const G4ElectronOccupancy& right) const 114 { 107 { 115 return !(*this == right); 108 return !(*this == right); 116 } 109 } 117 110 >> 111 118 void G4ElectronOccupancy::DumpInfo() const 112 void G4ElectronOccupancy::DumpInfo() const 119 { 113 { 120 G4cout << " -- Electron Occupancy -- " << G 114 G4cout << " -- Electron Occupancy -- " << G4endl; 121 for (G4int index = 0; index < theSizeOfOrbit << 115 G4int index; 122 G4cout << " " << index << "-th orbit << 116 for (index = 0; index < theSizeOfOrbit; index++) { >> 117 G4cout << " " << index << "-th orbit " >> 118 << theOccupancies[index] << G4endl; 123 } 119 } 124 } << 125 << 126 G4int G4ElectronOccupancy::AddElectron(G4int o << 127 { << 128 G4int value = 0; << 129 if (orbit >= theSizeOfOrbit) { << 130 std::ostringstream smsg; << 131 smsg << "Orbit (" << orbit << ") exceeds t << 132 G4String msg = smsg.str(); << 133 G4Exception("G4ElectronOccupancy::AddElect << 134 } << 135 else if (orbit >= 0) { << 136 theOccupancies[orbit] += number; << 137 theTotalOccupancy += number; << 138 value = number; << 139 } << 140 return value; << 141 } << 142 << 143 G4int G4ElectronOccupancy::RemoveElectron(G4in << 144 { << 145 G4int value = 0; << 146 if (orbit >= theSizeOfOrbit) { << 147 std::ostringstream smsg; << 148 smsg << "Orbit (" << orbit << ") exceeds t << 149 G4String msg = smsg.str(); << 150 G4Exception("G4ElectronOccupancy::RemoveEl << 151 } << 152 else if (orbit >= 0) { << 153 if (theOccupancies[orbit] < number) number << 154 theOccupancies[orbit] -= number; << 155 theTotalOccupancy -= number; << 156 value = number; << 157 } << 158 return value; << 159 } 120 } 160 121