G4NeutronHPPolynomExpansion.hh

Go to the documentation of this file.
00001 //
00002 // ********************************************************************
00003 // * License and Disclaimer                                           *
00004 // *                                                                  *
00005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
00006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
00007 // * conditions of the Geant4 Software License,  included in the file *
00008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
00009 // * include a list of copyright holders.                             *
00010 // *                                                                  *
00011 // * Neither the authors of this software system, nor their employing *
00012 // * institutes,nor the agencies providing financial support for this *
00013 // * work  make  any representation or  warranty, express or implied, *
00014 // * regarding  this  software system or assume any liability for its *
00015 // * use.  Please see the license in the file  LICENSE  and URL above *
00016 // * for the full disclaimer and the limitation of liability.         *
00017 // *                                                                  *
00018 // * This  code  implementation is the result of  the  scientific and *
00019 // * technical work of the GEANT4 collaboration.                      *
00020 // * By using,  copying,  modifying or  distributing the software (or *
00021 // * any work based  on the software)  you  agree  to acknowledge its *
00022 // * use  in  resulting  scientific  publications,  and indicate your *
00023 // * acceptance of all terms of the Geant4 Software license.          *
00024 // ********************************************************************
00025 //
00026 //
00027 // $Id$
00028 //
00029 #ifndef G4NeutronHPPolynomExpansion_h
00030 #define G4NeutronHPPolynomExpansion_h 1
00031 
00032 #include <fstream>
00033 #include <CLHEP/Units/SystemOfUnits.h>
00034 
00035 #include "globals.hh"
00036 #include "G4ios.hh"
00037 
00038 class G4NeutronHPPolynomExpansion
00039 {
00040   public:
00041   G4NeutronHPPolynomExpansion()
00042   {
00043     theCoeff = 0;
00044     nPoly=0;
00045   }
00046   ~G4NeutronHPPolynomExpansion()
00047   {
00048     if(theCoeff!=0) delete [] theCoeff;
00049   }
00050   
00051   inline void Init(std::ifstream & theData)
00052   {
00053     theData >> nPoly;
00054     theCoeff = new G4double[nPoly];
00055     G4int i;
00056     for(i=0;i<nPoly;i++)
00057     {
00058       theData >> theCoeff[i];
00059     }
00060   }
00061   
00062   inline G4double GetValue(G4double anEnergy) 
00063   {
00064     G4int i;
00065     G4double result=0;
00066     G4double base = anEnergy/CLHEP::eV;
00067     G4double running = 1;
00068     for(i=0; i<nPoly; i++)
00069     {
00070       result+=theCoeff[i]*running;
00071       running *= base;
00072     }
00073     return result;
00074   }
00075   private:
00076   
00077   G4int nPoly;
00078   G4double * theCoeff;
00079 };
00080 
00081 #endif

Generated on Mon May 27 17:49:04 2013 for Geant4 by  doxygen 1.4.7