| 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 //
27 // 070606 fix for Valgrind by T. Koi
28 //
29 // P. Arce, June-2014 Conversion neutron_hp to particle_hp
30 //
31 #ifndef G4ParticleHPFissionERelease_h
32 #define G4ParticleHPFissionERelease_h 1
33
34 #include "G4ios.hh"
35 #include "globals.hh"
36
37 #include <CLHEP/Units/SystemOfUnits.h>
38
39 #include <fstream>
40
41 class G4ParticleHPFissionERelease
42 {
43 public:
44 G4ParticleHPFissionERelease()
45
46 = default;
47 ~G4ParticleHPFissionERelease() = default;
48
49 inline void Init(std::istream& aDataFile)
50 {
51 G4double dummy;
52
53 aDataFile >> dummy >> fragmentKinetic >> promptNeutronKinetic >> delayedNeutronKinetic
54 >> promptGammaEnergy >> delayedGammaEnergy >> delayedBetaEnergy >> neutrinoEnergy
55 >> reducedTotalEnergy >> totalEnergy;
56
57 fragmentKinetic *= CLHEP::eV;
58 promptNeutronKinetic *= CLHEP::eV;
59 delayedNeutronKinetic *= CLHEP::eV;
60 promptGammaEnergy *= CLHEP::eV;
61 delayedGammaEnergy *= CLHEP::eV;
62 delayedBetaEnergy *= CLHEP::eV;
63 neutrinoEnergy *= CLHEP::eV;
64 reducedTotalEnergy *= CLHEP::eV;
65 totalEnergy *= CLHEP::eV;
66 }
67
68 inline G4double GetTotalEnergy(G4double deltaNNeu, G4double anEnergy)
69 {
70 G4double result, delta, energy;
71 energy = anEnergy / CLHEP::eV;
72 delta = -(1.057 * energy - 8.07 * deltaNNeu);
73 result = totalEnergy - delta * CLHEP::eV;
74 return result;
75 }
76 inline G4double GetFragmentKinetic() { return fragmentKinetic; }
77 inline G4double GetPromptNeutronKinetic(G4double deltaNNeu, G4double anEnergy)
78 {
79 G4double result, delta, energy;
80 energy = anEnergy / CLHEP::eV;
81 delta = -(1.307 * energy - 8.07 * deltaNNeu);
82 result = totalEnergy - delta * CLHEP::eV;
83 return result;
84 }
85 inline G4double GetDelayedNeutronKinetic() { return delayedNeutronKinetic; }
86 inline G4double GetPromptGammaEnergy() { return promptGammaEnergy; }
87 inline G4double GetDelayedGammaEnergy(G4double anEnergy)
88 {
89 G4double delta = 0.075 * anEnergy;
90 G4double result = delayedGammaEnergy - delta;
91 return result;
92 }
93 inline G4double GetDelayedBetaEnergy(G4double anEnergy)
94 {
95 G4double delta = 0.075 * anEnergy;
96 G4double result = delayedBetaEnergy - delta;
97 return result;
98 }
99 inline G4double GetNeutrinoEnergy(G4double anEnergy)
100 {
101 G4double delta = 0.1 * anEnergy;
102 G4double result = neutrinoEnergy - delta;
103 return result;
104 }
105 inline G4double GetReducedTotal(G4double deltaNNeu, G4double anEnergy)
106 {
107 return GetTotalEnergy(deltaNNeu, anEnergy) - GetNeutrinoEnergy(anEnergy);
108 }
109
110 private:
111 G4double totalEnergy{0.0};
112 G4double fragmentKinetic{0.0};
113 G4double promptNeutronKinetic{0.0};
114 G4double delayedNeutronKinetic{0.0};
115 G4double promptGammaEnergy{0.0};
116 G4double delayedGammaEnergy{0.0};
117 G4double delayedBetaEnergy{0.0};
118 G4double neutrinoEnergy{0.0};
119 G4double reducedTotalEnergy{0.0}; // total - neutrino
120 };
121
122 #endif
123