Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 // 26 // 27 #ifndef G4DNAEventScheduler_hh 28 #define G4DNAEventScheduler_hh 1 29 30 #include <G4VScheduler.hh> 31 #include <vector> 32 #include <map> 33 #include <memory> 34 #include "globals.hh" 35 #include "G4ITModelHandler.hh" 36 #include "G4ITStepStatus.hh" 37 #include "G4ITTrackHolder.hh" 38 #include "G4VStateDependent.hh" 39 #include "G4ITReaction.hh" 40 #include "G4DNAEventSet.hh" 41 #include "G4DNAUpdateSystemModel.hh" 42 #include "G4H2O2.hh" 43 #include <CLHEP/Units/SystemOfUnits.h> 44 class G4VITStepModel; 45 class G4DNAGillespieDirectMethod; 46 class G4UserMeshAction; 47 class IEventScheduler 48 { 49 public: 50 IEventScheduler() = default; 51 virtual ~IEventScheduler() = default; 52 }; 53 54 class G4DNAEventScheduler : public IEventScheduler 55 { 56 public: 57 using MolType = const G4MolecularConfiguration*; 58 using MapList = std::map<MolType, size_t>; 59 using MapCounter = std::map<MolType, G4int>; 60 G4DNAEventScheduler(); 61 ~G4DNAEventScheduler() override; 62 G4DNAEventScheduler(const G4DNAEventScheduler&) = delete; 63 G4DNAEventScheduler& operator=(const G4DNAEventScheduler& right) = delete; 64 void Initialize(const G4DNABoundingBox& boundingBox, G4int pixel); 65 void InitializeInMesh(); 66 void Voxelizing(); 67 void ReVoxelizing(G4int); 68 void SetEndTime(const G4double&); 69 G4double GetStartTime() const; 70 G4double GetEndTime() const; 71 G4double GetTimeStep() const; 72 void SetStartTime(G4double time); 73 74 inline void SetVerbose(G4int verbose) { fVerbose = verbose; } 75 inline G4int GetVerbose() const; 76 void Stepping(); 77 void SetChangeMesh(G4bool change) { fSetChangeMesh = change; } 78 79 void Reset(); 80 void ResetInMesh(); 81 void RunInMesh(); 82 void Run(); 83 84 void AddTimeToRecord(const G4double& time); 85 86 void RecordTime(); 87 void ClearAndReChargeCounter(); 88 void PrintRecordTime(); 89 void Stop(); 90 void SetMaxNbSteps(G4int); 91 std::map<G4double /*time*/, MapCounter> GetCounterMap() const; 92 G4DNAMesh* GetMesh() const; 93 G4int GetPixels() const; 94 void SetUserMeshAction(std::unique_ptr<G4UserMeshAction>); 95 static G4bool CheckingReactionRadius(G4double resolution); 96 97 private: 98 G4int fVerbose = 0; 99 G4bool fInitialized = false; 100 G4double fStartTime = 1 * CLHEP::picosecond; 101 G4double fEndTime = 10000 * CLHEP::second; 102 G4int fStepNumber = 0; 103 G4int fMaxStep = INT_MAX; 104 G4bool fRunning = true; 105 G4double fTimeStep = DBL_MAX; 106 G4double fGlobalTime = 1 * CLHEP::picosecond; 107 G4double fJumpingNumber = 0; 108 G4double fReactionNumber = 0; 109 G4int fPixel = 0; 110 G4bool fIsChangeMesh = false; 111 G4bool fSetChangeMesh = true; 112 G4int fStepNumberInMesh = 0; 113 G4double fTransferTime = 0.; 114 const G4double C = 20; 115 const G4double D = G4H2O2::Definition()->GetDiffusionCoefficient(); // this is the biggest D 116 117 std::unique_ptr<G4DNAMesh> fpMesh; 118 std::unique_ptr<G4DNAGillespieDirectMethod> fpGillespieReaction; 119 std::unique_ptr<G4DNAEventSet> fpEventSet; 120 std::unique_ptr<G4DNAUpdateSystemModel> fpUpdateSystem; 121 std::unique_ptr<G4UserMeshAction> fpUserMeshAction; 122 // an aternative Counter 123 std::map<G4double /*time*/, MapCounter> fCounterMap; 124 std::set<G4double> fTimeToRecord; 125 std::set<G4double>::iterator fLastRecoredTime; 126 void ResetEventSet(); 127 void LastRegisterForCounter(); 128 }; 129 #endif 130