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