| 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 // neutron_hp -- source file
27 // J.P. Wellisch, Nov-1996
28 // A prototype of the low energy neutron trans
29 //
30 // 080801 Prohibit level transition to oneself
31 //
32 // P. Arce, June-2014 Conversion neutron_hp to
33 //
34 #include "G4ParticleHPDeExGammas.hh"
35
36 #include "G4RandomDirection.hh"
37 #include "G4ReactionProduct.hh"
38 #include "G4SystemOfUnits.hh"
39
40 G4ParticleHPDeExGammas::G4ParticleHPDeExGammas
41
42 G4ParticleHPDeExGammas::~G4ParticleHPDeExGamma
43 {
44 for (auto& ptr : theLevels) {
45 delete ptr;
46 }
47 }
48
49 void G4ParticleHPDeExGammas::Init(std::istream
50 {
51 // G4cout << "### G4ParticleHPDeExGammas::In
52 // ground state
53 auto level = new G4ParticleHPNucLevel(0.0);
54
55 G4double elevel0 = 0.0;
56 G4double elevel = 0.0;
57 G4double egamma = 0.0;
58 G4double prob = 0.0;
59 constexpr G4double eps = 1 * CLHEP::eV;
60 for (;;) {
61 if (aDataFile >> elevel) {
62 // next line
63 aDataFile >> egamma >> prob;
64 egamma *= CLHEP::keV;
65 elevel *= CLHEP::keV;
66 prob = std::max(prob, 1.e-6);
67 // G4cout << " El0=" << elevel0 << "
68 // << " Eg=" << egamma << " w=" << prob
69
70 // save previous level and start a new l
71 if (std::abs(elevel - elevel0) > eps) {
72 level->Normalize();
73 theLevels.push_back(level);
74 ++nLevels;
75 level = new G4ParticleHPNucLevel(eleve
76 elevel0 = elevel;
77 // G4cout << " New level " << nLevels
78 }
79
80 // find the next level
81 G4double e = elevel - egamma;
82 G4int next = -1;
83 G4double del = DBL_MAX;
84 for (G4int i = 0; i < nLevels; ++i) {
85 G4double de = std::abs(theLevels[i]->G
86 if (de < del) {
87 next = i;
88 del = de;
89 }
90 }
91 // save level data
92 if (next >= 0) {
93 level->AddGamma(egamma, prob, next);
94 // G4cout << " NLevel=" << nLevel
95 // << " Egamma=" << egamma << " W
96 }
97 }
98 else {
99 // end of file - save recent level
100 level->Normalize();
101 theLevels.push_back(level);
102 ++nLevels;
103 // G4cout << "### End of file Nlevels="
104 break;
105 }
106 }
107 }
108
109 G4ReactionProductVector* G4ParticleHPDeExGamma
110 {
111 G4int idx = i;
112 if (idx >= nLevels || idx <= 0) return nullp
113 auto result = new G4ReactionProductVector();
114
115 for (;;) {
116 if (idx <= 0) {
117 break;
118 }
119 auto ptr = theLevels[idx]->GetDecayGamma(i
120 if (nullptr != ptr) {
121 result->push_back(ptr);
122 }
123 }
124 return result;
125 }
126