Geant4
v4-10.4-release
메인 페이지
관련된 페이지
모듈
네임스페이스
클래스
파일들
파일 목록
파일 멤버
모두
클래스
네임스페이스들
파일들
함수
변수
타입정의
열거형 타입
열거형 멤버
Friends
매크로
그룹들
페이지들
examples
advanced
composite_calorimeter
include
CCalDataSet.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
#ifndef CCalDataSet_h
27
#define CCalDataSet_h
28
29
#include "
G4ParticleHPInterpolator.hh
"
30
#include <vector>
31
32
// by JPW, working, but to be cleaned up. @@@@
33
34
class
CCalDataSet
35
{
36
public
:
37
38
void
insert
(
double
anEnergy,
double
aXsection)
39
{
40
pair<double, double> aPoint(anEnergy, aXsection);
41
theData
.push_back(aPoint);
42
}
43
44
double
getCrossSection
(
double
anEnergy)
45
{
46
double
result
;
47
if
(anEnergy <
theData
[0].
first
)
48
{
49
result = 0;
50
}
51
else
52
{
53
double
x1
,
x2
;
54
double
y1
,
y2
;
55
if
(anEnergy >
theData
[
theData
.size()-1].second)
56
{
57
// extrapolation
58
int
n
=
theData
.size();
59
x1 =
theData
[n-1].first;
60
y1 =
theData
[n-1].second;
61
x2 =
theData
[n-2].first;
62
y2 =
theData
[n-2].second;
63
}
64
else
65
{
66
// linear search and interpolation
67
unsigned
int
i;
68
for
(i=0; i<
theData
.size(); i++)
69
{
70
if
(
theData
[i].first>anEnergy)
break
;
71
}
72
x1 =
theData
[i-1].first;
73
y1 =
theData
[i-1].second;
74
x2 =
theData
[i].first;
75
y2 =
theData
[i].second;
76
}
77
// interpolate logarithmic in energy and linear in cross-section
78
result =
theInterpolator
.
Interpolate
(
LOGLIN
, anEnergy, x1,x2,y1,y2);
79
}
80
return
result*
barn
;
81
}
82
83
private
:
84
85
G4ParticleHPInterpolator
theInterpolator
;
86
87
private
:
88
89
vector<pair<double, double> >
theData
;
90
};
91
#endif
CCalDataSet::insert
void insert(double anEnergy, double aXsection)
Definition:
CCalDataSet.hh:38
first
G4int first
Definition:
G4ParticleHPJENDLHEData.cc:280
y1
Float_t y1[n_points_granero]
Definition:
compare.C:5
G4ParticleHPInterpolator::Interpolate
G4double Interpolate(G4InterpolationScheme aScheme, G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const
Definition:
G4ParticleHPInterpolator.hh:91
x1
Float_t x1[n_points_granero]
Definition:
compare.C:5
G4ParticleHPInterpolator
Definition:
G4ParticleHPInterpolator.hh:44
CCalDataSet::getCrossSection
double getCrossSection(double anEnergy)
Definition:
CCalDataSet.hh:44
y2
Float_t y2[n_points_geant4]
Definition:
compare.C:26
CCalDataSet::theInterpolator
G4ParticleHPInterpolator theInterpolator
Definition:
CCalDataSet.hh:85
CCalDataSet
Definition:
CCalDataSet.hh:34
LOGLIN
Definition:
G4InterpolationScheme.hh:39
G4ParticleHPInterpolator.hh
result
G4double G4ParticleHPJENDLHEData::G4double result
Definition:
G4ParticleHPJENDLHEData.cc:257
CCalDataSet::theData
vector< pair< double, double > > theData
Definition:
CCalDataSet.hh:89
barn
static constexpr double barn
Definition:
G4SIunits.hh:105
n
Char_t n[5]
Definition:
comparison_ascii.C:55
x2
Float_t x2[n_points_geant4]
Definition:
compare.C:26
다음에 의해 생성됨 :
1.8.5