Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/dna/management/include/G4Scheduler.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 /processes/electromagnetic/dna/management/include/G4Scheduler.hh (Version 11.3.0) and /processes/electromagnetic/dna/management/include/G4Scheduler.hh (Version 11.1)


  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 //                                                 26 //
 27 // Author: Mathieu Karamitros                      27 // Author: Mathieu Karamitros
 28                                                    28 
 29 // The code is developed in the framework of t     29 // The code is developed in the framework of the ESA AO7146
 30 //                                                 30 //
 31 // We would be very happy hearing from you, se     31 // We would be very happy hearing from you, send us your feedback! :)
 32 //                                                 32 //
 33 // In order for Geant4-DNA to be maintained an     33 // In order for Geant4-DNA to be maintained and still open-source,
 34 // article citations are crucial.              <<  34 // article citations are crucial. 
 35 // If you use Geant4-DNA chemistry and you pub <<  35 // If you use Geant4-DNA chemistry and you publish papers about your software, 
 36 // in addition to the general paper on Geant4-     36 // in addition to the general paper on Geant4-DNA:
 37 //                                                 37 //
 38 // Int. J. Model. Simul. Sci. Comput. 1 (2010)     38 // Int. J. Model. Simul. Sci. Comput. 1 (2010) 157–178
 39 //                                                 39 //
 40 // we would be very happy if you could please      40 // we would be very happy if you could please also cite the following
 41 // reference papers on chemistry:                  41 // reference papers on chemistry:
 42 //                                                 42 //
 43 // J. Comput. Phys. 274 (2014) 841-882             43 // J. Comput. Phys. 274 (2014) 841-882
 44 // Prog. Nucl. Sci. Tec. 2 (2011) 503-508      <<  44 // Prog. Nucl. Sci. Tec. 2 (2011) 503-508 
                                                   >>  45 
                                                   >>  46 
 45                                                    47 
 46 #ifndef G4Scheduler_h                              48 #ifndef G4Scheduler_h
 47 #define G4Scheduler_h                              49 #define G4Scheduler_h
 48                                                    50 
                                                   >>  51 #include <G4VScheduler.hh>
                                                   >>  52 #include <vector>
                                                   >>  53 #include <map>
                                                   >>  54 #include <memory>
                                                   >>  55 
                                                   >>  56 #include "globals.hh"
                                                   >>  57 
 49 #include "G4ITModelHandler.hh"                     58 #include "G4ITModelHandler.hh"
 50 #include "G4ITReaction.hh"                     << 
 51 #include "G4ITStepStatus.hh"                       59 #include "G4ITStepStatus.hh"
 52 #include "G4ITTrackHolder.hh"                      60 #include "G4ITTrackHolder.hh"
 53 #include "G4VScavengerMaterial.hh"             << 
 54 #include "G4VStateDependent.hh"                    61 #include "G4VStateDependent.hh"
 55 #include "globals.hh"                          <<  62 #include "G4ITReaction.hh"
 56                                                <<  63 #include "G4VScavengerMaterial.hh"
 57 #include <G4VScheduler.hh>                     << 
 58                                                << 
 59 #include <map>                                 << 
 60 #include <memory>                              << 
 61 #include <vector>                              << 
 62                                                    64 
 63 class G4ITTrackingManager;                         65 class G4ITTrackingManager;
 64 class G4ITModelProcessor;                          66 class G4ITModelProcessor;
 65 class G4ITStepProcessor;                           67 class G4ITStepProcessor;
 66 class G4Track;                                     68 class G4Track;
 67 class G4UserTimeStepAction;                        69 class G4UserTimeStepAction;
 68 class G4SchedulerMessenger;                        70 class G4SchedulerMessenger;
 69 class G4ITTrackingInteractivity;                   71 class G4ITTrackingInteractivity;
 70 class G4ITGun;                                     72 class G4ITGun;
 71                                                    73 
 72 #ifndef compTrackPerID__                           74 #ifndef compTrackPerID__
 73 #  define compTrackPerID__                     <<  75 #define compTrackPerID__
 74 struct compTrackPerID                          <<  76   struct compTrackPerID
 75 {                                              <<  77   {
 76     G4bool operator()(G4Track* rhs, G4Track* l     78     G4bool operator()(G4Track* rhs, G4Track* lhs) const
 77     {                                              79     {
 78       return rhs->GetTrackID() < lhs->GetTrack     80       return rhs->GetTrackID() < lhs->GetTrackID();
 79     }                                              81     }
 80 };                                             <<  82   };
 81 #endif                                             83 #endif
 82                                                    84 
 83 /**                                                85 /**
 84  * G4Scheduler synchronizes (in time) track st     86  * G4Scheduler synchronizes (in time) track stepping
 85  */                                                87  */
 86 class G4Scheduler : public G4VScheduler, publi <<  88 class G4Scheduler :
 87 {                                              <<  89     public G4VScheduler,
 88   protected:                                   <<  90     public G4VStateDependent
 89     ~G4Scheduler() override;                   <<  91 {
 90                                                <<  92 protected:
 91   public:                                      <<  93   virtual ~G4Scheduler();
 92     G4Scheduler(const G4Scheduler&) = delete;  <<  94 
 93     G4Scheduler& operator=(const G4Scheduler&) <<  95 public:
 94                                                <<  96   static G4Scheduler* Instance();
 95     static G4Scheduler* Instance();            <<  97   /** DeleteInstance should be used instead
 96     /** DeleteInstance should be used instead  <<  98    * of the destructor
 97      * of the destructor                       <<  99    */
 98      */                                        << 100   static void DeleteInstance();
 99     static void DeleteInstance();              << 101   virtual G4bool Notify(G4ApplicationState requestedState);
100     G4bool Notify(G4ApplicationState requested << 102 
101                                                << 103   virtual void RegisterModel(G4VITStepModel*, G4double);
102     void RegisterModel(G4VITStepModel*, G4doub << 104 
103                                                << 105   void Initialize();
104     void Initialize() override;                << 106   void ForceReinitialization();
105     void ForceReinitialization();              << 107   inline G4bool IsInitialized();
106     inline G4bool IsInitialized();             << 108   inline G4bool IsRunning(){return fRunning;}
107     inline G4bool IsRunning() override { retur << 109   void Reset();
108     void Reset() override;                     << 110   void Process();
109     void Process() override;                   << 111   void ClearList();
110     void ClearList();                          << 112 
111                                                << 113   inline void SetGun(G4ITGun*);
112     inline void SetGun(G4ITGun*) override;     << 114   inline G4ITGun* GetGun();
113     inline G4ITGun* GetGun();                  << 115 
114                                                << 116   inline void Stop();
115     inline void Stop();                        << 117   void Clear();
116     void Clear();                              << 118 
117                                                << 119   // To be called only in UserReactionAction::EndProcessing()
118     // To be called only in UserReactionAction << 120   // after fRunning flag has been turned off.
119     // after fRunning flag has been turned off << 121   // This is not done automatically before UserReactionAction::EndProcessing()
120     // This is not done automatically before U << 122   // is called in case one would like to access some track information
121     // is called in case one would like to acc << 123   void EndTracking();
122     void EndTracking();                        << 124 
123                                                << 125   void SetEndTime(const G4double);
124     void SetEndTime(const G4double) override;  << 126 
125                                                << 127   /* Two tracks below the time tolerance are supposed to be
126     /* Two tracks below the time tolerance are << 128    * in the same time slice
127      * in the same time slice                  << 129    */
128      */                                        << 130   inline void SetTimeTolerance(G4double);
129     inline void SetTimeTolerance(G4double) ove << 131   inline G4double GetTimeTolerance() const;
130     inline G4double GetTimeTolerance() const o << 132 
131                                                << 133   inline void SetMaxZeroTimeAllowed(G4int);
132     inline void SetMaxZeroTimeAllowed(G4int) o << 134   inline G4int GetMaxZeroTimeAllowed() const;
133     inline G4int GetMaxZeroTimeAllowed() const << 135 
134                                                << 136   inline G4ITModelHandler* GetModelHandler();
135     inline G4ITModelHandler* GetModelHandler() << 137 
136                                                << 138   inline void SetTimeSteps(std::map<G4double, G4double>*);
137     inline void SetTimeSteps(std::map<G4double << 139   inline void AddTimeStep(G4double, G4double);
138     inline void AddTimeStep(G4double, G4double << 140   inline void SetDefaultTimeStep(G4double);
139     inline void SetDefaultTimeStep(G4double) o << 141   G4double GetLimitingTimeStep() const;
140     G4double GetLimitingTimeStep() const overr << 142   inline G4int GetNbSteps() const;
141     inline G4int GetNbSteps() const override;  << 143   inline void SetMaxNbSteps(G4int);
142     inline void SetMaxNbSteps(G4int) override; << 144   inline G4int GetMaxNbSteps() const;
143     inline G4int GetMaxNbSteps() const overrid << 145   inline G4double GetStartTime() const;
144     inline G4double GetStartTime() const overr << 146   inline G4double GetEndTime() const;
145     inline G4double GetEndTime() const overrid << 147   virtual inline G4double GetTimeStep() const;
146     inline G4double GetTimeStep() const overri << 148   inline G4double GetPreviousTimeStep() const;
147     inline G4double GetPreviousTimeStep() cons << 149   inline G4double GetGlobalTime() const;
148     inline G4double GetGlobalTime() const over << 150   inline void SetUserAction(G4UserTimeStepAction*);
149     inline void SetUserAction(G4UserTimeStepAc << 151   inline G4UserTimeStepAction* GetUserTimeStepAction() const;
150     inline G4UserTimeStepAction* GetUserTimeSt << 152 
151                                                << 153   // To use with transportation only, no reactions
152     // To use with transportation only, no rea << 154   inline void UseDefaultTimeSteps(G4bool);
153     inline void UseDefaultTimeSteps(G4bool);   << 155   inline G4bool AreDefaultTimeStepsUsed();
154     inline G4bool AreDefaultTimeStepsUsed();   << 156 
155                                                << 157   inline G4ITStepStatus GetStatus() const;
156     inline G4ITStepStatus GetStatus() const;   << 158 
157                                                << 159   /* 1 : Reaction information
158     /* 1 : Reaction information                << 160    * 2 : (1) + time step information
159      * 2 : (1) + time step information         << 161    * 3 : (2) + step info for individual tracks
160      * 3 : (2) + step info for individual trac << 162    * 4 : (2) + trackList processing info + pushed and killed track info
161      * 4 : (2) + trackList processing info + p << 163    */
162      */                                        << 164   inline void SetVerbose(G4int);
163     inline void SetVerbose(G4int) override;    << 165   
164                                                << 166   inline G4int GetVerbose() const;
165     inline G4int GetVerbose() const;           << 167 
166                                                << 168   inline void WhyDoYouStop();
167     inline void WhyDoYouStop();                << 169 
168                                                << 170   void SetInteractivity(G4ITTrackingInteractivity*);
169     void SetInteractivity(G4ITTrackingInteract << 171   inline G4ITTrackingInteractivity* GetInteractivity();
170     inline G4ITTrackingInteractivity* GetInter << 172 
171                                                << 173   virtual size_t GetNTracks();
172     virtual size_t GetNTracks();               << 174 
173                                                << 175   void GetCollisionType(G4String& interactionType);
174     void GetCollisionType(G4String& interactio << 176 
175                                                << 177   void AddWatchedTime(G4double time)
176     void AddWatchedTime(G4double time) { fWatc << 178   {
                                                   >> 179     fWatchedTimes.insert(time);
                                                   >> 180   }
177                                                   181 
178     G4double GetNextWatchedTime() const;       << 182   G4double GetNextWatchedTime() const;
179                                                   183 
180     inline void SetMaxTimeStep(G4double maxTim << 184   inline void SetMaxTimeStep(G4double maxTimeStep)
                                                   >> 185   {
                                                   >> 186     fMaxTimeStep = maxTimeStep;
                                                   >> 187   }
181                                                   188 
182     inline G4double GetMaxTimeStep() const { r << 189   inline G4double GetMaxTimeStep() const
                                                   >> 190   {
                                                   >> 191     return fMaxTimeStep;
                                                   >> 192   }
183                                                   193 
184     inline G4VScavengerMaterial* GetScavengerM << 194   inline G4VScavengerMaterial* GetScavengerMaterial() const
185     inline void SetScavengerMaterial(std::uniq << 195   {
186     {                                          << 196       return fpUserScavenger.get();
                                                   >> 197   }
                                                   >> 198   inline void SetScavengerMaterial(std::unique_ptr<G4VScavengerMaterial> scavengerMaterial)
                                                   >> 199   {
187       fpUserScavenger = std::move(scavengerMat    200       fpUserScavenger = std::move(scavengerMaterial);
188     }                                          << 201   }
189                                                   202 
190   protected:                                   << 203 protected:
191     void DoProcess();                          << 
192     void SynchronizeTracks();                  << 
193     void Stepping();                           << 
194                                                << 
195     void FindUserPreDefinedTimeStep();         << 
196                                                << 
197     G4bool CanICarryOn();                      << 
198                                                << 
199     void PrintWhyDoYouStop();                  << 
200                                                << 
201   private:                                     << 
202     G4Scheduler();                             << 
203     void Create();                             << 
204                                                << 
205     G4SchedulerMessenger* fpMessenger = nullpt << 
206                                                << 
207     static G4ThreadLocal G4Scheduler* fgSchedu << 
208     G4int fVerbose;                            << 
209     G4bool fWhyDoYouStop;                      << 
210     G4bool fInitialized;                       << 
211     G4bool fRunning;                           << 
212     G4bool fContinue;                          << 
213                                                << 
214     G4int fNbSteps;                            << 
215     G4int fMaxSteps;                           << 
216                                                << 
217     G4ITStepStatus fITStepStatus;              << 
218                                                << 
219     // Time members                            << 
220     G4bool fUseDefaultTimeSteps;               << 
221     G4double fTimeTolerance;                   << 
222     G4double fGlobalTime;                      << 
223     G4double fStartTime;                       << 
224     G4double fStopTime;                        << 
225     G4double fEndTime;                         << 
226     G4double fPreviousTimeStep;                << 
227     G4int fZeroTimeCount;                      << 
228     G4int fMaxNZeroTimeStepsAllowed;           << 
229                                                << 
230     G4double fTimeStep;  // The selected minim << 
231     G4double fMaxTimeStep;                     << 
232                                                << 
233     // User steps                              << 
234     G4bool fUsePreDefinedTimeSteps;            << 
235     G4double fDefaultMinTimeStep;              << 
236     std::map<G4double, G4double>* fpUserTimeSt << 
237     // One can give time steps in respect to t << 
238     mutable G4double fUserUpperTimeLimit;      << 
239     G4double fDefinedMinTimeStep;              << 
240     // selected user time step in respect to t << 
241     G4bool fReachedUserTimeLimit;  // if fMinT << 
242                                                << 
243     std::set<G4double> fWatchedTimes;          << 
244                                                << 
245     G4UserTimeStepAction* fpUserTimeStepAction << 
246                                                << 
247     std::unique_ptr<G4VScavengerMaterial> fpUs << 
248                                                << 
249     // ======================================= << 
250     // TO BE REMOVED                           << 
251     G4ITStepProcessor* fpStepProcessor = nullp << 
252     G4ITModelProcessor* fpModelProcessor = nul << 
253     G4ITTrackingManager* fpTrackingManager = n << 
254     G4ITTrackingInteractivity* fpTrackingInter << 
255     G4ITReactionSet* fReactionSet = nullptr;   << 
256     G4ITTrackHolder& fTrackContainer;          << 
257     G4ITModelHandler* fpModelHandler = nullptr << 
258     // ======================================= << 
259                                                << 
260     G4double fTSTimeStep;                      << 
261     // Time calculated by the time stepper in  << 
262     G4double fILTimeStep;                      << 
263     // Time calculated by the interaction leng << 
264     // in ComputeInteractionLength()           << 
265                                                << 
266     G4bool fInteractionStep;                   << 
267     // Flag : if the step is driven by the int << 
268     // NOT by the reaction between tracks      << 
269                                                << 
270     G4ITGun* fpGun;                            << 
271                                                << 
272     // ======================================= << 
273     // Hoang                                   << 
274     G4bool fResetScavenger;                    << 
275                                                   204 
276   public:                                      << 205   void DoProcess();
277     void ResetScavenger(bool);                 << 206   void SynchronizeTracks();
                                                   >> 207   void Stepping();
                                                   >> 208 
                                                   >> 209   void FindUserPreDefinedTimeStep();
                                                   >> 210 
                                                   >> 211   G4bool CanICarryOn();
                                                   >> 212 
                                                   >> 213   void PrintWhyDoYouStop();
                                                   >> 214 
                                                   >> 215 private:
                                                   >> 216   G4Scheduler();
                                                   >> 217   void Create();
                                                   >> 218   G4Scheduler(const G4Scheduler&);
                                                   >> 219   G4Scheduler& operator=(const G4Scheduler&);
                                                   >> 220 
                                                   >> 221   G4SchedulerMessenger* fpMessenger;
                                                   >> 222 
                                                   >> 223   static G4ThreadLocal G4Scheduler* fgScheduler;
                                                   >> 224   G4int fVerbose;
                                                   >> 225   G4bool fWhyDoYouStop;
                                                   >> 226   G4bool fInitialized;
                                                   >> 227   G4bool fRunning;
                                                   >> 228   G4bool fContinue;
                                                   >> 229 
                                                   >> 230   G4int fNbSteps;
                                                   >> 231   G4int fMaxSteps;
                                                   >> 232 
                                                   >> 233   G4ITStepStatus fITStepStatus;
                                                   >> 234 
                                                   >> 235   // Time members
                                                   >> 236   G4bool fUseDefaultTimeSteps;
                                                   >> 237   G4double fTimeTolerance;
                                                   >> 238   G4double fGlobalTime;
                                                   >> 239   G4double fTmpGlobalTime;
                                                   >> 240   G4double fStartTime;
                                                   >> 241   G4double fStopTime;
                                                   >> 242   G4double fEndTime;
                                                   >> 243   G4double fPreviousTimeStep;
                                                   >> 244   G4int fZeroTimeCount;
                                                   >> 245   G4int fMaxNZeroTimeStepsAllowed;
                                                   >> 246 
                                                   >> 247   G4double fTimeStep; // The selected minimum time step
                                                   >> 248   G4double fMaxTimeStep;
                                                   >> 249 
                                                   >> 250   // User steps
                                                   >> 251   G4bool fUsePreDefinedTimeSteps;
                                                   >> 252   G4double fDefaultMinTimeStep;
                                                   >> 253   std::map<G4double, G4double>* fpUserTimeSteps;
                                                   >> 254   // One can give time steps in respect to the global time
                                                   >> 255   mutable G4double fUserUpperTimeLimit;
                                                   >> 256   G4double fDefinedMinTimeStep;
                                                   >> 257   // selected user time step in respect to the global time
                                                   >> 258   G4bool fReachedUserTimeLimit; // if fMinTimeStep == the user time step
                                                   >> 259 
                                                   >> 260   std::set<G4double> fWatchedTimes;
                                                   >> 261 
                                                   >> 262   G4UserTimeStepAction* fpUserTimeStepAction;
                                                   >> 263 
                                                   >> 264   std::unique_ptr<G4VScavengerMaterial> fpUserScavenger;
                                                   >> 265 
                                                   >> 266   // ==========================================
                                                   >> 267   // TO BE REMOVED
                                                   >> 268   G4ITStepProcessor* fpStepProcessor;
                                                   >> 269   G4ITModelProcessor* fpModelProcessor;
                                                   >> 270   G4ITTrackingManager* fpTrackingManager;
                                                   >> 271   G4ITTrackingInteractivity* fpTrackingInteractivity;
                                                   >> 272   G4ITReactionSet* fReactionSet;
                                                   >> 273   G4ITTrackHolder& fTrackContainer;
                                                   >> 274   G4ITModelHandler* fpModelHandler;
                                                   >> 275   // ==========================================
                                                   >> 276 
                                                   >> 277   G4double fTSTimeStep;
                                                   >> 278   // Time calculated by the time stepper in CalculateMinTimeStep()
                                                   >> 279   G4double fILTimeStep;
                                                   >> 280   // Time calculated by the interaction length methods
                                                   >> 281   // in ComputeInteractionLength()
                                                   >> 282 
                                                   >> 283   G4bool fInteractionStep;
                                                   >> 284   // Flag : if the step is driven by the interaction with the matter and
                                                   >> 285   // NOT by the reaction between tracks
                                                   >> 286 
                                                   >> 287   G4ITGun* fpGun;
                                                   >> 288 
                                                   >> 289   // ==========================================
                                                   >> 290   //Hoang
                                                   >> 291   bool fResetScavenger;
                                                   >> 292 public:
                                                   >> 293   void ResetScavenger(bool);
278 };                                                294 };
279                                                   295 
280 inline G4bool G4Scheduler::IsInitialized()        296 inline G4bool G4Scheduler::IsInitialized()
281 {                                                 297 {
282   return fInitialized;                            298   return fInitialized;
283 }                                                 299 }
284                                                   300 
285 inline G4ITModelHandler* G4Scheduler::GetModel    301 inline G4ITModelHandler* G4Scheduler::GetModelHandler()
286 {                                                 302 {
287   return fpModelHandler;                          303   return fpModelHandler;
288 }                                                 304 }
289                                                   305 
290 inline void G4Scheduler::SetEndTime(const G4do    306 inline void G4Scheduler::SetEndTime(const G4double __endtime)
291 {                                                 307 {
292   fEndTime = __endtime;                           308   fEndTime = __endtime;
293 }                                                 309 }
294                                                   310 
295 inline void G4Scheduler::SetTimeSteps(std::map << 311 inline
                                                   >> 312 void G4Scheduler::SetTimeSteps(std::map<G4double, G4double>* steps)
296 {                                                 313 {
297   fUsePreDefinedTimeSteps = true;                 314   fUsePreDefinedTimeSteps = true;
298   fpUserTimeSteps = steps;                        315   fpUserTimeSteps = steps;
299 }                                                 316 }
300                                                   317 
301 inline void G4Scheduler::AddTimeStep(G4double     318 inline void G4Scheduler::AddTimeStep(G4double startingTime, G4double timeStep)
302 {                                                 319 {
303   if (fpUserTimeSteps == nullptr) {            << 320   if (fpUserTimeSteps == 0)
                                                   >> 321   {
304     fpUserTimeSteps = new std::map<G4double, G    322     fpUserTimeSteps = new std::map<G4double, G4double>();
305     fUsePreDefinedTimeSteps = true;               323     fUsePreDefinedTimeSteps = true;
306   }                                               324   }
307                                                   325 
308   (*fpUserTimeSteps)[startingTime] = timeStep;    326   (*fpUserTimeSteps)[startingTime] = timeStep;
309 }                                                 327 }
310                                                   328 
311 inline G4int G4Scheduler::GetNbSteps() const      329 inline G4int G4Scheduler::GetNbSteps() const
312 {                                                 330 {
313   return fNbSteps;                                331   return fNbSteps;
314 }                                                 332 }
315                                                   333 
316 inline void G4Scheduler::SetMaxNbSteps(G4int m    334 inline void G4Scheduler::SetMaxNbSteps(G4int maxSteps)
317 {                                                 335 {
318   fMaxSteps = maxSteps;                           336   fMaxSteps = maxSteps;
319 }                                                 337 }
320                                                   338 
321 inline G4int G4Scheduler::GetMaxNbSteps() cons    339 inline G4int G4Scheduler::GetMaxNbSteps() const
322 {                                                 340 {
323   return fMaxSteps;                               341   return fMaxSteps;
324 }                                                 342 }
325                                                   343 
326 inline G4double G4Scheduler::GetStartTime() co    344 inline G4double G4Scheduler::GetStartTime() const
327 {                                                 345 {
328   return fStartTime;                              346   return fStartTime;
329 }                                                 347 }
330                                                   348 
331 inline G4double G4Scheduler::GetEndTime() cons    349 inline G4double G4Scheduler::GetEndTime() const
332 {                                                 350 {
333   return fEndTime;                                351   return fEndTime;
334 }                                                 352 }
335                                                   353 
336 inline G4double G4Scheduler::GetTimeStep() con    354 inline G4double G4Scheduler::GetTimeStep() const
337 {                                                 355 {
338   return fTimeStep;                               356   return fTimeStep;
339 }                                                 357 }
340                                                   358 
341 inline void G4Scheduler::SetDefaultTimeStep(G4    359 inline void G4Scheduler::SetDefaultTimeStep(G4double timeStep)
342 {                                                 360 {
343   fDefaultMinTimeStep = timeStep;                 361   fDefaultMinTimeStep = timeStep;
344 }                                                 362 }
345                                                   363 
346 inline G4double G4Scheduler::GetGlobalTime() c    364 inline G4double G4Scheduler::GetGlobalTime() const
347 {                                                 365 {
348   return fGlobalTime;                             366   return fGlobalTime;
349 }                                                 367 }
350                                                   368 
351 inline void G4Scheduler::SetUserAction(G4UserT << 369 inline
                                                   >> 370 void G4Scheduler::SetUserAction(G4UserTimeStepAction* userITAction)
352 {                                                 371 {
353   fpUserTimeStepAction = userITAction;            372   fpUserTimeStepAction = userITAction;
354 }                                                 373 }
355                                                   374 
356 inline G4UserTimeStepAction* G4Scheduler::GetU    375 inline G4UserTimeStepAction* G4Scheduler::GetUserTimeStepAction() const
357 {                                                 376 {
358   return fpUserTimeStepAction;                    377   return fpUserTimeStepAction;
359 }                                                 378 }
360                                                   379 
361 inline void G4Scheduler::SetVerbose(G4int verb    380 inline void G4Scheduler::SetVerbose(G4int verbose)
362 {                                                 381 {
363   fVerbose = verbose;                             382   fVerbose = verbose;
364 }                                                 383 }
365                                                   384 
366 inline G4int G4Scheduler::GetVerbose() const      385 inline G4int G4Scheduler::GetVerbose() const
367 {                                                 386 {
368   return fVerbose;                                387   return fVerbose;
369 }                                                 388 }
370                                                   389 
371 inline void G4Scheduler::SetMaxZeroTimeAllowed << 390 inline
                                                   >> 391 void G4Scheduler::SetMaxZeroTimeAllowed(G4int maxTimeStepAllowed)
372 {                                                 392 {
373   fMaxNZeroTimeStepsAllowed = maxTimeStepAllow    393   fMaxNZeroTimeStepsAllowed = maxTimeStepAllowed;
374 }                                                 394 }
375                                                   395 
376 inline G4int G4Scheduler::GetMaxZeroTimeAllowe    396 inline G4int G4Scheduler::GetMaxZeroTimeAllowed() const
377 {                                                 397 {
378   return fMaxNZeroTimeStepsAllowed;               398   return fMaxNZeroTimeStepsAllowed;
379 }                                                 399 }
380                                                   400 
381 inline void G4Scheduler::SetTimeTolerance(G4do    401 inline void G4Scheduler::SetTimeTolerance(G4double time)
382 {                                                 402 {
383   fTimeTolerance = time;                          403   fTimeTolerance = time;
384 }                                                 404 }
385                                                   405 
386 inline G4double G4Scheduler::GetTimeTolerance(    406 inline G4double G4Scheduler::GetTimeTolerance() const
387 {                                                 407 {
388   return fTimeTolerance;                          408   return fTimeTolerance;
389 }                                                 409 }
390                                                   410 
391 inline G4double G4Scheduler::GetPreviousTimeSt    411 inline G4double G4Scheduler::GetPreviousTimeStep() const
392 {                                                 412 {
393   return fPreviousTimeStep;                       413   return fPreviousTimeStep;
394 }                                                 414 }
395                                                   415 
396 inline G4ITStepStatus G4Scheduler::GetStatus()    416 inline G4ITStepStatus G4Scheduler::GetStatus() const
397 {                                                 417 {
398   return fITStepStatus;                           418   return fITStepStatus;
399 }                                                 419 }
400                                                   420 
401 inline void G4Scheduler::Stop()                   421 inline void G4Scheduler::Stop()
402 {                                                 422 {
403   fContinue = false;                              423   fContinue = false;
404 }                                                 424 }
405                                                   425 
406 inline G4ITTrackingInteractivity* G4Scheduler:    426 inline G4ITTrackingInteractivity* G4Scheduler::GetInteractivity()
407 {                                                 427 {
408   return fpTrackingInteractivity;                 428   return fpTrackingInteractivity;
409 }                                                 429 }
410                                                   430 
411 inline void G4Scheduler::SetGun(G4ITGun* gun)     431 inline void G4Scheduler::SetGun(G4ITGun* gun)
412 {                                                 432 {
413   fpGun = gun;                                    433   fpGun = gun;
414 }                                                 434 }
415                                                   435 
416 inline G4ITGun* G4Scheduler::GetGun()             436 inline G4ITGun* G4Scheduler::GetGun()
417 {                                                 437 {
418   return fpGun;                                   438   return fpGun;
419 }                                                 439 }
420                                                   440 
421 inline void G4Scheduler::WhyDoYouStop()           441 inline void G4Scheduler::WhyDoYouStop()
422 {                                                 442 {
423   fWhyDoYouStop = true;                           443   fWhyDoYouStop = true;
424 }                                                 444 }
425                                                   445 
426 inline void G4Scheduler::UseDefaultTimeSteps(G    446 inline void G4Scheduler::UseDefaultTimeSteps(G4bool flag)
427 {                                                 447 {
428   fUseDefaultTimeSteps = flag;                    448   fUseDefaultTimeSteps = flag;
429 }                                                 449 }
430                                                   450 
431 inline G4bool G4Scheduler::AreDefaultTimeSteps    451 inline G4bool G4Scheduler::AreDefaultTimeStepsUsed()
432 {                                                 452 {
433   return (!fUseDefaultTimeSteps && !fUsePreDef << 453   return (fUseDefaultTimeSteps == false && fUsePreDefinedTimeSteps == false);
434 }                                                 454 }
435                                                   455 
436 inline void G4Scheduler::ResetScavenger(bool v    456 inline void G4Scheduler::ResetScavenger(bool value)
437 {                                                 457 {
438   fResetScavenger = value;                     << 458     fResetScavenger = value;
439 }                                                 459 }
440                                                   460 
441 #endif                                            461 #endif
442                                                   462