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 5.2.p2)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                    <<   3 // * DISCLAIMER                                                       *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th <<   5 // * The following disclaimer summarizes all the specific disclaimers *
  6 // * the Geant4 Collaboration.  It is provided <<   6 // * of contributors to this software. The specific disclaimers,which *
  7 // * conditions of the Geant4 Software License <<   7 // * govern, are listed with their locations in:                      *
  8 // * LICENSE and available at  http://cern.ch/ <<   8 // *   http://cern.ch/geant4/license                                  *
  9 // * include a list of copyright holders.      << 
 10 // *                                                9 // *                                                                  *
 11 // * Neither the authors of this software syst     10 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     11 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     12 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     13 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file  <<  14 // * use.                                                             *
 16 // * for the full disclaimer and the limitatio << 
 17 // *                                               15 // *                                                                  *
 18 // * This  code  implementation is the result  <<  16 // * This  code  implementation is the  intellectual property  of the *
 19 // * technical work of the GEANT4 collaboratio <<  17 // * GEANT4 collaboration.                                            *
 20 // * By using,  copying,  modifying or  distri <<  18 // * By copying,  distributing  or modifying the Program (or any work *
 21 // * any work based  on the software)  you  ag <<  19 // * based  on  the Program)  you indicate  your  acceptance of  this *
 22 // * use  in  resulting  scientific  publicati <<  20 // * statement, and all its terms.                                    *
 23 // * acceptance of all terms of the Geant4 Sof << 
 24 // *******************************************     21 // ********************************************************************
 25 //                                                 22 //
 26 // G4VUserPhysicsList                          << 
 27 //                                                 23 //
 28 // Class description:                          <<  24 // $Id: G4VUserPhysicsList.hh,v 1.24 2003/05/21 21:06:00 asaim Exp $
                                                   >>  25 // GEANT4 tag $Name: geant4-05-02-patch-01 $
 29 //                                                 26 //
 30 // This class is an abstract class for constru <<  27 // 
 31 // User must implement the following two pure  <<  28 // ------------------------------------------------------------
 32 // class derived from this class:              <<  29 //  GEANT 4 class header file 
 33 // - G4VUserPhysicsList::ConstructParticle()   <<  30 // Class Description:
 34 //     Construct particles                     <<  31 //      This class is an abstruct class for
 35 // - G4VUserPhysicsList::ConstructProcess()    <<  32 //      constructing particles and processes.
 36 //     Construct procesess and register them t <<  33 //      User must implement following four virtual methods
 37                                                <<  34 //      in his own concrete class derived from this class. 
 38 // Original author: H.Kurashige (Kobe Universi <<  35 //        G4VUserPhysicsList::ConstructParticle() 
 39 // ------------------------------------------- <<  36 //           Construct particles
 40 #ifndef G4VUserPhysicsList_hh                  <<  37 //        G4VUserPhysicsList::constructPhysics() 
 41 #define G4VUserPhysicsList_hh 1                <<  38 //           Construct procesess and register them to particles
                                                   >>  39 //        G4VUserPhysicsList::SetCuts()
                                                   >>  40 //           set cut values in range to all particles
                                                   >>  41 //           (and rebuilding physics table will be invoked )
                                                   >>  42 //
                                                   >>  43 // ------------------------------------------- 
                                                   >>  44 //  History
                                                   >>  45 //        first version                   09 Jan. 1998 by H.Kurashige 
                                                   >>  46 //        modified                        24 Jan. 1998 by H.Kurashige
                                                   >>  47 //          rename DumpCutValues/DumpCutValuesTable 
                                                   >>  48 //          change SetCuts method
                                                   >>  49 //          add    SetCutsWithDefault method
                                                   >>  50 //        modified                       06 June 1998 by H.Kurashige
                                                   >>  51 //          add    AddProcessManager
                                                   >>  52 //          add    BuildPhysicsTable
                                                   >>  53 //       modified                        29  June 1998 by H.Kurashige
                                                   >>  54 //          add    AddProcessManager
                                                   >>  55 //       modified                        05 Dec. 1998 by H.Kurashige
                                                   >>  56 //          add    ConstructAllParticles()
                                                   >>  57 //        modified                        14, Apr 1999 by H.Kurashige
                                                   >>  58 //          change BuildPhysicsTable as public
                                                   >>  59 //          removed ConstructAllParticles() and related methods  
                                                   >>  60 //          changed SetCuts method argument
                                                   >>  61 //       modified                           08, Nov 2000 by H.Kurashige
                                                   >>  62 //          added   Retrieve/StorePhysicsTable and related methods
                                                   >>  63 //       modified                           08, Mar 2001 by H.Kurashige
                                                   >>  64 //          added   binary mode for Retrieve/StorePhysicsTable
                                                   >>  65 //          added   RetrieveCutValues and related
                                                   >>  66 //          added   Set/ResetStoredInAscii() to switch on ascii mode 
                                                   >>  67 //                  for Retrieve/StorePhysicsTable
                                                   >>  68 //       modified for CUTS per REGION      10, Oct 2002 by H.Kurashige 
                                                   >>  69 //          removed following methods 
                                                   >>  70 //           void ReCalcCutValue() 
                                                   >>  71 //           void SetCutValueForOthers()
                                                   >>  72 //           void SetCutValueForOtherThan() 
                                                   >>  73 //           void ReCalcCutValueForOthers()
                                                   >>  74 //           virtual G4bool  StoreMaterialInfo()
                                                   >>  75 //           virtual G4bool  StoreCutValues()
                                                   >>  76 //           virtual G4bool  RetrieveCutValues()
                                                   >>  77 //           virtual G4bool  CheckForRetrievePhysicsTable()
                                                   >>  78 //           virtual G4bool  CheckMaterialInfo()
                                                   >>  79 //          added    void BuildPhysicsTable()    
                                                   >>  80 // ------------------------------------------------------------
                                                   >>  81 #ifndef G4VUserPhysicsList_h
                                                   >>  82 #define G4VUserPhysicsList_h 1
                                                   >>  83 #include "globals.hh"
                                                   >>  84 #include "G4ios.hh"
 42                                                    85 
 43 #include "G4ParticleDefinition.hh"             << 
 44 #include "G4ParticleTable.hh"                      86 #include "G4ParticleTable.hh"
 45 #include "G4PhysicsModelCatalog.hh"            <<  87 #include "G4ParticleDefinition.hh" 
 46 #include "G4ProductionCutsTable.hh"                88 #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                                                    89 
 55 class G4UserPhysicsListMessenger;                  90 class G4UserPhysicsListMessenger;
 56 class G4PhysicsListHelper;                     << 
 57 class G4VProcess;                                  91 class G4VProcess;
 58                                                    92 
 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                           93 class G4VUserPhysicsList
104 {                                                  94 {
105   public:                                      <<  95   public: 
106     G4VUserPhysicsList();                          96     G4VUserPhysicsList();
107     virtual ~G4VUserPhysicsList();                 97     virtual ~G4VUserPhysicsList();
108     G4VUserPhysicsList(const G4VUserPhysicsLis << 
109     G4VUserPhysicsList& operator=(const G4VUse << 
110                                                << 
111     // Each particle type will be instantiated << 
112     // This method is invoked by the RunManger << 
113     virtual void ConstructParticle() = 0;      << 
114                                                    98 
115     // By calling the "Construct" method,      <<  99   public:  // with description
116     // process manager and processes are creat << 100     // By calling the "Construct" method, 
                                                   >> 101     // particles and processes are created    
117     void Construct();                             102     void Construct();
                                                   >> 103  
118                                                   104 
119     // Each physics process will be instantiat << 105   protected: // with description
120     // registered to the process manager of ea << 106    // These two methods of  ConstructParticle() and ConstructProcess()
121     // Invoked in the Construct() method.      << 107    // will be invoked in the Construct() method. 
122     virtual void ConstructProcess() = 0;       << 108 
123                                                << 109    // each particle type will be instantiated
124     // Sets a cut value for all particle types << 110    virtual void ConstructParticle() = 0;
125     virtual void SetCuts();                    << 111  
126                                                << 112    // each physics process will be instantiated and
127     // Set/get the default cut value. Calling  << 113    // registered to the process manager of each particle type 
128     // re-calcuration of cut values and physic << 114    virtual void ConstructProcess() = 0;
129     // next event loop.                        << 115 
130     void SetDefaultCutValue(G4double newCutVal << 116   protected: // with description
131     G4double GetDefaultCutValue() const;       << 117    //  User must invoke this method in his ConstructProcess() 
132                                                << 118    //  implementation in order to insures particle transportation.
133     // Invoke BuildPhysicsTable for all proces << 119    //  !! Caution: this class must not be overriden !!
                                                   >> 120    void AddTransportation();
                                                   >> 121 
                                                   >> 122   /////////////////////////////////////////////////////////////////
                                                   >> 123   public: // with description 
                                                   >> 124    //  "SetCuts" method sets a cut value for all particle types 
                                                   >> 125    //   in the particle table
                                                   >> 126    virtual void SetCuts() = 0; 
                                                   >> 127 
                                                   >> 128   public:  // with description
                                                   >> 129    //  set/get the default cut value
                                                   >> 130    //  Calling SetDefaultCutValue causes re-calcuration of cut values
                                                   >> 131    //  and physics tables just before the next event loop
                                                   >> 132    void     SetDefaultCutValue(G4double newCutValue);
                                                   >> 133    G4double GetDefaultCutValue() const;
                                                   >> 134 
                                                   >> 135   /////////////////////////////////////////////////////////////////////  
                                                   >> 136   public: // with description
                                                   >> 137     // Invoke BuildPhysicsTable for all processes for all particles
134     // In case of "Retrieve" flag is ON, Physi    138     // In case of "Retrieve" flag is ON, PhysicsTable will be
135     // retrieved from files.                   << 139     // retrieved from files
136     void BuildPhysicsTable();                  << 140     void BuildPhysicsTable();    
137                                                << 141   
138     // Prepare the PhysicsTable for specified  << 142    // do BuildPhysicsTable for specified particle type
139     void PreparePhysicsTable(G4ParticleDefinit << 143     void BuildPhysicsTable(G4ParticleDefinition* );    
140                                                   144 
141     // Build the PhysicsTable for specified pa << 145      // Store PhysicsTable together with both material and cut value 
142     void BuildPhysicsTable(G4ParticleDefinitio << 
143                                                << 
144     // Store PhysicsTable together with both m << 
145     // information in files under the specifie    146     // information in files under the specified directory.
146     // Returns "true" if files are successfull << 147     //  (return true if files are sucessfully created)
147     G4bool StorePhysicsTable(const G4String& d << 148     G4bool  StorePhysicsTable(const G4String& directory = ".");
148                                                << 149  
149     // Return true if "Retrieve" flag is ON.   << 150     // Return true if "Retrieve" flag is ON. 
150     // (i.e. PhysicsTable will be retrieved fr << 151     // (i.e. PhysicsTable will be retrieved from files)
151     G4bool IsPhysicsTableRetrieved() const;    << 152     G4bool  IsPhysicsTableRetrieved() const;
152     G4bool IsStoredInAscii() const;            << 153     G4bool  IsStoredInAscii() const;
153                                                   154 
154     // Get directory path for physics table fi    155     // Get directory path for physics table files.
155     const G4String& GetPhysicsTableDirectory()    156     const G4String& GetPhysicsTableDirectory() const;
156                                                   157 
157     // Set "Retrieve" flag. Directory path can << 158     // Set "Retrieve" flag
158     // Null string (default) means directory i << 159     // Directory path can be set together.
159     // from the current value.                 << 160     // Null string (default) means directory is not changed 
160     void SetPhysicsTableRetrieved(const G4Stri << 161     // from the current value 
161     void SetStoredInAscii();                   << 162     void    SetPhysicsTableRetrieved(const G4String& directory = "");
162                                                << 163     void    SetStoredInAscii();
163     // Reset "Retrieve" flag.                  << 164   
164     void ResetPhysicsTableRetrieved();         << 165     // Reset "Retrieve" flag
165     void ResetStoredInAscii();                 << 166     void    ResetPhysicsTableRetrieved();
166                                                << 167     void    ResetStoredInAscii();
167     // Print out the List of registered partic << 168 
                                                   >> 169  ///////////////////////////////////////////////////////////////////////
                                                   >> 170   public: // with description
                                                   >> 171     // Print out the List of registered particles types
168     void DumpList() const;                        172     void DumpList() const;
169                                                   173 
170     // Request to print out information of cut << 174   public: // with description
171     // Printing will be performed when all tab << 175     // Request to print out information of cut values
172     void DumpCutValuesTable(G4int flag = 1);   << 176     // Printing will be performed when all tables are made
173                                                << 177     void DumpCutValuesTable(G4int nParticles=3);
174     // Triggers the print-out requested by the << 178 
175     // This method must be invoked by RunManag << 179     // The following method actually trigger the print-out requested
                                                   >> 180     // by the above method. This method must be invoked by RunManager
                                                   >> 181     // at the proper moment.
176     void DumpCutValuesTableIfRequested();         182     void DumpCutValuesTableIfRequested();
177                                                   183 
178     // Set/get control flag for output message << 184     void DumpCutValues(const G4String &particle_name = "ALL")
                                                   >> 185    {
                                                   >> 186     G4cerr << "WARNING !" << G4endl;
                                                   >> 187     G4cerr << " Obsolete DumpCutValues() method is invoked for " << particle_name << G4endl;
                                                   >> 188     G4cerr << " Please use DumpCutValuesTable() instead." << G4endl;
                                                   >> 189     G4cerr << " This dummy method implementation will be removed soon." << G4endl;
                                                   >> 190     DumpCutValuesTable();
                                                   >> 191    }
                                                   >> 192     void DumpCutValues(G4ParticleDefinition* )
                                                   >> 193    {
                                                   >> 194     G4cerr << "WARNING !" << G4endl;
                                                   >> 195     G4cerr << " DumpCutValues() became obsolete." << G4endl;
                                                   >> 196     G4cerr << " Please use DumpCutValuesTable() instead." << G4endl;
                                                   >> 197     G4cerr << " This dummy method implementation will be removed soon." << G4endl;
                                                   >> 198     DumpCutValuesTable();
                                                   >> 199    }
                                                   >> 200 
                                                   >> 201   public: // with description
                                                   >> 202     void  SetVerboseLevel(G4int value);
                                                   >> 203     G4int GetVerboseLevel() const;
                                                   >> 204     // set/get controle flag for output message
179     //  0: Silent                                 205     //  0: Silent
180     //  1: Warning message                        206     //  1: Warning message
181     //  2: More                                   207     //  2: More
182     void SetVerboseLevel(G4int value);         << 
183     G4int GetVerboseLevel() const;             << 
184                                                << 
185     void UseCoupledTransportation(G4bool vl =  << 
186                                                   208 
187     // Invokes default SetCuts() method.       << 209   ///////////////////////////////////////////////////////////////////////////
188     // Note: cut values will not be overwritte << 210   public: // with description
189     // Use of default SetCuts() method is reco << 211    //  "SetCutsWithDefault" method sets the default cut value
190     void SetCutsWithDefault();                 << 212    //   for all particles for the default region.
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    // SetCutValue sets a cut value for a particle type for a region
199     void SetCutValue(G4double aCut, const G4St << 221    void SetCutValue(G4double aCut, const G4String& pname, const G4String& rname); 
200                                                << 222 
201     // Invoke SetCuts for specified particle f << 223    // Invoke SetCuts for specified particle for a region
202     // If the pointer to the region is NULL, t << 224    // If the pointer to the region is NULL, the default region is used
203     // In case of "Retrieve" flag is ON, cut v << 225    // In case of "Retrieve" flag is ON, 
204     // from files.                             << 226    // Cut values will be retrieved from files
205     void SetParticleCuts(G4double cut, G4Parti << 227    void SetParticleCuts(G4double cut,G4ParticleDefinition* particle,G4Region* region=0);
206     void SetParticleCuts(G4double cut, const G << 228 
207                                                << 229    // Invoke SetCuts for all particles in a region
208     // Invoke SetCuts() for all particles in a << 230    void SetCutsForRegion(G4double aCut, const G4String& rname);
209     void SetCutsForRegion(G4double aCut, const << 231 
210                                                << 232    // Following are utility methods are obsolete
211     // Gets/sets the flag for ApplyCuts().     << 233    void ResetCuts();
212     void SetApplyCuts(G4bool value, const G4St << 234 
213     G4bool GetApplyCuts(const G4String& name)  << 235 ///////////////////////////////////////////////////////////////////
                                                   >> 236   public:   
                                                   >> 237    // Get/SetApplyCuts gets/sets the flag for ApplyCuts
                                                   >> 238    void SetApplyCuts(G4bool value, const G4String& name); 
                                                   >> 239    G4bool GetApplyCuts(const G4String& name) const; 
                                                   >> 240 
                                                   >> 241 ///////////////////////////////////////////////////////////////////////////////
                                                   >> 242   protected:  
                                                   >> 243     // do BuildPhysicsTable for make the integral schema
                                                   >> 244     void BuildIntegralPhysicsTable(G4VProcess* ,G4ParticleDefinition*  );   
                                                   >> 245 
                                                   >> 246 
                                                   >> 247   protected: 
                                                   >> 248     // Retrieve PhysicsTable from files for proccess belongng the particle.
                                                   >> 249     // Normal BuildPhysics procedure of processes will be invoked, 
                                                   >> 250     // if it fails (in case of Process's RetrievePhysicsTable returns false)
                                                   >> 251     virtual void  RetrievePhysicsTable(G4ParticleDefinition* ,  
                                                   >> 252                const G4String& directory, 
                                                   >> 253                G4bool          ascii = false);
                                                   >> 254 
                                                   >> 255    /////////////////////////////////////////////////////////////////
                                                   >> 256   protected: 
                                                   >> 257     // adds new ProcessManager to all particles in the Particle Table
                                                   >> 258     //   this routine is used in Construct()
                                                   >> 259     void InitializeProcessManager();
214                                                   260 
215     // Remove and delete ProcessManagers for a << 261   public: // with description
216     // Particle Table.                         << 262     // remove and delete ProcessManagers for all particles in tha Particle Table
                                                   >> 263     //    this routine is invoked from RunManager 
217     void RemoveProcessManager();                  264     void RemoveProcessManager();
218                                                   265 
219     // Remove and delete TrackingManagers for  << 266   public: // with description
220     // Particle Table.                         << 267     // add process manager for particles created on-the-fly 
221     void RemoveTrackingManager();              << 
222                                                << 
223     // Add process manager for particles creat << 
224     void AddProcessManager(G4ParticleDefinitio    268     void AddProcessManager(G4ParticleDefinition* newParticle,
225                            G4ProcessManager* n << 269          G4ProcessManager*    newManager = 0 );
226                                                << 270  
227     // Check consistencies of list of particle << 271  
228     void CheckParticleList();                  << 272   ////////////////////////////////////////////////////////////////////////
229                                                << 
230     void DisableCheckParticleList();           << 
231                                                << 
232     inline G4int GetInstanceID() const;        << 
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:                                      273   protected:
245     // User must invoke this method in his Con << 274     // the particle table has the complete List of existing particle types
246     // implementation in order to enable parti << 275     G4ParticleTable* theParticleTable;
247     void AddTransportation();                  << 276     G4ParticleTable::G4PTblDicIterator* theParticleIterator;
248                                                << 277 
249     // Register a process to the particle type << 278   protected: 
250     // according to the ordering parameter tab << 279   // pointer to G4UserPhysicsListMessenger
251     // 'true' is returned if the process is re << 280     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                                                << 
263     // Adds new ProcessManager to all particle << 
264     // This function is used in Construct().   << 
265     void InitializeProcessManager();           << 
266                                                << 
267     G4ParticleTable::G4PTblDicIterator* GetPar << 
268                                                   281 
269   protected:                                      282   protected:
270     // The particle table has the complete Lis << 283     G4int verboseLevel;
271     G4ParticleTable* theParticleTable = nullpt << 
272                                                << 
273     G4int verboseLevel = 1;                    << 
274                                                   284 
275     // Default cut value for all particles     << 285   protected:
276     G4double defaultCutValue = 1.0;            << 286     // this is the default cut value for all particles
277     G4bool isSetDefaultCutValue = false;       << 287     G4double defaultCutValue;
278                                                << 
279     // Pointer to ProductionCutsTable.         << 
280     G4ProductionCutsTable* fCutsTable = nullpt << 
281                                                << 
282     // Flag to determine if physics table will << 
283     G4bool fRetrievePhysicsTable = false;      << 
284     G4bool fStoredInAscii = true;              << 
285                                                << 
286     G4bool fIsCheckedForRetrievePhysicsTable = << 
287     G4bool fIsRestoredCutValues = false;       << 
288                                                   288 
289     // Directory name for physics table files. << 289   protected:
290     G4String directoryPhysicsTable = ".";      << 290    // pointer to ProductionCutsTable
                                                   >> 291    G4ProductionCutsTable* fCutsTable;
291                                                   292 
292     // Flag for CheckParticleList().           << 293    // flag to determine physics table will be build from file or not
293     G4bool fDisableCheckParticleList = false;  << 294    G4bool fRetrievePhysicsTable;  
                                                   >> 295    G4bool fStoredInAscii;
                                                   >> 296  
                                                   >> 297    G4bool fIsCheckedForRetrievePhysicsTable;
                                                   >> 298    G4bool fIsRestoredCutValues;
                                                   >> 299 
                                                   >> 300    // directory name for physics table files 
                                                   >> 301    G4String directoryPhysicsTable;   
                                                   >> 302 
                                                   >> 303    // flag for displaying the range cuts & energy thresholds
                                                   >> 304    G4int fDisplayThreshold;
                                                   >> 305 
                                                   >> 306   private:
                                                   >> 307    enum { FixedStringLengthForStore = 32 }; 
                                                   >> 308 
                                                   >> 309 ////////////////////////////////////////////////////////////////////////////
                                                   >> 310 // Following method is for backward compatibility and removed soon
                                                   >> 311 ////////////////////////////////////////////////////////////////////////////
                                                   >> 312   protected:
                                                   >> 313    void SetCutValueForOthers(G4double)
                                                   >> 314    {
                                                   >> 315     G4cerr << "WARNING !" << G4endl;
                                                   >> 316     G4cerr << " SetCutValueForOthers became obsolete." << G4endl;
                                                   >> 317     G4cerr << " It is harmless to remove this invokation without any side effects." << G4endl;
                                                   >> 318     G4cerr << " This dummy method implementation will be removed soon." << G4endl;
                                                   >> 319    }
294                                                   320 
295     // MT data                                 << 
296     G4int g4vuplInstanceID = 0;                << 
297     G4RUN_DLL static G4VUPLManager subInstance << 
298 };                                                321 };
299                                                   322 
300 // Inline methods implementations              << 
301                                                   323 
302 inline void G4VUserPhysicsList::Construct()       324 inline void G4VUserPhysicsList::Construct()
303 {                                                 325 {
304 #ifdef G4VERBOSE                               << 326 #ifdef G4VERBOSE  
305   if (verboseLevel > 1) G4cout << "G4VUserPhys << 327   if (verboseLevel >1) G4cout << "G4VUserPhysicsList::Construct()" << G4endl;  
306 #endif                                         << 
307                                                   328 
308   if (G4Threading::IsMasterThread()) G4Physics << 329   if (verboseLevel >1) G4cout << "Construct particles " << G4endl;  
                                                   >> 330 #endif
                                                   >> 331   ConstructParticle();
309                                                   332 
310   InitializeProcessManager();                     333   InitializeProcessManager();
311                                                   334 
312 #ifdef G4VERBOSE                               << 335 #ifdef G4VERBOSE  
313   if (verboseLevel > 1) G4cout << "Construct p << 336   if (verboseLevel >1) G4cout << "Construct processes " << G4endl;  
314 #endif                                            337 #endif
315   ConstructProcess();                             338   ConstructProcess();
316 }                                                 339 }
317                                                   340 
318 inline G4double G4VUserPhysicsList::GetDefault    341 inline G4double G4VUserPhysicsList::GetDefaultCutValue() const
319 {                                                 342 {
320   return defaultCutValue;                         343   return defaultCutValue;
321 }                                                 344 }
322                                                   345 
323 inline G4int G4VUserPhysicsList::GetVerboseLev << 346 inline void G4VUserPhysicsList::SetVerboseLevel(G4int value)
                                                   >> 347 {
                                                   >> 348   verboseLevel = value;
                                                   >> 349 #ifdef G4VERBOSE
                                                   >> 350   if (verboseLevel >1){
                                                   >> 351     G4cout << "G4VUserPhysicsList::SetVerboseLevel  :";
                                                   >> 352     G4cout << " Verbose level is set to " << verboseLevel << G4endl;
                                                   >> 353   }
                                                   >> 354 #endif
                                                   >> 355 }
                                                   >> 356 
                                                   >> 357 inline  G4int G4VUserPhysicsList::GetVerboseLevel() const
324 {                                                 358 {
325   return verboseLevel;                         << 359   return  verboseLevel;
326 }                                                 360 }
327                                                   361 
328 inline G4bool G4VUserPhysicsList::IsPhysicsTab << 362 inline  
                                                   >> 363  G4bool  G4VUserPhysicsList::IsPhysicsTableRetrieved() const
329 {                                                 364 {
330   return fRetrievePhysicsTable;                << 365   return fRetrievePhysicsTable;  
331 }                                                 366 }
332                                                   367 
333 inline G4bool G4VUserPhysicsList::IsStoredInAs << 368 inline  
                                                   >> 369  G4bool  G4VUserPhysicsList::IsStoredInAscii() const
334 {                                                 370 {
335   return fStoredInAscii;                          371   return fStoredInAscii;
336 }                                                 372 }
337                                                   373 
338 inline const G4String& G4VUserPhysicsList::Get << 374 inline 
                                                   >> 375   const G4String& G4VUserPhysicsList::GetPhysicsTableDirectory() const
339 {                                                 376 {
340   return directoryPhysicsTable;                << 377   return directoryPhysicsTable;  
341 }                                                 378 }
342                                                   379 
343 inline void G4VUserPhysicsList::SetStoredInAsc << 380 inline 
                                                   >> 381  void  G4VUserPhysicsList::SetStoredInAscii()
344 {                                                 382 {
345   fStoredInAscii = true;                          383   fStoredInAscii = true;
346 }                                                 384 }
347                                                << 385     
348 inline void G4VUserPhysicsList::ResetPhysicsTa << 386     
                                                   >> 387 inline 
                                                   >> 388  void  G4VUserPhysicsList::ResetPhysicsTableRetrieved()
349 {                                                 389 {
350   fRetrievePhysicsTable = false;                  390   fRetrievePhysicsTable = false;
351   fIsRestoredCutValues = false;                   391   fIsRestoredCutValues = false;
352   fIsCheckedForRetrievePhysicsTable = false;   << 392   fIsCheckedForRetrievePhysicsTable=false;
353 }                                              << 
354                                                << 
355 inline void G4VUserPhysicsList::ResetStoredInA << 
356 {                                              << 
357   fStoredInAscii = false;                      << 
358 }                                                 393 }
359                                                   394 
360 inline void G4VUserPhysicsList::DisableCheckPa << 
361 {                                              << 
362   fDisableCheckParticleList = true;            << 
363 }                                              << 
364                                                << 
365 inline G4int G4VUserPhysicsList::GetInstanceID << 
366 {                                              << 
367   return g4vuplInstanceID;                     << 
368 }                                              << 
369                                                   395 
370 inline const G4VUPLManager& G4VUserPhysicsList << 396 inline 
                                                   >> 397  void  G4VUserPhysicsList::ResetStoredInAscii()
371 {                                                 398 {
372   return subInstanceManager;                   << 399   fStoredInAscii = false;
373 }                                                 400 }
374                                                << 
375 #endif                                            401 #endif
                                                   >> 402 
376                                                   403