| 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: G4ParticleHPElasticDataPT.c
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 elastic channel.
41 //
42 // Modifications:
43 //
44 // -------------------------------------------
45 //
46 //
47
48 #include "G4ParticleHPElasticDataPT.hh"
49 #include "G4DynamicParticle.hh"
50 #include "G4ParticleDefinition.hh"
51 #include "G4Neutron.hh"
52 #include "G4Element.hh"
53 #include "G4ParticleHPManager.hh"
54 #include "G4ParticleHPProbabilityTablesStore.h
55 #include "G4HadronicException.hh"
56
57
58 G4ParticleHPElasticDataPT::G4ParticleHPElastic
59 // minimum and maximum energy limit for URR
60 SetMinKinEnergy( 1.0*CLHEP::eV );
61 SetMaxKinEnergy( 1.2*CLHEP::MeV );
62 URRlimits = nullptr;
63 }
64
65
66 G4ParticleHPElasticDataPT::~G4ParticleHPElasti
67
68
69 G4bool G4ParticleHPElasticDataPT::IsIsoApplica
70
71 // checks applicability for the element
72 if ( dp->GetDefinition() != G4Neutron::Neutr
73 return false;
74 } else {
75 std::size_t elementI = elm->GetIndex();
76 G4double eKin = dp->GetKineticEnergy();
77 if ( eKin < (*URRlimits).at(elementI).firs
78 return false;
79 } else if ( eKin > (*URRlimits).at(element
80 return false;
81 }
82 }
83 return true;
84 }
85
86
87 G4double G4ParticleHPElasticDataPT::GetIsoCros
88
89 return G4ParticleHPProbabilityTablesStore::G
90 }
91
92
93 void G4ParticleHPElasticDataPT::BuildPhysicsTa
94 G4cout << "BuildPhysicsTable in G4ParticleHP
95 if ( &aP != G4Neutron::Neutron() ) {
96 throw G4HadronicException( __FILE__, __LIN
97 }
98 URRlimits = G4ParticleHPManager::GetInstance
99 if ( G4Threading::IsWorkerThread() ) {
100 // sets the overall limits of the URR, whi
101 // defines URR model energy range
102 SetMinKinEnergy( (*URRlimits).back().first
103 SetMaxKinEnergy( (*URRlimits).back().secon
104 } else {
105 if ( G4ParticleHPManager::GetInstance()->G
106 G4ParticleHPProbabilityTablesStore::GetI
107 G4ParticleHPManager::GetInstance()->Regi
108 }
109 if ( URRlimits == nullptr ) {
110 G4ParticleHPProbabilityTablesStore::GetI
111 URRlimits = G4ParticleHPProbabilityTable
112 G4ParticleHPManager::GetInstance()->Regi
113 }
114 // sets the overall limits of the URR, whi
115 // defines URR model energy range
116 SetMinKinEnergy( (*URRlimits).back().first
117 SetMaxKinEnergy( (*URRlimits).back().secon
118 }
119 }
120
121
122 G4int G4ParticleHPElasticDataPT::GetVerboseLev
123 return G4ParticleHPManager::GetInstance()->
124 }
125
126
127 void G4ParticleHPElasticDataPT::SetVerboseLeve
128 G4ParticleHPManager::GetInstance()->SetVerb
129 }
130
131
132 void G4ParticleHPElasticDataPT::CrossSectionDe
133 outFile << "Elastic probability tables." ;
134 }
135