| 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 // G4VUserPhysicsList <<
27 // 26 //
28 // Class description: << 27 // $Id: G4VUserPhysicsList.hh 73123 2013-08-19 07:52:00Z gcosmo $
29 // 28 //
30 // This class is an abstract class for constru << 29 //
31 // User must implement the following two pure << 30 // ------------------------------------------------------------
32 // class derived from this class: << 31 // GEANT 4 class header file
33 // - G4VUserPhysicsList::ConstructParticle() << 32 // Class Description:
34 // Construct particles << 33 // This class is an abstract class for
35 // - G4VUserPhysicsList::ConstructProcess() << 34 // constructing particles and processes.
36 // Construct procesess and register them t << 35 // User must implement following three virtual methods
37 << 36 // in his/her own concrete class derived from this class.
38 // Original author: H.Kurashige (Kobe Universi << 37 // G4VUserPhysicsList::ConstructParticle()
39 // ------------------------------------------- << 38 // Construct particles
40 #ifndef G4VUserPhysicsList_hh << 39 // G4VUserPhysicsList::ConstructProcess()
41 #define G4VUserPhysicsList_hh 1 << 40 // Construct procesess and register them to particles
>> 41 //
>> 42 // -------------------------------------------
>> 43 // History
>> 44 // first version 09 Jan. 1998 by H.Kurashige
>> 45 // modified 24 Jan. 1998 by H.Kurashige
>> 46 // rename DumpCutValues/DumpCutValuesTable
>> 47 // change SetCuts method
>> 48 // add SetCutsWithDefault method
>> 49 // modified 06 June 1998 by H.Kurashige
>> 50 // add AddProcessManager
>> 51 // add BuildPhysicsTable
>> 52 // modified 29 June 1998 by H.Kurashige
>> 53 // add AddProcessManager
>> 54 // modified 05 Dec. 1998 by H.Kurashige
>> 55 // add ConstructAllParticles()
>> 56 // modified 14, Apr 1999 by H.Kurashige
>> 57 // change BuildPhysicsTable as public
>> 58 // removed ConstructAllParticles() and related methods
>> 59 // changed SetCuts method argument
>> 60 // modified 08, Nov 2000 by H.Kurashige
>> 61 // added Retrieve/StorePhysicsTable and related methods
>> 62 // modified 08, Mar 2001 by H.Kurashige
>> 63 // added binary mode for Retrieve/StorePhysicsTable
>> 64 // added RetrieveCutValues and related
>> 65 // added Set/ResetStoredInAscii() to switch on ascii mode
>> 66 // for Retrieve/StorePhysicsTable
>> 67 // modified for CUTS per REGION 10, Oct 2002 by H.Kurashige
>> 68 // removed following methods
>> 69 // void ReCalcCutValue()
>> 70 // void SetCutValueForOthers()
>> 71 // void SetCutValueForOtherThan()
>> 72 // void ReCalcCutValueForOthers()
>> 73 // virtual G4bool StoreMaterialInfo()
>> 74 // virtual G4bool StoreCutValues()
>> 75 // virtual G4bool RetrieveCutValues()
>> 76 // virtual G4bool CheckForRetrievePhysicsTable()
>> 77 // virtual G4bool CheckMaterialInfo()
>> 78 // added void BuildPhysicsTable()
>> 79 // Added PhysicsListHelper 29 Apr. 2011 H.Kurashige
>> 80 // Added default impelmentation of SetCuts 10 June 2011 H.Kurashige
>> 81 // SetCuts is not 'pure virtual' any more
>> 82 // Trasnformations for multi-threading 26 Mar. 2013 A. Dotti
>> 83 // ------------------------------------------------------------
>> 84 #ifndef G4VUserPhysicsList_h
>> 85 #define G4VUserPhysicsList_h 1
>> 86
>> 87 #include "globals.hh"
>> 88 #include "tls.hh"
>> 89 #include "rundefs.hh"
>> 90 #include "G4ios.hh"
42 91
43 #include "G4ParticleDefinition.hh" <<
44 #include "G4ParticleTable.hh" 92 #include "G4ParticleTable.hh"
45 #include "G4PhysicsModelCatalog.hh" << 93 #include "G4ParticleDefinition.hh"
46 #include "G4ProductionCutsTable.hh" 94 #include "G4ProductionCutsTable.hh"
47 #include "G4Threading.hh" <<
48 #include "G4VUPLSplitter.hh" 95 #include "G4VUPLSplitter.hh"
49 #include "G4ios.hh" <<
50 #include "globals.hh" <<
51 96
52 #include "rundefs.hh" << 97 #include "G4Threading.hh"
53 #include "tls.hh" <<
54 98
55 class G4UserPhysicsListMessenger; 99 class G4UserPhysicsListMessenger;
56 class G4PhysicsListHelper; 100 class G4PhysicsListHelper;
57 class G4VProcess; 101 class G4VProcess;
58 102
59 // Encapsulate the fields of class G4VUserPhys <<
60 class G4VUPLData 103 class G4VUPLData
61 { 104 {
62 public: << 105 //Encapsulate the fields of class G4VUserPhysicsList
>> 106 //that are per-thread.
>> 107 public:
63 void initialize(); 108 void initialize();
64 << 109 G4ParticleTable::G4PTblDicIterator* _theParticleIterator;
65 G4ParticleTable::G4PTblDicIterator* _thePa << 110 G4UserPhysicsListMessenger* _theMessenger;
66 G4UserPhysicsListMessenger* _theMessenger << 111 G4PhysicsListHelper* _thePLHelper;
67 G4PhysicsListHelper* _thePLHelper = nullpt << 112 G4bool _fIsPhysicsTableBuilt;
68 G4bool _fIsPhysicsTableBuilt = false; << 113 G4int _fDisplayThreshold;
69 G4int _fDisplayThreshold = 0; <<
70 }; 114 };
71 <<
72 // The type G4VUPLManager is introduced to enc 115 // The type G4VUPLManager is introduced to encapsulate the methods used by
73 // both the master thread and worker threads t 116 // both the master thread and worker threads to allocate memory space for
74 // the fields encapsulated by the class G4VUPL 117 // the fields encapsulated by the class G4VUPLData. When each thread
75 // changes the value for these fields, it refe 118 // changes the value for these fields, it refers to them using a macro
76 // definition defined below. For every G4VUser 119 // definition defined below. For every G4VUserPhysicsList instance,
77 // there is a corresponding G4VUPLData instanc 120 // there is a corresponding G4VUPLData instance. All G4VUPLData instances
78 // are organized by the class G4VUPLManager as 121 // are organized by the class G4VUPLManager as an array.
79 // The field "int g4vuplInstanceID" is added t 122 // The field "int g4vuplInstanceID" is added to the class G4VUserPhysicsList.
80 // The value of this field in each G4VUserPhys 123 // The value of this field in each G4VUserPhysicsList instance is the
81 // subscript of the corresponding G44VUPLData 124 // subscript of the corresponding G44VUPLData instance.
82 // In order to use the class G44VUPLManager, w 125 // In order to use the class G44VUPLManager, we add a static member in the class
83 // G4VUserPhysicsList as follows: "static G4VU 126 // G4VUserPhysicsList as follows: "static G4VUPLManager subInstanceManager".
84 // Both the master thread and worker threads c 127 // Both the master thread and worker threads change the length of the array
85 // for G44VUPLData instances mutually along wi 128 // for G44VUPLData instances mutually along with G4VUserPhysicsList
86 // instances are created. For each worker thre 129 // instances are created. For each worker thread, it dynamically creates ions.
87 // Consider any thread A, if there is any othe 130 // Consider any thread A, if there is any other thread which creates an ion.
88 // This ion is shared by the thread A. So the 131 // This ion is shared by the thread A. So the thread A leaves an empty space
89 // in the array of G4PDefData instances for th 132 // in the array of G4PDefData instances for the ion.
90 // 133 //
91 // Important Note: you may wonder why we are i 134 // Important Note: you may wonder why we are introducing this mechanism
92 // since there is only one PL 135 // since there is only one PL for each application.
93 // This is true, in the sense 136 // This is true, in the sense that only one PL is allowed
94 // to be associated to a G4Run << 137 // to be associated to a G4RunManager, however user can instantiate
95 // instantiate as many PLs are << 138 // as many PLs are needed and at run-time select one of the PLs to be used
96 // of the PLs to be used we th << 139 // we thus need this mechanism to guarantee that the system works without
97 // guarantee that the system w << 140 // problems in case of this (unusual) case. This may be reviewed in the future
98 // this (unusual) case. This m << 141 typedef G4VUPLSplitter<G4VUPLData> G4VUPLManager;
>> 142 typedef G4VUPLManager G4VUserPhysicsListSubInstanceManager;
>> 143
>> 144 // This macros change the references to fields that are now encapsulated
>> 145 // in the class G4VUPLData.
99 // 146 //
100 using G4VUPLManager = G4VUPLSplitter<G4VUPLDat << 147 // Note1: the use of this-> this is needed to avoid compilation errors
101 using G4VUserPhysicsListSubInstanceManager = G << 148 // when using templated class with T=G4VUserPhysicsList. Don't know why.
>> 149 // Note2: the name of the first #define is different, because otherwise
>> 150 // we need to change its use in all classes that inherits from
>> 151 // this base class (all examples). However one should note comment
>> 152 // on JIRA task: http://jira-geant4.kek.jp/browse/DEV-27
>> 153 #define theParticleIterator ((this->subInstanceManager.offset[this->g4vuplInstanceID])._theParticleIterator)
>> 154 #define G4MT_theMessenger ((this->subInstanceManager.offset[this->g4vuplInstanceID])._theMessenger)
>> 155 #define G4MT_thePLHelper ((this->subInstanceManager.offset[this->g4vuplInstanceID])._thePLHelper)
>> 156 #define fIsPhysicsTableBuilt ((this->subInstanceManager.offset[this->g4vuplInstanceID])._fIsPhysicsTableBuilt)
>> 157 #define fDisplayThreshold ((this->subInstanceManager.offset[this->g4vuplInstanceID])._fDisplayThreshold)
102 158
103 class G4VUserPhysicsList 159 class G4VUserPhysicsList
104 { 160 {
105 public: << 161 public:
106 G4VUserPhysicsList(); 162 G4VUserPhysicsList();
107 virtual ~G4VUserPhysicsList(); 163 virtual ~G4VUserPhysicsList();
>> 164
>> 165 // copy constructor and assignment operator
108 G4VUserPhysicsList(const G4VUserPhysicsLis 166 G4VUserPhysicsList(const G4VUserPhysicsList&);
109 G4VUserPhysicsList& operator=(const G4VUse << 167 G4VUserPhysicsList & operator=(const G4VUserPhysicsList&);
110 168
111 // Each particle type will be instantiated << 169 public: // with description
112 // This method is invoked by the RunManger << 170 // Each particle type will be instantiated
113 virtual void ConstructParticle() = 0; << 171 // This method is invoked by the RunManger
114 << 172 virtual void ConstructParticle() = 0;
115 // By calling the "Construct" method, << 173
116 // process manager and processes are creat << 174 // By calling the "Construct" method,
117 void Construct(); << 175 // process manager and processes are created.
118 << 176 void Construct();
119 // Each physics process will be instantiat << 177
120 // registered to the process manager of ea << 178 // Each physics process will be instantiated and
121 // Invoked in the Construct() method. << 179 // registered to the process manager of each particle type
122 virtual void ConstructProcess() = 0; << 180 // This method is invoked in Construct method
123 << 181 virtual void ConstructProcess() = 0;
124 // Sets a cut value for all particle types << 182
125 virtual void SetCuts(); << 183 protected: // with description
126 << 184 // User must invoke this method in his ConstructProcess()
127 // Set/get the default cut value. Calling << 185 // implementation in order to insures particle transportation.
128 // re-calcuration of cut values and physic << 186 void AddTransportation();
129 // next event loop. << 187
130 void SetDefaultCutValue(G4double newCutVal << 188 //Register a process to the particle type
131 G4double GetDefaultCutValue() const; << 189 // according to the ordering parameter table
>> 190 // 'true' is returned if the process is registerd successfully
>> 191 G4bool RegisterProcess(G4VProcess* process,
>> 192 G4ParticleDefinition* particle);
132 193
133 // Invoke BuildPhysicsTable for all proces << 194 public:
134 // In case of "Retrieve" flag is ON, Physi << 195 void UseCoupledTransportation(G4bool vl=true);
135 // retrieved from files. <<
136 void BuildPhysicsTable(); <<
137 196
138 // Prepare the PhysicsTable for specified << 197 /////////////////////////////////////////////////////////////////
139 void PreparePhysicsTable(G4ParticleDefinit << 198 public: // with description
>> 199 // "SetCuts" method sets a cut value for all particle types
>> 200 // in the particle table
>> 201 virtual void SetCuts();
>> 202
>> 203 public: // with description
>> 204 // set/get the default cut value
>> 205 // Calling SetDefaultCutValue causes re-calcuration of cut values
>> 206 // and physics tables just before the next event loop
>> 207 void SetDefaultCutValue(G4double newCutValue);
>> 208 G4double GetDefaultCutValue() const;
>> 209
>> 210 /////////////////////////////////////////////////////////////////////
>> 211 public: // with description
>> 212 // Invoke BuildPhysicsTable for all processes for all particles
>> 213 // In case of "Retrieve" flag is ON, PhysicsTable will be
>> 214 // retrieved from files
>> 215 void BuildPhysicsTable();
>> 216
>> 217 // do PreparePhysicsTable for specified particle type
>> 218 void PreparePhysicsTable(G4ParticleDefinition* );
140 219
141 // Build the PhysicsTable for specified pa << 220 // do BuildPhysicsTable for specified particle type
142 void BuildPhysicsTable(G4ParticleDefinitio << 221 void BuildPhysicsTable(G4ParticleDefinition* );
143 222
144 // Store PhysicsTable together with both m << 223 // Store PhysicsTable together with both material and cut value
145 // information in files under the specifie 224 // information in files under the specified directory.
146 // Returns "true" if files are successfull << 225 // (return true if files are sucessfully created)
147 G4bool StorePhysicsTable(const G4String& d << 226 G4bool StorePhysicsTable(const G4String& directory = ".");
148 << 227
149 // Return true if "Retrieve" flag is ON. << 228 // Return true if "Retrieve" flag is ON.
150 // (i.e. PhysicsTable will be retrieved fr << 229 // (i.e. PhysicsTable will be retrieved from files)
151 G4bool IsPhysicsTableRetrieved() const; << 230 G4bool IsPhysicsTableRetrieved() const;
152 G4bool IsStoredInAscii() const; << 231 G4bool IsStoredInAscii() const;
153 232
154 // Get directory path for physics table fi 233 // Get directory path for physics table files.
155 const G4String& GetPhysicsTableDirectory() 234 const G4String& GetPhysicsTableDirectory() const;
156 235
157 // Set "Retrieve" flag. Directory path can << 236 // Set "Retrieve" flag
158 // Null string (default) means directory i << 237 // Directory path can be set together.
159 // from the current value. << 238 // Null string (default) means directory is not changed
160 void SetPhysicsTableRetrieved(const G4Stri << 239 // from the current value
161 void SetStoredInAscii(); << 240 void SetPhysicsTableRetrieved(const G4String& directory = "");
162 << 241 void SetStoredInAscii();
163 // Reset "Retrieve" flag. << 242
164 void ResetPhysicsTableRetrieved(); << 243 // Reset "Retrieve" flag
165 void ResetStoredInAscii(); << 244 void ResetPhysicsTableRetrieved();
166 << 245 void ResetStoredInAscii();
167 // Print out the List of registered partic << 246
>> 247 ///////////////////////////////////////////////////////////////////////
>> 248 public: // with description
>> 249 // Print out the List of registered particles types
168 void DumpList() const; 250 void DumpList() const;
169 251
170 // Request to print out information of cut << 252 public: // with description
171 // Printing will be performed when all tab << 253 // Request to print out information of cut values
172 void DumpCutValuesTable(G4int flag = 1); << 254 // Printing will be performed when all tables are made
173 << 255 void DumpCutValuesTable(G4int flag =1);
174 // Triggers the print-out requested by the << 256
175 // This method must be invoked by RunManag << 257 // The following method actually trigger the print-out requested
>> 258 // by the above method. This method must be invoked by RunManager
>> 259 // at the proper moment.
176 void DumpCutValuesTableIfRequested(); 260 void DumpCutValuesTableIfRequested();
177 261
178 // Set/get control flag for output message << 262 public: // with description
>> 263 void SetVerboseLevel(G4int value);
>> 264 G4int GetVerboseLevel() const;
>> 265 // set/get controle flag for output message
179 // 0: Silent 266 // 0: Silent
180 // 1: Warning message 267 // 1: Warning message
181 // 2: More 268 // 2: More
182 void SetVerboseLevel(G4int value); <<
183 G4int GetVerboseLevel() const; <<
184 <<
185 void UseCoupledTransportation(G4bool vl = <<
186 269
187 // Invokes default SetCuts() method. << 270 ///////////////////////////////////////////////////////////////////////////
188 // Note: cut values will not be overwritte << 271 public: // with description
189 // Use of default SetCuts() method is reco << 272 // "SetCutsWithDefault" method invokes default SetCuts method
190 void SetCutsWithDefault(); << 273 // Note: Cut values will not be overwriten with this method
191 << 274 // Using default SetCuts method is recommended
192 // Sets a cut value for a particle type fo << 275 // (i.e You do not need to implement SetCuts method)
193 void SetCutValue(G4double aCut, const G4St << 276 void SetCutsWithDefault();
194 << 277
195 // Gets a cut value for a particle type fo << 278 // Following are utility methods for SetCuts
196 G4double GetCutValue(const G4String& pname << 279
197 << 280 // SetCutValue sets a cut value for a particle type for the default region
198 // Sets a cut value for a particle type fo << 281 void SetCutValue(G4double aCut, const G4String& pname);
199 void SetCutValue(G4double aCut, const G4St << 282
200 << 283 // GetCutValue sets a cut value for a particle type for the default region
201 // Invoke SetCuts for specified particle f << 284 G4double GetCutValue(const G4String& pname) const;
202 // If the pointer to the region is NULL, t << 285
203 // In case of "Retrieve" flag is ON, cut v << 286 // SetCutValue sets a cut value for a particle type for a region
204 // from files. << 287 void SetCutValue(G4double aCut, const G4String& pname, const G4String& rname);
205 void SetParticleCuts(G4double cut, G4Parti << 288
206 void SetParticleCuts(G4double cut, const G << 289 // Invoke SetCuts for specified particle for a region
207 << 290 // If the pointer to the region is NULL, the default region is used
208 // Invoke SetCuts() for all particles in a << 291 // In case of "Retrieve" flag is ON,
209 void SetCutsForRegion(G4double aCut, const << 292 // Cut values will be retrieved from files
210 << 293 void SetParticleCuts(G4double cut,G4ParticleDefinition* particle,G4Region* region=0);
211 // Gets/sets the flag for ApplyCuts(). << 294 void SetParticleCuts( G4double cut, const G4String& particleName, G4Region* region=0);
212 void SetApplyCuts(G4bool value, const G4St << 295
213 G4bool GetApplyCuts(const G4String& name) << 296 // Invoke SetCuts for all particles in a region
>> 297 void SetCutsForRegion(G4double aCut, const G4String& rname);
>> 298
>> 299 // Following are utility methods are obsolete
>> 300 void ResetCuts();
>> 301
>> 302 ///////////////////////////////////////////////////////////////////
>> 303 public:
>> 304 // Get/SetApplyCuts gets/sets the flag for ApplyCuts
>> 305 void SetApplyCuts(G4bool value, const G4String& name);
>> 306 G4bool GetApplyCuts(const G4String& name) const;
>> 307
>> 308 ///////////////////////////////////////////////////////////////////////////////
>> 309 protected:
>> 310 // do BuildPhysicsTable for make the integral schema
>> 311 void BuildIntegralPhysicsTable(G4VProcess* ,G4ParticleDefinition* );
>> 312
>> 313
>> 314 protected:
>> 315 // Retrieve PhysicsTable from files for proccess belongng the particle.
>> 316 // Normal BuildPhysics procedure of processes will be invoked,
>> 317 // if it fails (in case of Process's RetrievePhysicsTable returns false)
>> 318 virtual void RetrievePhysicsTable(G4ParticleDefinition* ,
>> 319 const G4String& directory,
>> 320 G4bool ascii = false);
>> 321
>> 322 /////////////////////////////////////////////////////////////////
>> 323 protected:
>> 324 // adds new ProcessManager to all particles in the Particle Table
>> 325 // this routine is used in Construct()
>> 326 void InitializeProcessManager();
214 327
215 // Remove and delete ProcessManagers for a << 328 public: // with description
216 // Particle Table. << 329 // remove and delete ProcessManagers for all particles in tha Particle Table
>> 330 // this routine is invoked from RunManager
217 void RemoveProcessManager(); 331 void RemoveProcessManager();
218 332
219 // Remove and delete TrackingManagers for << 333 public: // with description
220 // Particle Table. << 334 // add process manager for particles created on-the-fly
221 void RemoveTrackingManager(); <<
222 <<
223 // Add process manager for particles creat <<
224 void AddProcessManager(G4ParticleDefinitio 335 void AddProcessManager(G4ParticleDefinition* newParticle,
225 G4ProcessManager* n << 336 G4ProcessManager* newManager = 0 );
>> 337
>> 338 /////////////////////////////////////////////////////////////////
>> 339 public:
>> 340 // check consistencies of list of particles
226 341
227 // Check consistencies of list of particle <<
228 void CheckParticleList(); 342 void CheckParticleList();
229 343
230 void DisableCheckParticleList(); 344 void DisableCheckParticleList();
231 << 345
232 inline G4int GetInstanceID() const; << 346 ////////////////////////////////////////////////////////////////////////
233 static const G4VUPLManager& GetSubInstance <<
234 <<
235 // Used by Worker threads on the shared in <<
236 // to initialise workers. Derived class re <<
237 // must also call this base class method. <<
238 virtual void InitializeWorker(); <<
239 <<
240 // Destroy thread-local data. Note that de <<
241 // implementing this method should still c <<
242 virtual void TerminateWorker(); <<
243 <<
244 protected: 347 protected:
245 // User must invoke this method in his Con << 348 // the particle table has the complete List of existing particle types
246 // implementation in order to enable parti << 349 G4ParticleTable* theParticleTable;
247 void AddTransportation(); << 350 //G4ParticleTable::G4PTblDicIterator* theParticleIterator; //AND
248 << 351
249 // Register a process to the particle type << 352 protected:
250 // according to the ordering parameter tab << 353 // pointer to G4UserPhysicsListMessenger
251 // 'true' is returned if the process is re << 354 //G4UserPhysicsListMessenger* theMessenger;
252 G4bool RegisterProcess(G4VProcess* process <<
253 <<
254 // Build PhysicsTable for making the integ <<
255 void BuildIntegralPhysicsTable(G4VProcess* <<
256 <<
257 // Retrieve PhysicsTable from files for pr <<
258 // Normal BuildPhysics procedure of proces <<
259 // fails (in case of process's RetrievePhy <<
260 virtual void RetrievePhysicsTable(G4Partic <<
261 G4bool a <<
262 355
263 // Adds new ProcessManager to all particle << 356 protected:
264 // This function is used in Construct(). << 357 G4int verboseLevel;
265 void InitializeProcessManager(); <<
266 358
267 G4ParticleTable::G4PTblDicIterator* GetPar << 359 protected:
>> 360 // this is the default cut value for all particles
>> 361 G4double defaultCutValue;
>> 362 G4bool isSetDefaultCutValue;
268 363
269 protected: 364 protected:
270 // The particle table has the complete Lis << 365 // pointer to ProductionCutsTable
271 G4ParticleTable* theParticleTable = nullpt << 366 G4ProductionCutsTable* fCutsTable;
272 367
273 G4int verboseLevel = 1; << 368 // flag to determine physics table will be build from file or not
>> 369 G4bool fRetrievePhysicsTable;
>> 370 G4bool fStoredInAscii;
>> 371
>> 372 G4bool fIsCheckedForRetrievePhysicsTable;
>> 373 G4bool fIsRestoredCutValues;
274 374
275 // Default cut value for all particles << 375 // directory name for physics table files
276 G4double defaultCutValue = 1.0; << 376 G4String directoryPhysicsTable;
277 G4bool isSetDefaultCutValue = false; <<
278 377
279 // Pointer to ProductionCutsTable. << 378 // flag for displaying the range cuts & energy thresholds
280 G4ProductionCutsTable* fCutsTable = nullpt << 379 //G4int fDisplayThreshold;
281 380
282 // Flag to determine if physics table will << 381 // flag for Physics Table has been built
283 G4bool fRetrievePhysicsTable = false; << 382 //G4bool fIsPhysicsTableBuilt;
284 G4bool fStoredInAscii = true; <<
285 383
286 G4bool fIsCheckedForRetrievePhysicsTable = << 384 // flag for CheckParticleList
287 G4bool fIsRestoredCutValues = false; << 385 G4bool fDisableCheckParticleList;
288 386
289 // Directory name for physics table files. << 387 // PhysicsListHelper
290 G4String directoryPhysicsTable = "."; << 388 //G4PhysicsListHelper* thePLHelper;
291 389
292 // Flag for CheckParticleList(). << 390 private:
293 G4bool fDisableCheckParticleList = false; << 391 enum { FixedStringLengthForStore = 32 };
294 392
295 // MT data << 393 //Changes for MT
296 G4int g4vuplInstanceID = 0; << 394 protected:
>> 395 G4int g4vuplInstanceID;
297 G4RUN_DLL static G4VUPLManager subInstance 396 G4RUN_DLL static G4VUPLManager subInstanceManager;
>> 397 public:
>> 398 inline G4int GetInstanceID() const;
>> 399 static const G4VUPLManager& GetSubInstanceManager();
>> 400 //Used by Worker threads on the shared instance of
>> 401 // PL to initialize workers
>> 402 void InitializeWorker();
298 }; 403 };
299 404
300 // Inline methods implementations <<
301 <<
302 inline void G4VUserPhysicsList::Construct() 405 inline void G4VUserPhysicsList::Construct()
303 { 406 {
304 #ifdef G4VERBOSE << 407 #ifdef G4VERBOSE
305 if (verboseLevel > 1) G4cout << "G4VUserPhys << 408 if (verboseLevel >1) G4cout << "G4VUserPhysicsList::Construct()" << G4endl;
306 #endif 409 #endif
307 410
308 if (G4Threading::IsMasterThread()) G4Physics <<
309 <<
310 InitializeProcessManager(); 411 InitializeProcessManager();
311 412
312 #ifdef G4VERBOSE << 413 #ifdef G4VERBOSE
313 if (verboseLevel > 1) G4cout << "Construct p << 414 if (verboseLevel >1) G4cout << "Construct processes " << G4endl;
314 #endif 415 #endif
315 ConstructProcess(); 416 ConstructProcess();
>> 417
316 } 418 }
317 419
318 inline G4double G4VUserPhysicsList::GetDefault 420 inline G4double G4VUserPhysicsList::GetDefaultCutValue() const
319 { 421 {
320 return defaultCutValue; 422 return defaultCutValue;
321 } 423 }
322 424
323 inline G4int G4VUserPhysicsList::GetVerboseLev << 425
>> 426 inline G4int G4VUserPhysicsList::GetVerboseLevel() const
324 { 427 {
325 return verboseLevel; << 428 return verboseLevel;
326 } 429 }
327 430
328 inline G4bool G4VUserPhysicsList::IsPhysicsTab << 431 inline
>> 432 G4bool G4VUserPhysicsList::IsPhysicsTableRetrieved() const
329 { 433 {
330 return fRetrievePhysicsTable; << 434 return fRetrievePhysicsTable;
331 } 435 }
332 436
333 inline G4bool G4VUserPhysicsList::IsStoredInAs << 437 inline
>> 438 G4bool G4VUserPhysicsList::IsStoredInAscii() const
334 { 439 {
335 return fStoredInAscii; 440 return fStoredInAscii;
336 } 441 }
337 442
338 inline const G4String& G4VUserPhysicsList::Get << 443 inline
>> 444 const G4String& G4VUserPhysicsList::GetPhysicsTableDirectory() const
339 { 445 {
340 return directoryPhysicsTable; << 446 return directoryPhysicsTable;
341 } 447 }
342 448
343 inline void G4VUserPhysicsList::SetStoredInAsc << 449 inline
>> 450 void G4VUserPhysicsList::SetStoredInAscii()
344 { 451 {
345 fStoredInAscii = true; 452 fStoredInAscii = true;
346 } 453 }
347 << 454
348 inline void G4VUserPhysicsList::ResetPhysicsTa << 455
>> 456 inline
>> 457 void G4VUserPhysicsList::ResetPhysicsTableRetrieved()
349 { 458 {
350 fRetrievePhysicsTable = false; 459 fRetrievePhysicsTable = false;
351 fIsRestoredCutValues = false; 460 fIsRestoredCutValues = false;
352 fIsCheckedForRetrievePhysicsTable = false; << 461 fIsCheckedForRetrievePhysicsTable=false;
353 } 462 }
354 463
355 inline void G4VUserPhysicsList::ResetStoredInA << 464
>> 465 inline
>> 466 void G4VUserPhysicsList::ResetStoredInAscii()
356 { 467 {
357 fStoredInAscii = false; 468 fStoredInAscii = false;
358 } 469 }
359 470
360 inline void G4VUserPhysicsList::DisableCheckPa << 471 inline
>> 472 void G4VUserPhysicsList::DisableCheckParticleList()
361 { 473 {
362 fDisableCheckParticleList = true; 474 fDisableCheckParticleList = true;
363 } 475 }
364 476
365 inline G4int G4VUserPhysicsList::GetInstanceID << 477 inline
>> 478 G4int G4VUserPhysicsList::GetInstanceID() const
366 { 479 {
367 return g4vuplInstanceID; << 480 return g4vuplInstanceID;
368 } 481 }
369 482
370 inline const G4VUPLManager& G4VUserPhysicsList << 483 inline
>> 484 const G4VUPLManager& G4VUserPhysicsList::GetSubInstanceManager()
371 { 485 {
372 return subInstanceManager; << 486 return subInstanceManager;
373 } 487 }
374 <<
375 #endif 488 #endif
>> 489
376 490