| 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 // P. Arce, June-2014 Conversion neutron_hp to particle_hp
28 // V. Ivanchenko, July-2023 Basic revision of particle HP classes
29 //
30 #ifndef G4ParticleHPProduct_h
31 #define G4ParticleHPProduct_h 1
32
33 #include "G4Cache.hh"
34 #include "G4ReactionProductVector.hh"
35 #include "G4VParticleHPEnergyAngular.hh"
36 #include "G4ParticleHPVector.hh"
37 #include "G4ios.hh"
38 #include "globals.hh"
39
40 #include <fstream>
41
42 class G4ParticleDefinition;
43
44 enum G4HPMultiMethod
45 {
46 G4HPMultiPoisson,
47 G4HPMultiBetweenInts
48 };
49
50 class G4ParticleHPProduct
51 {
52 struct toBeCached
53 {
54 G4ReactionProduct* theProjectileRP{nullptr};
55 G4ReactionProduct* theTarget{nullptr};
56 G4int theCurrentMultiplicity{-1};
57 toBeCached() = default;
58 };
59
60 public:
61 G4ParticleHPProduct();
62 ~G4ParticleHPProduct();
63
64 void Init(std::istream& aDataFile, const G4ParticleDefinition* projectile);
65
66 G4int GetMultiplicity(G4double anEnergy);
67 G4ReactionProductVector* Sample(G4double anEnergy, G4int nParticles);
68
69 G4double GetMeanYield(G4double anEnergy) { return theYield.GetY(anEnergy); }
70
71 void SetProjectileRP(G4ReactionProduct* aIncidentPart)
72 {
73 fCache.Get().theProjectileRP = aIncidentPart;
74 }
75
76 void SetTarget(G4ReactionProduct* aTarget) { fCache.Get().theTarget = aTarget; }
77
78 inline G4ReactionProduct* GetTarget() { return fCache.Get().theTarget; }
79
80 inline G4ReactionProduct* GetProjectileRP() { return fCache.Get().theProjectileRP; }
81
82 inline G4double MeanEnergyOfThisInteraction()
83 {
84 G4double result = 0.0;
85 if (theDist != nullptr) {
86 result = theDist->MeanEnergyOfThisInteraction();
87 result *= fCache.Get().theCurrentMultiplicity;
88 }
89 return result;
90 }
91
92 inline G4double GetQValue() { return theActualStateQValue; }
93
94 // TK120515 For migration of frameFlag (MF6 LCT) = 3 in
95 // G4ParticleHPEnAngCorrelation
96 G4double GetMassCode() { return theMassCode; }
97 G4double GetMass() { return theMass; }
98
99 private:
100 G4double theMassCode{0.0};
101 G4double theMass{0.0};
102 G4double theGroundStateQValue{0.0};
103 G4double theActualStateQValue{0.0};
104 G4int theIsomerFlag{0};
105 G4int theDistLaw{-1}; // redundant
106 G4VParticleHPEnergyAngular* theDist{nullptr};
107
108 // cashed values
109 //
110 G4Cache<toBeCached> fCache;
111
112 G4HPMultiMethod theMultiplicityMethod;
113 G4ParticleHPVector theYield;
114 };
115
116 #endif
117