| 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 /// \file radiobiology/src/IonLet.cc
28 /// \brief Implementation of the RadioBio::Ion
29
30 #include "IonLet.hh"
31
32 namespace RadioBio
33 {
34
35 //....oooOO0OOooo........oooOO0OOooo........oo
36
37 IonLet::IonLet(G4int trackID, G4int PDG, G4Str
38 G4int voxNumber)
39 : fIsPrimary(trackID == 1), fPDGencoding(PDG
40 {
41 fLETDN = array_type(0.0, voxNumber);
42 fLETDD = array_type(0.0, voxNumber);
43 fLETTN = array_type(0.0, voxNumber);
44 fLETTD = array_type(0.0, voxNumber);
45
46 fLETD = array_type(0.0, voxNumber);
47 fLETT = array_type(0.0, voxNumber);
48 }
49
50 //....oooOO0OOooo........oooOO0OOooo........oo
51
52 void IonLet::Update(G4int voxel, G4double DE,
53 {
54 // ions dose LET Numerator, including second
55 fLETDN[voxel] += (DE + DEELETrons) * Lsn;
56 // ions dose LET Denominator, including seco
57 fLETDD[voxel] += DE + DEELETrons;
58 // ions track LET Numerator
59 fLETTN[voxel] += DX * Lsn;
60 // ions track LET Denominator
61 fLETTD[voxel] += DX;
62 }
63
64 //....oooOO0OOooo........oooOO0OOooo........oo
65
66 void IonLet::Merge(const IonLet* rhs)
67 {
68 // If programmed correctly, this exception s
69 if (rhs->GetPDGencoding() != fPDGencoding ||
70 G4Exception("IonLet::merge", "mergingdiffe
71 "Cannotmerge ions, probably me
72
73 fLETDN += rhs->GetLETDN();
74 fLETDD += rhs->GetLETDD();
75 fLETTN += rhs->GetLETTN();
76 fLETTD += rhs->GetLETTD();
77 }
78
79 //....oooOO0OOooo........oooOO0OOooo........oo
80
81 void IonLet::Calculate()
82 {
83 for (unsigned int v = 0; v < fLETD.size(); v
84 if (fLETDD[v] > 0.) fLETD[v] = fLETDN[v] /
85 if (fLETTD[v] > 0.) fLETT[v] = fLETTN[v] /
86 }
87 }
88
89 //....oooOO0OOooo........oooOO0OOooo........oo
90
91 } // namespace RadioBio