ESCrossSec.hh

 
#ifndef __RAT_ESCrossSec__
#define __RAT_ESCrossSec__
 
// G4 headers
#include <globals.hh>
// RAT headers
#include <RAT/GLG4StringUtil.hh>
#include <RAT/LinearInterp.hh>
// ROOT headers
#include <TGraph.h>
 
namespace RAT {
 
// Forward declarations within the namespace
class ESCrossSecMessenger;
 
class ESCrossSec {
 public:
  enum NuEType { nue, nuebar, numu, numubar };
 
  ESCrossSec(const char *flavor = "nue");
 
  ~ESCrossSec();
 
  double Sigma(const double Enu) const;
 
  double dSigmadT(const double Enu, const double Te) const;
 
  double IntegraldSigmadT(const double Enu, const double T1, const double T2) const;
 
  double dSigmadCosTh(const double Enu, const double CosTh) const;
 
  //-----------------------------------------------------------------------
  void SetRadiativeCorrection(const int ii);
 
  inline int GetRadiativeCorrection() const { return fRadiativeCorrection; };
 
  void SetSinThetaW(const double &sintw);
 
  double GetSinThetaW() const { return fsinthetaW2; };
 
  void SetReaction(const std::string &reaction);
 
  const std::string &GetReaction() const { return fReactionStr; };
 
  TGraph *DrawdSigmadT(const double Enu) const;
 
  double CrossSecNorm() const { return 1e-42; };
 
 protected:
  void LoadTablesDB();
 
  void CalcG();
 
 private:
  // Sets the defaults for the calculation
  void Defaults();
 
  // add by y.t. 16-JAN-2003
  double fL(const double x) const;
 
  NuEType fReaction;         
  std::string fReactionStr;  
 
  double fEmin, fEmax;  
 
  // Some constants
  static const double fGf;      
  static const double fhbarc;   
  static const double fhbarc2;  
  static const double falpha;   
 
  static double fsinthetaW2;
 
  double fMe;           
  double fSigmaOverMe;  
 
  double fgL;
  double fgR;
 
  //-----------------------------------------------------------------------
  int fRadiativeCorrection;  // flag to use radiative correction or not
 
  static const double fTableTeMin;
  static const double fTableTeMax;
 
  // Total cross section table
  bool fExistTableTot;
  int fNDataTot;
  double fEnuStepTot;
  LinearInterp<double> fTableTot_e;
 
  // Differential cross section table
  bool fExistTableDif;
  int fNDataDif;
  double fEnuStepDif;
  LinearInterp<double> fTableDif_e;
  std::vector<double> fTableDif;
 
  ESCrossSecMessenger *fMessenger;
};
 
//
// Some inline definitions to make the code a bit faster
//
 
inline void ESCrossSec::CalcG() {
  // calculate the gL and gR for
  // the reaction type
  if (fReaction == nue) {
    fgL = 0.5 + fsinthetaW2;
    fgR = fsinthetaW2;
  } else if (fReaction == nuebar) {
    fgL = fsinthetaW2;
    fgR = 0.5 + fsinthetaW2;
  } else if (fReaction == numu) {
    fgL = -0.5 + fsinthetaW2;
    fgR = fsinthetaW2;
  } else if (fReaction == numubar) {
    fgL = fsinthetaW2;
    fgR = -0.5 + fsinthetaW2;
  } else {
    throw std::invalid_argument("Unknown reaction " + fReactionStr);
  }
}
}  // namespace RAT
 
#endif