Geant4.10
Main Page
Related Pages
Modules
Namespaces
Data Structures
Files
File List
Globals
All
Data Structures
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Groups
Pages
geant4.10.00.p01
source
processes
hadronic
cross_sections
src
G4IonsSihverCrossSection.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
// 18-Sep-2003 First version is written by T. Koi
27
// 23-Dec-2006 Isotope dependence added by D. Wright
28
// 14-Mar-2011 Moved constructor, destructor and virtual methods to source by V.Ivanchenko
29
// 19-Aug-2011 V.Ivanchenko move to new design and make x-section per element
30
//
31
32
#include "
G4IonsSihverCrossSection.hh
"
33
#include "
G4PhysicalConstants.hh
"
34
#include "
G4SystemOfUnits.hh
"
35
#include "
G4DynamicParticle.hh
"
36
#include "
G4HadTmpUtil.hh
"
37
#include "
G4NistManager.hh
"
38
#include "
G4Pow.hh
"
39
40
G4IonsSihverCrossSection::G4IonsSihverCrossSection
()
41
:
G4VCrossSectionDataSet
(
"IonsSihver"
), square_r0 ( (1.36*
fermi
) * (1.36*
fermi
) )
42
{}
43
44
G4IonsSihverCrossSection::~G4IonsSihverCrossSection
()
45
{}
46
47
void
48
G4IonsSihverCrossSection::CrossSectionDescription
(std::ostream&
outFile
)
const
49
{
50
outFile <<
"G4IonsSihverCrossSection calculates the total reaction cross\n"
51
<<
"section for nucleus-nucleus scattering using the Sihver\n"
52
<<
"parameterization. It is valid for projectiles and targets of\n"
53
<<
"all Z, and for all projectile energies above 100 MeV/n.\n"
;
54
}
55
56
G4bool
57
G4IonsSihverCrossSection::IsElementApplicable
(
const
G4DynamicParticle
* aDP,
58
G4int
,
const
G4Material
*)
59
{
60
G4int
BaryonNumber = aDP->
GetDefinition
()->
GetBaryonNumber
();
61
G4double
KineticEnergy = aDP->
GetKineticEnergy
();
62
if
( KineticEnergy / BaryonNumber >= 100*
MeV
&& BaryonNumber > 1 ) {
return
true
; }
63
return
false
;
64
}
65
66
G4double
67
G4IonsSihverCrossSection::GetElementCrossSection
(
68
const
G4DynamicParticle
* aParticle,
G4int
Z,
const
G4Material
*)
69
{
70
G4double
xsection = 0.0;
71
G4int
At =
G4lrint
(
G4NistManager::Instance
()->GetAtomicMassAmu(Z));
72
73
G4int
Ap = aParticle->
GetDefinition
()->
GetBaryonNumber
();
74
75
G4Pow
* g4pow =
G4Pow::GetInstance
();
76
77
G4double
cubicrAt = g4pow->
Z13
(At);
78
G4double
cubicrAp = g4pow->
Z13
(Ap);
79
80
G4double
b0 = 1.581 - 0.876 * (1.0/cubicrAp + 1.0/cubicrAt);
81
82
G4double
xr = cubicrAp + cubicrAt - b0 * (1.0/cubicrAp + 1.0/cubicrAt);
83
xsection =
pi
* square_r0 * xr * xr;
84
85
return
xsection;
86
}
87
G4Pow::GetInstance
static G4Pow * GetInstance()
Definition:
G4Pow.cc:53
python.hepunit.pi
float pi
Definition:
hepunit.py:247
G4IonsSihverCrossSection.hh
G4DynamicParticle::GetKineticEnergy
G4double GetKineticEnergy() const
G4Pow
Definition:
G4Pow.hh:56
outFile
std::ofstream outFile
Definition:
GammaRayTel.cc:68
G4NistManager.hh
python.hepunit.MeV
MeV
Definition:
hepunit.py:117
G4IonsSihverCrossSection::CrossSectionDescription
virtual void CrossSectionDescription(std::ostream &) const
Definition:
G4IonsSihverCrossSection.cc:48
G4Material
Definition:
G4Material.hh:118
G4DynamicParticle
Definition:
G4DynamicParticle.hh:73
G4IonsSihverCrossSection::~G4IonsSihverCrossSection
virtual ~G4IonsSihverCrossSection()
Definition:
G4IonsSihverCrossSection.cc:44
G4DynamicParticle.hh
G4DynamicParticle::GetDefinition
G4ParticleDefinition * GetDefinition() const
G4VCrossSectionDataSet
Definition:
G4VCrossSectionDataSet.hh:70
G4int
int G4int
Definition:
G4Types.hh:78
G4NistManager::Instance
static G4NistManager * Instance()
Definition:
G4NistManager.cc:68
python.hepunit.fermi
int fermi
Definition:
hepunit.py:37
G4Pow::Z13
G4double Z13(G4int Z) const
Definition:
G4Pow.hh:129
G4IonsSihverCrossSection::G4IonsSihverCrossSection
G4IonsSihverCrossSection()
Definition:
G4IonsSihverCrossSection.cc:40
G4bool
bool G4bool
Definition:
G4Types.hh:79
G4PhysicalConstants.hh
G4IonsSihverCrossSection::GetElementCrossSection
virtual G4double GetElementCrossSection(const G4DynamicParticle *, G4int Z, const G4Material *)
Definition:
G4IonsSihverCrossSection.cc:67
G4lrint
int G4lrint(double ad)
Definition:
templates.hh:163
G4IonsSihverCrossSection::IsElementApplicable
virtual G4bool IsElementApplicable(const G4DynamicParticle *aDP, G4int Z, const G4Material *)
Definition:
G4IonsSihverCrossSection.cc:57
G4double
double G4double
Definition:
G4Types.hh:76
G4SystemOfUnits.hh
G4Pow.hh
G4ParticleDefinition::GetBaryonNumber
G4int GetBaryonNumber() const
Definition:
G4ParticleDefinition.hh:183
G4HadTmpUtil.hh
Generated on Wed Apr 30 2014 15:54:51 for Geant4.10 by
1.8.7