Geant4
10.02.p01
Main Page
Related Pages
Modules
Namespaces
Classes
Files
File List
File Members
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:55
MeV
static const double MeV
Definition:
G4SIunits.hh:211
G4IonsSihverCrossSection.hh
G4DynamicParticle::GetKineticEnergy
G4double GetKineticEnergy() const
G4Pow
Definition:
G4Pow.hh:56
G4NistManager.hh
G4IonsSihverCrossSection::CrossSectionDescription
virtual void CrossSectionDescription(std::ostream &) const
Definition:
G4IonsSihverCrossSection.cc:48
G4Material
Definition:
G4Material.hh:120
G4DynamicParticle
Definition:
G4DynamicParticle.hh:73
G4IonsSihverCrossSection::~G4IonsSihverCrossSection
virtual ~G4IonsSihverCrossSection()
Definition:
G4IonsSihverCrossSection.cc:44
G4DynamicParticle.hh
G4DynamicParticle::GetDefinition
G4ParticleDefinition * GetDefinition() const
G4VCrossSectionDataSet
Definition:
G4VCrossSectionDataSet.hh:71
G4int
int G4int
Definition:
G4Types.hh:78
G4NistManager::Instance
static G4NistManager * Instance()
Definition:
G4NistManager.cc:68
G4Pow::Z13
G4double Z13(G4int Z) const
Definition:
G4Pow.hh:127
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
pi
static const double pi
Definition:
G4SIunits.hh:74
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
G4IonsSihverCrossSection::square_r0
const G4double square_r0
Definition:
G4IonsSihverCrossSection.hh:65
G4double
double G4double
Definition:
G4Types.hh:76
G4SystemOfUnits.hh
G4Pow.hh
b0
static const G4double b0
Definition:
G4BetheHeitlerModel.cc:76
fermi
static const double fermi
Definition:
G4SIunits.hh:102
G4ParticleDefinition::GetBaryonNumber
G4int GetBaryonNumber() const
Definition:
G4ParticleDefinition.hh:183
G4HadTmpUtil.hh
geant4.10.02.p01
source
processes
hadronic
cross_sections
src
G4IonsSihverCrossSection.cc
Generated on Thu Mar 3 2016 14:41:51 for Geant4 by
1.8.8