CCCrossSec.hh

 
#ifndef __RAT_CCCrossSec__
#define __RAT_CCCrossSec__
 
// G4 headers
#include <globals.hh>
// RAT headers
#include <RAT/GLG4StringUtil.hh>
#include <RAT/LinearInterp.hh>
// ROOT headers
#include <TGraph.h>
 
#include <vector>
 
namespace RAT {
 
// Forward declarations within the namespace
class CCCrossSecMessenger;
 
class CCCrossSec {
 public:
  enum NuEType { nue, nuebar, numu, numubar };
 
  CCCrossSec(const char *flavor = "nue");
 
  ~CCCrossSec();
 
  // Set's the defaults for the calculation
  void Defaults();
 
  double Sigma(const double Enu) const;
 
  std::vector<double> CalcAllowedElectronKE(const double Enu) const;
 
  std::vector<double> CalcAllowedNuclearEx(const double Enu) const;
 
  std::vector<double> GetAllowedTransitionTypes(const double Enu) const;
 
  std::vector<double> CalcdSigmadTNorms(const double Enu) const;
 
  double CrossSecNorm() const { return 1e-42; };
 
 private:
  NuEType fReaction;         
  std::string fReactionStr;  
 
  // Some constants
  static const double fGf;      
  static const double fhbarc;   
  static const double fhbarc2;  
  static const double falpha;   
 
  static double fsinthetaW2;
 
  double fMe;  
 
  std::vector<double> fLevels;      // Energy of level transitions (mass change + excitation energy)
  std::vector<double> fLevelTypes;  // Type of transition (F = 0, GT = 1)
  std::vector<double> fNorms;       // E_e indepedent normalization factors for each transition
 
  CCCrossSecMessenger *fMessenger;
};
 
}  // namespace RAT
 
#endif