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 * G4DNAIRT.hh 28 * 29 * Created on: Jul 23, 2019 30 * Author: W. G. Shin 31 * J. Ramos-Mendez and B. Faddegon 32 */ 33 34 #ifndef G4DNAIRT_HH_ 35 #define G4DNAIRT_HH_ 36 37 38 #include "G4DNAMolecularReaction.hh" 39 #include "G4DNAMolecularReactionTable.hh" 40 #include "G4ITReaction.hh" 41 #include "G4ITReactionTable.hh" 42 #include "G4ITTrackHolder.hh" 43 #include "G4Molecule.hh" 44 #include "G4MoleculeTable.hh" 45 #include "G4ParticleChange.hh" 46 #include "G4ThreeVector.hh" 47 #include "G4VDNAReactionModel.hh" 48 #include "G4VITReactionProcess.hh" 49 #include "globals.hh" 50 51 #include "AddClone_def.hh" 52 53 #include <map> 54 #include <vector> 55 56 class G4DNAMolecularReactionTable; 57 class G4VDNAReactionModel; 58 class G4ErrorFunction; 59 60 class G4DNAIRT : public G4VITReactionProcess 61 { 62 public: 63 64 G4DNAIRT(); 65 explicit G4DNAIRT(G4VDNAReactionModel*); 66 ~G4DNAIRT() override; 67 G4DNAIRT(const G4DNAIRT& other) = delete; 68 G4DNAIRT& operator=(const G4DNAIRT& other) = delete; 69 70 G4bool TestReactibility(const G4Track&, 71 const G4Track&, 72 G4double , 73 G4bool ) override; 74 std::vector<std::unique_ptr<G4ITReactionChange>> FindReaction(G4ITReactionSet*, const G4double, const G4double, const G4bool) override; 75 std::unique_ptr<G4ITReactionChange> MakeReaction(const G4Track&, const G4Track&) override; 76 77 void SetReactionModel(G4VDNAReactionModel*); 78 79 void Initialize() override; 80 void SpaceBinning(); 81 void IRTSampling(); 82 void Sampling(G4Track*); 83 84 G4double GetIndependentReactionTime(const G4MolecularConfiguration*, const G4MolecularConfiguration*, G4double); 85 G4int FindBin(G4int, G4double, G4double, G4double); 86 G4double SamplePDC(G4double , G4double ); 87 88 protected: 89 const G4DNAMolecularReactionTable*& fMolReactionTable; 90 G4VDNAReactionModel* fpReactionModel; 91 92 private: 93 G4ITTrackHolder* fTrackHolder; 94 G4ITReactionSet* fReactionSet; 95 G4ErrorFunction* erfc; 96 97 std::map<G4int,std::map<G4int,std::map<G4int,std::vector<G4Track*>>>> spaceBinned; 98 99 G4double fRCutOff; 100 G4double timeMin; 101 G4double timeMax; 102 103 G4double fXMin, fYMin, fZMin; 104 G4double fXMax, fYMax, fZMax; 105 G4int fNx, fNy, fNz; 106 G4int xiniIndex, yiniIndex, ziniIndex; 107 G4int xendIndex, yendIndex, zendIndex; 108 }; 109 110 #endif /* G4DNAIRT_HH_ */ 111