| 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 //
27 #include "G4DNADamage.hh"
28 #include "G4UnitsTable.hh"
29
30 G4ThreadLocal G4DNADamage* G4DNADamage::fpInst
31
32 G4DNAIndirectHit::G4DNAIndirectHit(const G4Str
33 const G4Mol
34 const G4Thr
35 G4double ti
36 fpMolecule(molecule)
37 {
38 fBaseName = baseName;
39 fPosition = position;
40 fTime = time;
41 }
42
43 G4DNAIndirectHit::~G4DNAIndirectHit()
44 {
45 delete fpMolecule;
46 fpMolecule = nullptr;
47 }
48
49 void G4DNAIndirectHit::Print()
50 {
51 G4cout << "Reaction : " << fpMolecule->GetNa
52 << " at position : " << G4BestUnit(fP
53 << " and time : " << G4BestUnit(fTime
54 }
55
56 G4DNADamage* G4DNADamage::Instance()
57 {
58 if (fpInstance == nullptr) new G4DNADamage()
59
60 return fpInstance;
61 }
62
63 G4DNADamage::G4DNADamage()
64 {
65 fJustCountDamage = false;
66 fNIndirectDamage = 0;
67 fpInstance = this;
68 }
69
70 G4DNADamage::~G4DNADamage()
71 {
72 for (auto & fIndirectHit : fIndirectHits)
73 {
74 delete fIndirectHit;
75 }
76 fIndirectHits.clear();
77 }
78
79 void G4DNADamage::DeleteInstance()
80 {
81 delete fpInstance;
82 fpInstance = nullptr;
83 }
84
85 void G4DNADamage::Reset()
86 {
87 fNIndirectDamage = 0;
88 for (auto & fIndirectHit : fIndirectHits)
89 {
90 delete fIndirectHit;
91 }
92 fIndirectHits.clear();
93 }
94
95 void G4DNADamage::AddIndirectDamage(const G4St
96 const G4Mo
97 const G4Th
98 G4double t
99 {
100 if (fJustCountDamage)
101 {
102 fNIndirectDamage++;
103 return;
104 }
105
106 G4DNAIndirectHit* indirectHit = nullptr;
107 auto it = fMolMap.find(*molecule);
108
109 if (it == fMolMap.cend())
110 {
111 G4Molecule* mol(nullptr);
112 fMolMap[*molecule] = (mol = new G4Molecule
113 indirectHit = new G4DNAIndirectHit(baseNam
114 }
115 else
116 {
117 indirectHit = new G4DNAIndirectHit(baseNam
118 }
119 fIndirectHits.push_back(indirectHit);
120 }
121