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