VertexGen_FastNeutron.hh

// RAT::VertexGen_FastNeutron
//
 
#ifndef __RAT_VertexGen_FastNeutron__
#define __RAT_VertexGen_FastNeutron__
 
#include <TF1.h>
#include <TF2.h>
 
#include <G4Event.hh>
#include <G4ThreeVector.hh>
#include <RAT/DB.hh>
#include <RAT/FastNeutronMessenger.hh>
#include <RAT/GLG4VertexGen.hh>
#include <globals.hh>
#include <vector>
 
namespace RAT {
class VertexGen_FastNeutron : public GLG4VertexGen {
 public:
  VertexGen_FastNeutron(const char *arg_dbname = "fastneutron");
  virtual ~VertexGen_FastNeutron();
  virtual void GeneratePrimaryVertex(G4Event *argEvent, G4ThreeVector &dx, G4double dt);
  virtual void SetState(G4String newValues);
  virtual G4String GetState();
 
  void SetDepth(double DAm);
  void SetEnThreshold(double EAm);
  void SetSideBool(double SBAm);
 
  inline double GetDepth() { return valueD; };
  inline double GetEnThreshold() { return valueE; };
  inline double GetSideBool() { return valueSB; };
 
  void GetMeiHimeParameters(double depth, double emin, double &cosTheta, double &neutronEnergy);
  void LoadTangHortonSmithCosTheta(double depth);
 
  double GetRandomMuonCosTheta();
 
  double evalIntegral(TF2 *func1, double x, double e_tmp);
 
 private:
  G4String _particle;       // name of the particle type
  G4ParticleDefinition *n;  // particle definition
  G4String _FastNeutron;    // name of the FastNeutron to use
 
  double valueD, valueE, valueSB;
  bool loadedTangHortonSmith;
 
  FastNeutronMessenger *messenger;
 
  G4ThreeVector nu_dir;
 
  TF1 *funcMuonCosTheta;
};
 
}  // namespace RAT
 
#endif