Geant4 11.1.1
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4LinLogInterpolation Class Reference

#include <G4LinLogInterpolation.hh>

+ Inheritance diagram for G4LinLogInterpolation:

Public Member Functions

 G4LinLogInterpolation ()
 
 ~G4LinLogInterpolation ()
 
G4double Calculate (G4double point, G4int bin, const G4DataVector &energies, const G4DataVector &data) const override
 
G4double Calculate (G4double point, G4int bin, const G4DataVector &energies, const G4DataVector &data, const G4DataVector &log_energies, const G4DataVector &log_data) const override
 
G4VDataSetAlgorithmClone () const override
 
- Public Member Functions inherited from G4VDataSetAlgorithm
 G4VDataSetAlgorithm ()
 
virtual ~G4VDataSetAlgorithm ()
 
virtual G4double Calculate (G4double point, G4int bin, const G4DataVector &energies, const G4DataVector &data) const =0
 
virtual G4double Calculate (G4double point, G4int bin, const G4DataVector &energies, const G4DataVector &data, const G4DataVector &log_energies, const G4DataVector &log_data) const =0
 
virtual G4VDataSetAlgorithmClone () const =0
 
 G4VDataSetAlgorithm (const G4VDataSetAlgorithm &)=delete
 
G4VDataSetAlgorithmoperator= (const G4VDataSetAlgorithm &right)=delete
 

Detailed Description

Definition at line 49 of file G4LinLogInterpolation.hh.

Constructor & Destructor Documentation

◆ G4LinLogInterpolation()

G4LinLogInterpolation::G4LinLogInterpolation ( )
explicit

Definition at line 41 of file G4LinLogInterpolation.cc.

42{ }

Referenced by Clone().

◆ ~G4LinLogInterpolation()

G4LinLogInterpolation::~G4LinLogInterpolation ( )

Definition at line 46 of file G4LinLogInterpolation.cc.

47{ }

Member Function Documentation

◆ Calculate() [1/2]

G4double G4LinLogInterpolation::Calculate ( G4double  point,
G4int  bin,
const G4DataVector energies,
const G4DataVector data 
) const
overridevirtual

Implements G4VDataSetAlgorithm.

Definition at line 50 of file G4LinLogInterpolation.cc.

53{
54 //G4cout << "G4LinLogInterpolation is performed on dataset (2 arguments) " << G4endl;
55 G4int nBins = G4int(data.size() - 1);
56 G4double value = 0.;
57 if (x < points[0])
58 {
59 value = 0.;
60 }
61 else if (bin < nBins)
62 {
63 G4double e1 = points[bin];
64 G4double e2 = points[bin+1];
65 G4double d1 = std::log(data[bin]);
66 G4double d2 = std::log(data[bin+1]);
67 value = G4Exp(d1 + (d2 - d1)*(x - e1)/ (e2 - e1));
68 }
69 else
70 {
71 value = data[nBins];
72 }
73 return value;
74}
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
Definition: G4Exp.hh:180
double G4double
Definition: G4Types.hh:83
int G4int
Definition: G4Types.hh:85

◆ Calculate() [2/2]

G4double G4LinLogInterpolation::Calculate ( G4double  point,
G4int  bin,
const G4DataVector energies,
const G4DataVector data,
const G4DataVector log_energies,
const G4DataVector log_data 
) const
overridevirtual

Implements G4VDataSetAlgorithm.

Definition at line 78 of file G4LinLogInterpolation.cc.

83{
84 //G4cout << "G4LinLogInterpolation is performed on dataset (4 arguments) " << G4endl;
85 G4int nBins = G4int(data.size() - 1);
86 G4double value = 0.;
87 //G4double log_x = std::log10(x);
88 if (x < points[0])
89 {
90 value = 0.;
91 }
92 else if (bin < nBins)
93 {
94 G4double e1 = points[bin];
95 G4double e2 = points[bin+1];
96 G4double d1 = data[bin];
97 G4double d2 = data[bin+1];
98 //G4double log_e1 = log_points[bin];
99 //G4double log_e2 = log_points[bin+1];
100 G4double log_d1 = log_data[bin];
101 G4double log_d2 = log_data[bin+1];
102 if (d1 > 0.0 && d2 > 0.0)
103 {
104 // Values e1, e2, d1 and d2 are the log values of the corresponding
105 // original energy and data values.
106 // Simple linear interpolation performed
107 // on logarithmic data should be equivalent to log-log interpolation
108 //value = log_d1 + (log_d2 - log_d1)*(log_x - log_e1)/(log_e2 - log_e1);
109 value = G4Exp(log_d1 + (log_d2 - log_d1)*(x - e1)/(e2 - e1));
110 }
111 else
112 {
113 if (d1 == 0.0) log_d1 = -300;
114 if (d2 == 0.0) log_d2 = -300;
115 value = G4Exp(log_d1 + (log_d2 - log_d1)*(x - e1)/(e2 - e1));
116 }
117 }
118 else
119 {
120 value = data[nBins];
121 }
122 return value;
123}

◆ Clone()

G4VDataSetAlgorithm * G4LinLogInterpolation::Clone ( ) const
inlineoverridevirtual

Implements G4VDataSetAlgorithm.

Definition at line 65 of file G4LinLogInterpolation.hh.


The documentation for this class was generated from the following files: