Geant4  9.6.p02
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4NeutronHPLegendreTable.hh
Go to the documentation of this file.
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$
28 //
29 #ifndef G4NeutronHPLegendreTable_h
30 #define G4NeutronHPLegendreTable_h 1
31 
32 #include <fstream>
34 
35 #include "globals.hh"
36 #include "G4ios.hh"
38 
40 {
41  public:
43  {
44  nCoeff=0;
45  theCoeff = 0;
46  }
47  ~G4NeutronHPLegendreTable(){if(theCoeff!=0) delete [] theCoeff;}
48 
50  {
51  if(&aSet!=this)
52  {
53  theRep = aSet.theRep;
54  theEnergy = aSet.theEnergy;
55  theTemp = aSet.theTemp;
56  theManager = aSet.theManager;
57  nCoeff = aSet.nCoeff;
58  if(theCoeff!=0) delete [] theCoeff;
59  theCoeff = new G4double[nCoeff];
60  for(G4int i=0; i<nCoeff; i++)
61  {
62  theCoeff[i] = aSet.theCoeff[i];
63  }
64  }
65  }
66 
67  inline void Init(std::ifstream & aDataFile)
68  {
69  G4double eNeu, coeff;
70  G4int nPoly;
71  aDataFile >> eNeu >> nPoly;
72  eNeu *= CLHEP::eV;
73  Init(eNeu, nPoly);
74  for(G4int l=0; l<nPoly; l++)
75  {
76  aDataFile >> coeff;
77  SetCoeff(l+1, coeff);
78  }
79  }
80 
81  inline void Init(G4double e, G4int n)
82  {
83  nCoeff = n+1;
84  theCoeff = new G4double[nCoeff];
85  for(G4int i=0; i<nCoeff; i++) theCoeff[i] = 0;
86  theCoeff[0]=1.;
87  theEnergy = e;
88 // G4cout << "G4NeutronHPLegendreTable::Init called "<<e<<" "<<n<<G4endl;
89  }
90  inline void SetEnergy(G4double energy){ theEnergy = energy; }
91  inline void SetTemperature(G4double temp){ theTemp = temp; }
92  inline void SetCoeff(G4int l, G4double coeff) {theCoeff[l]=coeff;}
93  inline void SetRepresentation(G4int aRep) {theRep = aRep;}
94 
95  inline G4double GetCoeff(G4int l) {return theCoeff[l];}
96  inline G4double GetEnergy(){return theEnergy;}
97  inline G4double GetTemperature(){return theTemp;}
98  inline G4int GetNumberOfPoly() {return nCoeff;}
99  inline G4int GetRepresentation() {return theRep;}
100  inline const G4InterpolationManager & GetManager() {return theManager;}
101  private:
102 
103  G4int theRep;
104  G4double theEnergy;
105  G4double theTemp;
106  G4int nCoeff;
107  G4InterpolationManager theManager; // knows the interpolation between stores
108  G4double * theCoeff;
109 };
110 
111 #endif