G4NeutronHPEvapSpectrum.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 G4NeutronHPEvapSpectrum_h
00030 #define G4NeutronHPEvapSpectrum_h 1
00031 
00032 #include <fstream>
00033 #include <CLHEP/Units/SystemOfUnits.h>
00034 
00035 #include "globals.hh"
00036 #include "G4ios.hh"
00037 #include "Randomize.hh"
00038 #include "G4NeutronHPVector.hh"
00039 #include "G4VNeutronHPEDis.hh"
00040 
00041 // we will need a List of these .... one per term.
00042 
00043 class G4NeutronHPEvapSpectrum : public G4VNeutronHPEDis
00044 {
00045   public:
00046   G4NeutronHPEvapSpectrum()
00047   {
00048   }
00049   ~G4NeutronHPEvapSpectrum()
00050   {
00051   }
00052   
00053   inline void Init(std::ifstream & aDataFile)
00054   {
00055     theFractionalProb.Init(aDataFile);
00056     theThetaDist.Init(aDataFile);
00057     theXDist.Init(aDataFile);
00058   }
00059   
00060   inline G4double GetFractionalProbability(G4double anEnergy)
00061   {
00062     return theFractionalProb.GetY(anEnergy);
00063   }
00064   
00065   inline G4double Sample(G4double anEnergy) 
00066   {
00067     // when this is called, theFractionalProb was used, and 'k' is sorted out already.
00068     G4double x = theXDist.Sample();
00069     G4double theta = theThetaDist.GetY(anEnergy);
00070     G4double result = x*theta;
00071     return result*CLHEP::eV;
00072   }
00073   
00074   private:
00075   
00076   G4NeutronHPVector theFractionalProb;
00077   
00078   G4NeutronHPVector theThetaDist;
00079   G4NeutronHPVector theXDist;
00080   
00081 };
00082 
00083 #endif

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