Geant4
v4-10.4-release
메인 페이지
관련된 페이지
모듈
네임스페이스
클래스
파일들
파일 목록
파일 멤버
모두
클래스
네임스페이스들
파일들
함수
변수
타입정의
열거형 타입
열거형 멤버
Friends
매크로
그룹들
페이지들
source
processes
hadronic
models
inclxx
incl_physics
src
G4INCLNuclearPotentialEnergyIsospinSmooth.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
// INCL++ intra-nuclear cascade model
27
// Alain Boudard, CEA-Saclay, France
28
// Joseph Cugnon, University of Liege, Belgium
29
// Jean-Christophe David, CEA-Saclay, France
30
// Pekka Kaitaniemi, CEA-Saclay, France, and Helsinki Institute of Physics, Finland
31
// Sylvie Leray, CEA-Saclay, France
32
// Davide Mancusi, CEA-Saclay, France
33
//
34
#define INCLXX_IN_GEANT4_MODE 1
35
36
#include "
globals.hh
"
37
48
#include "
G4INCLNuclearPotentialEnergyIsospinSmooth.hh
"
49
#include "
G4INCLParticleTable.hh
"
50
#include "
G4INCLGlobals.hh
"
51
52
namespace
G4INCL {
53
54
namespace
NuclearPotential {
55
56
const
G4double
NuclearPotentialEnergyIsospinSmooth::alpha
= 0.223;
57
const
G4double
NuclearPotentialEnergyIsospinSmooth::deltaE
= 25.;
58
59
// Constructors
60
NuclearPotentialEnergyIsospinSmooth::NuclearPotentialEnergyIsospinSmooth
(
const
G4int
A
,
const
G4int
Z
,
const
G4bool
aPionPotential)
61
:
NuclearPotentialIsospin
(A,Z,aPionPotential)
62
{}
63
64
// Destructor
65
NuclearPotentialEnergyIsospinSmooth::~NuclearPotentialEnergyIsospinSmooth
() {}
66
67
G4double
NuclearPotentialEnergyIsospinSmooth::computePotentialEnergy
(
const
Particle
*particle)
const
{
68
69
const
G4double
v0 =
NuclearPotentialIsospin::computePotentialEnergy
(particle);
70
71
if
(particle->
isNucleon
()) {
72
const
G4double
t = particle->
getKineticEnergy
();
73
const
G4double
tf =
getFermiEnergy
(particle);
74
// Constant potential for T<Tf
75
if
(t < tf)
76
return
v0;
77
78
// Linear function for Tf<T<T0, exponential function for T>T0
79
const
G4double
t0 = tf + v0*(1.-
alpha
)/
alpha
-
deltaE
;
// deltaE before the linear potential vanishes
80
G4double
v;
81
if
(t<t0) {
82
v = v0 -
alpha
*(t-tf)/(1.-
alpha
);
83
}
else
{
84
const
G4double
v_at_t0 = v0 -
alpha
*(t0-tf)/(1.-
alpha
);
85
const
G4double
kappa = alpha / (v_at_t0 * (1.-
alpha
));
86
v = v_at_t0 * std::exp(kappa * (t0-t));
87
}
88
return
(v>0.0) ? v : 0.0;
89
}
else
90
return
v0;
91
}
92
93
}
94
}
95
G4INCL::NuclearPotential::NuclearPotentialEnergyIsospinSmooth::alpha
static const G4double alpha
Slope of the V(T) curve.
Definition:
G4INCLNuclearPotentialEnergyIsospinSmooth.hh:67
G4INCL::Particle::getKineticEnergy
G4double getKineticEnergy() const
Get the particle kinetic energy.
Definition:
G4INCLParticle.hh:669
G4INCL::NuclearPotential::NuclearPotentialEnergyIsospinSmooth::deltaE
static const G4double deltaE
Distance from the cusp where the exponential kicks in.
Definition:
G4INCLNuclearPotentialEnergyIsospinSmooth.hh:70
G4INCL::NuclearPotential::NuclearPotentialIsospin::computePotentialEnergy
virtual G4double computePotentialEnergy(const Particle *const p) const
Definition:
G4INCLNuclearPotentialIsospin.cc:186
G4INCLNuclearPotentialEnergyIsospinSmooth.hh
Isospin- and energy-dependent nuclear potential.
G4INCLParticleTable.hh
Z
Float_t Z
Definition:
advanced/microbeam/plot.C:38
G4double
double G4double
Definition:
G4Types.hh:76
G4bool
bool G4bool
Definition:
G4Types.hh:79
G4INCL::NuclearPotential::NuclearPotentialEnergyIsospinSmooth::NuclearPotentialEnergyIsospinSmooth
NuclearPotentialEnergyIsospinSmooth(const G4int A, const G4int Z, const G4bool pionPotential)
Definition:
G4INCLNuclearPotentialEnergyIsospinSmooth.cc:60
A
double A(double temperature)
Definition:
G4DNAElectronHoleRecombination.cc:60
globals.hh
G4INCLGlobals.hh
G4int
int G4int
Definition:
G4Types.hh:78
G4INCL::Particle
Definition:
G4INCLParticle.hh:75
G4INCL::NuclearPotential::NuclearPotentialEnergyIsospinSmooth::~NuclearPotentialEnergyIsospinSmooth
virtual ~NuclearPotentialEnergyIsospinSmooth()
Definition:
G4INCLNuclearPotentialEnergyIsospinSmooth.cc:65
G4INCL::NuclearPotential::INuclearPotential::getFermiEnergy
G4double getFermiEnergy(const Particle *const p) const
Return the Fermi energy for a particle.
Definition:
G4INCLINuclearPotential.hh:105
G4INCL::NuclearPotential::NuclearPotentialIsospin
Definition:
G4INCLNuclearPotentialIsospin.hh:56
G4INCL::Particle::isNucleon
G4bool isNucleon() const
Definition:
G4INCLParticle.hh:296
G4INCL::NuclearPotential::NuclearPotentialEnergyIsospinSmooth::computePotentialEnergy
virtual G4double computePotentialEnergy(const Particle *const p) const
Definition:
G4INCLNuclearPotentialEnergyIsospinSmooth.cc:67
다음에 의해 생성됨 :
1.8.5