source: trunk/source/processes/hadronic/models/neutron_hp/include/G4NeutronHPInterpolator.hh@ 1199

Last change on this file since 1199 was 1196, checked in by garnier, 16 years ago

update CVS release candidate geant4.9.3.01
File size: 6.3 KB
Line 
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//
27// $Id: G4NeutronHPInterpolator.hh,v 1.20 2008/08/12 00:42:31 tkoi Exp $
28// GEANT4 tag $Name: geant4-09-03-cand-01 $
29//
30// 080809 Change interpolation scheme of "histogram", now using LinearLinear
31// For multidimensional interpolations By T. Koi
32//
33#ifndef G4NeutronHPInterpolator_h
34#define G4NeutronHPInterpolator_h 1
35
36#include "globals.hh"
37#include "G4InterpolationScheme.hh"
38#include "Randomize.hh"
39#include "G4ios.hh"
40#include "G4HadronicException.hh"
41
42
43class G4NeutronHPInterpolator
44{
45 public:
46
47 G4NeutronHPInterpolator(){}
48 ~G4NeutronHPInterpolator()
49 {
50 // G4cout <<"deleted the interpolator"<<G4endl;
51 }
52
53 inline G4double Lin(G4double x,G4double x1,G4double x2,G4double y1,G4double y2)
54 {
55 G4double slope=0, off=0;
56 if(x2-x1==0) return (y2+y1)/2.;
57 slope = (y2-y1)/(x2-x1);
58 off = y2-x2*slope;
59 G4double y = x*slope+off;
60 return y;
61 }
62
63 inline G4double Interpolate(G4InterpolationScheme aScheme,
64 G4double x, G4double x1, G4double x2,
65 G4double y1, G4double y2) const;
66
67 G4double
68 GetBinIntegral(const G4InterpolationScheme & aScheme,
69 const G4double x1,const G4double x2,const G4double y1,const G4double y2);
70
71 G4double
72 GetWeightedBinIntegral(const G4InterpolationScheme & aScheme,
73 const G4double x1,const G4double x2,const G4double y1,const G4double y2);
74
75 private:
76
77 inline G4double Histogram(G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const;
78 inline G4double LinearLinear(G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const;
79 inline G4double LinearLogarithmic(G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const;
80 inline G4double LogarithmicLinear(G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const;
81 inline G4double LogarithmicLogarithmic(G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const;
82 inline G4double Random(G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const;
83
84};
85
86inline G4double G4NeutronHPInterpolator::
87Interpolate(G4InterpolationScheme aScheme,
88 G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const
89{
90 G4double result(0);
91 G4int theScheme = aScheme;
92 theScheme = theScheme%CSTART_;
93 switch(theScheme)
94 {
95 case 1:
96 //080809
97 //result = Histogram(x, x1, x2, y1, y2);
98 result = LinearLinear(x, x1, x2, y1, y2);
99 break;
100 case 2:
101 result = LinearLinear(x, x1, x2, y1, y2);
102 break;
103 case 3:
104 result = LinearLogarithmic(x, x1, x2, y1, y2);
105 break;
106 case 4:
107 result = LogarithmicLinear(x, x1, x2, y1, y2);
108 break;
109 case 5:
110 result = LogarithmicLogarithmic(x, x1, x2, y1, y2);
111 break;
112 case 6:
113 result = Random(x, x1, x2, y1, y2);
114 break;
115 default:
116 G4cout << "theScheme = "<<theScheme<<G4endl;
117 throw G4HadronicException(__FILE__, __LINE__, "G4NeutronHPInterpolator::Carthesian Invalid InterpolationScheme");
118 break;
119 }
120 return result;
121}
122
123inline G4double G4NeutronHPInterpolator::
124Histogram(G4double , G4double , G4double , G4double y1, G4double ) const
125{
126 G4double result;
127 result = y1;
128 return result;
129}
130
131inline G4double G4NeutronHPInterpolator::
132LinearLinear(G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const
133{
134 G4double slope=0, off=0;
135 if(x2-x1==0) return (y2+y1)/2.;
136 slope = (y2-y1)/(x2-x1);
137 off = y2-x2*slope;
138 G4double y = x*slope+off;
139 return y;
140}
141
142inline G4double G4NeutronHPInterpolator::
143LinearLogarithmic(G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const
144{
145 G4double result;
146 if(x==0) result = y1+y2/2.;
147 else if(x1==0) result = y1;
148 else if(x2==0) result = y2;
149 else result = LinearLinear(std::log(x), std::log(x1), std::log(x2), y1, y2);
150 return result;
151}
152
153inline G4double G4NeutronHPInterpolator::
154LogarithmicLinear(G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const
155{
156 G4double result;
157 if(y1==0||y2==0) result = 0;
158 else
159 {
160 result = LinearLinear(x, x1, x2, std::log(y1), std::log(y2));
161 result = std::exp(result);
162 }
163 return result;
164}
165
166inline G4double G4NeutronHPInterpolator::
167LogarithmicLogarithmic(G4double x, G4double x1, G4double x2, G4double y1, G4double y2) const
168{
169 G4double result;
170 if(x==0) result = y1+y2/2.;
171 else if(x1==0) result = y1;
172 else if(x2==0) result = y2;
173 if(y1==0||y2==0) result = 0;
174 else
175 {
176 result = LinearLinear(std::log(x), std::log(x1), std::log(x2), std::log(y1), std::log(y2));
177 result = std::exp(result);
178 }
179 return result;
180}
181
182inline G4double G4NeutronHPInterpolator::
183Random(G4double , G4double , G4double , G4double y1, G4double y2) const
184{
185 G4double result;
186 result = y1+G4UniformRand()*(y2-y1);
187 return result;
188}
189
190#endif
Note: See TracBrowser for help on using the repository browser.