Loading...
Searching...
No Matches
G4SemiLogInterpolation Class Reference
#include <G4SemiLogInterpolation.hh>
Inheritance diagram for G4SemiLogInterpolation:
Public Member Functions | |
| G4SemiLogInterpolation () | |
| ~G4SemiLogInterpolation () | |
| G4double | Calculate (G4double point, G4int bin, const G4DataVector &energies, const G4DataVector &data) const |
| G4double | Calculate (G4double point, G4int bin, const G4DataVector &energies, const G4DataVector &data, const G4DataVector &log_energies, const G4DataVector &log_data) const |
| virtual G4VDataSetAlgorithm * | Clone () const |
| 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 G4VDataSetAlgorithm * | Clone () const =0 |
Detailed Description
Definition at line 53 of file G4SemiLogInterpolation.hh.
Constructor & Destructor Documentation
◆ G4SemiLogInterpolation()
| G4SemiLogInterpolation::G4SemiLogInterpolation | ( | ) |
◆ ~G4SemiLogInterpolation()
| G4SemiLogInterpolation::~G4SemiLogInterpolation | ( | ) |
Definition at line 47 of file G4SemiLogInterpolation.cc.
48{ }
Member Function Documentation
◆ Calculate() [1/2]
|
virtual |
Implements G4VDataSetAlgorithm.
Definition at line 54 of file G4SemiLogInterpolation.cc.
57{
58 //G4cout << "G4SemiLogInterpolation is performed(2 arguments) " << G4endl;
59 G4int nBins = data.size() - 1;
60 G4double value = 0.;
61 if (x < points[0])
62 {
63 value = 0.;
64 }
65 else if (bin < nBins)
66 {
67 G4double e1 = points[bin];
68 G4double e2 = points[bin+1];
69 G4double d1 = data[bin];
70 G4double d2 = data[bin+1];
71 value = (d1*std::log10(e2/x) + d2*std::log10(x/e1)) / std::log10(e2/e1);
72 }
73 else
74 {
75 value = data[nBins];
76 }
77 return value;
78}
◆ Calculate() [2/2]
|
virtual |
Implements G4VDataSetAlgorithm.
Definition at line 80 of file G4SemiLogInterpolation.cc.
85{
86//A combination of logarithmic interpolation on energy set and
87//linear Interpolation on data set
88 //G4cout << "G4SemiLogInterpolation is performed (4 arguments)" << G4endl;
89 G4int nBins = data.size() - 1;
90 G4double value = 0.;
91 G4double log_x = std::log10(x);
92 if (x < points[0])
93 {
94 value = 0.;
95 }
96 else if (bin < nBins)
97 {
98 G4double e1 = points[bin];
99 G4double e2 = points[bin+1];
100 G4double d1 = data[bin];
101 G4double d2 = data[bin+1];
102 G4double log_e1 = log_points[bin];
103 G4double log_e2 = log_points[bin+1];
104 //G4double log_d1 = log_data[bin];
105 //G4double log_d2 = log_data[bin+1];
106 //G4cout << "x = " << x << " , log_x = " << log_x << " , bin = " << bin << G4endl;
107 //G4cout << "e1 = " << e1 << " , d1 = " << d1 << G4endl;
108 //G4cout << "e2 = " << e2 << " , d2 = " << d2 << G4endl;
109// Values log_e1 and log_e2 are the log values of the corresponding
110// original energy actual values. Original d1 and d2 values are used.
111// Simple linear interpolation performed on loagarithmic data
112// should be equivalent to semi log-log interpolation
113 if (e1 == 0.0) log_e1 = -300;
114 if (e2 == 0.0) log_e2 = -300;
115 value = d1 + (d2 - d1)*(log_x - log_e1)/(log_e2 - log_e1);
116 //G4cout << "G4SemiLogInterpolation - Final Interpolated Value: " << value << G4endl << G4endl;
117 }
118 else
119 {
120 value = data[nBins];
121 }
122 return value;
123}
◆ Clone()
|
virtual |
Implements G4VDataSetAlgorithm.
Definition at line 50 of file G4SemiLogInterpolation.cc.
G4SemiLogInterpolation()
Definition: G4SemiLogInterpolation.cc:41
The documentation for this class was generated from the following files:
