Geant4
10.02.p02
Main Page
Related Pages
Modules
Namespaces
Classes
Files
File List
File Members
CCalDataSet.hh
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
#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
G4ParticleHPInterpolator
Definition:
G4ParticleHPInterpolator.hh:44
CCalDataSet::theData
vector< pair< double, double > > theData
Definition:
CCalDataSet.hh:89
G4ParticleHPInterpolator.hh
CCalDataSet
Definition:
CCalDataSet.hh:34
n
const G4int n
Definition:
G4UrQMD1_3Interface.hh:144
CCalDataSet::theInterpolator
G4ParticleHPInterpolator theInterpolator
Definition:
CCalDataSet.hh:85
G4ParticleHPInterpolator::Interpolate
G4double Interpolate(G4InterpolationScheme aScheme, G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const
Definition:
G4ParticleHPInterpolator.hh:91
CCalDataSet::getCrossSection
double getCrossSection(double anEnergy)
Definition:
CCalDataSet.hh:44
barn
static const double barn
Definition:
G4SIunits.hh:104
LOGLIN
Definition:
G4InterpolationScheme.hh:39
CCalDataSet::insert
void insert(double anEnergy, double aXsection)
Definition:
CCalDataSet.hh:38
geant4.10.02.p02
examples
advanced
composite_calorimeter
include
CCalDataSet.hh
Generated on Fri Jul 1 2016 14:10:32 for Geant4 by
1.8.8