CCgen.hh

#ifndef __RAT_CCgen__
#define __RAT_CCgen__
 
 
#include <CLHEP/Random/RandGeneral.h>
#include <CLHEP/Vector/LorentzVector.h>
#include <TF1.h>
 
#include <G4LorentzVector.hh>
#include <G4ThreeVector.hh>
#include <RAT/LinearInterp.hh>
 
class TGraph;
class RandGeneral;
 
namespace RAT {
 
class CCCrossSec;
 
class CCgen {
 public:
  CCgen();
  ~CCgen();
 
  // Generate random event vectors
  //    Pass in the neutrino direction (unit vector)
  //    Returns 4-momentum vectors for resulting electron.
  void GenerateEvent(const G4ThreeVector &nu_dir, G4LorentzVector &neutrino, G4LorentzVector &electron,
                     double &e_nucleus);
 
  void SetNuType(const G4String &nutype);
 
  inline const G4String &GetNuType() const { return fNuFlavor; };
 
  void SetNuFlavor(const G4String &nuflavor);
  inline const G4String &GetNuFlavor() const { return fNuFlavor; };
 
  inline G4double GetTotalFlux() {
    // Return the neutrino flux as the value loaded from the database
    return fTotalFlux;
  }
 
  G4double GetRatePerTarget();
 
  const G4String GetDBName() const { return fDBName; }
 
  void SetDBName(const G4String name);
 
 private:
  void LoadGenerator();
 
  void GenInteraction(double &Enu, double &CosThetaLab);
 
  void Reset();
  void Show();
 
  G4double SampleNuEnergy();
 
  G4double SampleRecoilEnergy(G4double Enu, int &Transition, double &Enucleus);
 
  G4double SampleRecoilAngle(G4double Enu, G4double Te, int Transition);
 
 protected:
  inline G4bool GetGenLoaded() { return fGenLoaded; };
 
  G4String fGenType;
  G4String fNuType;
 
  G4String fNuFlavor;
 
  CCCrossSec *fXS;
 
  TGraph *fNuSpectrum;
 
  TF1 *fFermiAngle;
 
  TF1 *fGTAngle;
 
  std::vector<double> fEnuTbl;
 
  std::vector<double> fFluxTbl;
 
  G4double fEmax;
  G4double fEmin;
 
  G4double fEnuMax;
  G4double fEnuMin;
 
  G4double fFluxMax;
 
  G4bool fGenLoaded;
 
  G4double fMassElectron;
 
  G4double fTotalFlux;
 
  CLHEP::RandGeneral *fSpectrumRndm;
 
  G4String fDBName;
};
 
}  // namespace RAT
 
#endif