Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/run/include/G4VUserPhysicsList.hh

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /run/include/G4VUserPhysicsList.hh (Version 11.3.0) and /run/include/G4VUserPhysicsList.hh (Version 9.6.p3)


  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$
 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 // ------------------------------------------------------------
                                                   >>  83 #ifndef G4VUserPhysicsList_h
                                                   >>  84 #define G4VUserPhysicsList_h 1
                                                   >>  85 #include "globals.hh"
                                                   >>  86 #include "G4ios.hh"
 42                                                    87 
 43 #include "G4ParticleDefinition.hh"             << 
 44 #include "G4ParticleTable.hh"                      88 #include "G4ParticleTable.hh"
 45 #include "G4PhysicsModelCatalog.hh"            <<  89 #include "G4ParticleDefinition.hh" 
 46 #include "G4ProductionCutsTable.hh"                90 #include "G4ProductionCutsTable.hh"
 47 #include "G4Threading.hh"                      << 
 48 #include "G4VUPLSplitter.hh"                   << 
 49 #include "G4ios.hh"                            << 
 50 #include "globals.hh"                          << 
 51                                                << 
 52 #include "rundefs.hh"                          << 
 53 #include "tls.hh"                              << 
 54                                                    91 
 55 class G4UserPhysicsListMessenger;                  92 class G4UserPhysicsListMessenger;
 56 class G4PhysicsListHelper;                         93 class G4PhysicsListHelper;
 57 class G4VProcess;                                  94 class G4VProcess;
 58                                                    95 
 59 // Encapsulate the fields of class G4VUserPhys << 
 60 class G4VUPLData                               << 
 61 {                                              << 
 62   public:                                      << 
 63     void initialize();                         << 
 64                                                << 
 65     G4ParticleTable::G4PTblDicIterator* _thePa << 
 66     G4UserPhysicsListMessenger* _theMessenger  << 
 67     G4PhysicsListHelper* _thePLHelper = nullpt << 
 68     G4bool _fIsPhysicsTableBuilt = false;      << 
 69     G4int _fDisplayThreshold = 0;              << 
 70 };                                             << 
 71                                                << 
 72 // The type G4VUPLManager is introduced to enc << 
 73 // both the master thread and worker threads t << 
 74 // the fields encapsulated by the class G4VUPL << 
 75 // changes the value for these fields, it refe << 
 76 // definition defined below. For every G4VUser << 
 77 // there is a corresponding G4VUPLData instanc << 
 78 // are organized by the class G4VUPLManager as << 
 79 // The field "int g4vuplInstanceID" is added t << 
 80 // The value of this field in each G4VUserPhys << 
 81 // subscript of the corresponding G44VUPLData  << 
 82 // In order to use the class G44VUPLManager, w << 
 83 // G4VUserPhysicsList as follows: "static G4VU << 
 84 // Both the master thread and worker threads c << 
 85 // for G44VUPLData instances mutually along wi << 
 86 // instances are created. For each worker thre << 
 87 // Consider any thread A, if there is any othe << 
 88 // This ion is shared by the thread A. So the  << 
 89 // in the array of G4PDefData instances for th << 
 90 //                                             << 
 91 // Important Note: you may wonder why we are i << 
 92 //                 since there is only one PL  << 
 93 //                 This is true, in the sense  << 
 94 //                 to be associated to a G4Run << 
 95 //                 instantiate as many PLs are << 
 96 //                 of the PLs to be used we th << 
 97 //                 guarantee that the system w << 
 98 //                 this (unusual) case. This m << 
 99 //                                             << 
100 using G4VUPLManager = G4VUPLSplitter<G4VUPLDat << 
101 using G4VUserPhysicsListSubInstanceManager = G << 
102                                                << 
103 class G4VUserPhysicsList                           96 class G4VUserPhysicsList
104 {                                                  97 {
105   public:                                      <<  98   public: 
106     G4VUserPhysicsList();                          99     G4VUserPhysicsList();
107     virtual ~G4VUserPhysicsList();                100     virtual ~G4VUserPhysicsList();
                                                   >> 101 
                                                   >> 102   // copy constructor and assignment operator
108     G4VUserPhysicsList(const G4VUserPhysicsLis    103     G4VUserPhysicsList(const G4VUserPhysicsList&);
109     G4VUserPhysicsList& operator=(const G4VUse << 104     G4VUserPhysicsList & operator=(const G4VUserPhysicsList&);
110                                                   105 
111     // Each particle type will be instantiated << 106   public:  // with description
112     // This method is invoked by the RunManger << 107    // Each particle type will be instantiated
113     virtual void ConstructParticle() = 0;      << 108    // This method is invoked by the RunManger 
114                                                << 109    virtual void ConstructParticle() = 0;
115     // By calling the "Construct" method,      << 110 
116     // process manager and processes are creat << 111    // By calling the "Construct" method, 
117     void Construct();                          << 112    // process manager and processes are created. 
118                                                << 113    void Construct();
119     // Each physics process will be instantiat << 114  
120     // registered to the process manager of ea << 115    // Each physics process will be instantiated and
121     // Invoked in the Construct() method.      << 116    // registered to the process manager of each particle type 
122     virtual void ConstructProcess() = 0;       << 117    // This method is invoked in Construct method 
123                                                << 118    virtual void ConstructProcess() = 0;
124     // Sets a cut value for all particle types << 119 
125     virtual void SetCuts();                    << 120   protected: // with description
126                                                << 121    //  User must invoke this method in his ConstructProcess() 
127     // Set/get the default cut value. Calling  << 122    //  implementation in order to insures particle transportation.
128     // re-calcuration of cut values and physic << 123    void AddTransportation();
129     // next event loop.                        << 124 
130     void SetDefaultCutValue(G4double newCutVal << 125    //Register a process to the particle type 
131     G4double GetDefaultCutValue() const;       << 126    // according to the ordering parameter table
                                                   >> 127    //  'true' is returned if the process is registerd successfully
                                                   >> 128    G4bool RegisterProcess(G4VProcess*            process,
                                                   >> 129         G4ParticleDefinition*  particle);
132                                                   130 
133     // Invoke BuildPhysicsTable for all proces << 131   public:
134     // In case of "Retrieve" flag is ON, Physi << 132    void UseCoupledTransportation(G4bool vl=true);
135     // retrieved from files.                   << 
136     void BuildPhysicsTable();                  << 
137                                                   133 
138     // Prepare the PhysicsTable for specified  << 134   /////////////////////////////////////////////////////////////////
139     void PreparePhysicsTable(G4ParticleDefinit << 135   public: // with description 
                                                   >> 136    //  "SetCuts" method sets a cut value for all particle types 
                                                   >> 137    //   in the particle table
                                                   >> 138    virtual void SetCuts(); 
                                                   >> 139 
                                                   >> 140   public:  // with description
                                                   >> 141    //  set/get the default cut value
                                                   >> 142    //  Calling SetDefaultCutValue causes re-calcuration of cut values
                                                   >> 143    //  and physics tables just before the next event loop
                                                   >> 144    void     SetDefaultCutValue(G4double newCutValue);
                                                   >> 145    G4double GetDefaultCutValue() const;
                                                   >> 146 
                                                   >> 147   /////////////////////////////////////////////////////////////////////  
                                                   >> 148   public: // with description
                                                   >> 149     // Invoke BuildPhysicsTable for all processes for all particles
                                                   >> 150     // In case of "Retrieve" flag is ON, PhysicsTable will be
                                                   >> 151     // retrieved from files
                                                   >> 152     void BuildPhysicsTable();    
                                                   >> 153   
                                                   >> 154    // do PreparePhysicsTable for specified particle type
                                                   >> 155     void PreparePhysicsTable(G4ParticleDefinition* );    
140                                                   156 
141     // Build the PhysicsTable for specified pa << 157    // do BuildPhysicsTable for specified particle type
142     void BuildPhysicsTable(G4ParticleDefinitio << 158     void BuildPhysicsTable(G4ParticleDefinition* );    
143                                                   159 
144     // Store PhysicsTable together with both m << 160      // Store PhysicsTable together with both material and cut value 
145     // information in files under the specifie    161     // information in files under the specified directory.
146     // Returns "true" if files are successfull << 162     //  (return true if files are sucessfully created)
147     G4bool StorePhysicsTable(const G4String& d << 163     G4bool  StorePhysicsTable(const G4String& directory = ".");
148                                                << 164  
149     // Return true if "Retrieve" flag is ON.   << 165     // Return true if "Retrieve" flag is ON. 
150     // (i.e. PhysicsTable will be retrieved fr << 166     // (i.e. PhysicsTable will be retrieved from files)
151     G4bool IsPhysicsTableRetrieved() const;    << 167     G4bool  IsPhysicsTableRetrieved() const;
152     G4bool IsStoredInAscii() const;            << 168     G4bool  IsStoredInAscii() const;
153                                                   169 
154     // Get directory path for physics table fi    170     // Get directory path for physics table files.
155     const G4String& GetPhysicsTableDirectory()    171     const G4String& GetPhysicsTableDirectory() const;
156                                                   172 
157     // Set "Retrieve" flag. Directory path can << 173     // Set "Retrieve" flag
158     // Null string (default) means directory i << 174     // Directory path can be set together.
159     // from the current value.                 << 175     // Null string (default) means directory is not changed 
160     void SetPhysicsTableRetrieved(const G4Stri << 176     // from the current value 
161     void SetStoredInAscii();                   << 177     void    SetPhysicsTableRetrieved(const G4String& directory = "");
162                                                << 178     void    SetStoredInAscii();
163     // Reset "Retrieve" flag.                  << 179   
164     void ResetPhysicsTableRetrieved();         << 180     // Reset "Retrieve" flag
165     void ResetStoredInAscii();                 << 181     void    ResetPhysicsTableRetrieved();
166                                                << 182     void    ResetStoredInAscii();
167     // Print out the List of registered partic << 183 
                                                   >> 184  ///////////////////////////////////////////////////////////////////////
                                                   >> 185   public: // with description
                                                   >> 186     // Print out the List of registered particles types
168     void DumpList() const;                        187     void DumpList() const;
169                                                   188 
170     // Request to print out information of cut << 189   public: // with description
171     // Printing will be performed when all tab << 190     // Request to print out information of cut values
172     void DumpCutValuesTable(G4int flag = 1);   << 191     // Printing will be performed when all tables are made
173                                                << 192     void DumpCutValuesTable(G4int flag =1);
174     // Triggers the print-out requested by the << 193 
175     // This method must be invoked by RunManag << 194     // The following method actually trigger the print-out requested
                                                   >> 195     // by the above method. This method must be invoked by RunManager
                                                   >> 196     // at the proper moment.
176     void DumpCutValuesTableIfRequested();         197     void DumpCutValuesTableIfRequested();
177                                                   198 
178     // Set/get control flag for output message << 199   public: // with description
                                                   >> 200     void  SetVerboseLevel(G4int value);
                                                   >> 201     G4int GetVerboseLevel() const;
                                                   >> 202     // set/get controle flag for output message
179     //  0: Silent                                 203     //  0: Silent
180     //  1: Warning message                        204     //  1: Warning message
181     //  2: More                                   205     //  2: More
182     void SetVerboseLevel(G4int value);         << 
183     G4int GetVerboseLevel() const;             << 
184                                                   206 
185     void UseCoupledTransportation(G4bool vl =  << 207   ///////////////////////////////////////////////////////////////////////////
186                                                << 208   public: // with description
187     // Invokes default SetCuts() method.       << 209    //  "SetCutsWithDefault" method invokes default SetCuts method
188     // Note: cut values will not be overwritte << 210    //   Note: Cut values will not be overwriten with this method 
189     // Use of default SetCuts() method is reco << 211    //   Using default SetCuts method is recommended
190     void SetCutsWithDefault();                 << 212    //  (i.e You do not need to implement SetCuts method) 
191                                                << 213    void SetCutsWithDefault();   
192     // Sets a cut value for a particle type fo << 214 
193     void SetCutValue(G4double aCut, const G4St << 215    // Following are utility methods for SetCuts
194                                                << 216   
195     // Gets a cut value for a particle type fo << 217    // SetCutValue sets a cut value for a particle type for the default region
196     G4double GetCutValue(const G4String& pname << 218    void SetCutValue(G4double aCut, const G4String& pname); 
197                                                << 219 
198     // Sets a cut value for a particle type fo << 220    // GetCutValue sets a cut value for a particle type for the default region
199     void SetCutValue(G4double aCut, const G4St << 221    G4double GetCutValue(const G4String& pname) const; 
200                                                << 222 
201     // Invoke SetCuts for specified particle f << 223    // SetCutValue sets a cut value for a particle type for a region
202     // If the pointer to the region is NULL, t << 224    void SetCutValue(G4double aCut, const G4String& pname, const G4String& rname); 
203     // In case of "Retrieve" flag is ON, cut v << 225 
204     // from files.                             << 226    // Invoke SetCuts for specified particle for a region
205     void SetParticleCuts(G4double cut, G4Parti << 227    // If the pointer to the region is NULL, the default region is used
206     void SetParticleCuts(G4double cut, const G << 228    // In case of "Retrieve" flag is ON, 
207                                                << 229    // Cut values will be retrieved from files
208     // Invoke SetCuts() for all particles in a << 230    void SetParticleCuts(G4double cut,G4ParticleDefinition* particle,G4Region* region=0);
209     void SetCutsForRegion(G4double aCut, const << 231   void SetParticleCuts( G4double cut, const G4String& particleName, G4Region* region=0);
210                                                << 232 
211     // Gets/sets the flag for ApplyCuts().     << 233    // Invoke SetCuts for all particles in a region
212     void SetApplyCuts(G4bool value, const G4St << 234    void SetCutsForRegion(G4double aCut, const G4String& rname);
213     G4bool GetApplyCuts(const G4String& name)  << 235 
                                                   >> 236    // Following are utility methods are obsolete
                                                   >> 237    void ResetCuts();
                                                   >> 238 
                                                   >> 239 ///////////////////////////////////////////////////////////////////
                                                   >> 240   public:   
                                                   >> 241    // Get/SetApplyCuts gets/sets the flag for ApplyCuts
                                                   >> 242    void SetApplyCuts(G4bool value, const G4String& name); 
                                                   >> 243    G4bool GetApplyCuts(const G4String& name) const; 
                                                   >> 244 
                                                   >> 245 ///////////////////////////////////////////////////////////////////////////////
                                                   >> 246   protected:  
                                                   >> 247     // do BuildPhysicsTable for make the integral schema
                                                   >> 248     void BuildIntegralPhysicsTable(G4VProcess* ,G4ParticleDefinition*  );   
                                                   >> 249 
                                                   >> 250 
                                                   >> 251   protected: 
                                                   >> 252     // Retrieve PhysicsTable from files for proccess belongng the particle.
                                                   >> 253     // Normal BuildPhysics procedure of processes will be invoked, 
                                                   >> 254     // if it fails (in case of Process's RetrievePhysicsTable returns false)
                                                   >> 255     virtual void  RetrievePhysicsTable(G4ParticleDefinition* ,  
                                                   >> 256                const G4String& directory, 
                                                   >> 257                G4bool          ascii = false);
                                                   >> 258 
                                                   >> 259    /////////////////////////////////////////////////////////////////
                                                   >> 260   protected: 
                                                   >> 261     // adds new ProcessManager to all particles in the Particle Table
                                                   >> 262     //   this routine is used in Construct()
                                                   >> 263     void InitializeProcessManager();
214                                                   264 
215     // Remove and delete ProcessManagers for a << 265   public: // with description
216     // Particle Table.                         << 266     // remove and delete ProcessManagers for all particles in tha Particle Table
                                                   >> 267     //    this routine is invoked from RunManager 
217     void RemoveProcessManager();                  268     void RemoveProcessManager();
218                                                   269 
219     // Remove and delete TrackingManagers for  << 270   public: // with description
220     // Particle Table.                         << 271     // add process manager for particles created on-the-fly 
221     void RemoveTrackingManager();              << 
222                                                << 
223     // Add process manager for particles creat << 
224     void AddProcessManager(G4ParticleDefinitio    272     void AddProcessManager(G4ParticleDefinition* newParticle,
225                            G4ProcessManager* n << 273          G4ProcessManager*    newManager = 0 );
                                                   >> 274  
                                                   >> 275    /////////////////////////////////////////////////////////////////
                                                   >> 276   public:
                                                   >> 277     // check consistencies of list of particles 
226                                                   278 
227     // Check consistencies of list of particle << 
228     void CheckParticleList();                     279     void CheckParticleList();
229                                                   280 
230     void DisableCheckParticleList();              281     void DisableCheckParticleList();
231                                                << 282  
232     inline G4int GetInstanceID() const;        << 283   ////////////////////////////////////////////////////////////////////////
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:                                      284   protected:
245     // User must invoke this method in his Con << 285    // the particle table has the complete List of existing particle types
246     // implementation in order to enable parti << 286    G4ParticleTable* theParticleTable;
247     void AddTransportation();                  << 287    G4ParticleTable::G4PTblDicIterator* theParticleIterator;
248                                                << 288 
249     // Register a process to the particle type << 289   protected: 
250     // according to the ordering parameter tab << 290    // pointer to G4UserPhysicsListMessenger
251     // 'true' is returned if the process is re << 291    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                                                   292 
263     // Adds new ProcessManager to all particle << 293   protected:
264     // This function is used in Construct().   << 294    G4int verboseLevel;
265     void InitializeProcessManager();           << 
266                                                   295 
267     G4ParticleTable::G4PTblDicIterator* GetPar << 296   protected:
                                                   >> 297    // this is the default cut value for all particles
                                                   >> 298    G4double defaultCutValue;
                                                   >> 299   G4bool   isSetDefaultCutValue;
268                                                   300 
269   protected:                                      301   protected:
270     // The particle table has the complete Lis << 302    // pointer to ProductionCutsTable
271     G4ParticleTable* theParticleTable = nullpt << 303    G4ProductionCutsTable* fCutsTable;
272                                                   304 
273     G4int verboseLevel = 1;                    << 305    // flag to determine physics table will be build from file or not
                                                   >> 306    G4bool fRetrievePhysicsTable;  
                                                   >> 307    G4bool fStoredInAscii;
                                                   >> 308  
                                                   >> 309    G4bool fIsCheckedForRetrievePhysicsTable;
                                                   >> 310    G4bool fIsRestoredCutValues;
274                                                   311 
275     // Default cut value for all particles     << 312    // directory name for physics table files 
276     G4double defaultCutValue = 1.0;            << 313    G4String directoryPhysicsTable;   
277     G4bool isSetDefaultCutValue = false;       << 
278                                                   314 
279     // Pointer to ProductionCutsTable.         << 315    // flag for displaying the range cuts & energy thresholds
280     G4ProductionCutsTable* fCutsTable = nullpt << 316    G4int fDisplayThreshold;
281                                                   317 
282     // Flag to determine if physics table will << 318   // flag for Physics Table has been built 
283     G4bool fRetrievePhysicsTable = false;      << 319    G4bool fIsPhysicsTableBuilt;
284     G4bool fStoredInAscii = true;              << 
285                                                   320 
286     G4bool fIsCheckedForRetrievePhysicsTable = << 321   // flag for CheckParticleList 
287     G4bool fIsRestoredCutValues = false;       << 322   G4bool fDisableCheckParticleList; 
288                                                   323 
289     // Directory name for physics table files. << 324   // PhysicsListHelper
290     G4String directoryPhysicsTable = ".";      << 325   G4PhysicsListHelper* thePLHelper;
291                                                   326 
292     // Flag for CheckParticleList().           << 327   private:
293     G4bool fDisableCheckParticleList = false;  << 328    enum { FixedStringLengthForStore = 32 }; 
294                                                   329 
295     // MT data                                 << 
296     G4int g4vuplInstanceID = 0;                << 
297     G4RUN_DLL static G4VUPLManager subInstance << 
298 };                                             << 
299                                                   330 
300 // Inline methods implementations              << 331 };
301                                                   332 
302 inline void G4VUserPhysicsList::Construct()       333 inline void G4VUserPhysicsList::Construct()
303 {                                                 334 {
304 #ifdef G4VERBOSE                               << 335 #ifdef G4VERBOSE  
305   if (verboseLevel > 1) G4cout << "G4VUserPhys << 336   if (verboseLevel >1) G4cout << "G4VUserPhysicsList::Construct()" << G4endl;  
306 #endif                                            337 #endif
307                                                   338 
308   if (G4Threading::IsMasterThread()) G4Physics << 
309                                                << 
310   InitializeProcessManager();                     339   InitializeProcessManager();
311                                                   340 
312 #ifdef G4VERBOSE                               << 341 #ifdef G4VERBOSE  
313   if (verboseLevel > 1) G4cout << "Construct p << 342   if (verboseLevel >1) G4cout << "Construct processes " << G4endl;  
314 #endif                                            343 #endif
315   ConstructProcess();                             344   ConstructProcess();
                                                   >> 345 
316 }                                                 346 }
317                                                   347 
318 inline G4double G4VUserPhysicsList::GetDefault    348 inline G4double G4VUserPhysicsList::GetDefaultCutValue() const
319 {                                                 349 {
320   return defaultCutValue;                         350   return defaultCutValue;
321 }                                                 351 }
322                                                   352 
323 inline G4int G4VUserPhysicsList::GetVerboseLev << 353 
                                                   >> 354 inline  G4int G4VUserPhysicsList::GetVerboseLevel() const
324 {                                                 355 {
325   return verboseLevel;                         << 356   return  verboseLevel;
326 }                                                 357 }
327                                                   358 
328 inline G4bool G4VUserPhysicsList::IsPhysicsTab << 359 inline  
                                                   >> 360  G4bool  G4VUserPhysicsList::IsPhysicsTableRetrieved() const
329 {                                                 361 {
330   return fRetrievePhysicsTable;                << 362   return fRetrievePhysicsTable;  
331 }                                                 363 }
332                                                   364 
333 inline G4bool G4VUserPhysicsList::IsStoredInAs << 365 inline  
                                                   >> 366  G4bool  G4VUserPhysicsList::IsStoredInAscii() const
334 {                                                 367 {
335   return fStoredInAscii;                          368   return fStoredInAscii;
336 }                                                 369 }
337                                                   370 
338 inline const G4String& G4VUserPhysicsList::Get << 371 inline 
                                                   >> 372   const G4String& G4VUserPhysicsList::GetPhysicsTableDirectory() const
339 {                                                 373 {
340   return directoryPhysicsTable;                << 374   return directoryPhysicsTable;  
341 }                                                 375 }
342                                                   376 
343 inline void G4VUserPhysicsList::SetStoredInAsc << 377 inline 
                                                   >> 378  void  G4VUserPhysicsList::SetStoredInAscii()
344 {                                                 379 {
345   fStoredInAscii = true;                          380   fStoredInAscii = true;
346 }                                                 381 }
347                                                << 382     
348 inline void G4VUserPhysicsList::ResetPhysicsTa << 383     
                                                   >> 384 inline 
                                                   >> 385  void  G4VUserPhysicsList::ResetPhysicsTableRetrieved()
349 {                                                 386 {
350   fRetrievePhysicsTable = false;                  387   fRetrievePhysicsTable = false;
351   fIsRestoredCutValues = false;                   388   fIsRestoredCutValues = false;
352   fIsCheckedForRetrievePhysicsTable = false;   << 389   fIsCheckedForRetrievePhysicsTable=false;
353 }                                                 390 }
354                                                   391 
355 inline void G4VUserPhysicsList::ResetStoredInA << 392 
                                                   >> 393 inline 
                                                   >> 394  void  G4VUserPhysicsList::ResetStoredInAscii()
356 {                                                 395 {
357   fStoredInAscii = false;                         396   fStoredInAscii = false;
358 }                                                 397 }
359                                                   398 
360 inline void G4VUserPhysicsList::DisableCheckPa << 399 inline
                                                   >> 400  void G4VUserPhysicsList::DisableCheckParticleList()
361 {                                                 401 {
362   fDisableCheckParticleList = true;               402   fDisableCheckParticleList = true;
363 }                                                 403 }
364                                                   404 
365 inline G4int G4VUserPhysicsList::GetInstanceID << 
366 {                                              << 
367   return g4vuplInstanceID;                     << 
368 }                                              << 
369                                                << 
370 inline const G4VUPLManager& G4VUserPhysicsList << 
371 {                                              << 
372   return subInstanceManager;                   << 
373 }                                              << 
374                                                << 
375 #endif                                            405 #endif
                                                   >> 406 
376                                                   407