Geant4
10.03.p01
Main Page
Related Pages
Namespaces
Classes
Files
File List
File Members
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Pages
G4INCLNuclearPotentialEnergyIsospinSmooth.cc
Go to the documentation of this file.
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::NuclearPotentialIsospin::computePotentialEnergy
virtual G4double computePotentialEnergy(const Particle *const p) const
Definition:
G4INCLNuclearPotentialIsospin.cc:148
G4INCL::NuclearPotential::NuclearPotentialEnergyIsospinSmooth::computePotentialEnergy
virtual G4double computePotentialEnergy(const Particle *const p) const
Definition:
G4INCLNuclearPotentialEnergyIsospinSmooth.cc:67
G4INCL::NuclearPotential::NuclearPotentialEnergyIsospinSmooth::~NuclearPotentialEnergyIsospinSmooth
virtual ~NuclearPotentialEnergyIsospinSmooth()
Definition:
G4INCLNuclearPotentialEnergyIsospinSmooth.cc:65
G4INCLGlobals.hh
G4INCLNuclearPotentialEnergyIsospinSmooth.hh
Isospin- and energy-dependent nuclear potential.
G4int
int G4int
Definition:
G4Types.hh:78
G4INCL::NuclearPotential::INuclearPotential::getFermiEnergy
G4double getFermiEnergy(const Particle *const p) const
Return the Fermi energy for a particle.
Definition:
G4INCLINuclearPotential.hh:98
G4INCL::NuclearPotential::NuclearPotentialIsospin
Definition:
G4INCLNuclearPotentialIsospin.hh:56
G4INCLParticleTable.hh
A
double A(double temperature)
Definition:
G4DNAElectronHoleRecombination.cc:60
G4bool
bool G4bool
Definition:
G4Types.hh:79
globals.hh
G4INCL::Particle
Definition:
G4INCLParticle.hh:73
G4INCL::NuclearPotential::NuclearPotentialEnergyIsospinSmooth::NuclearPotentialEnergyIsospinSmooth
NuclearPotentialEnergyIsospinSmooth(const G4int A, const G4int Z, const G4bool pionPotential)
Definition:
G4INCLNuclearPotentialEnergyIsospinSmooth.cc:60
G4INCL::Particle::isNucleon
G4bool isNucleon() const
Definition:
G4INCLParticle.hh:232
G4INCL::Particle::getKineticEnergy
G4double getKineticEnergy() const
Get the particle kinetic energy.
Definition:
G4INCLParticle.hh:543
G4double
double G4double
Definition:
G4Types.hh:76
Z
G4int Z
Definition:
G4ParticleHPJENDLHEData.cc:262
geant4.10.03.p01
source
processes
hadronic
models
inclxx
incl_physics
src
G4INCLNuclearPotentialEnergyIsospinSmooth.cc
Generated on Thu Mar 16 2017 22:38:10 for Geant4 by
1.8.5