| 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 //
27 // -------------------------------------------
28 //
29 // Geant4 source file
30 //
31 // File name: G4ParticleHPInelasticDataPT
32 //
33 // Authors: Marek Zmeskal (CTU, Czech Tec
34 // Loic Thulliez (CEA France)
35 //
36 // Creation date: 4 June 2024
37 //
38 // Description: Class for utilization of
39 // probability tables in the
40 // for inelastic channel.
41 //
42 // Modifications:
43 //
44 // -------------------------------------------
45 //
46 //
47
48 #include "G4ParticleHPInelasticDataPT.hh"
49 #include "G4DynamicParticle.hh"
50 #include "G4ParticleDefinition.hh"
51 #include "G4Neutron.hh"
52 #include "G4Element.hh"
53 #include "G4ParticleHPManager.hh"
54 #include "G4HadronicParameters.hh"
55 #include "G4ParticleHPProbabilityTablesStore.h
56 #include "G4HadronicException.hh"
57
58
59 G4ParticleHPInelasticDataPT::G4ParticleHPInela
60 // minimum and maximum energy limit for URR
61 SetMinKinEnergy( 1.0 * CLHEP::eV );
62 SetMaxKinEnergy( 1.2 * CLHEP::MeV );
63 URRlimits = nullptr;
64 }
65
66
67 G4ParticleHPInelasticDataPT::~G4ParticleHPInel
68
69
70 G4bool G4ParticleHPInelasticDataPT::IsIsoAppli
71
72 // checks applicability for the element
73 if ( doNOTusePTforInelastic ) {
74 // do not use for njoy
75 return false;
76 }
77 if ( dp->GetDefinition() != G4Neutron::Neutr
78 return false;
79 } else {
80 std::size_t elementI = elm->GetIndex();
81 G4double eKin = dp->GetKineticEnergy();
82 if ( eKin < (*URRlimits).at(elementI).fi
83 return false;
84 } else if ( eKin > (*URRlimits).at(eleme
85 return false;
86 }
87 return true;
88 }
89 return false;
90 }
91
92
93 G4double G4ParticleHPInelasticDataPT::GetIsoCr
94
95 return G4ParticleHPProbabilityTablesStore::G
96 }
97
98
99 void G4ParticleHPInelasticDataPT::BuildPhysics
100 if ( G4HadronicParameters::Instance()->GetTy
101 SetMinKinEnergy( DBL_MAX );
102 SetMaxKinEnergy( 0.0 );
103 doNOTusePTforInelastic = true;
104 } else if ( G4HadronicParameters::Instance()
105 doNOTusePTforInelastic = false;
106 G4cout << "BuildPhysicsTable in G4Particle
107 if ( &aP != G4Neutron::Neutron() ) {
108 throw G4HadronicException( __FILE__, _
109 }
110 URRlimits = G4ParticleHPManager::GetInstan
111 if ( G4Threading::IsWorkerThread() ) {
112 // sets the overall limits of the URR, w
113 // defines URR model energy range
114 SetMinKinEnergy( (*URRlimits).back().fir
115 SetMaxKinEnergy( (*URRlimits).back().sec
116 } else {
117 if (G4ParticleHPManager::GetInstance()->
118 G4ParticleHPProbabilityTablesStore::Ge
119 G4ParticleHPManager::GetInstance()->Re
120 }
121 if ( URRlimits == nullptr ) {
122 G4ParticleHPProbabilityTablesStore::Ge
123 URRlimits = G4ParticleHPProbabilityTab
124 G4ParticleHPManager::GetInstance()->Re
125 }
126 // sets the overall limits of the URR, w
127 // defines URR model energy range
128 SetMinKinEnergy( (*URRlimits).back().fir
129 SetMaxKinEnergy( (*URRlimits).back().sec
130 }
131 }
132 }
133
134
135 G4int G4ParticleHPInelasticDataPT::GetVerboseL
136 return G4ParticleHPManager::GetInstance()->G
137 }
138
139
140 void G4ParticleHPInelasticDataPT::SetVerboseLe
141 G4ParticleHPManager::GetInstance()->SetVerb
142 }
143
144
145 void G4ParticleHPInelasticDataPT::CrossSection
146 outFile << "Inelastic probability tables."
147 }
148