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