Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 // This example is provided by the Geant4-DNA 27 // Any report or published results obtained us 28 // shall cite the following Geant4-DNA collabo 29 // Med. Phys. 37 (2010) 4692-4708 30 // J. Comput. Phys. 274 (2014) 841-882 31 // The Geant4-DNA web site is available at htt 32 // 33 // $Id$ 34 // 35 /// \file TimeStepAction.hh 36 /// \brief Implementation of the TimeStepActio 37 38 #include "TimeStepAction.hh" 39 40 #include "DetectorConstruction.hh" 41 42 #include "G4AnalysisManager.hh" 43 #include "G4IT.hh" 44 #include "G4ITTrackHolder.hh" 45 #include "G4Molecule.hh" 46 #include "G4RunManager.hh" 47 #include "G4Scheduler.hh" 48 #include "G4SystemOfUnits.hh" 49 #include "G4UnitsTable.hh" 50 // using namespace G4DNAPARSER; 51 52 TimeStepAction::TimeStepAction() : G4UserTimeS 53 { 54 fpDetector = dynamic_cast<const DetectorCons 55 G4RunManager::GetRunManager()->GetUserDete 56 57 AddTimeStep(1 * picosecond, 0.1 * picosecond 58 AddTimeStep(10 * picosecond, 1 * picosecond) 59 AddTimeStep(100 * picosecond, 3 * picosecond 60 AddTimeStep(1000 * picosecond, 10 * picoseco 61 AddTimeStep(10000 * picosecond, 100 * picose 62 } 63 64 //....oooOO0OOooo........oooOO0OOooo........oo 65 66 TimeStepAction::~TimeStepAction() 67 { 68 // dtor 69 } 70 71 //....oooOO0OOooo........oooOO0OOooo........oo 72 73 TimeStepAction::TimeStepAction(const TimeStepA 74 { 75 // copy ctor 76 } 77 78 //....oooOO0OOooo........oooOO0OOooo........oo 79 80 TimeStepAction& TimeStepAction::operator=(cons 81 { 82 if (this == &rhs) return *this; 83 return *this; 84 } 85 86 //....oooOO0OOooo........oooOO0OOooo........oo 87 88 void TimeStepAction::Save(G4MolecularConfigura 89 { 90 G4int moleculeID = molconf->GetMoleculeID(); 91 const G4String& moleculeName = molconf->GetF 92 G4TrackList* trackList = G4ITTrackHolder::In 93 94 if (trackList == 0) return; 95 96 G4TrackList::iterator it = trackList->begin( 97 G4TrackList::iterator end = trackList->end() 98 99 for (; it != end; ++it) { 100 G4Track* track = *it; 101 SaveMoleculeInfo(track, moleculeID, molecu 102 } 103 } 104 105 //....oooOO0OOooo........oooOO0OOooo........oo 106 107 void TimeStepAction::SaveMoleculeInfo(G4Track* 108 { 109 G4AnalysisManager* analysisManager = G4Analy 110 if (!analysisManager->IsActive()) { 111 return; 112 } 113 114 const G4ThreeVector& position = track->GetPo 115 // H_{3}O^{1} ID = 0 116 // OH^{-1} ID = 1 117 // OH^{0} ID = 2 118 // e_{aq}^{1} ID = 3 119 // H0 ID = 4 120 // H_{2}^{0} ID = 5 121 // H2O2 ID = 6 122 // H_{2}O^{0} or H_{2}O^{1} ID=7-17 123 124 G4double xp = position.x(); 125 G4double yp = position.y(); 126 G4double zp = position.z(); 127 G4double R = std::sqrt(xp * xp + yp * yp + z 128 // G4double RNP = fpDetector->GetNPRadius( 129 130 G4int offset = 10; 131 132 if (molID < 7) { 133 analysisManager->FillH1(molID + offset, R) 134 } 135 else { 136 analysisManager->FillH1(17, R); 137 } 138 139 G4Scheduler::Instance()->Stop(); 140 } 141 //....oooOO0OOooo........oooOO0OOooo........oo 142 143 void TimeStepAction::UserPreTimeStepAction() 144 { 145 // Loop over defined molecules 146 G4ConfigurationIterator it = G4MoleculeTable 147 148 G4double time = G4Scheduler::Instance()->Get 149 150 if (time == 1. * picosecond) { 151 while (it()) { 152 G4MolecularConfiguration* molconf = it.v 153 Save(molconf); 154 } 155 } 156 } 157 158 //....oooOO0OOooo........oooOO0OOooo........oo 159 160 void TimeStepAction::UserReactionAction(const 161 const 162 {} 163