Geant4
9.6.p02
Main Page
Related Pages
Modules
Namespaces
Classes
Files
File List
File Members
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Groups
Pages
geant4_9_6_p02
source
processes
hadronic
models
inclxx
utils
include
G4INCLInverseInterpolationTable.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
// INCL++ intra-nuclear cascade model
27
// Pekka Kaitaniemi, CEA and Helsinki Institute of Physics
28
// Davide Mancusi, CEA
29
// Alain Boudard, CEA
30
// Sylvie Leray, CEA
31
// Joseph Cugnon, University of Liege
32
//
33
#define INCLXX_IN_GEANT4_MODE 1
34
35
#include "
globals.hh
"
36
44
#ifndef G4INCLINVERSEINTERPOLATIONTABLE_HH_
45
#define G4INCLINVERSEINTERPOLATIONTABLE_HH_
46
47
#include "
G4INCLIFunction1D.hh
"
48
#include <algorithm>
49
#include <functional>
50
#include <sstream>
51
52
namespace
G4INCL {
53
54
// Forward declaration
55
class
InverseInterpolationTable;
56
58
class
InterpolationNode
{
59
public
:
60
InterpolationNode
(
const
G4double
x0,
const
G4double
y0,
const
G4double
yPrime0) :
61
x
(x0),
62
y
(y0),
63
yPrime
(yPrime0)
64
{}
65
66
virtual
~InterpolationNode
() {}
67
68
G4bool
operator<
(
const
InterpolationNode
&rhs)
const
{
69
return
(
x
< rhs.
x
);
70
}
71
72
G4bool
operator<=
(
const
InterpolationNode
&rhs)
const
{
73
return
(
x
<= rhs.
x
);
74
}
75
76
G4bool
operator>
(
const
InterpolationNode
&rhs)
const
{
77
return
(
x
> rhs.
x
);
78
}
79
80
G4bool
operator>=
(
const
InterpolationNode
&rhs)
const
{
81
return
(
x
>= rhs.
x
);
82
}
83
85
friend
G4bool
operator<
(
const
InterpolationNode
&lhs,
const
G4double
rhs) {
86
return
lhs.
x
< rhs;
87
}
88
89
G4double
getX
()
const
{
return
x
; }
90
G4double
getY
()
const
{
return
y
; }
91
G4double
getYPrime
()
const
{
return
yPrime
; }
92
93
void
setX
(
const
G4double
x0) {
x
=x0; }
94
void
setY
(
const
G4double
y0) {
y
=y0; }
95
void
setYPrime
(
const
G4double
yPrime0) {
yPrime
=yPrime0; }
96
97
std::string
print
()
const
{
98
std::stringstream message;
99
message <<
"x, y, yPrime: "
<<
x
<<
'\t'
<<
y
<<
'\t'
<<
yPrime
<< std::endl;
100
return
message.str();
101
}
102
103
protected
:
105
G4double
x
;
107
G4double
y
;
109
G4double
yPrime
;
110
};
111
113
class
InverseInterpolationTable
:
public
IFunction1D
{
114
public
:
115
InverseInterpolationTable
(
IFunction1D
const
&
f
,
const
unsigned
int
nNodes=30);
116
InverseInterpolationTable
(std::vector<G4double>
const
&
x
, std::vector<G4double>
const
&
y
);
117
virtual
~InverseInterpolationTable
() {}
118
119
unsigned
int
getNumberOfNodes
()
const
{
return
nodes.size(); }
120
121
std::vector<G4double>
getNodeAbscissae
()
const
{
122
std::vector<G4double>
x
(nodes.size());
123
std::transform(nodes.begin(), nodes.end(), x.begin(),
124
std::mem_fun_ref(&
InterpolationNode::getX
));
125
return
x
;
126
}
127
128
std::vector<G4double>
getNodeValues
()
const
{
129
std::vector<G4double>
y
(nodes.size());
130
std::transform(nodes.begin(), nodes.end(), y.begin(),
131
std::mem_fun_ref(&
InterpolationNode::getY
));
132
return
y
;
133
}
134
135
G4double
operator()
(
const
G4double
x)
const
;
136
137
std::string
print
()
const
;
138
139
private
:
141
void
initDerivatives();
142
144
void
setFunctionDomain();
145
147
std::vector<InterpolationNode> nodes;
148
149
};
150
151
}
152
153
#endif // G4INCLINVERSEINTERPOLATIONTABLE_HH_
Generated on Sat May 25 2013 14:33:59 for Geant4 by
1.8.4