| 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 #include "G4DNAPTBExcitationStructure.hh"
27
28 #include "G4Material.hh"
29 #include "G4SystemOfUnits.hh"
30
31 G4DNAPTBExcitationStructure::G4DNAPTBExcitatio
32 {
33 fpN2 = G4Material::GetMaterial("N2", false);
34
35 // taken directly from PTra code by MPietrz
36 if (fpN2 != nullptr) {
37 auto index = fpN2->GetIndex();
38 energyConstant[index].push_back(1.85 * eV)
39 energyConstant[index].push_back(2.15 * eV)
40 energyConstant[index].push_back(8.00 * eV)
41 energyConstant[index].push_back(8.50 * eV)
42 energyConstant[index].push_back(8.60 * eV)
43 energyConstant[index].push_back(11.05 * eV
44 energyConstant[index].push_back(11.79 * eV
45 energyConstant[index].push_back(11.90 * eV
46 energyConstant[index].push_back(12.25 * eV
47 energyConstant[index].push_back(12.50 * eV
48 energyConstant[index].push_back(13.01 * eV
49 energyConstant[index].push_back(13.19 * eV
50 energyConstant[index].push_back(13.30 * eV
51 energyConstant[index].push_back(14.33 * eV
52 energyConstant[index].push_back(14.84 * eV
53 energyConstant[index].push_back(15.18 * eV
54 energyConstant[index].push_back(15.70 * eV
55 energyConstant[index].push_back(15.75 * eV
56 energyConstant[index].push_back(15.86 * eV
57 energyConstant[index].push_back(17.36 * eV
58 energyConstant[index].push_back(17.95 * eV
59 energyConstant[index].push_back(19.77 * eV
60 energyConstant[index].push_back(20.79 * eV
61 energyConstant[index].push_back(20.87 * eV
62 energyConstant[index].push_back(22.27 * eV
63 energyConstant[index].push_back(22.83 * eV
64 energyConstant[index].push_back(37.19 * eV
65 energyConstant[index].push_back(38.67 * eV
66 energyConstant[index].push_back(39.23 * eV
67 }
68
69 for (const auto& [index, levels] : energyCon
70 nExcLevels[index] = (G4int)levels.size();
71 }
72 }
73
74 G4double G4DNAPTBExcitationStructure::Excitati
75 const G4int& ExcLevel, const size_t& materia
76 {
77 size_t matNameModif = ReplaceMaterial(materi
78
79 // check if the material exist in the map
80 if (energyConstant.find(matNameModif) == ene
81 std::ostringstream oss;
82 oss << "Material name was not found in ene
83 << matNameModif;
84 G4Exception(
85 "G4DNAPTBExcitationStructure::Excitation
86 }
87
88 G4double excitation = 0.;
89
90 if (ExcLevel >= 0 && ExcLevel < nExcLevels[m
91 excitation = energyConstant[matNameModif][
92
93 return excitation;
94 }
95
96 G4int G4DNAPTBExcitationStructure::NumberOfExc
97 {
98 auto matNameModif = ReplaceMaterial(matID);
99
100 // check if the material exist in the map
101 if (nExcLevels.find(matNameModif) == nExcLev
102 std::ostringstream oss;
103 oss << "Material name was not found in ene
104 << matNameModif;
105 G4Exception("G4DNAPTBNDExcitationStructure
106 oss.str().c_str());
107 }
108
109 return nExcLevels[matNameModif];
110 }
111
112 size_t G4DNAPTBExcitationStructure::ReplaceMat
113 {
114 auto output = materialID;
115 auto G4_N2 = G4Material::GetMaterial("G4_N2"
116 if (materialID == G4_N2) {
117 output = fpN2->GetIndex();
118 }
119
120 return output;
121 }
122