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
G4PSTARStopping.hh
Go to the documentation of this file.
Or view
the newest version in sPHENIX GitHub for file G4PSTARStopping.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
#ifndef G4PSTARStopping_h
28
#define G4PSTARStopping_h 1
29
30
//---------------------------------------------------------------------------
31
//
32
// ClassName: G4PSTARStopping
33
//
34
// Description: Data on stopping power
35
//
36
// Author: Anton Ivantchenko 18.04.2006
37
//
38
// Organisation: QinetiQ Ltd, UK
39
// Customer: ESA/ESTEC, NOORDWIJK
40
// Contract: CSMAN-5288
41
//
42
// Modifications:
43
//
44
//----------------------------------------------------------------------------
45
//
46
// Class Description:
47
//
48
// Data on Stopping Powers from the NIST Data Base
49
// http://physics.nist.gov/PhysRefData/STAR/Text/PSTAR.html
50
//
51
// Current PSTAR database is available via url:
52
// http://physics.nist.gov/PhysRefData/Star/Text/PSTAR.html
53
//
54
55
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
56
57
#include "
globals.hh
"
58
#include "
G4Material.hh
"
59
#include "
G4LPhysicsFreeVector.hh
"
60
#include <vector>
61
62
class
G4PSTARStopping
63
{
64
public
:
65
66
explicit
G4PSTARStopping
();
67
68
~G4PSTARStopping
();
69
70
void
Initialise
();
71
72
inline
G4int
GetIndex
(
const
G4Material
*)
const
;
73
74
inline
G4int
GetIndex
(
const
G4String
&)
const
;
75
76
inline
G4double
GetElectronicDEDX
(
G4int
idx
,
G4double
energy
)
const
;
77
78
inline
G4double
GetElectronicDEDX
(
const
G4Material
*,
G4double
energy)
const
;
79
80
private
:
81
82
void
AddData
(
const
G4float
*
s
,
const
G4Material
*);
83
84
void
FindData
(
G4int
idx,
const
G4Material
*);
85
86
void
PrintWarning
(
G4int
idx)
const
;
87
88
// hide assignment operator
89
G4PSTARStopping
&
operator=
(
const
G4PSTARStopping
&
right
) =
delete
;
90
G4PSTARStopping
(
const
G4PSTARStopping
&) =
delete
;
91
92
G4int
nvectors
;
93
G4double
emin
;
94
std::vector<const G4Material*>
materials
;
95
std::vector<G4LPhysicsFreeVector*>
sdata
;
96
};
97
98
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
99
100
inline
G4int
G4PSTARStopping::GetIndex
(
const
G4Material
*
mat
)
const
101
{
102
G4int
idx
= -1;
103
for
(
G4int
i=0; i<
nvectors
; ++i){
104
if
(mat ==
materials
[i]){
105
idx = i;
106
break
;
107
}
108
}
109
return
idx
;
110
}
111
112
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
113
114
inline
G4int
G4PSTARStopping::GetIndex
(
const
G4String
& nam)
const
115
{
116
G4int
idx
= -1;
117
for
(
G4int
i=0; i<
nvectors
; ++i){
118
if
(nam ==
materials
[i]->GetName()){
119
idx = i;
120
break
;
121
}
122
}
123
return
idx
;
124
}
125
126
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
127
128
inline
G4double
129
G4PSTARStopping::GetElectronicDEDX
(
G4int
idx
,
G4double
energy
)
const
130
{
131
G4double
res = 0.0;
132
if
(idx<0 || idx>=
nvectors
) {
PrintWarning
(idx); }
133
else
if
(energy <
emin
) { res = (*(
sdata
[
idx
]))[0]*std::sqrt(energy/
emin
); }
134
else
{ res =
sdata
[
idx
]->Value(energy); }
135
return
res;
136
}
137
138
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
139
140
inline
G4double
141
G4PSTARStopping::GetElectronicDEDX
(
const
G4Material
*
mat
,
G4double
energy
)
const
142
{
143
return
GetElectronicDEDX
(
GetIndex
(mat), energy);
144
}
145
146
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
147
148
#endif
geant4
tree
geant4-10.6-release
source
processes
electromagnetic
standard
include
G4PSTARStopping.hh
Built by
Jin Huang
. updated:
Wed Jun 29 2022 17:25:34
using
1.8.2 with
ECCE GitHub integration