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