Geant4
10.02.p01
Main Page
Related Pages
Modules
Namespaces
Classes
Files
File List
File Members
G4IonPhysics.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
// $Id: G4IonPhysics.cc 80671 2014-05-06 13:59:16Z gcosmo $
27
// GEANT4 tag $Name: $
28
//
29
//---------------------------------------------------------------------------
30
//
31
// Class: G4IonPhysics
32
//
33
// Author: A.Ivanchenko 02.03.2011
34
//
35
// Modified:
36
// 16.10.2012 A.Ribon: renamed G4IonFTFPBinaryCascadePhysics as G4IonPhysics
37
//
38
//---------------------------------------------------------------------------
39
40
#include "
G4IonPhysics.hh
"
41
#include "
G4SystemOfUnits.hh
"
42
#include "
G4ParticleDefinition.hh
"
43
#include "
G4ProcessManager.hh
"
44
#include "
G4Deuteron.hh
"
45
#include "
G4Triton.hh
"
46
#include "
G4He3.hh
"
47
#include "
G4Alpha.hh
"
48
#include "
G4GenericIon.hh
"
49
#include "
G4IonConstructor.hh
"
50
51
#include "
G4HadronInelasticProcess.hh
"
52
#include "
G4BinaryLightIonReaction.hh
"
53
#include "
G4ComponentGGNuclNuclXsc.hh
"
54
#include "
G4CrossSectionInelastic.hh
"
55
56
#include "
G4PreCompoundModel.hh
"
57
#include "
G4ExcitationHandler.hh
"
58
#include "
G4FTFBuilder.hh
"
59
#include "
G4HadronicInteraction.hh
"
60
#include "
G4BuilderType.hh
"
61
#include "
G4HadronicInteractionRegistry.hh
"
62
63
using namespace
std
;
64
65
// factory
66
#include "
G4PhysicsConstructorFactory.hh
"
67
//
68
G4_DECLARE_PHYSCONSTR_FACTORY
(
G4IonPhysics
);
69
70
G4ThreadLocal
G4bool
G4IonPhysics::wasActivated
=
false
;
71
G4ThreadLocal
G4BinaryLightIonReaction
*
G4IonPhysics::theIonBC
= 0;
72
G4ThreadLocal
G4HadronicInteraction
*
G4IonPhysics::theFTFP
= 0;
73
G4ThreadLocal
G4VCrossSectionDataSet
*
G4IonPhysics::theNuclNuclData
= 0;
74
G4ThreadLocal
G4VComponentCrossSection
*
G4IonPhysics::theGGNuclNuclXS
= 0;
75
G4ThreadLocal
G4FTFBuilder
*
G4IonPhysics::theBuilder
;
76
77
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
78
79
G4IonPhysics::G4IonPhysics
(
G4int
ver)
80
:
G4VPhysicsConstructor
(
"ionInelasticFTFP_BIC"
),verbose(ver)
81
{
82
SetPhysicsType
(
bIons
);
83
if
(
verbose
> 1) {
G4cout
<<
"### G4IonPhysics"
<<
G4endl
; }
84
}
85
86
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
87
88
G4IonPhysics::G4IonPhysics
(
const
G4String
& nname)
89
:
G4VPhysicsConstructor
(nname),verbose(1)
90
{
91
SetPhysicsType
(
bIons
);
92
if
(
verbose
> 1) {
G4cout
<<
"### G4IonPhysics"
<<
G4endl
; }
93
}
94
95
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
96
97
G4IonPhysics::~G4IonPhysics
()
98
{
99
//Explictly setting pointers to zero is actually needed.
100
//These are static variables, in case we restart threads we need to re-create objects
101
delete
theBuilder
;
theBuilder
= 0;
102
delete
theGGNuclNuclXS
;
theGGNuclNuclXS
= 0;
103
delete
theNuclNuclData
;
theNuclNuclData
= 0;
104
delete
theIonBC
;
theIonBC
= 0;
105
delete
theFTFP
;
theFTFP
= 0;
106
}
107
108
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
109
110
void
G4IonPhysics::ConstructParticle
()
111
{
112
// Construct ions
113
G4IonConstructor
pConstructor;
114
pConstructor.
ConstructParticle
();
115
}
116
117
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
118
119
void
G4IonPhysics::ConstructProcess
()
120
{
121
if
(
wasActivated
) {
return
; }
122
wasActivated
=
true
;
123
124
G4double
emax
= 100.*
TeV
;
125
126
G4HadronicInteraction
* p =
127
G4HadronicInteractionRegistry::Instance
()->
FindModel
(
"PRECO"
);
128
G4PreCompoundModel
* thePreCompound =
static_cast<
G4PreCompoundModel
*
>
(p);
129
if
(!thePreCompound) { thePreCompound =
new
G4PreCompoundModel
; }
130
131
// Binary Cascade
132
theIonBC
=
new
G4BinaryLightIonReaction
(thePreCompound);
133
theIonBC
->
SetMinEnergy
(0.0);
134
theIonBC
->
SetMaxEnergy
(4*
GeV
);
135
136
// FTFP
137
theBuilder
=
new
G4FTFBuilder
(
"FTFP"
,thePreCompound);
138
theFTFP
=
theBuilder
->
GetModel
();
139
theFTFP
->
SetMinEnergy
(2*
GeV
);
140
theFTFP
->
SetMaxEnergy
(emax);
141
142
theNuclNuclData
=
143
new
G4CrossSectionInelastic
(
theGGNuclNuclXS
=
new
G4ComponentGGNuclNuclXsc
() );
144
145
AddProcess
(
"dInelastic"
,
G4Deuteron::Deuteron
(),
false
);
146
AddProcess
(
"tInelastic"
,
G4Triton::Triton
(),
false
);
147
AddProcess
(
"He3Inelastic"
,
G4He3::He3
(),
true
);
148
AddProcess
(
"alphaInelastic"
,
G4Alpha::Alpha
(),
true
);
149
AddProcess
(
"ionInelastic"
,
G4GenericIon::GenericIon
(),
true
);
150
151
if
(
verbose
> 1) {
152
G4cout
<<
"G4IonPhysics::ConstructProcess done! "
153
<<
G4endl
;
154
}
155
}
156
157
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
158
159
void
G4IonPhysics::AddProcess
(
const
G4String
&
name
,
160
G4ParticleDefinition
* part,
161
G4bool
)
//isIon)
162
{
163
G4HadronInelasticProcess
* hadi =
new
G4HadronInelasticProcess
(name, part);
164
G4ProcessManager
* pManager = part->
GetProcessManager
();
165
pManager->
AddDiscreteProcess
(hadi);
166
167
hadi->
AddDataSet
(
theNuclNuclData
);
168
169
hadi->
RegisterMe
(
theIonBC
);
170
hadi->
RegisterMe
(
theFTFP
);
171
}
172
173
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
G4BinaryLightIonReaction.hh
G4GenericIon.hh
G4Triton.hh
G4IonPhysics::theGGNuclNuclXS
static G4ThreadLocal G4VComponentCrossSection * theGGNuclNuclXS
Definition:
G4IonPhysics.hh:72
G4ExcitationHandler.hh
name
G4String name
Definition:
TRTMaterials.hh:40
G4IonPhysics.hh
G4ProcessManager::AddDiscreteProcess
G4int AddDiscreteProcess(G4VProcess *aProcess, G4int ord=ordDefault)
G4IonPhysics::G4IonPhysics
G4IonPhysics(G4int ver=0)
Definition:
G4IonPhysics.cc:79
std
G4ComponentGGNuclNuclXsc
Definition:
G4ComponentGGNuclNuclXsc.hh:50
G4IonPhysics::theIonBC
static G4ThreadLocal G4BinaryLightIonReaction * theIonBC
Definition:
G4IonPhysics.hh:74
G4HadronicInteractionRegistry.hh
G4IonConstructor.hh
G4BuilderType.hh
G4ParticleDefinition
Definition:
G4ParticleDefinition.hh:111
G4ThreadLocal
#define G4ThreadLocal
Definition:
tls.hh:89
G4IonPhysics::theBuilder
static G4ThreadLocal G4FTFBuilder * theBuilder
Definition:
G4IonPhysics.hh:76
G4VCrossSectionDataSet
Definition:
G4VCrossSectionDataSet.hh:71
G4IonPhysics
Definition:
G4IonPhysics.hh:53
G4ParticleDefinition::GetProcessManager
G4ProcessManager * GetProcessManager() const
G4int
int G4int
Definition:
G4Types.hh:78
G4HadronicInteraction
Definition:
G4HadronicInteraction.hh:64
G4FTFBuilder.hh
G4IonPhysics::~G4IonPhysics
virtual ~G4IonPhysics()
Definition:
G4IonPhysics.cc:97
G4VHadronModelBuilder::GetModel
G4HadronicInteraction * GetModel()
Definition:
G4VHadronModelBuilder.cc:48
G4IonConstructor::ConstructParticle
static void ConstructParticle()
Definition:
G4IonConstructor.cc:59
G4VComponentCrossSection
Definition:
G4VComponentCrossSection.hh:52
G4HadronicProcess::RegisterMe
void RegisterMe(G4HadronicInteraction *a)
Definition:
G4HadronicProcess.cc:149
G4FTFBuilder
Definition:
G4FTFBuilder.hh:50
G4HadronicInteraction::SetMinEnergy
void SetMinEnergy(G4double anEnergy)
Definition:
G4HadronicInteraction.hh:91
G4IonPhysics::theNuclNuclData
static G4ThreadLocal G4VCrossSectionDataSet * theNuclNuclData
Definition:
G4IonPhysics.hh:71
G4HadronicProcess::AddDataSet
void AddDataSet(G4VCrossSectionDataSet *aDataSet)
Definition:
G4HadronicProcess.hh:111
G4PreCompoundModel.hh
G4IonPhysics::wasActivated
static G4ThreadLocal G4bool wasActivated
Definition:
G4IonPhysics.hh:79
G4cout
G4GLOB_DLL std::ostream G4cout
G4PreCompoundModel
Definition:
G4PreCompoundModel.hh:64
G4bool
bool G4bool
Definition:
G4Types.hh:79
G4IonPhysics::theFTFP
static G4ThreadLocal G4HadronicInteraction * theFTFP
Definition:
G4IonPhysics.hh:75
G4CrossSectionInelastic
Definition:
G4CrossSectionInelastic.hh:61
G4Alpha.hh
G4Triton::Triton
static G4Triton * Triton()
Definition:
G4Triton.cc:95
GeV
static const double GeV
Definition:
G4SIunits.hh:214
G4ProcessManager.hh
G4ParticleDefinition.hh
G4BinaryLightIonReaction
Definition:
G4BinaryLightIonReaction.hh:34
G4IonConstructor
Definition:
G4IonConstructor.hh:39
G4HadronInelasticProcess.hh
G4Deuteron::Deuteron
static G4Deuteron * Deuteron()
Definition:
G4Deuteron.cc:94
G4Deuteron.hh
G4IonPhysics::ConstructProcess
void ConstructProcess()
Definition:
G4IonPhysics.cc:119
bIons
Definition:
G4BuilderType.hh:54
emax
static const G4double emax
Definition:
G4ChatterjeeCrossSection.hh:39
G4HadronicInteractionRegistry::FindModel
G4HadronicInteraction * FindModel(const G4String &name)
Definition:
G4HadronicInteractionRegistry.cc:108
G4HadronInelasticProcess
Definition:
G4HadronInelasticProcess.hh:48
G4GenericIon::GenericIon
static G4GenericIon * GenericIon()
Definition:
G4GenericIon.cc:93
G4HadronicInteraction.hh
G4IonPhysics::AddProcess
void AddProcess(const G4String &, G4ParticleDefinition *, G4bool isIon)
Definition:
G4IonPhysics.cc:159
G4He3.hh
G4PhysicsConstructorFactory.hh
G4HadronicInteractionRegistry::Instance
static G4HadronicInteractionRegistry * Instance()
Definition:
G4HadronicInteractionRegistry.cc:36
G4_DECLARE_PHYSCONSTR_FACTORY
G4_DECLARE_PHYSCONSTR_FACTORY(G4IonPhysics)
G4HadronicInteraction::SetMaxEnergy
void SetMaxEnergy(const G4double anEnergy)
Definition:
G4HadronicInteraction.hh:104
G4endl
#define G4endl
Definition:
G4ios.hh:61
G4ProcessManager
Definition:
G4ProcessManager.hh:106
G4IonPhysics::ConstructParticle
void ConstructParticle()
Definition:
G4IonPhysics.cc:110
TeV
static const double TeV
Definition:
G4SIunits.hh:215
G4Alpha::Alpha
static G4Alpha * Alpha()
Definition:
G4Alpha.cc:89
G4double
double G4double
Definition:
G4Types.hh:76
G4SystemOfUnits.hh
G4CrossSectionInelastic.hh
G4He3::He3
static G4He3 * He3()
Definition:
G4He3.cc:94
G4IonPhysics::verbose
G4int verbose
Definition:
G4IonPhysics.hh:78
G4VPhysicsConstructor::SetPhysicsType
void SetPhysicsType(G4int)
Definition:
G4VPhysicsConstructor.hh:199
G4VPhysicsConstructor
Definition:
G4VPhysicsConstructor.hh:121
G4ComponentGGNuclNuclXsc.hh
G4String
Definition:
G4String.hh:45
geant4.10.02.p01
source
physics_lists
constructors
ions
src
G4IonPhysics.cc
Generated on Thu Mar 3 2016 14:41:31 for Geant4 by
1.8.8