Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer 3 // * License and Disclaimer * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. 9 // * include a list of copyright holders. * 10 // * 10 // * * 11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. * 17 // * 17 // * * 18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************* 24 // ******************************************************************** 25 // 25 // 26 // 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 << 57 #include <G4VScheduler.hh> << 58 << 59 #include <map> << 60 #include <memory> << 61 #include <vector> << 62 63 63 class G4ITTrackingManager; 64 class G4ITTrackingManager; 64 class G4ITModelProcessor; 65 class G4ITModelProcessor; 65 class G4ITStepProcessor; 66 class G4ITStepProcessor; 66 class G4Track; 67 class G4Track; 67 class G4UserTimeStepAction; 68 class G4UserTimeStepAction; 68 class G4SchedulerMessenger; 69 class G4SchedulerMessenger; 69 class G4ITTrackingInteractivity; 70 class G4ITTrackingInteractivity; 70 class G4ITGun; 71 class G4ITGun; 71 72 72 #ifndef compTrackPerID__ 73 #ifndef compTrackPerID__ 73 # define compTrackPerID__ << 74 #define compTrackPerID__ 74 struct compTrackPerID << 75 struct compTrackPerID 75 { << 76 { 76 G4bool operator()(G4Track* rhs, G4Track* l << 77 bool operator()(G4Track* rhs, G4Track* lhs) const 77 { 78 { 78 return rhs->GetTrackID() < lhs->GetTrack 79 return rhs->GetTrackID() < lhs->GetTrackID(); 79 } 80 } 80 }; << 81 }; 81 #endif 82 #endif 82 83 83 /** 84 /** 84 * G4Scheduler synchronizes (in time) track st 85 * G4Scheduler synchronizes (in time) track stepping 85 */ 86 */ 86 class G4Scheduler : public G4VScheduler, publi << 87 class G4Scheduler : 87 { << 88 public G4VScheduler, 88 protected: << 89 public G4VStateDependent 89 ~G4Scheduler() override; << 90 { 90 << 91 protected: 91 public: << 92 virtual ~G4Scheduler(); 92 G4Scheduler(const G4Scheduler&) = delete; << 93 93 G4Scheduler& operator=(const G4Scheduler&) << 94 public: 94 << 95 static G4Scheduler* Instance(); 95 static G4Scheduler* Instance(); << 96 /** DeleteInstance should be used instead 96 /** DeleteInstance should be used instead << 97 * of the destructor 97 * of the destructor << 98 */ 98 */ << 99 static void DeleteInstance(); 99 static void DeleteInstance(); << 100 virtual G4bool Notify(G4ApplicationState requestedState); 100 G4bool Notify(G4ApplicationState requested << 101 101 << 102 virtual void RegisterModel(G4VITStepModel*, double); 102 void RegisterModel(G4VITStepModel*, G4doub << 103 103 << 104 void Initialize(); 104 void Initialize() override; << 105 void ForceReinitialization(); 105 void ForceReinitialization(); << 106 inline bool IsInitialized(); 106 inline G4bool IsInitialized(); << 107 inline bool IsRunning(){return fRunning;} 107 inline G4bool IsRunning() override { retur << 108 void Reset(); 108 void Reset() override; << 109 void Process(); 109 void Process() override; << 110 void ClearList(); 110 void ClearList(); << 111 111 << 112 inline void SetGun(G4ITGun*); 112 inline void SetGun(G4ITGun*) override; << 113 inline G4ITGun* GetGun(); 113 inline G4ITGun* GetGun(); << 114 114 << 115 inline void Stop(); 115 inline void Stop(); << 116 void Clear(); 116 void Clear(); << 117 117 << 118 // To be called only in UserReactionAction::EndProcessing() 118 // To be called only in UserReactionAction << 119 // after fRunning flag has been turned off. 119 // after fRunning flag has been turned off << 120 // This is not done automatically before UserReactionAction::EndProcessing() 120 // This is not done automatically before U << 121 // is called in case one would like to access some track information 121 // is called in case one would like to acc << 122 void EndTracking(); 122 void EndTracking(); << 123 123 << 124 void SetEndTime(const double); 124 void SetEndTime(const G4double) override; << 125 125 << 126 /* Two tracks below the time tolerance are supposed to be 126 /* Two tracks below the time tolerance are << 127 * in the same time slice 127 * in the same time slice << 128 */ 128 */ << 129 inline void SetTimeTolerance(double); 129 inline void SetTimeTolerance(G4double) ove << 130 inline double GetTimeTolerance() const; 130 inline G4double GetTimeTolerance() const o << 131 131 << 132 inline void SetMaxZeroTimeAllowed(int); 132 inline void SetMaxZeroTimeAllowed(G4int) o << 133 inline int GetMaxZeroTimeAllowed() const; 133 inline G4int GetMaxZeroTimeAllowed() const << 134 134 << 135 inline G4ITModelHandler* GetModelHandler(); 135 inline G4ITModelHandler* GetModelHandler() << 136 136 << 137 inline void SetTimeSteps(std::map<double, double>*); 137 inline void SetTimeSteps(std::map<G4double << 138 inline void AddTimeStep(double, double); 138 inline void AddTimeStep(G4double, G4double << 139 inline void SetDefaultTimeStep(double); 139 inline void SetDefaultTimeStep(G4double) o << 140 double GetLimitingTimeStep() const; 140 G4double GetLimitingTimeStep() const overr << 141 inline G4int GetNbSteps() const; 141 inline G4int GetNbSteps() const override; << 142 inline void SetMaxNbSteps(G4int); 142 inline void SetMaxNbSteps(G4int) override; << 143 inline G4int GetMaxNbSteps() const; 143 inline G4int GetMaxNbSteps() const overrid << 144 inline G4double GetStartTime() const; 144 inline G4double GetStartTime() const overr << 145 inline G4double GetEndTime() const; 145 inline G4double GetEndTime() const overrid << 146 virtual inline G4double GetTimeStep() const; 146 inline G4double GetTimeStep() const overri << 147 inline G4double GetPreviousTimeStep() const; 147 inline G4double GetPreviousTimeStep() cons << 148 inline G4double GetGlobalTime() const; 148 inline G4double GetGlobalTime() const over << 149 inline void SetUserAction(G4UserTimeStepAction*); 149 inline void SetUserAction(G4UserTimeStepAc << 150 inline G4UserTimeStepAction* GetUserTimeStepAction() const; 150 inline G4UserTimeStepAction* GetUserTimeSt << 151 151 << 152 // To use with transportation only, no reactions 152 // To use with transportation only, no rea << 153 inline void UseDefaultTimeSteps(G4bool); 153 inline void UseDefaultTimeSteps(G4bool); << 154 inline G4bool AreDefaultTimeStepsUsed(); 154 inline G4bool AreDefaultTimeStepsUsed(); << 155 155 << 156 inline G4ITStepStatus GetStatus() const; 156 inline G4ITStepStatus GetStatus() const; << 157 157 << 158 /* 1 : Reaction information 158 /* 1 : Reaction information << 159 * 2 : (1) + time step information 159 * 2 : (1) + time step information << 160 * 3 : (2) + step info for individual tracks 160 * 3 : (2) + step info for individual trac << 161 * 4 : (2) + trackList processing info + pushed and killed track info 161 * 4 : (2) + trackList processing info + p << 162 */ 162 */ << 163 inline void SetVerbose(int); 163 inline void SetVerbose(G4int) override; << 164 164 << 165 inline int GetVerbose() const; 165 inline G4int GetVerbose() const; << 166 166 << 167 inline void WhyDoYouStop(); 167 inline void WhyDoYouStop(); << 168 168 << 169 void SetInteractivity(G4ITTrackingInteractivity*); 169 void SetInteractivity(G4ITTrackingInteract << 170 inline G4ITTrackingInteractivity* GetInteractivity(); 170 inline G4ITTrackingInteractivity* GetInter << 171 171 << 172 virtual size_t GetNTracks(); 172 virtual size_t GetNTracks(); << 173 >> 174 void GetCollisionType(G4String& interactionType); >> 175 >> 176 void AddWatchedTime(double time) >> 177 { >> 178 fWatchedTimes.insert(time); >> 179 } 173 180 174 void GetCollisionType(G4String& interactio << 181 double GetNextWatchedTime() const; 175 182 176 void AddWatchedTime(G4double time) { fWatc << 183 inline void SetMaxTimeStep(double maxTimeStep) >> 184 { >> 185 fMaxTimeStep = maxTimeStep; >> 186 } 177 187 178 G4double GetNextWatchedTime() const; << 188 inline double GetMaxTimeStep() const >> 189 { >> 190 return fMaxTimeStep; >> 191 } 179 192 180 inline void SetMaxTimeStep(G4double maxTim << 193 protected: 181 194 182 inline G4double GetMaxTimeStep() const { r << 195 void DoProcess(); >> 196 void SynchronizeTracks(); >> 197 void Stepping(); >> 198 >> 199 void FindUserPreDefinedTimeStep(); >> 200 >> 201 bool CanICarryOn(); >> 202 >> 203 void PrintWhyDoYouStop(); >> 204 >> 205 private: >> 206 G4Scheduler(); >> 207 void Create(); >> 208 G4Scheduler(const G4Scheduler&); >> 209 G4Scheduler& operator=(const G4Scheduler&); >> 210 >> 211 G4SchedulerMessenger* fpMessenger; >> 212 >> 213 static G4ThreadLocal G4Scheduler* fgScheduler; >> 214 int fVerbose; >> 215 bool fWhyDoYouStop; >> 216 bool fInitialized; >> 217 bool fRunning; >> 218 G4bool fContinue; >> 219 >> 220 int fNbSteps; >> 221 int fMaxSteps; >> 222 >> 223 G4ITStepStatus fITStepStatus; >> 224 >> 225 // Time members >> 226 G4bool fUseDefaultTimeSteps; >> 227 double fTimeTolerance; >> 228 double fGlobalTime; >> 229 double fTmpGlobalTime; >> 230 double fStartTime; >> 231 double fStopTime; >> 232 double fEndTime; >> 233 double fPreviousTimeStep; >> 234 int fZeroTimeCount; >> 235 int fMaxNZeroTimeStepsAllowed; >> 236 >> 237 double fTimeStep; // The selected minimum time step >> 238 double fMaxTimeStep; >> 239 >> 240 // User steps >> 241 bool fUsePreDefinedTimeSteps; >> 242 double fDefaultMinTimeStep; >> 243 std::map<double, double>* fpUserTimeSteps; >> 244 // One can give time steps in respect to the global time >> 245 mutable double fUserUpperTimeLimit; >> 246 double fDefinedMinTimeStep; >> 247 // selected user time step in respect to the global time >> 248 bool fReachedUserTimeLimit; // if fMinTimeStep == the user time step >> 249 >> 250 std::set<double> fWatchedTimes; >> 251 >> 252 G4UserTimeStepAction* fpUserTimeStepAction; >> 253 >> 254 // ========================================== >> 255 // TO BE REMOVED >> 256 G4ITStepProcessor* fpStepProcessor; >> 257 G4ITModelProcessor* fpModelProcessor; >> 258 G4ITTrackingManager* fpTrackingManager; >> 259 G4ITTrackingInteractivity* fpTrackingInteractivity; >> 260 G4ITReactionSet* fReactionSet; >> 261 G4ITTrackHolder& fTrackContainer; >> 262 G4ITModelHandler* fpModelHandler; >> 263 // ========================================== >> 264 >> 265 double fTSTimeStep; >> 266 // Time calculated by the time stepper in CalculateMinTimeStep() >> 267 double fILTimeStep; >> 268 // Time calculated by the interaction length methods >> 269 // in ComputeInteractionLength() >> 270 >> 271 bool fInteractionStep; >> 272 // Flag : if the step is driven by the interaction with the matter and >> 273 // NOT by the reaction between tracks 183 274 184 inline G4VScavengerMaterial* GetScavengerM << 275 G4ITGun* fpGun; 185 inline void SetScavengerMaterial(std::uniq << 186 { << 187 fpUserScavenger = std::move(scavengerMat << 188 } << 189 276 190 protected: << 277 // ========================================== 191 void DoProcess(); << 278 // 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 279 276 public: << 277 void ResetScavenger(bool); << 278 }; 280 }; 279 281 280 inline G4bool G4Scheduler::IsInitialized() << 282 inline bool G4Scheduler::IsInitialized() 281 { 283 { 282 return fInitialized; 284 return fInitialized; 283 } 285 } 284 286 285 inline G4ITModelHandler* G4Scheduler::GetModel 287 inline G4ITModelHandler* G4Scheduler::GetModelHandler() 286 { 288 { 287 return fpModelHandler; 289 return fpModelHandler; 288 } 290 } 289 291 290 inline void G4Scheduler::SetEndTime(const G4do << 292 inline void G4Scheduler::SetEndTime(const double __endtime) 291 { 293 { 292 fEndTime = __endtime; 294 fEndTime = __endtime; 293 } 295 } 294 296 295 inline void G4Scheduler::SetTimeSteps(std::map << 297 inline >> 298 void G4Scheduler::SetTimeSteps(std::map<double, double>* steps) 296 { 299 { 297 fUsePreDefinedTimeSteps = true; 300 fUsePreDefinedTimeSteps = true; 298 fpUserTimeSteps = steps; 301 fpUserTimeSteps = steps; 299 } 302 } 300 303 301 inline void G4Scheduler::AddTimeStep(G4double << 304 inline void G4Scheduler::AddTimeStep(double startingTime, double timeStep) 302 { 305 { 303 if (fpUserTimeSteps == nullptr) { << 306 if (fpUserTimeSteps == 0) 304 fpUserTimeSteps = new std::map<G4double, G << 307 { >> 308 fpUserTimeSteps = new std::map<double, double>(); 305 fUsePreDefinedTimeSteps = true; 309 fUsePreDefinedTimeSteps = true; 306 } 310 } 307 311 308 (*fpUserTimeSteps)[startingTime] = timeStep; 312 (*fpUserTimeSteps)[startingTime] = timeStep; 309 } 313 } 310 314 311 inline G4int G4Scheduler::GetNbSteps() const 315 inline G4int G4Scheduler::GetNbSteps() const 312 { 316 { 313 return fNbSteps; 317 return fNbSteps; 314 } 318 } 315 319 316 inline void G4Scheduler::SetMaxNbSteps(G4int m 320 inline void G4Scheduler::SetMaxNbSteps(G4int maxSteps) 317 { 321 { 318 fMaxSteps = maxSteps; 322 fMaxSteps = maxSteps; 319 } 323 } 320 324 321 inline G4int G4Scheduler::GetMaxNbSteps() cons 325 inline G4int G4Scheduler::GetMaxNbSteps() const 322 { 326 { 323 return fMaxSteps; 327 return fMaxSteps; 324 } 328 } 325 329 326 inline G4double G4Scheduler::GetStartTime() co 330 inline G4double G4Scheduler::GetStartTime() const 327 { 331 { 328 return fStartTime; 332 return fStartTime; 329 } 333 } 330 334 331 inline G4double G4Scheduler::GetEndTime() cons 335 inline G4double G4Scheduler::GetEndTime() const 332 { 336 { 333 return fEndTime; 337 return fEndTime; 334 } 338 } 335 339 336 inline G4double G4Scheduler::GetTimeStep() con 340 inline G4double G4Scheduler::GetTimeStep() const 337 { 341 { 338 return fTimeStep; 342 return fTimeStep; 339 } 343 } 340 344 341 inline void G4Scheduler::SetDefaultTimeStep(G4 << 345 inline void G4Scheduler::SetDefaultTimeStep(double timeStep) 342 { 346 { 343 fDefaultMinTimeStep = timeStep; 347 fDefaultMinTimeStep = timeStep; 344 } 348 } 345 349 346 inline G4double G4Scheduler::GetGlobalTime() c 350 inline G4double G4Scheduler::GetGlobalTime() const 347 { 351 { 348 return fGlobalTime; 352 return fGlobalTime; 349 } 353 } 350 354 351 inline void G4Scheduler::SetUserAction(G4UserT << 355 inline >> 356 void G4Scheduler::SetUserAction(G4UserTimeStepAction* userITAction) 352 { 357 { 353 fpUserTimeStepAction = userITAction; 358 fpUserTimeStepAction = userITAction; 354 } 359 } 355 360 356 inline G4UserTimeStepAction* G4Scheduler::GetU 361 inline G4UserTimeStepAction* G4Scheduler::GetUserTimeStepAction() const 357 { 362 { 358 return fpUserTimeStepAction; 363 return fpUserTimeStepAction; 359 } 364 } 360 365 361 inline void G4Scheduler::SetVerbose(G4int verb << 366 inline void G4Scheduler::SetVerbose(int verbose) 362 { 367 { 363 fVerbose = verbose; 368 fVerbose = verbose; 364 } 369 } 365 370 366 inline G4int G4Scheduler::GetVerbose() const << 371 inline int G4Scheduler::GetVerbose() const 367 { 372 { 368 return fVerbose; 373 return fVerbose; 369 } 374 } 370 375 371 inline void G4Scheduler::SetMaxZeroTimeAllowed << 376 inline >> 377 void G4Scheduler::SetMaxZeroTimeAllowed(int maxTimeStepAllowed) 372 { 378 { 373 fMaxNZeroTimeStepsAllowed = maxTimeStepAllow 379 fMaxNZeroTimeStepsAllowed = maxTimeStepAllowed; 374 } 380 } 375 381 376 inline G4int G4Scheduler::GetMaxZeroTimeAllowe << 382 inline int G4Scheduler::GetMaxZeroTimeAllowed() const 377 { 383 { 378 return fMaxNZeroTimeStepsAllowed; 384 return fMaxNZeroTimeStepsAllowed; 379 } 385 } 380 386 381 inline void G4Scheduler::SetTimeTolerance(G4do << 387 inline void G4Scheduler::SetTimeTolerance(double time) 382 { 388 { 383 fTimeTolerance = time; 389 fTimeTolerance = time; 384 } 390 } 385 391 386 inline G4double G4Scheduler::GetTimeTolerance( << 392 inline double G4Scheduler::GetTimeTolerance() const 387 { 393 { 388 return fTimeTolerance; 394 return fTimeTolerance; 389 } 395 } 390 396 391 inline G4double G4Scheduler::GetPreviousTimeSt 397 inline G4double G4Scheduler::GetPreviousTimeStep() const 392 { 398 { 393 return fPreviousTimeStep; 399 return fPreviousTimeStep; 394 } 400 } 395 401 396 inline G4ITStepStatus G4Scheduler::GetStatus() 402 inline G4ITStepStatus G4Scheduler::GetStatus() const 397 { 403 { 398 return fITStepStatus; 404 return fITStepStatus; 399 } 405 } 400 406 401 inline void G4Scheduler::Stop() 407 inline void G4Scheduler::Stop() 402 { 408 { 403 fContinue = false; 409 fContinue = false; 404 } 410 } 405 411 406 inline G4ITTrackingInteractivity* G4Scheduler: 412 inline G4ITTrackingInteractivity* G4Scheduler::GetInteractivity() 407 { 413 { 408 return fpTrackingInteractivity; 414 return fpTrackingInteractivity; 409 } 415 } 410 416 411 inline void G4Scheduler::SetGun(G4ITGun* gun) 417 inline void G4Scheduler::SetGun(G4ITGun* gun) 412 { 418 { 413 fpGun = gun; 419 fpGun = gun; 414 } 420 } 415 421 416 inline G4ITGun* G4Scheduler::GetGun() 422 inline G4ITGun* G4Scheduler::GetGun() 417 { 423 { 418 return fpGun; 424 return fpGun; 419 } 425 } 420 426 421 inline void G4Scheduler::WhyDoYouStop() 427 inline void G4Scheduler::WhyDoYouStop() 422 { 428 { 423 fWhyDoYouStop = true; 429 fWhyDoYouStop = true; 424 } 430 } 425 431 426 inline void G4Scheduler::UseDefaultTimeSteps(G 432 inline void G4Scheduler::UseDefaultTimeSteps(G4bool flag) 427 { 433 { 428 fUseDefaultTimeSteps = flag; 434 fUseDefaultTimeSteps = flag; 429 } 435 } 430 436 431 inline G4bool G4Scheduler::AreDefaultTimeSteps 437 inline G4bool G4Scheduler::AreDefaultTimeStepsUsed() 432 { 438 { 433 return (!fUseDefaultTimeSteps && !fUsePreDef << 439 return (fUseDefaultTimeSteps == false && fUsePreDefinedTimeSteps == false); 434 } << 435 << 436 inline void G4Scheduler::ResetScavenger(bool v << 437 { << 438 fResetScavenger = value; << 439 } 440 } 440 441 441 #endif 442 #endif 442 443