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 /// \file TimeStepAction.hh 28 /// \brief Implementation of the TimeStepAction class 29 30 #include "TimeStepAction.hh" 31 32 #include "G4AnalysisManager.hh" 33 #include "G4DNAMolecule.hh" 34 #include "G4Event.hh" 35 #include "G4EventManager.hh" 36 #include "G4Molecule.hh" 37 #include "G4MoleculeTable.hh" 38 #include "G4OH.hh" 39 #include "G4SystemOfUnits.hh" 40 #include "G4UnitsTable.hh" 41 42 #include <G4Scheduler.hh> 43 44 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 45 46 TimeStepAction::TimeStepAction() : G4UserTimeStepAction() 47 { 48 AddTimeStep(1 * picosecond, 0.35 * picosecond); 49 AddTimeStep(10 * picosecond, 1 * picosecond); 50 AddTimeStep(100 * picosecond, 3 * picosecond); 51 AddTimeStep(1000 * picosecond, 10 * picosecond); 52 AddTimeStep(10000 * picosecond, 100 * picosecond); 53 } 54 55 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 56 57 TimeStepAction::~TimeStepAction() {} 58 59 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 60 61 TimeStepAction::TimeStepAction(const TimeStepAction& other) : G4UserTimeStepAction(other) {} 62 63 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 64 65 TimeStepAction& TimeStepAction::operator=(const TimeStepAction& rhs) 66 { 67 if (this == &rhs) return *this; 68 return *this; 69 } 70 71 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 72 73 void TimeStepAction::UserReactionAction(const G4Track& trackA, const G4Track& trackB, 74 const std::vector<G4Track*>* pProducts) 75 { 76 if (pProducts 77 && (GetMolecule((*pProducts)[0])->GetDefinition() != G4DamagedDeoxyribose::Definition())) 78 { 79 return; 80 } 81 82 // check for the case "no product" 83 if (!pProducts) { 84 return; 85 } 86 87 auto phyEventId = G4EventManager::GetEventManager()->GetConstCurrentEvent()->GetEventID(); 88 89 const G4Track* DNAElement = nullptr; 90 const G4Track* radical = nullptr; 91 if (GetMolecule(&trackA)->GetDefinition() == G4Deoxyribose::Definition()) { 92 DNAElement = &trackA; 93 radical = &trackB; 94 } 95 else { 96 DNAElement = &trackB; 97 radical = &trackA; 98 } 99 100 if (GetMolecule(radical)->GetDefinition() != G4OH::Definition()) { 101 return; 102 } 103 104 G4AnalysisManager* analysisManager = G4AnalysisManager::Instance(); 105 analysisManager->FillNtupleDColumn(2, 0, DNAElement->GetPosition().getX() / nm); 106 analysisManager->FillNtupleDColumn(2, 1, DNAElement->GetPosition().getY() / nm); 107 analysisManager->FillNtupleDColumn(2, 2, DNAElement->GetPosition().getZ() / nm); 108 analysisManager->FillNtupleSColumn(2, 3, GetMolecule(radical)->GetName()); 109 analysisManager->FillNtupleIColumn(2, 4, G4int(phyEventId)); 110 analysisManager->AddNtupleRow(2); 111 } 112