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 // G4ParticleTable << 27 // 26 // 28 // Class description: << 29 // 27 // 30 // G4ParticleTable is the table of pointers to << 28 // 31 // It is a "singleton" (only one static object << 29 // ------------------------------------------------------------ 32 // Each G4ParticleDefinition pointer is stored << 30 // GEANT 4 class header file 33 // to itself. So, each G4ParticleDefinition ob << 31 // 34 // name. << 32 // History: first implementation, based on object model of 35 << 33 // 27 June 1996, H.Kurashige 36 // Authors: G.Cosmo, 2 December 1995 - Design, << 34 // ------------------------------------------------------------ 37 // H.Kurashige, 27 June 1996 - First << 35 // added fParticleMessenger 14 Nov., 97 H.Kurashige 38 // History: << 36 // added Create/DeleteMessenger 06 Jul., 98 H.Kurashige 39 // - 14 Nov 1997, H.Kurashige - Added messenge << 37 // modified FindIon 02 Aug., 98 H.Kurashige 40 // - 24 Sep 1998, H.Kurashige - Added dictiona << 38 // added dictionary for encoding 24 Sep., 98 H.Kurashige 41 // - 28 Oct 1999, H.Kurashige - Migration to S << 39 // added RemoveAllParticles() 8 Nov., 98 H.Kurashige 42 // - 15 Sep 2017, K.L.Genser - Added support f << 40 // -------------------------------- 43 // ------------------------------------------- << 41 // fixed some improper codings 08 Apr., 99 H.Kurashige 44 #ifndef G4ParticleTable_hh << 42 // modified FindIon/GetIon methods 17 AUg., 99 H.Kurashige 45 #define G4ParticleTable_hh 1 << 43 // implement new version for using STL map instaed of RW PtrHashedDictionary >> 44 // 28 ct., 99 H.Kurashige >> 45 // modified implementation of Remove 21 Mar.,08 H.Kurashige >> 46 // remove G4ShortLivedTable 25 July, 13 H.Kurashige >> 47 // added support for MuonicAtom September, 17 K.L.Genser >> 48 // 46 49 47 #include "G4ParticleDefinition.hh" << 50 #ifndef G4ParticleTable_h 48 #include "G4ParticleTableIterator.hh" << 51 #define G4ParticleTable_h 1 49 #include "G4Threading.hh" << 50 #include "G4ios.hh" << 51 #include "globals.hh" << 52 52 53 #include <map> 53 #include <map> 54 54 >> 55 #include "G4ios.hh" >> 56 #include "globals.hh" >> 57 #include "G4Threading.hh" >> 58 #include "G4ParticleDefinition.hh" >> 59 #include "G4ParticleTableIterator.hh" >> 60 55 class G4UImessenger; 61 class G4UImessenger; 56 class G4ParticleMessenger; 62 class G4ParticleMessenger; 57 class G4IonTable; 63 class G4IonTable; 58 64 59 class G4ParticleTable 65 class G4ParticleTable 60 { 66 { 61 public: << 67 // Class Description 62 using G4PTblDictionary = G4ParticleTableIt << 68 // G4ParticleTable is the table of pointer to G4ParticleDefinition 63 using G4PTblDicIterator = G4ParticleTableI << 69 // G4ParticleTable is a "singleton" (only one and staic object) 64 using G4PTblEncodingDictionary = G4Particl << 70 // In G4ParticleTable, each G4ParticleDefinition pointer is stored 65 using G4PTblEncodingDicIterator = G4Partic << 71 // with its name as a key to itself. So, each G4ParticleDefinition 66 << 72 // object must have unique name for itself. 67 virtual ~G4ParticleTable(); << 73 // 68 << 74 69 // Copy constructor and assignment operato << 75 public: 70 G4ParticleTable(const G4ParticleTable&) = << 76 71 G4ParticleTable& operator=(const G4Particl << 77 typedef G4ParticleTableIterator<G4String, G4ParticleDefinition*>::Map G4PTblDictionary; 72 << 78 typedef G4ParticleTableIterator<G4String, G4ParticleDefinition*> G4PTblDicIterator; 73 // This method is similar to the construct << 79 typedef G4ParticleTableIterator<G4int, G4ParticleDefinition*>::Map G4PTblEncodingDictionary; 74 // thread to achieve the partial effect as << 80 typedef G4ParticleTableIterator<G4int, G4ParticleDefinition*> G4PTblEncodingDicIterator; 75 void WorkerG4ParticleTable(); << 81 76 << 82 protected: 77 // This method is similar to the destructo << 83 // default constructor 78 // thread to achieve the partial effect as << 84 G4ParticleTable(); 79 void DestroyWorkerG4ParticleTable(); << 85 // Copy constructor and assignment operator 80 << 86 G4ParticleTable(const G4ParticleTable &right); 81 // Return the pointer to the G4ParticleTab << 87 G4ParticleTable & operator=(const G4ParticleTable &); 82 // G4ParticleTable is a "singleton" and ca << 88 83 // function. At the first time of calling << 89 public: 84 // G4ParticleTable object is instantiated << 90 85 static G4ParticleTable* GetParticleTable() << 91 void SlaveG4ParticleTable(); 86 << 92 void WorkerG4ParticleTable(); 87 // Returns TRUE if the ParticleTable conta << 93 // This method is similar to the constructor. It is used by each worker 88 inline G4bool contains(const G4ParticleDef << 94 // thread to achieve the partial effect as that of the master thread. 89 G4bool contains(const G4String& particle_n << 95 90 << 96 virtual ~G4ParticleTable(); 91 // Returns the number of particles in the << 97 void DestroyWorkerG4ParticleTable(); 92 G4int entries() const; << 98 // This method is similar to the destructor. It is used by each worker 93 G4int size() const; << 99 // thread to achieve the partial effect as that of the master thread. 94 << 100 95 // Returns a pointer to the i-th particle << 101 public: // With Description 96 // 0 <= index < entries() << 102 static G4ParticleTable* GetParticleTable(); 97 G4ParticleDefinition* GetParticle(G4int in << 103 // return the pointer to G4ParticleTable object 98 << 104 // G4ParticleTable is a "singleton" and can get its pointer by this function 99 // Returns the name of i-th particle in th << 105 // At the first time of calling this function, the G4ParticleTable object 100 const G4String& GetParticleName(G4int inde << 106 // is instantiated 101 << 107 102 // Returns a pointer to the particle (0 if << 108 G4bool contains(const G4ParticleDefinition *particle) const; 103 G4ParticleDefinition* FindParticle(G4int P << 109 G4bool contains(const G4String &particle_name) const; 104 G4ParticleDefinition* FindParticle(const G << 110 // returns TRUE if the ParticleTable contains 105 G4ParticleDefinition* FindParticle(const G << 111 106 << 112 G4int entries() const; 107 // Returns a pointer to its anti-particle << 113 G4int size() const; 108 inline G4ParticleDefinition* FindAntiParti << 114 // returns the number of Particles in the ParticleTable 109 inline G4ParticleDefinition* FindAntiParti << 115 110 inline G4ParticleDefinition* FindAntiParti << 116 G4ParticleDefinition* GetParticle(G4int index) const; 111 << 117 // returns a pointer to i-th particles in the ParticleTable 112 // Returns the pointer to the Iterator << 118 // 0<= index < entries() 113 G4PTblDicIterator* GetIterator() const; << 119 114 << 120 const G4String& GetParticleName(G4int index) const; 115 // Dumps information of particles specifie << 121 // returns name of i-th particles in the ParticleTable 116 void DumpTable(const G4String& particle_na << 122 117 << 123 G4ParticleDefinition* FindParticle(G4int PDGEncoding ); 118 // Returns the pointer to the G4IonTable o << 124 G4ParticleDefinition* FindParticle(const G4String &particle_name); 119 G4IonTable* GetIonTable() const; << 125 G4ParticleDefinition* FindParticle(const G4ParticleDefinition *particle); 120 << 126 // returns a pointer to the particle (0 if not contained) 121 // Inserts the particle into ParticleTable << 127 122 // Returned value is the same as particle << 128 G4ParticleDefinition* FindAntiParticle(G4int PDGEncoding ); 123 // or the pointer to another G4ParticleD << 129 G4ParticleDefinition* FindAntiParticle(const G4String &particle_name); 124 // which has same particle name << 130 G4ParticleDefinition* FindAntiParticle(const G4ParticleDefinition *particle); 125 // or nullptr if failing to insert by ot << 131 // returns a pointer to its anti-particle (0 if not contained) 126 G4ParticleDefinition* Insert(G4ParticleDef << 132 127 << 133 G4PTblDicIterator* GetIterator() const; 128 // Removes the particle from the table (no << 134 // return the pointer of Iterator (RW compatible) 129 G4ParticleDefinition* Remove(G4ParticleDef << 135 130 << 136 void DumpTable(const G4String &particle_name = "ALL"); 131 // Removes all particles from G4ParticleTa << 137 // dump information of particles specified by name 132 void RemoveAllParticles(); << 138 133 << 139 public: //With Description 134 // Removes and deletes all particles from << 140 135 void DeleteAllParticles(); << 141 G4IonTable* GetIonTable() const; 136 << 142 // return the pointer to G4IonTable object 137 // Creates messenger << 143 138 G4UImessenger* CreateMessenger(); << 144 139 << 145 public: // With Description 140 void SelectParticle(const G4String& name); << 146 G4ParticleDefinition* Insert(G4ParticleDefinition *particle); 141 << 147 // insert the particle into ParticleTable 142 inline const G4ParticleDefinition* GetSele << 148 // return value is same as particle if successfully inserted 143 << 149 // or pointer to another G4ParticleDefinition object 144 inline void SetVerboseLevel(G4int value); << 150 // which has same name of particle 145 inline G4int GetVerboseLevel() const; << 151 // or 0 if fail to insert by another reason 146 << 152 147 inline void SetReadiness(G4bool val = true << 153 G4ParticleDefinition* Remove(G4ParticleDefinition *particle); 148 inline G4bool GetReadiness() const; << 154 // Remove the particle from the table (not delete) 149 << 155 150 inline G4ParticleDefinition* GetGenericIon << 156 void RemoveAllParticles(); 151 inline void SetGenericIon(G4ParticleDefini << 157 // remove all particles from G4ParticleTable 152 << 158 153 inline G4ParticleDefinition* GetGenericMuo << 159 void DeleteAllParticles(); 154 inline void SetGenericMuonicAtom(G4Particl << 160 // remove and delete all particles from G4ParticleTable 155 << 161 156 // Public data --------------------------- << 162 public: 157 << 163 G4UImessenger* CreateMessenger(); 158 // These fields should be thread local or << 164 void DeleteMessenger(); 159 // class, we can change any member field a << 165 // create/delete messenger for the particle table 160 // because there is only one instance. The << 166 // these methods are supposed to be invoked by G4RunManager only 161 // "G4ThreadLocal" << 167 162 G4ParticleMessenger* fParticleMessenger = << 163 static G4ThreadLocal G4PTblDictionary* fDi << 164 static G4ThreadLocal G4PTblDicIterator* fI << 165 static G4ThreadLocal G4PTblEncodingDiction << 166 << 167 // Particle table is being shared << 168 static G4ParticleTable* fgParticleTable; << 169 << 170 // This field should be thread private. Ho << 171 // of the ion table pointer. So we change << 172 // G4IonTable to be thread local << 173 G4IonTable* fIonTable = nullptr; << 174 << 175 // These shadow pointers are used by each << 176 // from the master thread << 177 // << 178 static G4PTblDictionary* fDictionaryShadow << 179 static G4PTblDicIterator* fIteratorShadow; << 180 static G4PTblEncodingDictionary* fEncoding << 181 << 182 #ifdef G4MULTITHREADED << 183 // Shared instance of a mutex << 184 static G4GLOB_DLL G4Mutex& particleTableMu << 185 static G4GLOB_DLL G4int& lockCount(); << 186 #endif << 187 << 188 protected: 168 protected: 189 const G4PTblDictionary* GetDictionary() co << 190 << 191 // Returns key value of the particle (i.e. << 192 inline const G4String& GetKey(const G4Part << 193 169 194 // Returns the pointer to EncodingDictiona << 170 const G4PTblDictionary* GetDictionary() const; 195 const G4PTblEncodingDictionary* GetEncodin << 196 171 197 private: << 172 const G4String& GetKey(const G4ParticleDefinition *particle) const; 198 // Provate default constructor << 173 // return key value of the particle (i.e. particle name) 199 G4ParticleTable(); << 200 174 201 void CheckReadiness() const; << 175 const G4PTblEncodingDictionary* GetEncodingDictionary() const; >> 176 // return the pointer to EncodingDictionary 202 177 203 // Private data -------------------------- << 178 private: >> 179 G4int verboseLevel; >> 180 // controle flag for output message >> 181 // 0: Silent >> 182 // 1: Warning message >> 183 // 2: More >> 184 >> 185 public: >> 186 void SetVerboseLevel(G4int value); >> 187 G4int GetVerboseLevel() const; >> 188 >> 189 static G4ThreadLocal G4ParticleMessenger* fParticleMessenger; >> 190 static G4ThreadLocal G4PTblDictionary* fDictionary; >> 191 static G4ThreadLocal G4PTblDicIterator* fIterator; >> 192 static G4ThreadLocal G4PTblEncodingDictionary* fEncodingDictionary; >> 193 // These fields should be thread local or thread private. For a singleton >> 194 // class, we can change any member field as static without any problem >> 195 // because there is only one instance. Then we are allowed to add >> 196 // "G4ThreadLocal". >> 197 >> 198 //01.25.2009 Xin Dong: Phase II change for Geant4 multi-threading. >> 199 //Phase I changes this member to be thread local >> 200 //,while each thread holds its own copy of particles. >> 201 //Phase II changes this member back in order to share particles. >> 202 static G4ParticleTable* fgParticleTable; >> 203 >> 204 static G4IonTable* fIonTable; >> 205 // This field should be thread private. However, we have to keep one copy >> 206 // of the ion table pointer. So we change all important fields of G4IonTable >> 207 // to the thread local variable. >> 208 >> 209 // These shadow pointers are used by each worker thread to copy the content >> 210 // from the master thread. >> 211 >> 212 static G4ParticleMessenger* fParticleMessengerShadow; >> 213 static G4PTblDictionary* fDictionaryShadow; >> 214 static G4PTblDicIterator* fIteratorShadow; >> 215 static G4PTblEncodingDictionary* fEncodingDictionaryShadow; >> 216 >> 217 private: >> 218 const G4String noName; >> 219 >> 220 G4bool readyToUse; >> 221 G4ParticleDefinition* genericIon; >> 222 G4ParticleDefinition* genericMuonicAtom; >> 223 >> 224 public: >> 225 void SetReadiness(G4bool val=true); >> 226 G4bool GetReadiness() const; >> 227 G4ParticleDefinition* GetGenericIon() const; >> 228 void SetGenericIon(G4ParticleDefinition*); >> 229 G4ParticleDefinition* GetGenericMuonicAtom() const; >> 230 void SetGenericMuonicAtom(G4ParticleDefinition*); >> 231 private: >> 232 void CheckReadiness() const; 204 233 205 G4ParticleDefinition* genericIon = nullptr << 206 G4ParticleDefinition* genericMuonicAtom = << 207 const G4ParticleDefinition* selectedPartic << 208 234 209 const G4String noName = " "; << 235 #ifdef G4MULTITHREADED 210 G4String selectedName = "undefined"; << 236 public: 211 << 237 //Andrea Dotti January 16. Shared instance of a mutex 212 // Control flag for output message << 238 static G4GLOB_DLL G4Mutex& particleTableMutex(); 213 // 0: Silent << 239 static G4GLOB_DLL G4int& lockCount(); 214 // 1: Warning message << 240 #endif 215 // 2: More << 216 G4int verboseLevel = 1; << 217 << 218 G4bool readyToUse = false; << 219 }; 241 }; 220 << 221 #include "G4ParticleTable.icc" 242 #include "G4ParticleTable.icc" 222 243 223 #endif 244 #endif 224 245