| 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: G4ParticleHPElasticURR.cc
32 //
33 // Authors: Marek Zmeskal (CTU, Czech Tec
34 // Loic Thulliez (CEA France)
35 //
36 // Creation date: 4 June 2024
37 //
38 // Description: Class to handle URR range
39 // proper isotope cross-sect
40 //
41 // Modifications:
42 //
43 // -------------------------------------------
44 //
45 //
46
47 #include "G4ParticleHPElasticURR.hh"
48 #include "G4ParticleHPManager.hh"
49 #include "G4ParticleHPChannel.hh"
50 #include "G4ParticleHPElastic.hh"
51 #include "G4ParticleHPProbabilityTablesStore.h
52 #include "G4SystemOfUnits.hh"
53 #include "G4Threading.hh"
54
55
56 G4ParticleHPElasticURR::G4ParticleHPElasticURR
57 G4HadronicInteraction( "NeutronHPElasticURR"
58 G4double minEnergy = 0.0;
59 if ( isThermalScatteringOn ) minEnergy = 4.0
60 SetMinEnergy( minEnergy );
61 SetMaxEnergy( 20.0 * CLHEP::MeV );
62 particleHPelastic = new G4ParticleHPElastic;
63 }
64
65
66 G4ParticleHPElasticURR::~G4ParticleHPElasticUR
67
68
69 G4HadFinalState* G4ParticleHPElasticURR::Apply
70 const G4Material* theMaterial = aTrack.GetMa
71 G4double kineticEnergy = aTrack.GetKineticEn
72 G4HadFinalState* theFinalState = nullptr;
73 if ( kineticEnergy < (*URRlimits).back().fir
74 return particleHPelastic->ApplyYourself( a
75 }
76 G4int elementI = -1;
77 G4int isotopeJ = -1;
78 G4int A = aNucleus.GetA_asInt();
79 G4int Z = aNucleus.GetZ_asInt();
80 // finds the element and isotope of the sele
81 for ( G4int i = 0; i < (G4int)theMaterial->G
82 if ( Z == theMaterial->GetElement(i)->GetZ
83 for ( G4int j = 0; j < (G4int)theMateria
84 if ( A == theMaterial->GetElement(i)->
85 isotopeJ = j;
86 break;
87 }
88 }
89 // the loop cannot be ended here because
90 if ( isotopeJ != -1 ) {
91 // isotope was found and for loop is ended
92 elementI = (G4int)theMaterial->GetElement(i)
93 break;
94 }
95 } // end if find element
96 } // end element loop
97 // Check whether the energy is out of the UR
98 if ( kineticEnergy < (*URRlimits).at(element
99 // Call elastic final state in G4ParicleHP
100 G4ParticleHPManager::GetInstance()->OpenRe
101 theFinalState = (*G4ParticleHPManager::Get
102 // Update target nucleus information accor
103 G4int selectedIsotope_A = G4ParticleHPMana
104 aNucleus.SetParameters( selectedIsotope_A,
105 const G4Element* target_element = (*G4Elem
106 const G4Isotope* target_isotope = nullptr;
107 // Find the selected isotope among in the
108 for ( G4int j = 0; j < (G4int)target_eleme
109 target_isotope = target_element->GetIsot
110 if ( target_isotope->GetN() == selectedI
111 }
112 aNucleus.SetIsotope( target_isotope );
113 G4ParticleHPManager::GetInstance()->CloseR
114 } else {
115 // the energy is inside the limits of the
116 theFinalState = ((*G4ParticleHPManager::Ge
117 }
118 return theFinalState;
119 }
120
121
122 void G4ParticleHPElasticURR::BuildPhysicsTable
123 particleHPelastic->BuildPhysicsTable( *(G4Ne
124 URRlimits = G4ParticleHPManager::GetInstance
125 if ( URRlimits == nullptr ) {
126 G4ParticleHPProbabilityTablesStore::GetIns
127 URRlimits = G4ParticleHPProbabilityTablesS
128 G4ParticleHPManager::GetInstance()->Regist
129 }
130 }
131
132
133 const std::pair< G4double, G4double > G4Partic
134 // max energy non-conservation is mass of he
135 return std::pair< G4double, G4double >( 10.0
136 }
137
138
139 G4int G4ParticleHPElasticURR::GetVerboseLevel(
140 return G4ParticleHPManager::GetInstance()->G
141 }
142
143
144 void G4ParticleHPElasticURR::SetVerboseLevel(
145 G4ParticleHPManager::GetInstance()->SetVerb
146 }
147
148
149 void G4ParticleHPElasticURR::ModelDescription(
150 outFile << "High Precision model based on E
151 }
152