source: trunk/source/processes/hadronic/models/neutron_hp/include/G4NeutronHPFastLegendre.hh @ 1340

Last change on this file since 1340 was 1340, checked in by garnier, 14 years ago

update ti head

File size: 6.5 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: G4NeutronHPFastLegendre.hh,v 1.11 2006/06/29 20:47:39 gunter Exp $
28// GEANT4 tag $Name: geant4-09-03-ref-09 $
29//
30#ifndef G4NeutronHPFastLegendre_h
31#define G4NeutronHPFastLegendre_h 1
32
33#include "globals.hh"
34
35class G4NeutronHPFastLegendre
36{
37  public:
38 
39  G4NeutronHPFastLegendre()
40  {
41    value = new G4double * [31];
42    value[0] = l0;
43    value[1] = l1;
44    value[2] = l2;
45    value[3] = l3;
46    value[4] = l4;
47    value[5] = l5;
48    value[6] = l6;
49    value[7] = l7;
50    value[8] = l8;
51    value[9] = l9;
52    value[10] = l10;
53    value[11] = l11;
54    value[12] = l12;
55    value[13] = l13;
56    value[14] = l14;
57    value[15] = l15;
58    value[16] = l16;
59    value[17] = l17;
60    value[18] = l18;
61    value[19] = l19;
62    value[20] = l20;
63    value[21] = l21;
64    value[22] = l22;
65    value[23] = l23;
66    value[24] = l24;
67    value[25] = l25;
68    value[26] = l26;
69    value[27] = l27;
70    value[28] = l28;
71    value[29] = l29;
72    value[30] = l30;
73    integral = new G4double * [31];
74    integral[0] = i0;
75    integral[1] = i1;
76    integral[2] = i2;
77    integral[3] = i3;
78    integral[4] = i4;
79    integral[5] = i5;
80    integral[6] = i6;
81    integral[7] = i7;
82    integral[8] = i8;
83    integral[9] = i9;
84    integral[10] = i10;
85    integral[11] = i11;
86    integral[12] = i12;
87    integral[13] = i13;
88    integral[14] = i14;
89    integral[15] = i15;
90    integral[16] = i16;
91    integral[17] = i17;
92    integral[18] = i18;
93    integral[19] = i19;
94    integral[20] = i20;
95    integral[21] = i21;
96    integral[22] = i22;
97    integral[23] = i23;
98    integral[24] = i24;
99    integral[25] = i25;
100    integral[26] = i26;
101    integral[27] = i27;
102    integral[28] = i28;
103    integral[29] = i29;
104    integral[30] = i30;
105   
106    G4int i;
107    for(i=0;i<31;i++) theNbin[i]=1+200*(i+1);
108  }
109 
110  ~G4NeutronHPFastLegendre()
111  {
112    delete [] value;
113    delete [] integral;
114  }
115 
116  G4double Integrate(G4int l, G4double costh)
117  {
118    G4int bin = GetBin(l, costh);
119    G4double y1, y2;
120//    G4cout <<"Testhpw G4NeutronHPFastLegendre::Integrate "<<l<<" "<<bin<<G4endl;
121    y1 = integral[l][bin];
122    y2 = integral[l][bin+1];
123//    G4cout <<"Testhpw G4NeutronHPFastLegendre::Integrate exit"<<G4endl;
124    return Interpolate(bin, l, y1, y2, costh);
125  }
126 
127  inline G4double Evaluate(G4int l, G4double costh)
128  {
129    G4double result;
130    G4int bin = GetBin(l, costh);
131    if(bin != theNbin[l]-1)
132    {
133      G4double y1, y2;
134      y1 = value[l][bin];
135      y2 = value[l][bin+1];
136      result = Interpolate(bin, l, y1, y2, costh);
137    }
138    else
139    {
140      result = value[l][bin];
141    }
142    return result;
143  }
144  private:
145 
146  inline G4int GetBin(G4int l, G4double costh)
147  {
148    G4int bin=0;
149    bin = G4int( (theNbin[l]-1)*(costh+1)/2. );
150    if(bin == theNbin[l]-1) bin--;
151    return bin;
152  }
153 
154  inline G4double Interpolate(G4int bin, G4int l, G4double y1, G4double y2, G4double x)
155  {
156    G4double slope = 0, off = 0, x2=0, x1mx2;
157    G4int half = (theNbin[l]-1)/2;
158//    x1 = (bin-half)/G4double(half);
159    x2 = (bin+1-half)/G4double(half);
160    x1mx2 = 1./G4double( (theNbin[l]-1)/2 );
161//    slope = (y2-y1)/(x2-x1);
162    slope = (y2-y1)/x1mx2;
163    off = y2-x2*slope;
164    return x*slope+off;
165  }
166 
167  G4double ** value;
168  G4double ** integral;
169  G4int theNbin[31];
170  static G4double l0[201];
171  static G4double i0[201];
172  static G4double l1[401];
173  static G4double i1[401];
174  static G4double l2[601];
175  static G4double i2[601];
176  static G4double l3[801];
177  static G4double i3[801];
178  static G4double l4[1001];
179  static G4double i4[1001];
180  static G4double l5[1201];
181  static G4double i5[1201];
182  static G4double l6[1401];
183  static G4double i6[1401];
184  static G4double l7[1601];
185  static G4double i7[1601];
186  static G4double l8[1801];
187  static G4double i8[1801];
188  static G4double l9[2001];
189  static G4double i9[2001];
190  static G4double l10[2201];
191  static G4double i10[2201];
192  static G4double l11[2401];
193  static G4double i11[2401];
194  static G4double l12[2601];
195  static G4double i12[2601];
196  static G4double l13[2801];
197  static G4double i13[2801];
198  static G4double l14[3001];
199  static G4double i14[3001];
200  static G4double l15[3201];
201  static G4double i15[3201];
202  static G4double l16[3401];
203  static G4double i16[3401];
204  static G4double l17[3601];
205  static G4double i17[3601];
206  static G4double l18[3801];
207  static G4double i18[3801];
208  static G4double l19[4001];
209  static G4double i19[4001];
210  static G4double l20[4201];
211  static G4double i20[4201];
212  static G4double l21[4401];
213  static G4double i21[4401];
214  static G4double l22[4601];
215  static G4double i22[4601];
216  static G4double l23[4801];
217  static G4double i23[4801];
218  static G4double l24[5001];
219  static G4double i24[5001];
220  static G4double l25[5201];
221  static G4double i25[5201];
222  static G4double l26[5401];
223  static G4double i26[5401];
224  static G4double l27[5601];
225  static G4double i27[5601];
226  static G4double l28[5801];
227  static G4double i28[5801];
228  static G4double l29[6001];
229  static G4double i29[6001];
230  static G4double l30[6201];
231  static G4double i30[6201];
232};
233#endif
Note: See TracBrowser for help on using the repository browser.