EnergeticDD.h

Go to the documentation of this file.

#include "NSDG4PrimaryGenerator.h"
namespace NSDG4
{
class EnergeticDD : public PrimaryGeneratorAction
{
public:
  //default constructor puts the class into a working state.
  EnergeticDD();
 
  virtual void GeneratePrimaries(G4Event* event);
  
private:
  // Gets the neutron energy as a function of angle for a given deueteron energy
  double getNeutronEnergyThetaEd(double a_Ed,double a_theta);
  //reads the required anisotropy file from disk
  void readAnisoFile();
 
  //linear array of cdf data for specified energies
  std::vector<double> m_cdfData;
  //energies assoicated with the above array
  std::vector<double> m_cdfEn;
  //angles associated with the cdf
  std::vector<double> m_theta;
 
  //integer width of a set of cdf values for a given energy
  int m_cdfDataWidth;
  //loads an interpolated cdf for a given energy
  void loadCDF(double a_Ed);
  //used to store a current estimated CDF for an interpolated energy
  std::vector<double> m_interpCDF;
  //returns a theta for a given energy
  double getTheta(double a_Ed);
 
 
  bool m_isReflected;
 
};
 
}