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
G4eeToHadronsMultiModel.hh
Go to the documentation of this file.
Or view
the newest version in sPHENIX GitHub for file G4eeToHadronsMultiModel.hh
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
// GEANT4 Class header file
30
//
31
//
32
// File name: G4eeToHadronsMultiModel
33
//
34
// Author: Vladimir Ivanchenko
35
//
36
// Creation date: 18.05.2005
37
//
38
// Modifications:
39
//
40
41
//
42
// Class Description: vector of e+e- -> hadrons models
43
//
44
45
// -------------------------------------------------------------------
46
//
47
48
#ifndef G4eeToHadronsMultiModel_h
49
#define G4eeToHadronsMultiModel_h 1
50
51
#include "
G4VEmModel.hh
"
52
#include "
G4eeToHadronsModel.hh
"
53
#include "
G4ParticleChangeForGamma.hh
"
54
#include "
G4TrackStatus.hh
"
55
#include "
Randomize.hh
"
56
#include <
CLHEP/Units/SystemOfUnits.h
>
57
#include <
CLHEP/Units/PhysicalConstants.h
>
58
#include <vector>
59
60
class
G4eeCrossSections
;
61
class
G4Vee2hadrons
;
62
63
class
G4eeToHadronsMultiModel
:
public
G4VEmModel
64
{
65
66
public
:
67
68
explicit
G4eeToHadronsMultiModel
(
G4int
ver=0,
69
const
G4String
& nam =
"eeToHadrons"
);
70
71
virtual
~G4eeToHadronsMultiModel
();
72
73
virtual
void
Initialise
(
const
G4ParticleDefinition
*,
74
const
G4DataVector
&)
override
;
75
76
virtual
G4double
CrossSectionPerVolume
(
const
G4Material
*,
77
const
G4ParticleDefinition
*,
78
G4double
kineticEnergy,
79
G4double
cutEnergy,
80
G4double
maxEnergy)
override
;
81
82
virtual
G4double
ComputeCrossSectionPerAtom
(
83
const
G4ParticleDefinition
*,
84
G4double
kineticEnergy,
85
G4double
Z
,
G4double
A
,
86
G4double
cutEnergy = 0.0,
87
G4double
maxEnergy =
DBL_MAX
)
override
;
88
89
virtual
void
SampleSecondaries
(std::vector<G4DynamicParticle*>*,
90
const
G4MaterialCutsCouple
*,
91
const
G4DynamicParticle
*,
92
G4double
tmin = 0.0,
93
G4double
maxEnergy =
DBL_MAX
)
override
;
94
95
virtual
void
ModelDescription
(std::ostream& outFile)
const override
;
96
97
// Set the factor to artificially increase the crossSection (default 1)
98
void
SetCrossSecFactor
(
G4double
fac
);
99
100
inline
G4double
ComputeCrossSectionPerElectron
(
101
const
G4ParticleDefinition
*,
102
G4double
kineticEnergy,
103
G4double
cutEnergy = 0.0,
104
G4double
maxEnergy =
DBL_MAX
);
105
106
private
:
107
108
void
AddEEModel
(
G4Vee2hadrons
*,
const
G4DataVector
&);
109
110
//change incident e+ kinetic energy into CM total energy(sum of e+ and e-)
111
inline
G4double
LabToCM
(
G4double
);
112
113
// hide assignment operator
114
G4eeToHadronsMultiModel
&
operator=
(
const
G4eeToHadronsMultiModel
&
right
);
115
G4eeToHadronsMultiModel
(
const
G4eeToHadronsMultiModel
&);
116
117
G4eeCrossSections
*
cross
;
118
G4ParticleChangeForGamma
*
fParticleChange
;
119
G4double
delta
;
120
121
std::vector<G4eeToHadronsModel*>
models
;
122
123
std::vector<G4double>
ekinMin
;
124
std::vector<G4double>
ekinPeak
;
125
std::vector<G4double>
ekinMax
;
126
std::vector<G4double>
cumSum
;
127
128
G4double
thKineticEnergy
;
129
G4double
maxKineticEnergy
;
130
G4double
csFactor
;
131
132
G4int
nModels
;
133
G4int
verbose
;
134
G4bool
isInitialised
;
135
};
136
137
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
138
139
//change incident e+ kinetic energy into CM total energy(sum of e+ and e-)
140
inline
G4double
G4eeToHadronsMultiModel::LabToCM
(
G4double
kinE_lab)
141
{
142
G4double
totE_CM = 0.0;
143
G4double
mass
=
CLHEP::electron_mass_c2
;
144
G4double
totE_lab = kinE_lab +
mass
;
145
totE_CM = std::sqrt(2*mass*(mass+totE_lab));
146
147
return
totE_CM;
148
}
149
150
inline
G4double
G4eeToHadronsMultiModel::ComputeCrossSectionPerElectron
(
151
const
G4ParticleDefinition
*,
152
G4double
kineticEnergy,
153
G4double
,
G4double
)
154
{
155
G4double
res = 0.0;
156
157
G4double
energy
=
LabToCM
(kineticEnergy);
158
159
if
(energy >
thKineticEnergy
) {
160
for
(
G4int
i=0; i<
nModels
; i++) {
161
if
(energy >=
ekinMin
[i] && energy <=
ekinMax
[i]){
162
res += (
models
[i])->
ComputeCrossSectionPerElectron
(0,energy);
163
}
164
cumSum
[i] = res;
165
}
166
}
167
return
res*
csFactor
;
168
}
169
170
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
171
172
#endif
geant4
tree
geant4-10.6-release
source
processes
electromagnetic
highenergy
include
G4eeToHadronsMultiModel.hh
Built by
Jin Huang
. updated:
Wed Jun 29 2022 17:25:31
using
1.8.2 with
ECCE GitHub integration