| Geant4 Cross Reference |
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 // neutron_hp -- source file
27 // J.P. Wellisch, Nov-1996
28 // A prototype of the low energy neutron transport model.
29 //
30 // 080801 Prohibit level transition to oneself by T. Koi
31 //
32 // P. Arce, June-2014 Conversion neutron_hp to particle_hp
33 //
34 #include "G4ParticleHPDeExGammas.hh"
35
36 #include "G4RandomDirection.hh"
37 #include "G4ReactionProduct.hh"
38 #include "G4SystemOfUnits.hh"
39
40 G4ParticleHPDeExGammas::G4ParticleHPDeExGammas() = default;
41
42 G4ParticleHPDeExGammas::~G4ParticleHPDeExGammas()
43 {
44 for (auto& ptr : theLevels) {
45 delete ptr;
46 }
47 }
48
49 void G4ParticleHPDeExGammas::Init(std::istream& aDataFile)
50 {
51 // G4cout << "### G4ParticleHPDeExGammas::Init new file " << G4endl;
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 << " El=" << elevel
68 // << " Eg=" << egamma << " w=" << prob << G4endl;
69
70 // save previous level and start a new level
71 if (std::abs(elevel - elevel0) > eps) {
72 level->Normalize();
73 theLevels.push_back(level);
74 ++nLevels;
75 level = new G4ParticleHPNucLevel(elevel);
76 elevel0 = elevel;
77 // G4cout << " New level " << nLevels << " E=" << elevel << G4endl;
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]->GetLevelEnergy() - e);
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=" << nLevels << " Elevel=" << elevel
95 // << " Egamma=" << egamma << " W=" << prob << " next=" << next << G4endl;
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=" << nLevels << G4endl;
104 break;
105 }
106 }
107 }
108
109 G4ReactionProductVector* G4ParticleHPDeExGammas::GetDecayGammas(G4int i) const
110 {
111 G4int idx = i;
112 if (idx >= nLevels || idx <= 0) return nullptr;
113 auto result = new G4ReactionProductVector();
114
115 for (;;) {
116 if (idx <= 0) {
117 break;
118 }
119 auto ptr = theLevels[idx]->GetDecayGamma(idx);
120 if (nullptr != ptr) {
121 result->push_back(ptr);
122 }
123 }
124 return result;
125 }
126