ECCE @ EIC Software
Reference for
ECCE @ EIC
simulation and reconstruction software on GitHub
Home page
Related Pages
Modules
Namespaces
Classes
Files
External Links
File List
File Members
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Groups
Pages
GammaRayTelIonPhysics.cc
Go to the documentation of this file.
Or view
the newest version in sPHENIX GitHub for file GammaRayTelIonPhysics.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
//
27
//
28
//
29
30
#include <iomanip>
31
32
#include "
GammaRayTelIonPhysics.hh
"
33
34
#include "
globals.hh
"
35
#include "
G4ios.hh
"
36
#include "
G4SystemOfUnits.hh
"
37
#include "
G4HadronicParameters.hh
"
38
39
40
GammaRayTelIonPhysics::GammaRayTelIonPhysics
(
const
G4String
&
name
)
41
:
G4VPhysicsConstructor
(name)
42
{;}
43
44
GammaRayTelIonPhysics::~GammaRayTelIonPhysics
()
45
{;}
46
47
void
GammaRayTelIonPhysics::ConstructParticle
()
48
{;}
49
50
51
#include "
G4ProcessManager.hh
"
52
53
54
void
GammaRayTelIonPhysics::ConstructProcess
()
55
{
56
G4ProcessManager
* pManager = 0;
57
58
const
G4double
theBERTMin = 0.0*
GeV
;
59
const
G4double
theBERTMax = 5.0*
GeV
;
60
const
G4double
theFTFMin = 4.0*
GeV
;
61
const
G4double
theFTFMax =
G4HadronicParameters::Instance
()->
GetMaxEnergy
();
62
63
G4FTFModel
* theStringModel =
new
G4FTFModel
;
64
G4ExcitedStringDecay
* theStringDecay =
new
G4ExcitedStringDecay
(
new
G4LundStringFragmentation
);
65
theStringModel->
SetFragmentationModel
( theStringDecay );
66
G4PreCompoundModel
* thePreEquilib =
new
G4PreCompoundModel
(
new
G4ExcitationHandler
);
67
G4GeneratorPrecompoundInterface
* theCascade =
new
G4GeneratorPrecompoundInterface
( thePreEquilib );
68
69
G4TheoFSGenerator
* theModel =
new
G4TheoFSGenerator
(
"FTFP"
);
70
theModel->
SetHighEnergyGenerator
( theStringModel );
71
theModel->
SetTransport
( theCascade );
72
theModel->
SetMinEnergy
( theFTFMin );
73
theModel->
SetMaxEnergy
( theFTFMax );
74
75
G4CascadeInterface
* theBERTModel =
new
G4CascadeInterface
;
76
theBERTModel->
SetMinEnergy
( theBERTMin );
77
theBERTModel->
SetMaxEnergy
( theBERTMax );
78
79
// Elastic Process
80
theElasticModel
=
new
G4HadronElastic
();
81
theElasticProcess
.
RegisterMe
(
theElasticModel
);
82
83
// Generic Ion
84
pManager =
G4GenericIon::GenericIon
()->
GetProcessManager
();
85
// add process
86
pManager->
AddDiscreteProcess
(&
theElasticProcess
);
87
88
pManager->
AddProcess
(&
fIonIonisation
,
ordInActive
, 2, 2);
89
90
pManager->
AddProcess
(&
fIonMultipleScattering
);
91
pManager->
SetProcessOrdering
(&
fIonMultipleScattering
,
idxAlongStep
, 1);
92
pManager->
SetProcessOrdering
(&
fIonMultipleScattering
,
idxPostStep
, 1);
93
94
G4ComponentGGNuclNuclXsc
* ggNuclNuclXsec =
new
G4ComponentGGNuclNuclXsc
();
95
G4VCrossSectionDataSet
* theGGNuclNuclData =
new
G4CrossSectionInelastic
(ggNuclNuclXsec);
96
97
// Deuteron
98
pManager =
G4Deuteron::Deuteron
()->
GetProcessManager
();
99
// add process
100
pManager->
AddDiscreteProcess
(&
theElasticProcess
);
101
fDeuteronProcess
.
AddDataSet
(theGGNuclNuclData);
102
fDeuteronProcess
.
RegisterMe
(theBERTModel);
103
fDeuteronProcess
.
RegisterMe
(theModel);
104
pManager->
AddDiscreteProcess
(&
fDeuteronProcess
);
105
106
pManager->
AddProcess
(&
fDeuteronIonisation
,
ordInActive
, 2, 2);
107
108
pManager->
AddProcess
(&
fDeuteronMultipleScattering
);
109
pManager->
SetProcessOrdering
(&
fDeuteronMultipleScattering
,
idxAlongStep
, 1);
110
pManager->
SetProcessOrdering
(&
fDeuteronMultipleScattering
,
idxPostStep
, 1);
111
112
// Triton
113
pManager =
G4Triton::Triton
()->
GetProcessManager
();
114
// add process
115
pManager->
AddDiscreteProcess
(&
theElasticProcess
);
116
fTritonProcess
.
AddDataSet
(theGGNuclNuclData);
117
fTritonProcess
.
RegisterMe
(theBERTModel);
118
fTritonProcess
.
RegisterMe
(theModel);
119
pManager->
AddDiscreteProcess
(&
fTritonProcess
);
120
121
pManager->
AddProcess
(&
fTritonIonisation
,
ordInActive
, 2, 2);
122
123
pManager->
AddProcess
(&
fTritonMultipleScattering
);
124
pManager->
SetProcessOrdering
(&
fTritonMultipleScattering
,
idxAlongStep
, 1);
125
pManager->
SetProcessOrdering
(&
fTritonMultipleScattering
,
idxPostStep
, 1);
126
127
// Alpha
128
pManager =
G4Alpha::Alpha
()->
GetProcessManager
();
129
// add process
130
pManager->
AddDiscreteProcess
(&
theElasticProcess
);
131
132
fAlphaProcess
.
AddDataSet
(theGGNuclNuclData);
133
fAlphaProcess
.
RegisterMe
(theBERTModel);
134
fAlphaProcess
.
RegisterMe
(theModel);
135
pManager->
AddDiscreteProcess
(&
fAlphaProcess
);
136
137
pManager->
AddProcess
(&
fAlphaIonisation
,
ordInActive
, 2, 2);
138
139
pManager->
AddProcess
(&
fAlphaMultipleScattering
);
140
pManager->
SetProcessOrdering
(&
fAlphaMultipleScattering
,
idxAlongStep
, 1);
141
pManager->
SetProcessOrdering
(&
fAlphaMultipleScattering
,
idxPostStep
, 1);
142
143
// He3
144
pManager =
G4He3::He3
()->
GetProcessManager
();
145
// add process
146
pManager->
AddDiscreteProcess
(&
theElasticProcess
);
147
148
pManager->
AddProcess
(&
fHe3Ionisation
,
ordInActive
, 2, 2);
149
150
pManager->
AddProcess
(&
fHe3MultipleScattering
);
151
pManager->
SetProcessOrdering
(&
fHe3MultipleScattering
,
idxAlongStep
, 1);
152
pManager->
SetProcessOrdering
(&
fHe3MultipleScattering
,
idxPostStep
, 1);
153
154
}
155
156
157
geant4
tree
geant4-10.6-release
examples
advanced
gammaray_telescope
src
GammaRayTelIonPhysics.cc
Built by
Jin Huang
. updated:
Wed Jun 29 2022 17:24:58
using
1.8.2 with
ECCE GitHub integration