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 // G4ElectronOccupancy << 27 // 26 // 28 // Class description: << 27 // $Id: G4ElectronOccupancy.hh 67971 2013-03-13 10:13:24Z gcosmo $ 29 // 28 // 30 // This class has information of occupation of << 29 // 31 // GetOccupancy(N) gives the number of electro << 30 // ------------------------------------------------------------ 32 // For example : Carbon atom should be << 31 // GEANT 4 class header file >> 32 // >> 33 // History: first implementation, based on object model of >> 34 // Hisaya Kurashige, 17 Aug 1999 >> 35 // ---------------------------------------------------------------- >> 36 // Class Description >> 37 // This class has information of occupation of electrons >> 38 // in atomic orbits >> 39 // - >> 40 // GetOccupancy(N) gives the number of electron >> 41 // in N-th orbit >> 42 // For example : Carbon atom should be 33 // GetOccupancy(0) --> 2 43 // GetOccupancy(0) --> 2 34 // GetOccupancy(1) --> 4 44 // GetOccupancy(1) --> 4 35 // GetOccupancy(2..7) --> 0 45 // GetOccupancy(2..7) --> 0 >> 46 // - >> 47 // GetTotalOccupancy() gives the total number of electrons 36 // 48 // 37 // GetTotalOccupancy() gives the total number << 49 // --------------------------------------------------------------- >> 50 38 51 39 // Author: Hisaya Kurashige, 17 Aug 1999 << 40 // ------------------------------------------- << 41 #ifndef G4ElectronOccupancy_hh << 42 #define G4ElectronOccupancy_hh 1 << 43 52 >> 53 #ifndef G4ElectronOccupancy_h >> 54 #define G4ElectronOccupancy_h 1 >> 55 >> 56 #include "globals.hh" 44 #include "G4Allocator.hh" 57 #include "G4Allocator.hh" 45 #include "G4ios.hh" 58 #include "G4ios.hh" 46 #include "globals.hh" << 47 59 48 #include "pwdefs.hh" 60 #include "pwdefs.hh" 49 61 50 class G4ElectronOccupancy << 62 class G4ElectronOccupancy 51 { 63 { 52 public: << 64 public: 53 enum << 65 enum { MaxSizeOfOrbit = 10}; 54 { << 66 55 MaxSizeOfOrbit = 20 << 67 public: // With Description 56 }; << 68 G4ElectronOccupancy( G4int sizeOrbit = MaxSizeOfOrbit ); 57 << 69 G4ElectronOccupancy( const G4ElectronOccupancy& right ); 58 G4ElectronOccupancy(G4int sizeOrbit = MaxS << 70 59 G4ElectronOccupancy(const G4ElectronOccupa << 71 public: 60 << 72 virtual ~G4ElectronOccupancy(); 61 virtual ~G4ElectronOccupancy(); << 73 62 << 74 // new/delete operators are oberloded to use G4Allocator 63 // new/delete operators are overloaded to << 75 inline void *operator new(size_t); 64 inline void* operator new(size_t); << 76 inline void operator delete(void *aElectronOccupancy); 65 inline void operator delete(void* aElectro << 77 66 << 78 67 G4ElectronOccupancy& operator=(const G4Ele << 79 //- operators 68 G4bool operator==(const G4ElectronOccupanc << 80 G4ElectronOccupancy & operator=(const G4ElectronOccupancy &right); 69 G4bool operator!=(const G4ElectronOccupanc << 81 G4int operator==(const G4ElectronOccupancy &right) const; 70 << 82 G4int operator!=(const G4ElectronOccupancy &right) const; 71 inline G4int GetTotalOccupancy() const; << 83 72 inline G4int GetOccupancy(G4int orbit) con << 84 public: // With Description 73 // The following methods returns << 85 // The following methods returns 74 // 0: if the orbit(atom) is vacant << 86 // 0: if the orbit(atom) is vacant 75 // >0: number of electrons in orbit << 87 // >0: number of electrons in orbit 76 << 88 G4int GetTotalOccupancy() const; 77 G4int AddElectron(G4int orbit, G4int numbe << 89 G4int GetOccupancy(G4int orbit) const; 78 G4int RemoveElectron(G4int orbit, G4int nu << 90 79 << 91 // 80 inline G4int GetSizeOfOrbit() const; << 92 G4int AddElectron(G4int orbit, G4int number = 1); 81 void DumpInfo() const; << 93 G4int RemoveElectron(G4int orbit, G4int number = 1); 82 << 94 83 private: << 95 G4int GetSizeOfOrbit() const; 84 G4int theSizeOfOrbit = 0; << 96 void DumpInfo() const; 85 G4int theTotalOccupancy = 0; << 97 86 G4int* theOccupancies = nullptr; << 98 private: >> 99 G4int theSizeOfOrbit; >> 100 G4int theTotalOccupancy; >> 101 G4int* theOccupancies; >> 102 87 }; 103 }; 88 104 89 extern G4PART_DLL G4Allocator<G4ElectronOccupa << 105 extern G4PART_DLL G4ThreadLocal G4Allocator<G4ElectronOccupancy> *aElectronOccupancyAllocator; 90 106 91 // ------------------------ 107 // ------------------------ 92 // Inline methods << 108 // Inlined operators 93 // ------------------------ 109 // ------------------------ 94 110 95 inline void* G4ElectronOccupancy::operator new << 111 inline void * G4ElectronOccupancy::operator new(size_t) 96 { 112 { 97 if (aElectronOccupancyAllocator() == nullptr << 113 if (!aElectronOccupancyAllocator) 98 aElectronOccupancyAllocator() = new G4Allo << 114 { >> 115 aElectronOccupancyAllocator = new G4Allocator<G4ElectronOccupancy>; 99 } 116 } 100 return (void*)aElectronOccupancyAllocator()- << 117 return (void *) aElectronOccupancyAllocator->MallocSingle(); 101 } 118 } 102 119 103 inline void G4ElectronOccupancy::operator dele << 120 inline void G4ElectronOccupancy::operator delete(void * aElectronOccupancy) 104 { 121 { 105 aElectronOccupancyAllocator()->FreeSingle((G << 122 aElectronOccupancyAllocator->FreeSingle((G4ElectronOccupancy *) aElectronOccupancy); 106 } 123 } 107 124 108 inline G4int G4ElectronOccupancy::GetSizeOfOrb << 125 inline >> 126 G4int G4ElectronOccupancy::GetSizeOfOrbit() const 109 { 127 { 110 return theSizeOfOrbit; << 128 return theSizeOfOrbit; 111 } 129 } 112 130 113 inline G4int G4ElectronOccupancy::GetTotalOccu << 131 inline >> 132 G4int G4ElectronOccupancy::GetTotalOccupancy() const 114 { 133 { 115 return theTotalOccupancy; << 134 return theTotalOccupancy; 116 } 135 } 117 136 118 inline G4int G4ElectronOccupancy::GetOccupancy << 137 inline >> 138 G4int G4ElectronOccupancy::GetOccupancy(G4int orbit) const 119 { 139 { 120 G4int value = 0; 140 G4int value = 0; 121 if ((orbit >= 0) && (orbit < theSizeOfOrbit) << 141 if ((orbit >=0)&&(orbit<theSizeOfOrbit)){ 122 value = theOccupancies[orbit]; 142 value = theOccupancies[orbit]; 123 } 143 } >> 144 return value; >> 145 } >> 146 >> 147 inline >> 148 G4int G4ElectronOccupancy::AddElectron(G4int orbit, G4int number) >> 149 { >> 150 G4int value =0; >> 151 if ((orbit >=0)&&(orbit<theSizeOfOrbit)){ >> 152 theOccupancies[orbit] += number; >> 153 theTotalOccupancy += number; >> 154 value = number; >> 155 } 124 return value; 156 return value; 125 } 157 } 126 158 >> 159 inline >> 160 G4int G4ElectronOccupancy::RemoveElectron(G4int orbit, G4int number) >> 161 { >> 162 G4int value =0; >> 163 if ((orbit >=0)&&(orbit<theSizeOfOrbit) ){ >> 164 if ( theOccupancies[orbit] < number ) number = theOccupancies[orbit]; >> 165 theOccupancies[orbit] -= number; >> 166 theTotalOccupancy -= number; >> 167 value = number; >> 168 } >> 169 return value; >> 170 } 127 #endif 171 #endif 128 172