Geant4
v4-10.4-release
메인 페이지
관련된 페이지
모듈
네임스페이스
클래스
파일들
파일 목록
파일 멤버
모두
클래스
네임스페이스들
파일들
함수
변수
타입정의
열거형 타입
열거형 멤버
Friends
매크로
그룹들
페이지들
source
processes
hadronic
models
pre_equilibrium
exciton_model
src
G4HETCNeutron.cc
이 파일의 문서화 페이지로 가기
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
// $Id: G4HETCNeutron.cc 90337 2015-05-26 08:34:27Z gcosmo $
27
//
28
// by V. Lara
29
//
30
// Modified:
31
// 23.08.2010 V.Ivanchenko general cleanup, move constructor and destructor
32
// the source, use G4Pow
33
34
#include "
G4HETCNeutron.hh
"
35
#include "
G4PhysicalConstants.hh
"
36
#include "
G4SystemOfUnits.hh
"
37
#include "
G4Neutron.hh
"
38
39
G4HETCNeutron::G4HETCNeutron
()
40
:
G4HETCFragment
(
G4Neutron
::
Neutron
(), &theNeutronCoulombBarrier)
41
{}
42
43
G4HETCNeutron::~G4HETCNeutron
()
44
{}
45
46
G4double
G4HETCNeutron::GetAlpha
()
const
47
{
48
return
0.76+2.2/
theResA13
;
49
}
50
51
G4double
G4HETCNeutron::GetBeta
()
const
52
{
53
return
(2.12/(
theResA13
*
theResA13
)-0.05)*
MeV
/
GetAlpha
();
54
}
55
56
G4double
G4HETCNeutron::GetSpinFactor
()
const
57
{
58
// (2s+1)
59
return
2.0;
60
}
61
62
G4double
G4HETCNeutron::K
(
const
G4Fragment
& aFragment)
63
{
64
// Number of protons in emitted fragment
65
G4int
Pa =
theZ
;
66
// Number of neutrons in emitted fragment
67
G4int
Na =
theA
- Pa;
68
69
G4double
r =
G4double
(
theResZ
)/
G4double
(
theResA
);
70
71
G4int
P
= aFragment.
GetNumberOfParticles
();
72
G4int
H = aFragment.
GetNumberOfHoles
();
73
74
G4double
result
= 0.0;
75
if
(P > 0)
76
{
77
result = (H + Na/(1.0-r))/
P
;
78
}
79
80
return
std::max
(0.0,result);
81
}
82
83
G4double
G4HETCNeutron::SampleKineticEnergy
(
const
G4Fragment
& aFragment)
84
{
85
G4int
H = aFragment.
GetNumberOfHoles
();
86
G4int
Pb = aFragment.
GetNumberOfParticles
();
87
G4int
Nb = Pb + H;
88
G4double
g0 = (6.0/
pi2
)*
theFragA
*
theParameters
->
GetLevelDensity
();
89
90
G4double
Ab =
std::max
(0.0,
G4double
(Pb*Pb+H*H+Pb-3*H)/(4.0*g0));
91
G4double
Emax
=
theMaxKinEnergy
- Ab;
92
93
G4double
cut =
GetBeta
() / (
GetBeta
()+Emax/
G4double
(Nb+1));
94
G4double
x
(0.0);
95
if
(
G4UniformRand
() <= cut)
96
{
97
x =
BetaRand
(Nb,1);
98
}
99
else
100
{
101
x =
BetaRand
(Nb,2);
102
}
103
104
return
Emax * (1.0 -
x
);
105
}
x
Float_t x
Definition:
compare.C:6
G4HETCNeutron::G4HETCNeutron
G4HETCNeutron()
Definition:
G4HETCNeutron.cc:39
G4HETCFragment::BetaRand
G4double BetaRand(G4int N, G4int L) const
Definition:
G4HETCFragment.hh:79
G4HETCNeutron.hh
G4INCL::Math::max
T max(const T t1, const T t2)
brief Return the largest of the two arguments
Definition:
G4INCLGlobals.hh:112
G4PhysicalConstants.hh
G4VPreCompoundFragment::theResA
G4int theResA
Definition:
G4VPreCompoundFragment.hh:146
G4HETCFragment
Definition:
G4HETCFragment.hh:40
G4VPreCompoundFragment::theParameters
G4DeexPrecoParameters * theParameters
Definition:
G4VPreCompoundFragment.hh:141
MeV
static constexpr double MeV
Definition:
G4SIunits.hh:214
G4VPreCompoundFragment::theFragA
G4int theFragA
Definition:
G4VPreCompoundFragment.hh:148
G4VPreCompoundFragment::theResA13
G4double theResA13
Definition:
G4VPreCompoundFragment.hh:151
G4HETCNeutron::SampleKineticEnergy
virtual G4double SampleKineticEnergy(const G4Fragment &aFragment)
Definition:
G4HETCNeutron.cc:83
G4Fragment::GetNumberOfHoles
G4int GetNumberOfHoles() const
Definition:
G4Fragment.hh:358
G4Fragment::GetNumberOfParticles
G4int GetNumberOfParticles() const
Definition:
G4Fragment.hh:338
G4double
double G4double
Definition:
G4Types.hh:76
Neutron
Definition:
G4RadioactiveDecayMode.hh:68
G4UniformRand
#define G4UniformRand()
Definition:
Randomize.hh:53
G4Fragment
Definition:
G4Fragment.hh:66
G4SystemOfUnits.hh
G4HETCNeutron::GetBeta
virtual G4double GetBeta() const
Definition:
G4HETCNeutron.cc:51
pi2
static constexpr double pi2
Definition:
G4SIunits.hh:78
result
G4double G4ParticleHPJENDLHEData::G4double result
Definition:
G4ParticleHPJENDLHEData.cc:257
P
static double P[]
Definition:
nf_gammaFunctions.cc:89
G4HETCNeutron::~G4HETCNeutron
~G4HETCNeutron()
Definition:
G4HETCNeutron.cc:43
G4DeexPrecoParameters::GetLevelDensity
G4double GetLevelDensity() const
Definition:
G4DeexPrecoParameters.hh:253
G4VPreCompoundFragment::theZ
G4int theZ
Definition:
G4VPreCompoundFragment.hh:145
G4int
int G4int
Definition:
G4Types.hh:78
G4Neutron.hh
Emax
static const G4double Emax
Definition:
G4PhotoNuclearCrossSection.cc:51
G4HETCNeutron::GetSpinFactor
virtual G4double GetSpinFactor() const
Definition:
G4HETCNeutron.cc:56
G4HETCNeutron::K
virtual G4double K(const G4Fragment &aFragment)
Definition:
G4HETCNeutron.cc:62
G4HETCNeutron::GetAlpha
virtual G4double GetAlpha() const
Definition:
G4HETCNeutron.cc:46
G4VPreCompoundFragment::theMaxKinEnergy
G4double theMaxKinEnergy
Definition:
G4VPreCompoundFragment.hh:154
G4VPreCompoundFragment::theA
G4int theA
Definition:
G4VPreCompoundFragment.hh:144
G4VPreCompoundFragment::theResZ
G4int theResZ
Definition:
G4VPreCompoundFragment.hh:147
G4Neutron
Definition:
G4Neutron.hh:51
다음에 의해 생성됨 :
1.8.5