genie::INukeHadroData2018 Class Reference

#include <INukeHadroData2018.h>

Collaboration diagram for genie::INukeHadroData2018:

Classes
struct	Cleaner

Public Member Functions
double	XSec (int hpdgc, int tgt, int nprod, INukeFateHN_t rxnType, double ke, double costh) const
double	XSec (int hpdgc, INukeFateHN_t fate, double ke, int targA, int targZ) const
double	FracADep (int hpdgc, INukeFateHA_t fate, double ke, int targA) const
double	FracAIndep (int hpdgc, INukeFateHA_t fate, double ke) const
double	Frac (int hpdgc, INukeFateHN_t fate, double ke, int targA=0, int targZ=0) const
double	IntBounce (const GHepParticle *p, int target, int s1, INukeFateHN_t fate)
const Spline *	XSecPipn_Tot (void) const
const Spline *	XSecPipn_CEx (void) const
const Spline *	XSecPipn_Elas (void) const
const Spline *	XSecPipn_Reac (void) const
const Spline *	XSecPipp_Tot (void) const
const Spline *	XSecPipp_CEx (void) const
const Spline *	XSecPipp_Elas (void) const
const Spline *	XSecPipp_Reac (void) const
const Spline *	XSecPipp_Abs (void) const
const Spline *	XSecPi0n_Tot (void) const
const Spline *	XSecPi0n_CEx (void) const
const Spline *	XSecPi0n_Elas (void) const
const Spline *	XSecPi0n_Reac (void) const
const Spline *	XSecPi0p_Tot (void) const
const Spline *	XSecPi0p_CEx (void) const
const Spline *	XSecPi0p_Elas (void) const
const Spline *	XSecPi0p_Reac (void) const
const Spline *	XSecPi0p_Abs (void) const
const Spline *	XSecPp_Tot (void) const
const Spline *	XSecPp_Elas (void) const
const Spline *	XSecPp_Reac (void) const
const Spline *	XSecPn_Tot (void) const
const Spline *	XSecPn_Elas (void) const
const Spline *	XSecPn_Reac (void) const
const Spline *	XSecNn_Tot (void) const
const Spline *	XSecNn_Elas (void) const
const Spline *	XSecNn_Reac (void) const
const Spline *	XSecKpn_Elas (void) const
const Spline *	XSecKpn_CEx (void) const
const Spline *	XSecKpp_Elas (void) const
const Spline *	XSecKpN_Abs (void) const
const Spline *	XSecKpN_Tot (void) const
const Spline *	XSecGamp_fs (void) const
const Spline *	XSecGamn_fs (void) const
const Spline *	XSecGamN_Tot (void) const
const Spline *	FracPA_Tot (void) const
const Spline *	FracPA_Elas (void) const
const Spline *	FracPA_Inel (void) const
const Spline *	FracPA_CEx (void) const
const Spline *	FracPA_Abs (void) const
const Spline *	FracPA_PiPro (void) const
const Spline *	FracNA_Tot (void) const
const Spline *	FracNA_Elas (void) const
const Spline *	FracNA_Inel (void) const
const Spline *	FracNA_CEx (void) const
const Spline *	FracNA_Abs (void) const
const Spline *	FracNA_PiPro (void) const
const Spline *	FracKA_Tot (void) const
const Spline *	FracKA_Elas (void) const
const Spline *	FracKA_Inel (void) const
const Spline *	FracKA_CEx (void) const
const Spline *	FracKA_Abs (void) const
const Spline *	FracPA_Cmp (void) const
const Spline *	FracNA_Cmp (void) const
const BLI2DNonUnifGrid *	hN2dXSecPP_Elas (void) const
const BLI2DNonUnifGrid *	hN2dXSecNP_Elas (void) const
const BLI2DNonUnifGrid *	hN2dXSecPipN_Elas (void) const
const BLI2DNonUnifGrid *	hN2dXSecPi0N_Elas (void) const
const BLI2DNonUnifGrid *	hN2dXSecPimN_Elas (void) const
const BLI2DNonUnifGrid *	hN2dXSecKpN_Elas (void) const
const BLI2DNonUnifGrid *	hN2dXSecKpP_Elas (void) const
const BLI2DNonUnifGrid *	hN2dXSecPiN_CEx (void) const
const BLI2DNonUnifGrid *	hN2dXSecPiN_Abs (void) const
const BLI2DNonUnifGrid *	hN2dXSecGamPi0P_Inelas (void) const
const BLI2DNonUnifGrid *	hN2dXSecGamPi0N_Inelas (void) const
const BLI2DNonUnifGrid *	hN2dXSecGamPipN_Inelas (void) const
const BLI2DNonUnifGrid *	hN2dXSecGamPimP_Inelas (void) const

Static Public Member Functions
static INukeHadroData2018 *	Instance (void)

Static Public Attributes
static double	fMinKinEnergy = 1.0
static double	fMaxKinEnergyHA = 999.0
static double	fMaxKinEnergyHN = 1799.0

Private Member Functions
	INukeHadroData2018 ()
	INukeHadroData2018 (const INukeHadroData2018 &shx)
	~INukeHadroData2018 ()
void	LoadCrossSections (void)
void	ReadhNFile (string filename, double ke, int npoints, int &curr_point, double *costh_array, double *xsec_array, int cols)

Private Attributes
Spline *	fXSecPipn_Tot
	pi+n hN x-section splines
Spline *	fXSecPipn_CEx
Spline *	fXSecPipn_Elas
Spline *	fXSecPipn_Reac
Spline *	fXSecPipp_Tot
	pi+p hN x-section splines
Spline *	fXSecPipp_CEx
Spline *	fXSecPipp_Elas
Spline *	fXSecPipp_Reac
Spline *	fXSecPipd_Abs
Spline *	fXSecPi0n_Tot
	pi0n hN x-section splines
Spline *	fXSecPi0n_CEx
Spline *	fXSecPi0n_Elas
Spline *	fXSecPi0n_Reac
Spline *	fXSecPi0p_Tot
	pi0p hN x-section splines
Spline *	fXSecPi0p_CEx
Spline *	fXSecPi0p_Elas
Spline *	fXSecPi0p_Reac
Spline *	fXSecPi0d_Abs
Spline *	fXSecPp_Tot
	p/nN x-section splines
Spline *	fXSecPp_Elas
Spline *	fXSecPp_Reac
Spline *	fXSecPn_Tot
Spline *	fXSecPn_Elas
Spline *	fXSecPn_Reac
Spline *	fXSecNn_Tot
Spline *	fXSecNn_Elas
Spline *	fXSecNn_Reac
Spline *	fXSecKpn_Elas
	K+N x-section splines.
Spline *	fXSecKpp_Elas
Spline *	fXSecKpn_CEx
Spline *	fXSecKpN_Abs
Spline *	fXSecKpN_Tot
Spline *	fFracPA_Tot
	N+A x-section splines.
Spline *	fFracPA_Elas
Spline *	fFracPA_Inel
Spline *	fFracPA_CEx
Spline *	fFracPA_Abs
Spline *	fFracPA_PiPro
Spline *	fFracNA_Tot
Spline *	fFracNA_Elas
Spline *	fFracNA_Inel
Spline *	fFracNA_CEx
Spline *	fFracNA_Abs
Spline *	fFracNA_PiPro
Spline *	fFracKA_Tot
	K+A x-section splines.
Spline *	fFracKA_Elas
Spline *	fFracKA_CEx
Spline *	fFracKA_Inel
Spline *	fFracKA_Abs
Spline *	fXSecGamp_fs
	gamma A x-section splines
Spline *	fXSecGamn_fs
Spline *	fXSecGamN_Tot
Spline *	fXSecPp_Cmp
	NN cmp (compound nucleus) fate.
Spline *	fXSecPn_Cmp
Spline *	fXSecNn_Cmp
Spline *	fFracPA_Cmp
Spline *	fFracNA_Cmp
TGraph2D *	TfracPipA_CEx
TGraph2D *	TfracPipA_Inelas
TGraph2D *	TfracPipA_Abs
TGraph2D *	TfracPipA_PiPro
BLI2DNonUnifGrid *	fhN2dXSecPP_Elas
BLI2DNonUnifGrid *	fhN2dXSecNP_Elas
BLI2DNonUnifGrid *	fhN2dXSecPipN_Elas
BLI2DNonUnifGrid *	fhN2dXSecPi0N_Elas
BLI2DNonUnifGrid *	fhN2dXSecPimN_Elas
BLI2DNonUnifGrid *	fhN2dXSecKpN_Elas
BLI2DNonUnifGrid *	fhN2dXSecKpP_Elas
BLI2DNonUnifGrid *	fhN2dXSecKpN_CEx
BLI2DNonUnifGrid *	fhN2dXSecKpN_Abs
BLI2DNonUnifGrid *	fhN2dXSecPiN_CEx
BLI2DNonUnifGrid *	fhN2dXSecPiN_Abs
BLI2DNonUnifGrid *	fhN2dXSecGamPi0P_Inelas
BLI2DNonUnifGrid *	fhN2dXSecGamPi0N_Inelas
BLI2DNonUnifGrid *	fhN2dXSecGamPipN_Inelas
BLI2DNonUnifGrid *	fhN2dXSecGamPimP_Inelas

Static Private Attributes
static INukeHadroData2018 *	fInstance = 0

Friends
struct	Cleaner

Detailed Description

Definition at line 51 of file INukeHadroData2018.h.

Constructor & Destructor Documentation

INukeHadroData2018() [1/2]

INukeHadroData2018::INukeHadroData2018 ( )

private

Definition at line 70 of file INukeHadroData2018.cxx.

{
  this->LoadCrossSections();
  fInstance = 0;
}

References fInstance, and LoadCrossSections().

Referenced by Instance(), and INukeHadroData2018().

INukeHadroData2018() [2/2]

genie::INukeHadroData2018::INukeHadroData2018 ( const INukeHadroData2018 & shx )

private

References INukeHadroData2018().

~INukeHadroData2018()

INukeHadroData2018::~INukeHadroData2018 ( )

private

Definition at line 76 of file INukeHadroData2018.cxx.

{
  
  // pi+n/p hA x-section splines
  delete fXSecPipn_Tot;
  delete fXSecPipn_CEx;
  delete fXSecPipn_Elas;
  delete fXSecPipn_Reac;
  delete fXSecPipp_Tot;
  delete fXSecPipp_CEx;
  delete fXSecPipp_Elas;
  delete fXSecPipp_Reac;
  delete fXSecPipd_Abs;
 
  delete fXSecPp_Cmp; //added to fix memory leak; no noticeable effect, but good convention.
  delete fXSecPn_Cmp;
  delete fXSecNn_Cmp;
  delete fXSecPp_Tot;
  delete fXSecPp_Elas;
  delete fXSecPp_Reac;
  delete fXSecPn_Tot;
  delete fXSecPn_Elas;
  delete fXSecPn_Reac;
  delete fXSecNn_Tot;
  delete fXSecNn_Elas;
  delete fXSecNn_Reac;
 
  // pi0n/p hA x-section splines
  delete fXSecPi0n_Tot;
  delete fXSecPi0n_CEx;
  delete fXSecPi0n_Elas;
  delete fXSecPi0n_Reac;
  delete fXSecPi0p_Tot;
  delete fXSecPi0p_CEx;
  delete fXSecPi0p_Elas;
  delete fXSecPi0p_Reac;
  delete fXSecPi0d_Abs;
 
  // K+N x-section splines
  delete fXSecKpn_Elas;
  delete fXSecKpp_Elas;
  delete fXSecKpn_CEx;
  delete fXSecKpN_Abs;
  delete fXSecKpN_Tot;
 
  // gamma x-section splines (inelastic only)
  delete fXSecGamp_fs;
  delete fXSecGamn_fs;
  delete fXSecGamN_Tot;
 
  // N+A x-section splines
  delete fFracPA_Tot;
  //  delete fFracPA_Elas;
  delete fFracPA_Inel;
  delete fFracPA_CEx;
  delete fFracPA_Abs;
  delete fFracPA_PiPro;
  delete fFracNA_Tot;
  //  delete fFracNA_Elas;
  delete fFracNA_Inel;
  delete fFracNA_CEx;
  delete fFracNA_Abs;
  delete fFracNA_PiPro;
 
  delete fFracPA_Cmp;  //  cmp - add support later
  delete fFracNA_Cmp;
 
  // hN data
  delete fhN2dXSecPP_Elas;
  delete fhN2dXSecNP_Elas;
  delete fhN2dXSecPipN_Elas;
  delete fhN2dXSecPi0N_Elas;
  delete fhN2dXSecPimN_Elas;
  delete fhN2dXSecKpN_Elas;
  delete fhN2dXSecKpP_Elas;
  delete fhN2dXSecPiN_CEx;
  delete fhN2dXSecPiN_Abs;
  delete fhN2dXSecGamPi0P_Inelas;
  delete fhN2dXSecGamPi0N_Inelas;
  delete fhN2dXSecGamPipN_Inelas;
  delete fhN2dXSecGamPimP_Inelas;
  delete fhN2dXSecKpN_CEx;
  
//  delete TPipA_Tot; 
  delete TfracPipA_Abs;
  delete TfracPipA_CEx;
//  delete TfracPipA_Elas;
  delete TfracPipA_Inelas;
  delete TfracPipA_PiPro;
  
  // K+A x-section fraction splines
  delete fFracKA_Tot;
  delete fFracKA_Elas;
  delete fFracKA_CEx;
  delete fFracKA_Inel;
  delete fFracKA_Abs;
  
}

References fFracKA_Abs, fFracKA_CEx, fFracKA_Elas, fFracKA_Inel, fFracKA_Tot, fFracNA_Abs, fFracNA_CEx, fFracNA_Cmp, fFracNA_Inel, fFracNA_PiPro, fFracNA_Tot, fFracPA_Abs, fFracPA_CEx, fFracPA_Cmp, fFracPA_Inel, fFracPA_PiPro, fFracPA_Tot, fhN2dXSecGamPi0N_Inelas, fhN2dXSecGamPi0P_Inelas, fhN2dXSecGamPimP_Inelas, fhN2dXSecGamPipN_Inelas, fhN2dXSecKpN_CEx, fhN2dXSecKpN_Elas, fhN2dXSecKpP_Elas, fhN2dXSecNP_Elas, fhN2dXSecPi0N_Elas, fhN2dXSecPimN_Elas, fhN2dXSecPiN_Abs, fhN2dXSecPiN_CEx, fhN2dXSecPipN_Elas, fhN2dXSecPP_Elas, fXSecGamn_fs, fXSecGamN_Tot, fXSecGamp_fs, fXSecKpN_Abs, fXSecKpn_CEx, fXSecKpn_Elas, fXSecKpN_Tot, fXSecKpp_Elas, fXSecNn_Cmp, fXSecNn_Elas, fXSecNn_Reac, fXSecNn_Tot, fXSecPi0d_Abs, fXSecPi0n_CEx, fXSecPi0n_Elas, fXSecPi0n_Reac, fXSecPi0n_Tot, fXSecPi0p_CEx, fXSecPi0p_Elas, fXSecPi0p_Reac, fXSecPi0p_Tot, fXSecPipd_Abs, fXSecPipn_CEx, fXSecPipn_Elas, fXSecPipn_Reac, fXSecPipn_Tot, fXSecPipp_CEx, fXSecPipp_Elas, fXSecPipp_Reac, fXSecPipp_Tot, fXSecPn_Cmp, fXSecPn_Elas, fXSecPn_Reac, fXSecPn_Tot, fXSecPp_Cmp, fXSecPp_Elas, fXSecPp_Reac, fXSecPp_Tot, TfracPipA_Abs, TfracPipA_CEx, TfracPipA_Inelas, and TfracPipA_PiPro.

Member Function Documentation

Frac()

double INukeHadroData2018::Frac	(	int	hpdgc,
		INukeFateHN_t	fate,
		double	ke,
		int	targA = 0,
		int	targZ = 0 ) const

Definition at line 1636 of file INukeHadroData2018.cxx.

{
// return the x-section fraction for the input fate for the particle with the
// input pdg code at the input kinetic energy
 
  ke = TMath::Max(fMinKinEnergy,   ke);
  ke = TMath::Min(fMaxKinEnergyHN, ke);
 
  // get x-section
  double xsec = this->XSec(hpdgc,fate,ke,targA,targZ);
 
  // get max x-section
  double xsec_tot = 0;
       if (hpdgc == kPdgPiP    ){xsec_tot = TMath::Max(0., fXSecPipp_Tot  -> Evaluate(ke)) *  targZ;
                                 xsec_tot+= TMath::Max(0., fXSecPipn_Tot  -> Evaluate(ke)) * (targA-targZ);}
  else if (hpdgc == kPdgPiM    ){xsec_tot = TMath::Max(0., fXSecPipn_Tot  -> Evaluate(ke)) *  targZ;
                                 xsec_tot+= TMath::Max(0., fXSecPipp_Tot  -> Evaluate(ke)) * (targA-targZ);}
  else if (hpdgc == kPdgPi0    ){xsec_tot = TMath::Max(0., fXSecPi0p_Tot  -> Evaluate(ke)) *  targZ;
                                 xsec_tot+= TMath::Max(0., fXSecPi0n_Tot  -> Evaluate(ke)) * (targA-targZ);}
  else if (hpdgc == kPdgProton ){xsec_tot = TMath::Max(0., fXSecPp_Tot    -> Evaluate(ke)) *  targZ;
                                 xsec_tot+= TMath::Max(0., fXSecPn_Tot    -> Evaluate(ke)) * (targA-targZ);}
  else if (hpdgc == kPdgNeutron){xsec_tot = TMath::Max(0., fXSecPn_Tot    -> Evaluate(ke)) *  targZ;
                                 xsec_tot+= TMath::Max(0., fXSecNn_Tot    -> Evaluate(ke)) * (targA-targZ);}
  else if (hpdgc == kPdgGamma  ) xsec_tot = TMath::Max(0., fXSecGamN_Tot  -> Evaluate(ke));
  else if (hpdgc == kPdgKP     ) xsec_tot = TMath::Max(0., fXSecKpN_Tot   -> Evaluate(ke));
 
  // compute fraction
  double frac = (xsec_tot>0) ? xsec/xsec_tot : 0.;
  return frac;
}

References fMaxKinEnergyHN, fMinKinEnergy, fXSecGamN_Tot, fXSecKpN_Tot, fXSecNn_Tot, fXSecPi0n_Tot, fXSecPi0p_Tot, fXSecPipn_Tot, fXSecPipp_Tot, fXSecPn_Tot, fXSecPp_Tot, genie::kPdgGamma, genie::kPdgKP, genie::kPdgNeutron, genie::kPdgPi0, genie::kPdgPiM, genie::kPdgPiP, genie::kPdgProton, and XSec().

FracADep()

double INukeHadroData2018::FracADep	(	int	hpdgc,
		INukeFateHA_t	fate,
		double	ke,
		int	targA ) const

Definition at line 1377 of file INukeHadroData2018.cxx.

{
  // return the x-section fraction for the input fate for the particle with the input pdg
  // code and the target with the input mass number at the input kinetic energy
 
  ke = TMath::Max(fMinKinEnergy,   ke);  // ke >= 1 MeV
  ke = TMath::Min(fMaxKinEnergyHA, ke);  // ke <= 999 MeV
 
  targA = TMath::Min(208, targA);  // A <= 208
 
  LOG("INukeData", pDEBUG)  << "Querying hA cross section at ke  = " << ke << " and target " << targA;
 
  // Handle pions (currently the same cross sections are used for pi+, pi-, and pi0)
  if ( hpdgc == kPdgPiP || hpdgc == kPdgPiM || hpdgc == kPdgPi0 ) {
 
    double frac_cex = TfracPipA_CEx->Interpolate(targA, ke);
    //double frac_elas = TfracPipA_Elas->Interpolate(targA, ke);
    double frac_inelas = TfracPipA_Inelas->Interpolate(targA, ke);
    double frac_abs = TfracPipA_Abs->Interpolate(targA, ke);
    double frac_pipro = TfracPipA_PiPro->Interpolate(targA, ke);
 
    // Protect against unitarity violation due to interpolation problems
    // by renormalizing all available fate fractions to unity.
    double total = frac_cex + frac_inelas + frac_abs + frac_pipro; // + frac_elas
 
    if ( fate == kIHAFtCEx ) return frac_cex / total;
  //else if ( fate == kIHAFtElas   ) return frac_elas / total;
    else if ( fate == kIHAFtInelas ) return frac_inelas / total;
    else if ( fate == kIHAFtAbs    ) return frac_abs / total;
    else if ( fate == kIHAFtPiProd ) return frac_pipro / total;
    else {
      std::string sign("+");
      if ( hpdgc == kPdgPiM ) sign = "-";
      else if ( hpdgc == kPdgPi0 ) sign = "0";
      LOG("INukeData", pWARN) << "Pi" << sign << "'s don't have this fate: " << INukeHadroFates::AsString(fate);
      return 0.;
    }
  }
 
  LOG("INukeData", pWARN) << "Can't handle particles with pdg code = " << hpdgc;
  return 0.;
}

References genie::INukeHadroFates::AsString(), fMaxKinEnergyHA, fMinKinEnergy, genie::kIHAFtAbs, genie::kIHAFtCEx, genie::kIHAFtInelas, genie::kIHAFtPiProd, genie::kPdgPi0, genie::kPdgPiM, genie::kPdgPiP, LOG, pDEBUG, pWARN, TfracPipA_Abs, TfracPipA_CEx, TfracPipA_Inelas, and TfracPipA_PiPro.

FracAIndep()

double INukeHadroData2018::FracAIndep	(	int	hpdgc,
		INukeFateHA_t	fate,
		double	ke ) const

Definition at line 1420 of file INukeHadroData2018.cxx.

{
  // return the x-section fraction for the input fate for the particle with the input pdg
  // code at the input kinetic energy
  ke = TMath::Max(fMinKinEnergy,   ke);
  ke = TMath::Min(fMaxKinEnergyHA, ke);
 
  LOG("INukeData", pDEBUG)  << "Querying hA cross section at ke = " << ke;
 
  // TODO: reduce code duplication here
  if (hpdgc == kPdgProton) {
    // handle protons
    double frac_cex = fFracPA_CEx->Evaluate(ke);
    double frac_inelas = fFracPA_Inel->Evaluate(ke);
    double frac_abs = fFracPA_Abs->Evaluate(ke);
    double frac_pipro = fFracPA_PiPro->Evaluate(ke);
    double frac_comp = fFracPA_Cmp->Evaluate(ke);
 
    // Protect against unitarity violation due to interpolation problems
    // by renormalizing all available fate fractions to unity.
    double total = frac_cex + frac_inelas + frac_abs + frac_pipro + frac_comp; // + frac_elas
 
    if ( fate == kIHAFtCEx ) return frac_cex / total;
  //else if ( fate == kIHAFtElas   ) return frac_elas / total;
    else if ( fate == kIHAFtInelas ) return frac_inelas / total;
    else if ( fate == kIHAFtAbs    ) return frac_abs / total;
    else if ( fate == kIHAFtPiProd ) return frac_pipro / total;
    else if ( fate == kIHAFtCmp    ) return frac_comp / total; // cmp - add support for this later
    else {
      LOG("INukeData", pWARN)
        << "Protons don't have this fate: " << INukeHadroFates::AsString(fate);
      return 0;
    }
  }
  else if (hpdgc == kPdgNeutron) {
    // handle neutrons
    double frac_cex = fFracNA_CEx->Evaluate(ke);
    double frac_inelas = fFracNA_Inel->Evaluate(ke);
    double frac_abs = fFracNA_Abs->Evaluate(ke);
    double frac_pipro = fFracNA_PiPro->Evaluate(ke);
    double frac_comp = fFracNA_Cmp->Evaluate(ke);
 
    // Protect against unitarity violation due to interpolation problems
    // by renormalizing all available fate fractions to unity.
    double total = frac_cex + frac_inelas + frac_abs + frac_pipro + frac_comp; // + frac_elas
 
    if ( fate == kIHAFtCEx ) return frac_cex / total;
  //else if ( fate == kIHAFtElas   ) return frac_elas / total;
    else if ( fate == kIHAFtInelas ) return frac_inelas / total;
    else if ( fate == kIHAFtAbs    ) return frac_abs / total;
    else if ( fate == kIHAFtPiProd ) return frac_pipro / total;
    else if ( fate == kIHAFtCmp    ) return frac_comp / total; // cmp - add support for this later
    else {
      LOG("INukeData", pWARN)
        << "Neutrons don't have this fate: " << INukeHadroFates::AsString(fate);
      return 0;
    }
  }
  else if (hpdgc == kPdgKP) {
    // handle K+
    double frac_inelas = fFracKA_Inel->Evaluate(ke);
  //double frac_elas = fFracKA_Elas->Evaluate(ke);
    double frac_abs = fFracKA_Abs->Evaluate(ke);
 
    // Protect against unitarity violation due to interpolation problems
    // by renormalizing all available fate fractions to unity.
    double total = frac_inelas + frac_abs; // + frac_elas
 
    if ( fate == kIHAFtInelas ) return frac_inelas / total;
    else if ( fate == kIHAFtAbs ) return frac_abs / total;
    else {
      LOG("INukeData", pWARN)
        << "K+'s don't have this fate: " << INukeHadroFates::AsString(fate);
      return 0.;
    }
  }
  LOG("INukeData", pWARN) << "Can't handle particles with pdg code = " << hpdgc;
  return 0.;
}

References genie::INukeHadroFates::AsString(), fFracKA_Abs, fFracKA_Inel, fFracNA_Abs, fFracNA_CEx, fFracNA_Cmp, fFracNA_Inel, fFracNA_PiPro, fFracPA_Abs, fFracPA_CEx, fFracPA_Cmp, fFracPA_Inel, fFracPA_PiPro, fMaxKinEnergyHA, fMinKinEnergy, genie::kIHAFtAbs, genie::kIHAFtCEx, genie::kIHAFtCmp, genie::kIHAFtInelas, genie::kIHAFtPiProd, genie::kPdgKP, genie::kPdgNeutron, genie::kPdgProton, LOG, pDEBUG, and pWARN.

FracKA_Abs()

const Spline * genie::INukeHadroData2018::FracKA_Abs ( void ) const

inline

Definition at line 125 of file INukeHadroData2018.h.

{ return fFracKA_Abs;       }

References fFracKA_Abs.

FracKA_CEx()

const Spline * genie::INukeHadroData2018::FracKA_CEx ( void ) const

inline

Definition at line 124 of file INukeHadroData2018.h.

{ return fFracKA_CEx;       }

References fFracKA_CEx.

FracKA_Elas()

const Spline * genie::INukeHadroData2018::FracKA_Elas ( void ) const

inline

Definition at line 122 of file INukeHadroData2018.h.

{ return fFracKA_Elas;      }

References fFracKA_Elas.

FracKA_Inel()

const Spline * genie::INukeHadroData2018::FracKA_Inel ( void ) const

inline

Definition at line 123 of file INukeHadroData2018.h.

{ return fFracKA_Inel;      }

References fFracKA_Inel.

FracKA_Tot()

const Spline * genie::INukeHadroData2018::FracKA_Tot ( void ) const

inline

Definition at line 121 of file INukeHadroData2018.h.

{ return fFracKA_Tot;       }

References fFracKA_Tot.

FracNA_Abs()

const Spline * genie::INukeHadroData2018::FracNA_Abs ( void ) const

inline

Definition at line 119 of file INukeHadroData2018.h.

{ return fFracNA_Abs;       }

References fFracNA_Abs.

FracNA_CEx()

const Spline * genie::INukeHadroData2018::FracNA_CEx ( void ) const

inline

Definition at line 118 of file INukeHadroData2018.h.

{ return fFracNA_CEx;       }

References fFracNA_CEx.

FracNA_Cmp()

const Spline * genie::INukeHadroData2018::FracNA_Cmp ( void ) const

inline

Definition at line 128 of file INukeHadroData2018.h.

{ return fFracNA_Cmp;      } //suarez

References fFracNA_Cmp.

FracNA_Elas()

const Spline * genie::INukeHadroData2018::FracNA_Elas ( void ) const

inline

Definition at line 116 of file INukeHadroData2018.h.

{ return fFracNA_Elas;      }

References fFracNA_Elas.

FracNA_Inel()

const Spline * genie::INukeHadroData2018::FracNA_Inel ( void ) const

inline

Definition at line 117 of file INukeHadroData2018.h.

{ return fFracNA_Inel;      }

References fFracNA_Inel.

FracNA_PiPro()

const Spline * genie::INukeHadroData2018::FracNA_PiPro ( void ) const

inline

Definition at line 120 of file INukeHadroData2018.h.

{ return fFracNA_PiPro;     }

References fFracNA_PiPro.

FracNA_Tot()

const Spline * genie::INukeHadroData2018::FracNA_Tot ( void ) const

inline

Definition at line 115 of file INukeHadroData2018.h.

{ return fFracNA_Tot;       }

References fFracNA_Tot.

FracPA_Abs()

const Spline * genie::INukeHadroData2018::FracPA_Abs ( void ) const

inline

Definition at line 113 of file INukeHadroData2018.h.

{ return fFracPA_Abs;       }

References fFracPA_Abs.

FracPA_CEx()

const Spline * genie::INukeHadroData2018::FracPA_CEx ( void ) const

inline

Definition at line 112 of file INukeHadroData2018.h.

{ return fFracPA_CEx;       }

References fFracPA_CEx.

FracPA_Cmp()

const Spline * genie::INukeHadroData2018::FracPA_Cmp ( void ) const

inline

Definition at line 127 of file INukeHadroData2018.h.

{ return fFracPA_Cmp;      }

References fFracPA_Cmp.

FracPA_Elas()

const Spline * genie::INukeHadroData2018::FracPA_Elas ( void ) const

inline

Definition at line 110 of file INukeHadroData2018.h.

{ return fFracPA_Elas;      }

References fFracPA_Elas.

FracPA_Inel()

const Spline * genie::INukeHadroData2018::FracPA_Inel ( void ) const

inline

Definition at line 111 of file INukeHadroData2018.h.

{ return fFracPA_Inel;      }

References fFracPA_Inel.

FracPA_PiPro()

const Spline * genie::INukeHadroData2018::FracPA_PiPro ( void ) const

inline

Definition at line 114 of file INukeHadroData2018.h.

{ return fFracPA_PiPro;     }

References fFracPA_PiPro.

FracPA_Tot()

const Spline * genie::INukeHadroData2018::FracPA_Tot ( void ) const

inline

Definition at line 109 of file INukeHadroData2018.h.

{ return fFracPA_Tot;       }

References fFracPA_Tot.

hN2dXSecGamPi0N_Inelas()

const BLI2DNonUnifGrid * genie::INukeHadroData2018::hN2dXSecGamPi0N_Inelas ( void ) const

inline

Definition at line 141 of file INukeHadroData2018.h.

{ return fhN2dXSecGamPi0N_Inelas;   }

References fhN2dXSecGamPi0N_Inelas.

hN2dXSecGamPi0P_Inelas()

const BLI2DNonUnifGrid * genie::INukeHadroData2018::hN2dXSecGamPi0P_Inelas ( void ) const

inline

Definition at line 140 of file INukeHadroData2018.h.

{ return fhN2dXSecGamPi0P_Inelas;   }

References fhN2dXSecGamPi0P_Inelas.

hN2dXSecGamPimP_Inelas()

const BLI2DNonUnifGrid * genie::INukeHadroData2018::hN2dXSecGamPimP_Inelas ( void ) const

inline

Definition at line 143 of file INukeHadroData2018.h.

{ return fhN2dXSecGamPimP_Inelas;   }

References fhN2dXSecGamPimP_Inelas.

hN2dXSecGamPipN_Inelas()

const BLI2DNonUnifGrid * genie::INukeHadroData2018::hN2dXSecGamPipN_Inelas ( void ) const

inline

Definition at line 142 of file INukeHadroData2018.h.

{ return fhN2dXSecGamPipN_Inelas;   }

References fhN2dXSecGamPipN_Inelas.

hN2dXSecKpN_Elas()

const BLI2DNonUnifGrid * genie::INukeHadroData2018::hN2dXSecKpN_Elas ( void ) const

inline

Definition at line 136 of file INukeHadroData2018.h.

{ return fhN2dXSecKpN_Elas;         }

References fhN2dXSecKpN_Elas.

hN2dXSecKpP_Elas()

const BLI2DNonUnifGrid * genie::INukeHadroData2018::hN2dXSecKpP_Elas ( void ) const

inline

Definition at line 137 of file INukeHadroData2018.h.

{ return fhN2dXSecKpP_Elas;         }

References fhN2dXSecKpP_Elas.

hN2dXSecNP_Elas()

const BLI2DNonUnifGrid * genie::INukeHadroData2018::hN2dXSecNP_Elas ( void ) const

inline

Definition at line 132 of file INukeHadroData2018.h.

{ return fhN2dXSecNP_Elas;          }

References fhN2dXSecNP_Elas.

hN2dXSecPi0N_Elas()

const BLI2DNonUnifGrid * genie::INukeHadroData2018::hN2dXSecPi0N_Elas ( void ) const

inline

Definition at line 134 of file INukeHadroData2018.h.

{ return fhN2dXSecPi0N_Elas;        }

References fhN2dXSecPi0N_Elas.

hN2dXSecPimN_Elas()

const BLI2DNonUnifGrid * genie::INukeHadroData2018::hN2dXSecPimN_Elas ( void ) const

inline

Definition at line 135 of file INukeHadroData2018.h.

{ return fhN2dXSecPimN_Elas;        }

References fhN2dXSecPimN_Elas.

hN2dXSecPiN_Abs()

const BLI2DNonUnifGrid * genie::INukeHadroData2018::hN2dXSecPiN_Abs ( void ) const

inline

Definition at line 139 of file INukeHadroData2018.h.

{ return fhN2dXSecPiN_Abs;          }

References fhN2dXSecPiN_Abs.

hN2dXSecPiN_CEx()

const BLI2DNonUnifGrid * genie::INukeHadroData2018::hN2dXSecPiN_CEx ( void ) const

inline

Definition at line 138 of file INukeHadroData2018.h.

{ return fhN2dXSecPiN_CEx;          }

References fhN2dXSecPiN_CEx.

hN2dXSecPipN_Elas()

const BLI2DNonUnifGrid * genie::INukeHadroData2018::hN2dXSecPipN_Elas ( void ) const

inline

Definition at line 133 of file INukeHadroData2018.h.

{ return fhN2dXSecPipN_Elas;        }

References fhN2dXSecPipN_Elas.

hN2dXSecPP_Elas()

const BLI2DNonUnifGrid * genie::INukeHadroData2018::hN2dXSecPP_Elas ( void ) const

inline

Definition at line 131 of file INukeHadroData2018.h.

{ return fhN2dXSecPP_Elas;          }

References fhN2dXSecPP_Elas.

Instance()

INukeHadroData2018 * INukeHadroData2018::Instance ( void )

static

Definition at line 175 of file INukeHadroData2018.cxx.

{
  if(fInstance == 0) {
    LOG("INukeData", pINFO) << "INukeHadroData2018 late initialization";
    static INukeHadroData2018::Cleaner cleaner;
    cleaner.DummyMethodAndSilentCompiler();
    fInstance = new INukeHadroData2018;
  }
  return fInstance;
}

References genie::INukeHadroData2018::Cleaner::DummyMethodAndSilentCompiler(), fInstance, INukeHadroData2018(), LOG, and pINFO.

Referenced by INukeNucleonCorr::AvgCorrection(), genie::HAIntranuke2018::LoadConfig(), genie::HNIntranuke2018::LoadConfig(), and genie::utils::intranuke2018::MeanFreePath().

IntBounce()

double INukeHadroData2018::IntBounce	(	const GHepParticle *	p,
		int	target,
		int	s1,
		INukeFateHN_t	fate )

Definition at line 1667 of file INukeHadroData2018.cxx.

{
  // This method returns a random cos(ang) according to a distribution
  // based upon the particle and fate. The sampling uses the
  // Accept/Reject method, whereby a distribution is bounded above by
  // an envelope, or in this case, a number of envelopes, which can be
  // easily sampled (here, we use uniform distributions).
  // To get a random value, first the envelope is sampled to
  // obtain an x-coordinate (cos(ang)), and then another random value
  // is obtained uniformally in the range [0,h(j,0)], where h(j,0)
  // is the height of the j-th envelope. If the point is beneath the
  // distribution, the x-coordinate is accepted, otherwise, we try
  // again.
 
  RandomGen * rnd = RandomGen::Instance();
 
  // numEnv is the number of envelopes in the total envelope,
  // that is, the number of seperate simple uniform distributions
  // that will be fit against the distribution in question in the
  // Accept/Reject process of sampling
  int    numEnv = 4;
  int    numPoints = 1000;            // The number of points to be evaluated
                                      //   for the purpose of finding the max
                                      //   value of the distribution
  assert((numPoints%numEnv)==0);      // numPoints/numEnv has to be an integer
  double sr = 2.0 / numEnv;           // Subrange, i.e., range of an envelope
  double cstep = 2.0 / (numPoints);   // Magnitude of the step between eval. points
 
  double ke = (p->E() - p->Mass()) * 1000.0; // ke in MeV
  if (TMath::Abs((int)ke-ke)<.01) ke+=.3;    // make sure ke isn't an integer,
                                             // otherwise sometimes gives weird results
                                             // due to ROOT's Interpolate() function
  double avg = 0.0; // average value in envelop
 
  // Matrices to hold data; buff holds the distribution
  //   data per envelope from which the max value is
  //   obtained. That value is then recorded in dist, where
  //   the integral of the envelope to that point is
  //   also recorded
 
  double * buff = new double[numPoints/numEnv + 1];
  double ** dist = new double*[numEnv];
  for(int ih=0;ih<numEnv;ih++)
  {
    dist[ih] = new double[3];
  }
 
  // Acc-Rej Sampling Method
  // -- Starting at the beginning of each envelope,
  // this loop evaluates (numPoints) amount of points
  // on the distribution and holds them in buff;
  // then takes the max and  places it in the first row
  // of h. The second row of h contains the interval of
  // the total envelope, up to the current envelope.
  // Thus, when properly normalized, the last value
  // in the second row of h should be 1.
  double totxsec = 0.0;
  for(int i=0;i<numEnv;i++)
    {
      double lbound = -1 + i*sr;
 
      for(int j=0;j<=numPoints / numEnv; j++)
        {
          buff[j] = this->XSec(p->Pdg(),target,scode,fate,ke,lbound+j*cstep);
          avg += buff[j];
        }
 
      totxsec+=avg;
      avg/= (double(numPoints)/double(numEnv));
      dist[i][0] = TMath::MaxElement(numPoints/numEnv+1,buff);
      dist[i][1] = avg;
      dist[i][2] = dist[i][1] + ((i==0)?0.0:dist[i-1][2]);
      avg=0.0;
    }
 
 
  delete [] buff;
 
  int iter=1;         // keep track of iterations
  int env=0;          // envelope index
  double rval = 0.0;  // random value
  double val = 0.0;   // angle value
 
  // Get a random point, see if its in the distribution, and if not
  // then try again.
 
  rval = rnd->RndFsi().Rndm()*dist[numEnv-1][2];
 
  env=0;
  // Check which envelope it's in, to
  // get proper height
  while(env<numEnv)
    {
      if(rval<=dist[env][2]) break;
      else env++;
    }
  if(env==numEnv) env=numEnv - 1;
 
while(iter)
    {
 
      // Obtain the correct x-coordinate from the random sample
      val = rnd->RndFsi().Rndm()*sr;
      val += sr*env-1;
      rval = rnd->RndFsi().Rndm()*dist[env][0];
 
      // Test to see if point is in distribution, if it is, stop and return
      if(rval < this->XSec(p->Pdg(),target,scode,fate,ke,val)) break;
 
      // Possibly an extremely long loop, don't want to
      // hold up the program
      if(iter==1000)
        {
          int NUM_POINTS=2000;
          int pvalues=0;
          double points[200]={0};
          for(int k=0;k<NUM_POINTS;k++)
          {
            points[int(k/10)]=this->XSec(p->Pdg(),target,scode,fate,ke,-1+(2.0/NUM_POINTS)*k);
            if(points[int(k/10)]>0) pvalues++;
          }
          if(pvalues<(.05*NUM_POINTS))
          {
            // if it reaches here, one more test...if momenta of particle is
            // extremely low, just give it an angle from a uniform distribution
            if(p->P4()->P()<.005) // 5 MeV
            {
              val = 2*rnd->RndFsi().Rndm()-1;
              break;
            }
            else
            {
              LOG("Intranuke", pWARN) << "Hung-up in IntBounce method - Exiting";
              LOG("Intranuke", pWARN) << (*p);
              LOG("Intranuke", pWARN) << "Target: " << target << ", Scode: " << scode << ", fate: " << INukeHadroFates::AsString(fate);
              for(int ie=0;ie<200;ie+=10) {
                LOG("Intranuke", pWARN)   << points[ie+0] << ", " << points[ie+1] << ", " << points[ie+2] << ", "
                   << points[ie+3] << ", " << points[ie+4] << ", " << points[ie+5] << ", " << points[ie+6] << ", "
                                           << points[ie+7] << ", " << points[ie+8] << ", " << points[ie+9];
              }
              for(int ih=0;ih<numEnv;ih++)
              {
                delete [] dist[ih];
              }
              delete [] dist;
 
              return -2.;
            }
          }
        }
      iter++;
    }
 
  for(int ih=0;ih<numEnv;ih++)
  {
    delete [] dist[ih];
  }
  delete [] dist;
 
  return val;
}

References genie::INukeHadroFates::AsString(), genie::GHepParticle::E(), genie::RandomGen::Instance(), LOG, genie::GHepParticle::Mass(), genie::GHepParticle::P4(), genie::GHepParticle::Pdg(), pWARN, genie::RandomGen::RndFsi(), and XSec().

Referenced by INukeNucleonCorr::AvgCorrection().

LoadCrossSections()

void INukeHadroData2018::LoadCrossSections ( void )

private

Definition at line 186 of file INukeHadroData2018.cxx.

{
// Loads hadronic x-section data
 
  //-- Get the top-level directory with input hadron cross-section data
  //   (search for $GINUKEHADRONDATA or use default location)
  string data_dir = (gSystem->Getenv("GINUKEHADRONDATA")) ?
             string(gSystem->Getenv("GINUKEHADRONDATA")) :
             string(gSystem->Getenv("GENIE")) + string("/data/evgen/intranuke");
 
  LOG("INukeData", pINFO)
      << "Loading INTRANUKE hadron data from: " << data_dir;
 
  //-- Build filenames
 
  string datafile_NN   = data_dir + "/tot_xsec/intranuke-xsections-NN2014.dat";
  string datafile_pipN = data_dir + "/tot_xsec/intranuke-xsections-pi+N.dat";
  string datafile_pi0N = data_dir + "/tot_xsec/intranuke-xsections-pi0N.dat";
  string datafile_NA   = data_dir + "/tot_xsec/intranuke-fractions-NA2016.dat";
  string datafile_KA   = data_dir + "/tot_xsec/intranuke-fractions-KA.dat";
  string datafile_gamN = data_dir + "/tot_xsec/intranuke-xsections-gamN.dat";
  string datafile_kN   = data_dir + "/tot_xsec/intranuke-xsections-kaonN2018.dat";
 
  //-- Make sure that all data files are available
 
  assert( ! gSystem->AccessPathName(datafile_NN.  c_str()) );
  assert( ! gSystem->AccessPathName(datafile_pipN.c_str()) );
  assert( ! gSystem->AccessPathName(datafile_pi0N.c_str()) );
  assert( ! gSystem->AccessPathName(datafile_NA.  c_str()) );
  assert( ! gSystem->AccessPathName(datafile_KA. c_str())  );
  assert( ! gSystem->AccessPathName(datafile_gamN.c_str())  );
  assert( ! gSystem->AccessPathName(datafile_kN.  c_str())  );
 
  LOG("INukeData", pINFO)  << "Found all necessary data files...";
 
  //-- Load data files
 
  TTree data_NN;
  TTree data_pipN;
  TTree data_pi0N;
  TTree data_NA;
  TTree data_KA;
  TTree data_gamN;
  TTree data_kN;
 
  data_NN.ReadFile(datafile_NN.c_str(),"ke/D:pp_tot/D:pp_elas/D:pp_reac/D:pn_tot/D:pn_elas/D:pn_reac/D:nn_tot/D:nn_elas/D:nn_reac/D:pp_cmp/D:pn_cmp/D:nn_cmp/D");
  data_pipN.ReadFile(datafile_pipN.c_str(),
     "ke/D:pipn_tot/D:pipn_cex/D:pipn_elas/D:pipn_reac/D:pipp_tot/D:pipp_cex/D:pipp_elas/D:pipp_reac/D:pipd_abs");
  data_pi0N.ReadFile(datafile_pi0N.c_str(),
     "ke/D:pi0n_tot/D:pi0n_cex/D:pi0n_elas/D:pi0n_reac/D:pi0p_tot/D:pi0p_cex/D:pi0p_elas/D:pi0p_reac/D:pi0d_abs");
  //data_NA.ReadFile(datafile_NA.c_str(),
  //"ke/D:pA_tot/D:pA_elas/D:pA_inel/D:pA_cex/D:pA_abs/D:pA_pipro/D");  // add support for cmp here (?)
  data_NA.ReadFile(datafile_NA.c_str(),
                   "ke/D:pA_tot/D:pA_inel/D:pA_cex/D:pA_abs/D:pA_pipro/D:pA_cmp/D");  // add support for cmp here (?)
  data_gamN.ReadFile(datafile_gamN.c_str(),
    "ke/D:pi0p_tot/D:pipn_tot/D:pimp_tot/D:pi0n_tot/D:gamp_fs/D:gamn_fs/D:gamN_tot/D");
  data_kN.ReadFile(datafile_kN.c_str(),
                   "ke/D:kpp_elas/D:kpn_elas/D:kpn_cex/D:kp_abs/D:kpN_tot/D");
  data_KA.ReadFile(datafile_KA.c_str(),
     "ke/D:KA_tot/D:KA_elas/D:KA_inel/D:KA_abs/D");
 
  LOG("INukeData", pDEBUG)  << "Number of data rows in NN : "   << data_NN.GetEntries();
  LOG("INukeData", pDEBUG)  << "Number of data rows in pipN : " << data_pipN.GetEntries();
  LOG("INukeData", pDEBUG)  << "Number of data rows in pi0N : " << data_pi0N.GetEntries();
  LOG("INukeData", pDEBUG)  << "Number of data rows in NA  : "  << data_NA.GetEntries();
  LOG("INukeData", pDEBUG)  << "Number of data rows in KA : "   << data_KA.GetEntries();
  LOG("INukeData", pDEBUG)  << "Number of data rows in gamN : " << data_gamN.GetEntries();
  LOG("INukeData", pDEBUG)  << "Number of data rows in kN  : "  << data_kN.GetEntries();
 
  LOG("INukeData", pINFO)  << "Done loading all x-section files...";
 
  //-- Build x-section splines
 
  // p/n+p/n hA x-section splines
  fXSecPp_Tot      = new Spline(&data_NN, "ke:pp_tot");
  fXSecPp_Elas     = new Spline(&data_NN, "ke:pp_elas");
  fXSecPp_Reac     = new Spline(&data_NN, "ke:pp_reac");
  fXSecPn_Tot      = new Spline(&data_NN, "ke:pn_tot");
  fXSecPn_Elas     = new Spline(&data_NN, "ke:pn_elas");
  fXSecPn_Reac     = new Spline(&data_NN, "ke:pn_reac");
  fXSecNn_Tot      = new Spline(&data_NN, "ke:nn_tot");
  fXSecNn_Elas     = new Spline(&data_NN, "ke:nn_elas");
  fXSecNn_Reac     = new Spline(&data_NN, "ke:nn_reac");
  fXSecPp_Cmp      = new Spline(&data_NN, "ke:pp_cmp"); //for compound nucleus fate
  fXSecPn_Cmp      = new Spline(&data_NN, "ke:pn_cmp");
  fXSecNn_Cmp      = new Spline(&data_NN, "ke:nn_cmp");
 
  // pi+n/p hA x-section splines
  fXSecPipn_Tot     = new Spline(&data_pipN, "ke:pipn_tot");
  fXSecPipn_CEx     = new Spline(&data_pipN, "ke:pipn_cex");
  fXSecPipn_Elas    = new Spline(&data_pipN, "ke:pipn_elas");
  fXSecPipn_Reac    = new Spline(&data_pipN, "ke:pipn_reac");
  fXSecPipp_Tot     = new Spline(&data_pipN, "ke:pipp_tot");
  fXSecPipp_CEx     = new Spline(&data_pipN, "ke:pipp_cex");
  fXSecPipp_Elas    = new Spline(&data_pipN, "ke:pipp_elas");
  fXSecPipp_Reac    = new Spline(&data_pipN, "ke:pipp_reac");
  fXSecPipd_Abs     = new Spline(&data_pipN, "ke:pipd_abs");
 
  // pi0n/p hA x-section splines
  fXSecPi0n_Tot     = new Spline(&data_pi0N, "ke:pi0n_tot");
  fXSecPi0n_CEx     = new Spline(&data_pi0N, "ke:pi0n_cex");
  fXSecPi0n_Elas    = new Spline(&data_pi0N, "ke:pi0n_elas");
  fXSecPi0n_Reac    = new Spline(&data_pi0N, "ke:pi0n_reac");
  fXSecPi0p_Tot     = new Spline(&data_pi0N, "ke:pi0p_tot");
  fXSecPi0p_CEx     = new Spline(&data_pi0N, "ke:pi0p_cex");
  fXSecPi0p_Elas    = new Spline(&data_pi0N, "ke:pi0p_elas");
  fXSecPi0p_Reac    = new Spline(&data_pi0N, "ke:pi0p_reac");
  fXSecPi0d_Abs     = new Spline(&data_pi0N, "ke:pi0d_abs");
 
   // K+N x-section splines
  fXSecKpn_Elas   = new Spline(&data_kN,  "ke:kpn_elas");
  fXSecKpp_Elas   = new Spline(&data_kN,  "ke:kpp_elas");
  fXSecKpn_CEx    = new Spline(&data_kN,  "ke:kpn_cex");
  fXSecKpN_Abs    = 0; // new Spline(&data_kN,  "ke:kp_abs");  why not used?
  fXSecKpN_Tot    = new Spline(&data_kN,  "ke:kpN_tot");
 
  // gamma x-section splines
  fXSecGamp_fs     = new Spline(&data_gamN, "ke:gamp_fs");
  fXSecGamn_fs     = new Spline(&data_gamN, "ke:gamn_fs");
  fXSecGamN_Tot    = new Spline(&data_gamN, "ke:gamN_tot");
 
  // N+A x-section fraction splines
  fFracPA_Tot      = new Spline(&data_NA, "ke:pA_tot");
  //  fFracPA_Elas     = new Spline(&data_NA, "ke:pA_elas");
  fFracPA_Inel     = new Spline(&data_NA, "ke:pA_inel");
  fFracPA_CEx      = new Spline(&data_NA, "ke:pA_cex");
  fFracPA_Abs      = new Spline(&data_NA, "ke:pA_abs");
  fFracPA_PiPro    = new Spline(&data_NA, "ke:pA_pipro");
  fFracNA_Tot      = new Spline(&data_NA, "ke:pA_tot");  // assuming nA same as pA
  //  fFracNA_Elas     = new Spline(&data_NA, "ke:pA_elas");
  fFracNA_Inel     = new Spline(&data_NA, "ke:pA_inel");
  fFracNA_CEx      = new Spline(&data_NA, "ke:pA_cex");
  fFracNA_Abs      = new Spline(&data_NA, "ke:pA_abs");
  fFracNA_PiPro    = new Spline(&data_NA, "ke:pA_pipro");
 
  fFracPA_Cmp      = new Spline(&data_NA, "ke:pA_cmp");  //cmp - add support later
  fFracNA_Cmp      = new Spline(&data_NA, "ke:pA_cmp");
 
  // K+A x-section fraction splines
  fFracKA_Tot      = new Spline(&data_KA, "ke:KA_tot");
  fFracKA_Elas     = new Spline(&data_KA, "ke:KA_elas");
  fFracKA_CEx      = 0; // new Spline(&data_KA, "ke:KA_cex"); //Added, but needs to be computed
  fFracKA_Inel     = new Spline(&data_KA, "ke:KA_inel");
  fFracKA_Abs      = new Spline(&data_KA, "ke:KA_abs");
  //
  // hN stuff
  //
 
 
  //    kIHNFtElas
  //      pp, nn --> read from pp/pp%.txt
  //      pn, np --> read from pp/pn%.txt
  //      pi+ N  --> read from pip/pip%.txt
  //      pi0 N  --> read from pip/pip%.txt
  //      pi- N  --> read from pim/pim%.txt
  //      K+  N  --> read from kpn/kpn%.txt
  //      K+  P  --> read from kpp/kpp%.txt
  //    kIHNFtCEx
  //      pi+, pi0, pi- --> read from pie/pie%.txt (using pip+n->pi0+p data)
  //    kIHNFtAbs
  //      pi+, pi0, pi- --> read from pid2p/pid2p%.txt (using pip+D->2p data)
  //    kIHNFtInelas
  //      gamma p -> p pi0 --> read from gampi0p/%-pi0p.txt
  //      gamma p -> n pi+ --> read from gampi+n/%-pi+n.txt
  //      gamma n -> n pi0 --> read from gampi0n/%-pi0n.txt
  //      gamma n -> p pi- --> read from gampi-p/%-pi-p.txt
 
 
  // kIHNFtElas, pp&nn :
  {
    const int hN_ppelas_nfiles = 20;
    const int hN_ppelas_points_per_file = 21;
    const int hN_ppelas_npoints = hN_ppelas_points_per_file * hN_ppelas_nfiles;
 
    double hN_ppelas_energies[hN_ppelas_nfiles] = {
        50, 100, 150, 200, 250, 300, 350, 400, 450, 500,
       550, 600, 650, 700, 750, 800, 850, 900, 950, 1000
    };
 
    double hN_ppelas_costh [hN_ppelas_points_per_file];
    double hN_ppelas_xsec  [hN_ppelas_npoints];
 
    int ipoint=0;
 
    for(int ifile = 0; ifile < hN_ppelas_nfiles; ifile++) {
     // build filename
     ostringstream hN_datafile;
     double ke = hN_ppelas_energies[ifile];
     hN_datafile << data_dir << "/diff_ang/pp/pp" << ke << ".txt";
     // read data
     ReadhNFile(
                hN_datafile.str(), ke, hN_ppelas_points_per_file,
                ipoint, hN_ppelas_costh, hN_ppelas_xsec,2);
    }//loop over files
 
    fhN2dXSecPP_Elas = new BLI2DNonUnifGrid(hN_ppelas_nfiles,hN_ppelas_points_per_file,
                           hN_ppelas_energies,hN_ppelas_costh,hN_ppelas_xsec);
  }
 
  // kIHNFtElas, pn&np :
  {
    const int hN_npelas_nfiles = 20;
    const int hN_npelas_points_per_file = 21;
    const int hN_npelas_npoints = hN_npelas_points_per_file * hN_npelas_nfiles;
 
    double hN_npelas_energies[hN_npelas_nfiles] = {
        50, 100, 150, 200, 250, 300, 350, 400, 450, 500,
       550, 600, 650, 700, 750, 800, 850, 900, 950, 1000
    };
 
    double hN_npelas_costh [hN_npelas_points_per_file];
    double hN_npelas_xsec  [hN_npelas_npoints];
 
    int ipoint=0;
 
    for(int ifile = 0; ifile < hN_npelas_nfiles; ifile++) {
     // build filename
     ostringstream hN_datafile;
     double ke = hN_npelas_energies[ifile];
     hN_datafile << data_dir << "/diff_ang/pn/pn" << ke << ".txt";
     // read data
     ReadhNFile(
       hN_datafile.str(), ke, hN_npelas_points_per_file,
       ipoint, hN_npelas_costh, hN_npelas_xsec,2);
    }//loop over files
 
    fhN2dXSecNP_Elas = new BLI2DNonUnifGrid(hN_npelas_nfiles,hN_npelas_points_per_file,
                           hN_npelas_energies,hN_npelas_costh,hN_npelas_xsec);
  }
 
  // kIHNFtElas, pipN :
  {
    const int hN_pipNelas_nfiles = 60;
    const int hN_pipNelas_points_per_file = 21;
    const int hN_pipNelas_npoints = hN_pipNelas_points_per_file * hN_pipNelas_nfiles;
 
    double hN_pipNelas_energies[hN_pipNelas_nfiles] = {
      10,  20,  30,  40,  50,  60,  70,  80,  90,
     100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
     200, 210, 220, 230, 240, 250, 260, 270, 280, 290,
     300, 340, 380, 420, 460, 500, 540, 580, 620, 660,
     700, 740, 780, 820, 860, 900, 940, 980,
     1020, 1060, 1100, 1140, 1180, 1220, 1260,
     1300, 1340, 1380, 1420, 1460, 1500
    };
 
    double hN_pipNelas_costh [hN_pipNelas_points_per_file];
    double hN_pipNelas_xsec  [hN_pipNelas_npoints];
 
    int ipoint=0;
 
    for(int ifile = 0; ifile < hN_pipNelas_nfiles; ifile++) {
     // build filename
     ostringstream hN_datafile;
     double ke = hN_pipNelas_energies[ifile];
     hN_datafile << data_dir << "/diff_ang/pip/pip" << ke << ".txt";
     // read data
     ReadhNFile(
       hN_datafile.str(), ke, hN_pipNelas_points_per_file,
       ipoint, hN_pipNelas_costh, hN_pipNelas_xsec,2);
    }//loop over files
 
    fhN2dXSecPipN_Elas = new BLI2DNonUnifGrid(hN_pipNelas_nfiles,hN_pipNelas_points_per_file,
                           hN_pipNelas_energies,hN_pipNelas_costh,hN_pipNelas_xsec);
  }
 
  // kIHNFtElas, pi0N :
  {
    const int hN_pi0Nelas_nfiles = 60;
    const int hN_pi0Nelas_points_per_file = 21;
    const int hN_pi0Nelas_npoints = hN_pi0Nelas_points_per_file * hN_pi0Nelas_nfiles;
 
    double hN_pi0Nelas_energies[hN_pi0Nelas_nfiles] = {
      10,  20,  30,  40,  50,  60,  70,  80,  90,
     100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
     200, 210, 220, 230, 240, 250, 260, 270, 280, 290,
     300, 340, 380, 420, 460, 500, 540, 580, 620, 660,
     700, 740, 780, 820, 860, 900, 940, 980,
     1020, 1060, 1100, 1140, 1180, 1220, 1260,
     1300, 1340, 1380, 1420, 1460, 1500
    };
 
    double hN_pi0Nelas_costh [hN_pi0Nelas_points_per_file];
    double hN_pi0Nelas_xsec  [hN_pi0Nelas_npoints];
 
    int ipoint=0;
 
    for(int ifile = 0; ifile < hN_pi0Nelas_nfiles; ifile++) {
     // build filename
     ostringstream hN_datafile;
     double ke = hN_pi0Nelas_energies[ifile];
     hN_datafile << data_dir << "/diff_ang/pip/pip" << ke << ".txt";
     // read data
     ReadhNFile(
       hN_datafile.str(), ke, hN_pi0Nelas_points_per_file,
       ipoint, hN_pi0Nelas_costh, hN_pi0Nelas_xsec,2);
    }//loop over files
 
    fhN2dXSecPi0N_Elas = new BLI2DNonUnifGrid(hN_pi0Nelas_nfiles,hN_pi0Nelas_points_per_file,
                           hN_pi0Nelas_energies,hN_pi0Nelas_costh,hN_pi0Nelas_xsec);
  }
 
  // kIHNFtElas, pimN :
  {
    const int hN_pimNelas_nfiles = 60;
    const int hN_pimNelas_points_per_file = 21;
    const int hN_pimNelas_npoints = hN_pimNelas_points_per_file * hN_pimNelas_nfiles;
 
    double hN_pimNelas_energies[hN_pimNelas_nfiles] = {
      10,  20,  30,  40,  50,  60,  70,  80,  90,
     100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
     200, 210, 220, 230, 240, 250, 260, 270, 280, 290,
     300, 340, 380, 420, 460, 500, 540, 580, 620, 660,
     700, 740, 780, 820, 860, 900, 940, 980,
     1020, 1060, 1100, 1140, 1180, 1220, 1260,
     1300, 1340, 1380, 1420, 1460, 1500
    };
 
    double hN_pimNelas_costh [hN_pimNelas_points_per_file];
    double hN_pimNelas_xsec  [hN_pimNelas_npoints];
 
    int ipoint=0;
 
    for(int ifile = 0; ifile < hN_pimNelas_nfiles; ifile++) {
     // build filename
     ostringstream hN_datafile;
     double ke = hN_pimNelas_energies[ifile];
     hN_datafile << data_dir << "/diff_ang/pim/pim" << ke << ".txt";
     // read data
     ReadhNFile(
       hN_datafile.str(), ke, hN_pimNelas_points_per_file,
       ipoint, hN_pimNelas_costh, hN_pimNelas_xsec,2);
    }//loop over files
 
    fhN2dXSecPimN_Elas = new BLI2DNonUnifGrid(hN_pimNelas_nfiles,hN_pimNelas_points_per_file,
                           hN_pimNelas_energies,hN_pimNelas_costh,hN_pimNelas_xsec);
  }
 
  // kIHNFtElas, kpn :
  {
    const int hN_kpNelas_nfiles = 18;
    const int hN_kpNelas_points_per_file = 37;
    const int hN_kpNelas_npoints = hN_kpNelas_points_per_file * hN_kpNelas_nfiles;
 
    double hN_kpNelas_energies[hN_kpNelas_nfiles] = {
     100, 200, 300, 400, 500, 600, 700, 800, 900, 1000,
     1100, 1200, 1300, 1400, 1500, 1600, 1700, 1800
    };
 
    double hN_kpNelas_costh [hN_kpNelas_points_per_file];
    double hN_kpNelas_xsec  [hN_kpNelas_npoints];
 
    int ipoint=0;
 
    for(int ifile = 0; ifile < hN_kpNelas_nfiles; ifile++) {
     // build filename
     ostringstream hN_datafile;
     double ke = hN_kpNelas_energies[ifile];
     hN_datafile << data_dir << "/diff_ang/kpn/kpn" << ke << ".txt";
     // read data
     ReadhNFile(
       hN_datafile.str(), ke, hN_kpNelas_points_per_file,
       ipoint, hN_kpNelas_costh, hN_kpNelas_xsec,2);
    }//loop over files
 
    fhN2dXSecKpN_Elas = new BLI2DNonUnifGrid(hN_kpNelas_nfiles,hN_kpNelas_points_per_file,
                           hN_kpNelas_energies,hN_kpNelas_costh,hN_kpNelas_xsec);
  }
  // kIHNFtCEx, kpn :
  {
    const int hN_kpNcex_nfiles = 18;
    const int hN_kpNcex_points_per_file = 37;
    const int hN_kpNcex_npoints = hN_kpNcex_points_per_file * hN_kpNcex_nfiles;
 
    double hN_kpNcex_energies[hN_kpNcex_nfiles] = {
     100, 200, 300, 400, 500, 600, 700, 800, 900, 1000,
     1100, 1200, 1300, 1400, 1500, 1600, 1700, 1800
    };
 
    double hN_kpNcex_costh [hN_kpNcex_points_per_file];
    double hN_kpNcex_xsec  [hN_kpNcex_npoints];
 
    int ipoint=0;
 
    for(int ifile = 0; ifile < hN_kpNcex_nfiles; ifile++) {
     // build filename
     ostringstream hN_datafile;
     double ke = hN_kpNcex_energies[ifile];
     hN_datafile << data_dir << "/diff_ang/kpncex/kpcex" << ke << ".txt";
     // read data
     ReadhNFile(
       hN_datafile.str(), ke, hN_kpNcex_points_per_file,
       ipoint, hN_kpNcex_costh, hN_kpNcex_xsec,2);
    }//loop over files
 
    /*double hN_kpNcex_costh_cond [hN_kpNcex_points_per_file];
    for (int ient = 0; ient < hN_kpNcex_points_per_file; ient++) {
      hN_kpNcex_costh_cond[ient] = hN_kpNcex_costh[ient];
      }*/
 
    fhN2dXSecKpN_CEx = new BLI2DNonUnifGrid(hN_kpNcex_nfiles,hN_kpNcex_points_per_file,
                           hN_kpNcex_energies,hN_kpNcex_costh,hN_kpNcex_xsec);
  }
//----------------------------------------------------------------------------------------
 
 
  // kIHNFtElas, kpp :
  {
    const int hN_kpPelas_nfiles = 18;
    const int hN_kpPelas_points_per_file = 37;
    const int hN_kpPelas_npoints = hN_kpPelas_points_per_file * hN_kpPelas_nfiles;
 
    double hN_kpPelas_energies[hN_kpPelas_nfiles] = {
     100, 200, 300, 400, 500, 600, 700, 800, 900, 1000,
     1100, 1200, 1300, 1400, 1500, 1600, 1700, 1800
    };
 
    double hN_kpPelas_costh [hN_kpPelas_points_per_file];
    double hN_kpPelas_xsec  [hN_kpPelas_npoints];
 
    int ipoint=0;
 
    for(int ifile = 0; ifile < hN_kpPelas_nfiles; ifile++) {
     // build filename
     ostringstream hN_datafile;
     double ke = hN_kpPelas_energies[ifile];
     hN_datafile << data_dir << "/diff_ang/kpp/kpp" << ke << ".txt";
     // read data
     ReadhNFile(
       hN_datafile.str(), ke, hN_kpPelas_points_per_file,
       ipoint, hN_kpPelas_costh, hN_kpPelas_xsec,2);
    }//loop over files
 
    fhN2dXSecKpP_Elas = new BLI2DNonUnifGrid(hN_kpPelas_nfiles,hN_kpPelas_points_per_file,
                           hN_kpPelas_energies,hN_kpPelas_costh,hN_kpPelas_xsec);
        }
 
  // kIHNFtCEx, (pi+, pi0, pi-) N
  {
    const int hN_piNcex_nfiles = 60;
    const int hN_piNcex_points_per_file = 21;
    const int hN_piNcex_npoints = hN_piNcex_points_per_file * hN_piNcex_nfiles;
 
    double hN_piNcex_energies[hN_piNcex_nfiles] = {
      10,  20,  30,  40,  50,  60,  70,  80,  90,
     100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
     200, 210, 220, 230, 240, 250, 260, 270, 280, 290,
     300, 340, 380, 420, 460, 500, 540, 580, 620, 660,
     700, 740, 780, 820, 860, 900, 940, 980,
     1020, 1060, 1100, 1140, 1180, 1220, 1260,
     1300, 1340, 1380, 1420, 1460, 1500
    };
 
    double hN_piNcex_costh [hN_piNcex_points_per_file];
    double hN_piNcex_xsec  [hN_piNcex_npoints];
 
    int ipoint=0;
 
    for(int ifile = 0; ifile < hN_piNcex_nfiles; ifile++) {
     // build filename
     ostringstream hN_datafile;
     double ke = hN_piNcex_energies[ifile];
     hN_datafile << data_dir << "/diff_ang/pie/pie" << ke << ".txt";
     // read data
     ReadhNFile(
       hN_datafile.str(), ke, hN_piNcex_points_per_file,
       ipoint, hN_piNcex_costh, hN_piNcex_xsec,2);
    }//loop over files
 
    fhN2dXSecPiN_CEx = new BLI2DNonUnifGrid(hN_piNcex_nfiles,hN_piNcex_points_per_file,
                           hN_piNcex_energies,hN_piNcex_costh,hN_piNcex_xsec);
  }
 
  // kIHNFtAbs, (pi+, pi0, pi-) N
  {
    const int hN_piNabs_nfiles = 19;
    const int hN_piNabs_points_per_file = 21;
    const int hN_piNabs_npoints = hN_piNabs_points_per_file * hN_piNabs_nfiles;
 
    double hN_piNabs_energies[hN_piNabs_nfiles] = {
      50,  75, 100, 125, 150, 175, 200, 225, 250, 275,
     300, 325, 350, 375, 400, 425, 450, 475, 500
    };
 
    double hN_piNabs_costh [hN_piNabs_points_per_file];
    double hN_piNabs_xsec  [hN_piNabs_npoints];
 
    int ipoint=0;
 
    for(int ifile = 0; ifile < hN_piNabs_nfiles; ifile++) {
     // build filename
     ostringstream hN_datafile;
     double ke = hN_piNabs_energies[ifile];
     hN_datafile << data_dir << "/diff_ang/pid2p/pid2p" << ke << ".txt";
     // read data
     ReadhNFile(
       hN_datafile.str(), ke, hN_piNabs_points_per_file,
       ipoint, hN_piNabs_costh, hN_piNabs_xsec,2);
    }//loop over files
 
    fhN2dXSecPiN_Abs = new BLI2DNonUnifGrid(hN_piNabs_nfiles,hN_piNabs_points_per_file,
                           hN_piNabs_energies,hN_piNabs_costh,hN_piNabs_xsec);
  }
 
  // kIHNFtInelas, gamma p -> p pi0
    {
    const int hN_gampi0pInelas_nfiles = 29;
    const int hN_gampi0pInelas_points_per_file = 37;
    const int hN_gampi0pInelas_npoints = hN_gampi0pInelas_points_per_file * hN_gampi0pInelas_nfiles;
 
    double hN_gampi0pInelas_energies[hN_gampi0pInelas_nfiles] = {
      160,  180,  200,  220,  240,  260,  280,  300,  320,  340,
      360,  380,  400,  450,  500,  550,  600,  650,  700,  750,
      800,  850,  900,  950,  1000, 1050, 1100, 1150, 1200
    };
 
    double hN_gampi0pInelas_costh [hN_gampi0pInelas_points_per_file];
    double hN_gampi0pInelas_xsec  [hN_gampi0pInelas_npoints];
 
    int ipoint=0;
 
    for(int ifile = 0; ifile < hN_gampi0pInelas_nfiles; ifile++) {
     // build filename
     ostringstream hN_datafile;
     double ke = hN_gampi0pInelas_energies[ifile];
     hN_datafile << data_dir << "/diff_ang/gampi0p/" << ke << "-pi0p.txt";
     // read data
     ReadhNFile(
       hN_datafile.str(), ke, hN_gampi0pInelas_points_per_file,
       ipoint, hN_gampi0pInelas_costh, hN_gampi0pInelas_xsec,3);
    }//loop over files
 
    fhN2dXSecGamPi0P_Inelas = new BLI2DNonUnifGrid(hN_gampi0pInelas_nfiles,hN_gampi0pInelas_points_per_file,
                           hN_gampi0pInelas_energies,hN_gampi0pInelas_costh,hN_gampi0pInelas_xsec);
  }
 
  // kIHNFtInelas, gamma n -> n pi0
  {
    const int hN_gampi0nInelas_nfiles = 29;
    const int hN_gampi0nInelas_points_per_file = 37;
    const int hN_gampi0nInelas_npoints = hN_gampi0nInelas_points_per_file * hN_gampi0nInelas_nfiles;
 
    double hN_gampi0nInelas_energies[hN_gampi0nInelas_nfiles] = {
      160,  180,  200,  220,  240,  260,  280,  300,  320,  340,
      360,  380,  400,  450,  500,  550,  600,  650,  700,  750,
      800,  850,  900,  950,  1000, 1050, 1100, 1150, 1200
    };
 
    double hN_gampi0nInelas_costh [hN_gampi0nInelas_points_per_file];
    double hN_gampi0nInelas_xsec  [hN_gampi0nInelas_npoints];
    int ipoint=0;
 
    for(int ifile = 0; ifile < hN_gampi0nInelas_nfiles; ifile++) {
     // build filename
     ostringstream hN_datafile;
     double ke = hN_gampi0nInelas_energies[ifile];
     hN_datafile << data_dir << "/diff_ang/gampi0n/" << ke << "-pi0n.txt";
     // read data
     ReadhNFile(
       hN_datafile.str(), ke, hN_gampi0nInelas_points_per_file,
       ipoint, hN_gampi0nInelas_costh, hN_gampi0nInelas_xsec,3);
    }//loop over files
 
    fhN2dXSecGamPi0N_Inelas = new BLI2DNonUnifGrid(hN_gampi0nInelas_nfiles,hN_gampi0nInelas_points_per_file,
                           hN_gampi0nInelas_energies,hN_gampi0nInelas_costh,hN_gampi0nInelas_xsec);
  }
 
  // kIHNFtInelas, gamma p -> n pi+
  {
    const int hN_gampipnInelas_nfiles = 29;
    const int hN_gampipnInelas_points_per_file = 37;
    const int hN_gampipnInelas_npoints = hN_gampipnInelas_points_per_file * hN_gampipnInelas_nfiles;
 
    double hN_gampipnInelas_energies[hN_gampipnInelas_nfiles] = {
      160,  180,  200,  220,  240,  260,  280,  300,  320,  340,
      360,  380,  400,  450,  500,  550,  600,  650,  700,  750,
      800,  850,  900,  950,  1000, 1050, 1100, 1150, 1200
    };
 
    double hN_gampipnInelas_costh [hN_gampipnInelas_points_per_file];
    double hN_gampipnInelas_xsec  [hN_gampipnInelas_npoints];
 
    int ipoint=0;
 
    for(int ifile = 0; ifile < hN_gampipnInelas_nfiles; ifile++) {
     // build filename
     ostringstream hN_datafile;
     double ke = hN_gampipnInelas_energies[ifile];
     hN_datafile << data_dir << "/diff_ang/gampi+n/" << ke << "-pi+n.txt";
     // read data
     ReadhNFile(
       hN_datafile.str(), ke, hN_gampipnInelas_points_per_file,
       ipoint, hN_gampipnInelas_costh, hN_gampipnInelas_xsec,3);
    }//loop over files
 
    fhN2dXSecGamPipN_Inelas = new BLI2DNonUnifGrid(hN_gampipnInelas_nfiles,hN_gampipnInelas_points_per_file,
                           hN_gampipnInelas_energies,hN_gampipnInelas_costh,hN_gampipnInelas_xsec);
  }
 
  // kIHNFtInelas, gamma n -> p pi-
  {
    const int hN_gampimpInelas_nfiles = 29;
    const int hN_gampimpInelas_points_per_file = 37;
    const int hN_gampimpInelas_npoints = hN_gampimpInelas_points_per_file * hN_gampimpInelas_nfiles;
 
    double hN_gampimpInelas_energies[hN_gampimpInelas_nfiles] = {
      160,  180,  200,  220,  240,  260,  280,  300,  320,  340,
      360,  380,  400,  450,  500,  550,  600,  650,  700,  750,
      800,  850,  900,  950,  1000, 1050, 1100, 1150, 1200
    };
 
    double hN_gampimpInelas_costh [hN_gampimpInelas_points_per_file];
    double hN_gampimpInelas_xsec  [hN_gampimpInelas_npoints];
 
    int ipoint=0;
 
    for(int ifile = 0; ifile < hN_gampimpInelas_nfiles; ifile++) {
     // build filename
     ostringstream hN_datafile;
     double ke = hN_gampimpInelas_energies[ifile];
     hN_datafile << data_dir << "/diff_ang/gampi-p/" << ke << "-pi-p.txt";
     // read data
     ReadhNFile(
       hN_datafile.str(), ke, hN_gampimpInelas_points_per_file,
       ipoint, hN_gampimpInelas_costh, hN_gampimpInelas_xsec,3);
    }//loop over files
 
    fhN2dXSecGamPimP_Inelas = new BLI2DNonUnifGrid(hN_gampimpInelas_nfiles,hN_gampimpInelas_points_per_file,
                           hN_gampimpInelas_energies,hN_gampimpInelas_costh,hN_gampimpInelas_xsec);
  }
 
 
 
 
  TFile TGraphs_file;
  bool saveTGraphsToFile = false;  //true;
 
  if (saveTGraphsToFile) {
    string filename = "TGraphs.root";
    LOG("INukeHadroData2018", pNOTICE) << "Saving INTRANUKE hadron x-section data to ROOT file: " << filename;
    TGraphs_file.Open(filename.c_str(), "RECREATE");
  }
 
  /*
  // kIHNFtTot,   pip + A                                            PipA_Tot
  {
    const int pipATot_nfiles = 22;
    const int pipATot_nuclei[pipATot_nfiles] = {1, 2, 3, 4, 6, 7, 9, 12, 16, 27, 28, 32, 40, 48, 56, 58, 63, 93, 120, 165, 181, 209};
    const int pipATot_npoints = 203;
 
    TPipA_Tot = new TGraph2D(pipATot_npoints);
    TPipA_Tot->SetNameTitle("TPipA_Tot","TPipA_Tot");
    TPipA_Tot->SetDirectory(0);
 
    int ipoint=0;
    double x, y;
 
    for(int ifile=0; ifile < pipATot_nfiles; ifile++) {
      ostringstream ADep_datafile;
      int nucleus = pipATot_nuclei[ifile];
      ADep_datafile << data_dir << "/tot_xsec/pipA_tot/pip" << nucleus << "_tot.txt";
      TGraph * buff = new TGraph(ADep_datafile.str().c_str());
      buff->SetNameTitle("buff","buff");
      for(int i=0; i < buff->GetN(); i++) {
        buff -> GetPoint(i,x,y);
        TPipA_Tot -> SetPoint(ipoint,(double)nucleus,x,y);
        ipoint++;
      }
      delete buff;
    }
 
    if (saveTGraphsToFile) {
      TPipA_Tot -> Write("TPipA_Tot"); // TPipA_Tot will be _key_ name
    }
  }
  */
 
  // kIHNFtAbs, pip + A                                                            PipA_Abs_frac
  {
    const int pipAAbs_f_nfiles = 18;
    const int pipAAbs_f_nuclei[pipAAbs_f_nfiles] = {1, 2, 3, 4, 7, 9, 12, 16, 27, 48, 56, 58, 63, 93, 120, 165, 181, 209};
    const int pipAAbs_f_npoints = 111;
 
    TfracPipA_Abs = new TGraph2D(pipAAbs_f_npoints);
    TfracPipA_Abs->SetNameTitle("TfracPipA_Abs","TfracPipA_Abs");
    TfracPipA_Abs->SetDirectory(0);
 
    int ipoint=0;
    double x, y;
    for(int ifile=0; ifile < pipAAbs_f_nfiles; ifile++) {
      ostringstream ADep_datafile;
      int nucleus = pipAAbs_f_nuclei[ifile];
      ADep_datafile << data_dir << "/tot_xsec/pipA_abs_frac/pip" << nucleus << "_abs_frac.txt";
      TGraph * buff = new TGraph(ADep_datafile.str().c_str());
      buff->SetNameTitle("buff","buff");
      for(int i=0; i < buff->GetN(); i++) {
        buff -> GetPoint(i,x,y);
        TfracPipA_Abs -> SetPoint(ipoint,(double)nucleus,x,y);
        ipoint++;
      }
      delete buff;
    }
    if (saveTGraphsToFile) {
      TfracPipA_Abs -> Write("TfracPipA_Abs");
    }
 
  }
 
 
  // kIHNFtCEx, pip + A      PipA_CEx_frac
  {
    const int pipACEx_f_nfiles = 18;
    const int pipACEx_f_nuclei[pipACEx_f_nfiles] = {1, 2, 3, 4, 7, 9, 12, 16, 27, 48, 56, 58, 63, 93, 120, 165, 181, 209};
    const int pipACEx_f_npoints = 129;
 
    TfracPipA_CEx = new TGraph2D(pipACEx_f_npoints);
    TfracPipA_CEx->SetNameTitle("TfracPipA_CEx","TfracPipA_CEx");
    TfracPipA_CEx->SetDirectory(0);
 
    int ipoint=0;
    double x, y;
 
    for(int ifile=0; ifile < pipACEx_f_nfiles; ifile++) {
      ostringstream ADep_datafile;
      int nucleus = pipACEx_f_nuclei[ifile];
      ADep_datafile << data_dir << "/tot_xsec/pipA_cex_frac/pip" << nucleus << "_cex_frac.txt";
      TGraph * buff = new TGraph(ADep_datafile.str().c_str());
      buff->SetNameTitle("buff","buff");
      for(int i=0; i < buff->GetN(); i++) {
        buff -> GetPoint(i,x,y);
        TfracPipA_CEx -> SetPoint(ipoint,(double)nucleus,x,y);
        ipoint++;
      }
      delete buff;
    }
 
    if (saveTGraphsToFile) {
      TfracPipA_CEx -> Write("TfracPipA_CEx");
    }
 
  }
 
 
 
  // kIHNFtCEx, pip + A                                                            PipA_CEx (just for developmental purposes)
  {
    TGraph2D * TPipA_CEx;
 
    const int pipACEx_nfiles = 18;
    const int pipACEx_nuclei[pipACEx_nfiles] = {1, 2, 3, 4, 7, 9, 12, 16, 27, 48, 56, 58, 63, 93, 120, 165, 181, 209};
    const int pipACEx_npoints = 129;
 
    TPipA_CEx = new TGraph2D(pipACEx_npoints);
    TPipA_CEx->SetNameTitle("TPipA_CEx","TPipA_CEx");
    TPipA_CEx->SetDirectory(0);
 
    int ipoint=0;
    double x, y;
 
    for(int ifile=0; ifile < pipACEx_nfiles; ifile++) {
      ostringstream ADep_datafile;
      int nucleus = pipACEx_nuclei[ifile];
      ADep_datafile << data_dir << "/tot_xsec/pipA_cex/pip" << nucleus << "_cex.txt";
      TGraph * buff = new TGraph(ADep_datafile.str().c_str());
      buff->SetNameTitle("buff","buff");
      for(int i=0; i < buff->GetN(); i++) {
        buff -> GetPoint(i,x,y);
        TPipA_CEx -> SetPoint(ipoint,(double)nucleus,x,y);
        ipoint++;
      }
      delete buff;
    }
 
    if (saveTGraphsToFile) {
      TPipA_CEx -> Write("TPipA_CEx");
    }
    delete TPipA_CEx;
   
  }
 
  // kIHNFtAbs, pip + A                                                            PipA_Abs (just for developmental purposes)
  {
    TGraph2D * TPipA_Abs;
 
    const int pipAAbs_nfiles = 18;
    const int pipAAbs_nuclei[pipAAbs_nfiles] = {1, 2, 3, 4, 7, 9, 12, 16, 27, 48, 56, 58, 63, 93, 120, 165, 181, 209};
    const int pipAAbs_npoints = 111;
 
    TPipA_Abs = new TGraph2D(pipAAbs_npoints);
    TPipA_Abs->SetNameTitle("TPipA_Abs","TPipA_Abs");
    TPipA_Abs->SetDirectory(0);
 
    int ipoint=0;
    double x, y;
 
    for(int ifile=0; ifile < pipAAbs_nfiles; ifile++) {
      ostringstream ADep_datafile;
      int nucleus = pipAAbs_nuclei[ifile];
      ADep_datafile << data_dir << "/tot_xsec/pipA_abs/pip" << nucleus << "_abs.txt";
      TGraph * buff = new TGraph(ADep_datafile.str().c_str());
      buff->SetNameTitle("buff","buff");
      for(int i=0; i < buff->GetN(); i++) {
        buff -> GetPoint(i,x,y);
        TPipA_Abs -> SetPoint(ipoint,(double)nucleus,x,y);
        ipoint++;
      }
      delete buff;
    }
 
    if (saveTGraphsToFile) {
      TPipA_Abs -> Write("TPipA_Abs");
    }
    delete TPipA_Abs;
   
  }
 
  // kIHNFtElas, pip + A                                                            PipA_Elas (just for developmental purposes)
  {
    TGraph2D * TPipA_Elas;
 
    const int pipAElas_nfiles = 18;
    const int pipAElas_nuclei[pipAElas_nfiles] = {1, 2, 3, 4, 7, 9, 12, 16, 27, 48, 56, 58, 63, 93, 120, 165, 181, 209};
    const int pipAElas_npoints = 125;
 
    TPipA_Elas = new TGraph2D(pipAElas_npoints);
    TPipA_Elas->SetNameTitle("TPipA_Elas","TPipA_Elas");
    TPipA_Elas->SetDirectory(0);
 
    int ipoint=0;
    double x, y;
 
    for(int ifile=0; ifile < pipAElas_nfiles; ifile++) {
      ostringstream ADep_datafile;
      int nucleus = pipAElas_nuclei[ifile];
      ADep_datafile << data_dir << "/tot_xsec/pipA_elas/pip" << nucleus << "_elas.txt";
      TGraph * buff = new TGraph(ADep_datafile.str().c_str());
      buff->SetNameTitle("buff","buff");
      for(int i=0; i < buff->GetN(); i++) {
        buff -> GetPoint(i,x,y);
        TPipA_Elas -> SetPoint(ipoint,(double)nucleus,x,y);
        ipoint++;
      }
      delete buff;
    }
 
    if (saveTGraphsToFile) {
      TPipA_Elas -> Write("TPipA_Elas");
    }
    delete TPipA_Elas; 
  }
 
  // kIHNFtInelas, pip + A                                                            PipA_Inelas (just for developmental purposes)
  {
    TGraph2D * TPipA_Inelas;
 
    const int pipAInelas_nfiles = 20;
    const int pipAInelas_nuclei[pipAInelas_nfiles] = {1, 2, 3, 4, 7, 9, 12, 16, 27, 40, 48, 56, 58, 63, 93, 120, 165, 181, 208, 209};
    const int pipAInelas_npoints = 118;
 
    TPipA_Inelas = new TGraph2D(pipAInelas_npoints);
    TPipA_Inelas->SetNameTitle("TPipA_Inelas","TPipA_Inelas");
    TPipA_Inelas->SetDirectory(0);
 
    int ipoint=0;
    double x, y;
 
    for(int ifile=0; ifile < pipAInelas_nfiles; ifile++) {
      ostringstream ADep_datafile;
      int nucleus = pipAInelas_nuclei[ifile];
      ADep_datafile << data_dir << "/tot_xsec/pipA_inelas/pip" << nucleus << "_inelas.txt";
      TGraph * buff = new TGraph(ADep_datafile.str().c_str());
      buff->SetNameTitle("buff","buff");
      for(int i=0; i < buff->GetN(); i++) {
        buff -> GetPoint(i,x,y);
        TPipA_Inelas -> SetPoint(ipoint,(double)nucleus,x,y);
        ipoint++;
      }
      delete buff;
    }
 
    if (saveTGraphsToFile) {
      TPipA_Inelas -> Write("TPipA_Inelas");
    }
    delete TPipA_Inelas;
  }
 
 
  /*
  // kIHNFtElas, pip + A                                                            PipA_Elas_frac
  {
    const int pipAElas_f_nfiles = 18;
    const int pipAElas_f_nuclei[pipAElas_f_nfiles] = {1, 2, 3, 4, 7, 9, 12, 16, 27, 48, 56, 58, 63, 93, 120, 165, 181, 209};
    const int pipAElas_f_npoints = 125;
 
    TfracPipA_Elas = new TGraph2D(pipAElas_f_npoints);
    TfracPipA_Elas->SetNameTitle("TfracPipA_Elas","TfracPipA_Elas");
    TfracPipA_Elas->SetDirectory(0);
 
    int ipoint=0;
    double x, y;
 
    for(int ifile=0; ifile < pipAElas_f_nfiles; ifile++) {
      ostringstream ADep_datafile;
      int nucleus = pipAElas_f_nuclei[ifile];
      ADep_datafile << data_dir << "/tot_xsec/pipA_elas_frac/pip" << nucleus << "_elas_frac.txt";
      TGraph * buff = new TGraph(ADep_datafile.str().c_str());
      buff->SetNameTitle("buff","buff");
      for(int i=0; i < buff->GetN(); i++) {
        buff -> GetPoint(i,x,y);
        TfracPipA_Elas -> SetPoint(ipoint,(double)nucleus,x,y);
        ipoint++;
      }
      delete buff;
    }
 
    if (saveTGraphsToFile) {
      TfracPipA_Elas -> Write("TfracPipA_Elas");
    }
 
  }
  */
 
  // kIHNFtInelas, pip + A                                                            PipA_Inelas_frac
  {
    const int pipAInelas_f_nfiles = 20;
    const int pipAInelas_f_nuclei[pipAInelas_f_nfiles] = {1, 2, 3, 4, 7, 9, 12, 16, 27, 40, 48, 56, 58, 63, 93, 120, 165, 181, 208, 209};
    const int pipAInelas_f_npoints = 118;
 
    TfracPipA_Inelas = new TGraph2D(pipAInelas_f_npoints);
    TfracPipA_Inelas->SetNameTitle("TfracPipA_Inelas","TfracPipA_Inelas");
    TfracPipA_Inelas->SetDirectory(0);
 
    int ipoint=0;
    double x, y;
 
    for(int ifile=0; ifile < pipAInelas_f_nfiles; ifile++) {
      ostringstream ADep_datafile;
      int nucleus = pipAInelas_f_nuclei[ifile];
      ADep_datafile << data_dir << "/tot_xsec/pipA_inelas_frac/pip" << nucleus << "_inelas_frac.txt";
      TGraph * buff = new TGraph(ADep_datafile.str().c_str());
      buff->SetNameTitle("buff","buff");
      for(int i=0; i < buff->GetN(); i++) {
        buff -> GetPoint(i,x,y);
        TfracPipA_Inelas -> SetPoint(ipoint,(double)nucleus,x,y);
        ipoint++;
      }
      delete buff;
    }
 
    if (saveTGraphsToFile) {
      TfracPipA_Inelas -> Write("TfracPipA_Inelas");
    }
 
  }
 
 
  // kIHNFtPiPro, pip + A                                                            PipA_PiPro_frac
   {
    const int pipAPiPro_f_nfiles = 17;
    const int pipAPiPro_f_nuclei[pipAPiPro_f_nfiles] = {1, 2, 3, 4, 7, 9, 12, 16, 48, 56, 58, 63, 93, 120, 165, 181, 209};
    const int pipAPiPro_f_npoints = 76;
 
    TfracPipA_PiPro = new TGraph2D(pipAPiPro_f_npoints);
    TfracPipA_PiPro->SetNameTitle("TfracPipA_PiPro","TfracPipA_PiPro");
    TfracPipA_PiPro->SetDirectory(0);
 
    int ipoint=0;
    double x, y;
 
    for(int ifile=0; ifile < pipAPiPro_f_nfiles; ifile++) {
      ostringstream ADep_datafile;
      int nucleus = pipAPiPro_f_nuclei[ifile];
      ADep_datafile << data_dir << "/tot_xsec/pipA_pipro_frac/pip" << nucleus << "_pipro_frac.txt";
      TGraph * buff = new TGraph(ADep_datafile.str().c_str());
      buff->SetNameTitle("buff","buff");
      for(int i=0; i < buff->GetN(); i++) {
        buff -> GetPoint(i,x,y);
        TfracPipA_PiPro -> SetPoint(ipoint,(double)nucleus,x,y);
        ipoint++;
      }
      delete buff;
    }
 
    if (saveTGraphsToFile) {
      TfracPipA_PiPro -> Write("TfracPipA_PiPro");
    }
   }
 
   TGraphs_file.Close();
 
   LOG("INukeData", pINFO)  << "Done building x-section splines...";
   
}

References fFracKA_Abs, fFracKA_CEx, fFracKA_Elas, fFracKA_Inel, fFracKA_Tot, fFracNA_Abs, fFracNA_CEx, fFracNA_Cmp, fFracNA_Inel, fFracNA_PiPro, fFracNA_Tot, fFracPA_Abs, fFracPA_CEx, fFracPA_Cmp, fFracPA_Inel, fFracPA_PiPro, fFracPA_Tot, fhN2dXSecGamPi0N_Inelas, fhN2dXSecGamPi0P_Inelas, fhN2dXSecGamPimP_Inelas, fhN2dXSecGamPipN_Inelas, fhN2dXSecKpN_CEx, fhN2dXSecKpN_Elas, fhN2dXSecKpP_Elas, fhN2dXSecNP_Elas, fhN2dXSecPi0N_Elas, fhN2dXSecPimN_Elas, fhN2dXSecPiN_Abs, fhN2dXSecPiN_CEx, fhN2dXSecPipN_Elas, fhN2dXSecPP_Elas, fXSecGamn_fs, fXSecGamN_Tot, fXSecGamp_fs, fXSecKpN_Abs, fXSecKpn_CEx, fXSecKpn_Elas, fXSecKpN_Tot, fXSecKpp_Elas, fXSecNn_Cmp, fXSecNn_Elas, fXSecNn_Reac, fXSecNn_Tot, fXSecPi0d_Abs, fXSecPi0n_CEx, fXSecPi0n_Elas, fXSecPi0n_Reac, fXSecPi0n_Tot, fXSecPi0p_CEx, fXSecPi0p_Elas, fXSecPi0p_Reac, fXSecPi0p_Tot, fXSecPipd_Abs, fXSecPipn_CEx, fXSecPipn_Elas, fXSecPipn_Reac, fXSecPipn_Tot, fXSecPipp_CEx, fXSecPipp_Elas, fXSecPipp_Reac, fXSecPipp_Tot, fXSecPn_Cmp, fXSecPn_Elas, fXSecPn_Reac, fXSecPn_Tot, fXSecPp_Cmp, fXSecPp_Elas, fXSecPp_Reac, fXSecPp_Tot, LOG, pDEBUG, pINFO, pNOTICE, ReadhNFile(), TfracPipA_Abs, TfracPipA_CEx, TfracPipA_Inelas, TfracPipA_PiPro, and Write().

Referenced by INukeHadroData2018().

ReadhNFile()

void INukeHadroData2018::ReadhNFile ( string filename, double ke, int npoints, int & curr_point, double * costh_array, double * xsec_array, int cols )

private

Definition at line 1180 of file INukeHadroData2018.cxx.

{
  // open
  std::ifstream hN_stream(filename.c_str(), ios::in);
  if(!hN_stream.good()) {
      LOG("INukeData", pERROR)
          << "Error reading INTRANUKE/hN data from: " << filename;
      return;
  }
 
  if(cols<2) {
    LOG("INukeData", pERROR)
      << "Error reading INTRANUKE/hN data from: " << filename;
    LOG("INukeData", pERROR)
      << "Too few columns: " << cols;
    return;
  }
 
  LOG("INukeData", pINFO)
     << "Reading INTRANUKE/hN data from: " << filename;
 
  // skip initial comments
  char cbuf[501];
  hN_stream.getline(cbuf,400);
  hN_stream.getline(cbuf,400);
  hN_stream.getline(cbuf,400);
 
  // read
  double angle = 0;
  double xsec  = 0;
  double trash = 0;
 
  for(int ip = 0; ip < npoints; ip++) {
     hN_stream >> angle >> xsec;
 
     for(int ic = 0; ic < (cols-2); ic++) {
       hN_stream >> trash;
     }
 
     LOG("INukeData", pDEBUG)
       << "Adding data point: (KE = " << ke << " MeV, angle = "
       << angle << ", sigma = " << xsec << " mbarn)";
     costh_array[ip] = TMath::Cos(angle*kPi/180.);
     xsec_array [curr_point] = xsec;
     curr_point++;
  }
}

References genie::constants::kPi, LOG, pDEBUG, pERROR, and pINFO.

Referenced by LoadCrossSections().

XSec() [1/2]

double INukeHadroData2018::XSec	(	int	hpdgc,
		int	tgt,
		int	nprod,
		INukeFateHN_t	rxnType,
		double	ke,
		double	costh ) const

Definition at line 1230 of file INukeHadroData2018.cxx.

{
// inputs
//      fate    : h+N fate code
//      hpdgc   : h PDG code
//      tgtpdgc : N PDG code
//      nppdgc  : product N PDG code
//      ke      : kinetic energy (MeV)
//      costh   : cos(scattering angle)
// returns
//      xsec    : mbarn
 
  double ke_eval    = ke;
  double costh_eval = costh;
 
  costh_eval = TMath::Min(costh,  1.);
  costh_eval = TMath::Max(costh_eval, -1.);
 
  if(fate==kIHNFtElas) {
 
     if( (hpdgc==kPdgProton  && tgtpdgc==kPdgProton) ||
         (hpdgc==kPdgNeutron && tgtpdgc==kPdgNeutron) )
     {
       ke_eval = TMath::Min(ke_eval, 999.);
       ke_eval = TMath::Max(ke_eval,  50.);
       return fhN2dXSecPP_Elas->Evaluate(ke_eval, costh_eval);
     }
     else
     if( (hpdgc==kPdgProton  && tgtpdgc==kPdgNeutron) ||
         (hpdgc==kPdgNeutron && tgtpdgc==kPdgProton) )
     {
       ke_eval = TMath::Min(ke_eval, 999.);
       ke_eval = TMath::Max(ke_eval,  50.);
       return fhN2dXSecNP_Elas->Evaluate(ke_eval, costh_eval);
     }
     else
     if(hpdgc==kPdgPiP)
     {
       ke_eval = TMath::Min(ke_eval, 1499.);
       ke_eval = TMath::Max(ke_eval,   10.);
       return fhN2dXSecPipN_Elas->Evaluate(ke_eval, costh_eval);
     }
     else
     if(hpdgc==kPdgPi0)
     {
       ke_eval = TMath::Min(ke_eval, 1499.);
       ke_eval = TMath::Max(ke_eval,   10.);
       return fhN2dXSecPi0N_Elas->Evaluate(ke_eval, costh_eval);
     }
     else
     if(hpdgc==kPdgPiM)
     {
       ke_eval = TMath::Min(ke_eval, 1499.);
       ke_eval = TMath::Max(ke_eval,   10.);
       return fhN2dXSecPimN_Elas->Evaluate(ke_eval, costh_eval);
     }
     else
     if(hpdgc==kPdgKP && tgtpdgc==kPdgNeutron)
     {
       ke_eval = TMath::Min(ke_eval, 1799.);
       ke_eval = TMath::Max(ke_eval,  100.);
       return fhN2dXSecKpN_Elas->Evaluate(ke_eval, costh_eval);
     }
     else
     if(hpdgc==kPdgKP && tgtpdgc==kPdgProton)
     {
       ke_eval = TMath::Min(ke_eval, 1799.);
       ke_eval = TMath::Max(ke_eval,  100.);
       return fhN2dXSecKpP_Elas->Evaluate(ke_eval, costh_eval);
     }
  }
 
  else if(fate == kIHNFtCEx) {
    if( (hpdgc==kPdgPiP || hpdgc==kPdgPi0 || hpdgc==kPdgPiM) &&
         (tgtpdgc==kPdgProton || tgtpdgc==kPdgNeutron) )
     {
        ke_eval = TMath::Min(ke_eval, 1499.);
        ke_eval = TMath::Max(ke_eval,   10.);
        return fhN2dXSecPiN_CEx->Evaluate(ke_eval, costh_eval);
     }
    else if( (hpdgc == kPdgProton && tgtpdgc == kPdgProton) ||
             (hpdgc == kPdgNeutron && tgtpdgc == kPdgNeutron) )
      {
        LOG("INukeData", pWARN)  << "Inelastic pp does not exist!";
        ke_eval = TMath::Min(ke_eval, 999.);
        ke_eval = TMath::Max(ke_eval,  50.);
        return fhN2dXSecPP_Elas->Evaluate(ke_eval, costh_eval);
      }
    else if( (hpdgc == kPdgProton && tgtpdgc == kPdgNeutron) ||
             (hpdgc == kPdgNeutron && tgtpdgc == kPdgProton) )
      {
        ke_eval = TMath::Min(ke_eval, 999.);
        ke_eval = TMath::Max(ke_eval,  50.);
        return fhN2dXSecNP_Elas->Evaluate(ke_eval, costh_eval);
      }
    else if(hpdgc == kPdgKP && tgtpdgc == kPdgNeutron) {
        ke_eval = TMath::Min(ke_eval, 1799.);
        ke_eval = TMath::Max(ke_eval,  100.);
        return fhN2dXSecKpN_CEx->Evaluate(ke_eval, costh_eval);
    }
  }
 
  else if(fate == kIHNFtAbs) {
    if( (hpdgc==kPdgPiP || hpdgc==kPdgPi0 || hpdgc==kPdgPiM) &&
         (tgtpdgc==kPdgProton || tgtpdgc==kPdgNeutron) )
     {
        ke_eval = TMath::Min(ke_eval, 499.);
        ke_eval = TMath::Max(ke_eval,  50.);
        return fhN2dXSecPiN_Abs->Evaluate(ke_eval, costh_eval);
     }
    if(hpdgc==kPdgKP) return 1.;  //isotropic since no data ???
  }
 
  else if(fate == kIHNFtInelas) {
    if( hpdgc==kPdgGamma && tgtpdgc==kPdgProton  &&nppdgc==kPdgProton  )
    {
       ke_eval = TMath::Min(ke_eval, 1199.);
       ke_eval = TMath::Max(ke_eval,  160.);
       return fhN2dXSecGamPi0P_Inelas->Evaluate(ke_eval, costh_eval);
    }
    else
    if( hpdgc==kPdgGamma && tgtpdgc==kPdgProton  && nppdgc==kPdgNeutron )
    {
       ke_eval = TMath::Min(ke_eval, 1199.);
       ke_eval = TMath::Max(ke_eval,  160.);
       return fhN2dXSecGamPipN_Inelas->Evaluate(ke_eval, costh_eval);
    }
    else
    if( hpdgc==kPdgGamma && tgtpdgc==kPdgNeutron && nppdgc==kPdgProton  )
    {
       ke_eval = TMath::Min(ke_eval, 1199.);
       ke_eval = TMath::Max(ke_eval,  160.);
       return fhN2dXSecGamPimP_Inelas->Evaluate(ke_eval, costh_eval);
    }
    else
    if( hpdgc==kPdgGamma && tgtpdgc==kPdgNeutron && nppdgc==kPdgNeutron )
    {
       ke_eval = TMath::Min(ke_eval, 1199.);
       ke_eval = TMath::Max(ke_eval,  160.);
      return fhN2dXSecGamPi0N_Inelas->Evaluate(ke_eval, costh_eval);
    }
  }
 
  return 0;
}

References fhN2dXSecGamPi0N_Inelas, fhN2dXSecGamPi0P_Inelas, fhN2dXSecGamPimP_Inelas, fhN2dXSecGamPipN_Inelas, fhN2dXSecKpN_CEx, fhN2dXSecKpN_Elas, fhN2dXSecKpP_Elas, fhN2dXSecNP_Elas, fhN2dXSecPi0N_Elas, fhN2dXSecPimN_Elas, fhN2dXSecPiN_Abs, fhN2dXSecPiN_CEx, fhN2dXSecPipN_Elas, fhN2dXSecPP_Elas, genie::kIHNFtAbs, genie::kIHNFtCEx, genie::kIHNFtElas, genie::kIHNFtInelas, genie::kPdgGamma, genie::kPdgKP, genie::kPdgNeutron, genie::kPdgPi0, genie::kPdgPiM, genie::kPdgPiP, genie::kPdgProton, LOG, and pWARN.

Referenced by Frac(), and IntBounce().

XSec() [2/2]

double INukeHadroData2018::XSec	(	int	hpdgc,
		INukeFateHN_t	fate,
		double	ke,
		int	targA,
		int	targZ ) const

Definition at line 1500 of file INukeHadroData2018.cxx.

{
// return the x-section for the input fate for the particle with the input pdg
// code at the input kinetic energy
//
  ke = TMath::Max(fMinKinEnergy,   ke);
  ke = TMath::Min(fMaxKinEnergyHN, ke);
 
  LOG("INukeData", pDEBUG)  << "Querying hN cross section at ke = " << ke;
 
  double xsec=0;
 
    if (hpdgc == kPdgPiP) {
    /* handle pi+ */
         if (fate == kIHNFtCEx   ) {xsec = TMath::Max(0., fXSecPipp_CEx  -> Evaluate(ke)) *  targZ;
                                    xsec+= TMath::Max(0., fXSecPipn_CEx  -> Evaluate(ke)) * (targA-targZ);
                                    return xsec;}
    else if (fate == kIHNFtElas  ) {xsec = TMath::Max(0., fXSecPipp_Elas -> Evaluate(ke)) *  targZ;
                                    xsec+= TMath::Max(0., fXSecPipn_Elas -> Evaluate(ke)) * (targA-targZ);
                                    return xsec;}
    else if (fate == kIHNFtInelas) {xsec = TMath::Max(0., fXSecPipp_Reac -> Evaluate(ke)) *  targZ;
                                    xsec+= TMath::Max(0., fXSecPipn_Reac -> Evaluate(ke)) * (targA-targZ);
                                    return xsec;}
    else if (fate == kIHNFtAbs   ) {xsec = TMath::Max(0., fXSecPipd_Abs  -> Evaluate(ke)) *  targA;
                                    return xsec;}
    else {
     LOG("INukeData", pWARN)
        << "Pi+'s don't have this fate: " << INukeHadroFates::AsString(fate);
     return 0;
    }
 
  } else if (hpdgc == kPdgPiM) {
    /* handle pi- */
         if (fate == kIHNFtCEx   ) {xsec = TMath::Max(0., fXSecPipn_CEx  -> Evaluate(ke)) *  targZ;
                                    xsec+= TMath::Max(0., fXSecPipp_CEx  -> Evaluate(ke)) * (targA-targZ);
                                    return xsec;}
    else if (fate == kIHNFtElas  ) {xsec = TMath::Max(0., fXSecPipn_Elas -> Evaluate(ke)) *  targZ;
                                    xsec+= TMath::Max(0., fXSecPipp_Elas -> Evaluate(ke)) * (targA-targZ);
                                    return xsec;}
    else if (fate == kIHNFtInelas) {xsec = TMath::Max(0., fXSecPipn_Reac -> Evaluate(ke)) *  targZ;
                                    xsec+= TMath::Max(0., fXSecPipp_Reac -> Evaluate(ke)) * (targA-targZ);
                                    return xsec;}
    else if (fate == kIHNFtAbs   ) {xsec = TMath::Max(0., fXSecPipd_Abs  -> Evaluate(ke)) *  targA;
                                    return xsec;}
    else {
     LOG("INukeData", pWARN)
        << "Pi-'s don't have this fate: " << INukeHadroFates::AsString(fate);
     return 0;
    }
 
  } else if (hpdgc == kPdgPi0) {
    /* handle pi0 */
         if (fate == kIHNFtCEx   ) {xsec = TMath::Max(0., fXSecPi0p_CEx  -> Evaluate(ke)) *  targZ;
                                    xsec+= TMath::Max(0., fXSecPi0n_CEx  -> Evaluate(ke)) * (targA-targZ);
                                    return xsec;}
    else if (fate == kIHNFtElas  ) {xsec = TMath::Max(0., fXSecPi0p_Elas -> Evaluate(ke)) *  targZ;
                                    xsec+= TMath::Max(0., fXSecPi0n_Elas -> Evaluate(ke)) * (targA-targZ);
                                    return xsec;}
    else if (fate == kIHNFtInelas) {xsec = TMath::Max(0., fXSecPi0p_Reac -> Evaluate(ke)) *  targZ;
                                    xsec+= TMath::Max(0., fXSecPi0n_Reac -> Evaluate(ke)) * (targA-targZ);
                                    return xsec;}
    else if (fate == kIHNFtAbs   ) {xsec = TMath::Max(0., fXSecPi0d_Abs  -> Evaluate(ke)) *  targA;
                                    return xsec;}
    else {
     LOG("INukeData", pWARN)
        << "Pi0's don't have this fate: " << INukeHadroFates::AsString(fate);
     return 0;
    }
 
  } else if (hpdgc == kPdgProton) {
    /* handle protons */
      if (fate == kIHNFtElas  ) {xsec = TMath::Max(0., fXSecPp_Elas -> Evaluate(ke)) *  targZ;
                                    xsec+= TMath::Max(0., fXSecPn_Elas -> Evaluate(ke)) * (targA-targZ);
                                    return xsec;}
    else if (fate == kIHNFtInelas) {xsec = TMath::Max(0., fXSecPp_Reac -> Evaluate(ke)) *  targZ;
                                    xsec+= TMath::Max(0., fXSecPn_Reac -> Evaluate(ke)) * (targA-targZ);
                                    return xsec;}
    else if (fate == kIHNFtCmp) {xsec = TMath::Max(0., fXSecPp_Cmp -> Evaluate(ke)) *  targZ;
                                    xsec+= TMath::Max(0., fXSecPn_Cmp -> Evaluate(ke)) * (targA-targZ);
                                    return xsec;}
    else {
     LOG("INukeData", pWARN)
        << "Protons don't have this fate: " << INukeHadroFates::AsString(fate);
     return 0;
    }
 
  } else if (hpdgc == kPdgNeutron) {
    /* handle protons */
         if (fate == kIHNFtElas  ) {xsec = TMath::Max(0., fXSecPn_Elas -> Evaluate(ke)) *  targZ;
                                    xsec+= TMath::Max(0., fXSecNn_Elas -> Evaluate(ke)) * (targA-targZ);
                                    return xsec;}
    else if (fate == kIHNFtInelas) {xsec = TMath::Max(0., fXSecPn_Reac -> Evaluate(ke)) *  targZ;
                                    xsec+= TMath::Max(0., fXSecNn_Reac -> Evaluate(ke)) * (targA-targZ);
                                    return xsec;}
    else if (fate == kIHNFtCmp) {xsec = TMath::Max(0., fXSecPp_Cmp -> Evaluate(ke)) *  targZ;
                                    xsec+= TMath::Max(0., fXSecPn_Cmp -> Evaluate(ke)) * (targA-targZ);
                                    return xsec;}
    else {
     LOG("INukeData", pWARN)
        << "Neutrons don't have this fate: " << INukeHadroFates::AsString(fate);
     return 0;
    }
    //Adding here kaons, why elastic only on protons? hA or hN? No _Reac for kaons...
    } else if (hpdgc == kPdgKP) {
    /* handle K+ */
        if (fate == kIHNFtCEx   ) {xsec = TMath::Max(0., fXSecKpn_CEx  -> Evaluate(ke)) *  targZ;
                                    xsec+= TMath::Max(0., fXSecKpn_CEx  -> Evaluate(ke)) * (targA-targZ);
                                    return xsec;}
        else if (fate == kIHNFtElas  ) {xsec = TMath::Max(0., fXSecKpn_Elas -> Evaluate(ke)) *  targZ;
                                    xsec+= TMath::Max(0., fXSecKpn_Elas -> Evaluate(ke)) * (targA-targZ);
                                    return xsec;}
        /*else if (fate == kIHNFtAbs   ) {xsec = TMath::Max(0., fXSecKpd_Abs  -> Evaluate(ke)) *  targA;
                                    return xsec;}*/
        else {
                LOG("INukeData", pWARN)
                << "K+'s don't have this fate: " << INukeHadroFates::AsString(fate);
        return 0;
    }
    //------------------------------------------------
         /*   }  else if (hpdgc == kPdgGamma) {
    / * handle gamma * /
         if (fate == kIHNFtInelas) {xsec = TMath::Max(0., fXSecGamp_fs   -> Evaluate(ke)) *  targZ;
                                    xsec+= TMath::Max(0., fXSecGamn_fs   -> Evaluate(ke)) * (targA-targZ);
                                    return xsec;}
    else {
     LOG("INukeData", pWARN)
        << "Gamma's don't have this fate: " << INukeHadroFates::AsString(fate);
     return 0;
     }*/
   }
  LOG("INukeData", pWARN)
      << "Can't handle particles with pdg code = " << hpdgc;
 
  return 0;
}

References genie::INukeHadroFates::AsString(), fMaxKinEnergyHN, fMinKinEnergy, fXSecKpn_CEx, fXSecKpn_Elas, fXSecNn_Elas, fXSecNn_Reac, fXSecPi0d_Abs, fXSecPi0n_CEx, fXSecPi0n_Elas, fXSecPi0n_Reac, fXSecPi0p_CEx, fXSecPi0p_Elas, fXSecPi0p_Reac, fXSecPipd_Abs, fXSecPipn_CEx, fXSecPipn_Elas, fXSecPipn_Reac, fXSecPipp_CEx, fXSecPipp_Elas, fXSecPipp_Reac, fXSecPn_Cmp, fXSecPn_Elas, fXSecPn_Reac, fXSecPp_Cmp, fXSecPp_Elas, fXSecPp_Reac, genie::kIHNFtAbs, genie::kIHNFtCEx, genie::kIHNFtCmp, genie::kIHNFtElas, genie::kIHNFtInelas, genie::kPdgKP, genie::kPdgNeutron, genie::kPdgPi0, genie::kPdgPiM, genie::kPdgPiP, genie::kPdgProton, LOG, pDEBUG, and pWARN.

XSecGamn_fs()

const Spline * genie::INukeHadroData2018::XSecGamn_fs ( void ) const

inline

Definition at line 104 of file INukeHadroData2018.h.

{ return fXSecGamn_fs;      }

References fXSecGamn_fs.

XSecGamN_Tot()

const Spline * genie::INukeHadroData2018::XSecGamN_Tot ( void ) const

inline

Definition at line 105 of file INukeHadroData2018.h.

{ return fXSecGamN_Tot;     }

References fXSecGamN_Tot.

XSecGamp_fs()

const Spline * genie::INukeHadroData2018::XSecGamp_fs ( void ) const

inline

Definition at line 103 of file INukeHadroData2018.h.

{ return fXSecGamp_fs;      }

References fXSecGamp_fs.

XSecKpN_Abs()

const Spline * genie::INukeHadroData2018::XSecKpN_Abs ( void ) const

inline

Definition at line 101 of file INukeHadroData2018.h.

{ return fXSecKpN_Abs;      } //not implemented

References fXSecKpN_Abs.

XSecKpn_CEx()

const Spline * genie::INukeHadroData2018::XSecKpn_CEx ( void ) const

inline

Definition at line 99 of file INukeHadroData2018.h.

{ return fXSecKpn_CEx;      }

References fXSecKpn_CEx.

XSecKpn_Elas()

const Spline * genie::INukeHadroData2018::XSecKpn_Elas ( void ) const

inline

Definition at line 98 of file INukeHadroData2018.h.

{ return fXSecKpn_Elas;     }

References fXSecKpn_Elas.

XSecKpN_Tot()

const Spline * genie::INukeHadroData2018::XSecKpN_Tot ( void ) const

inline

Definition at line 102 of file INukeHadroData2018.h.

{ return fXSecKpN_Tot;      }

References fXSecKpN_Tot.

XSecKpp_Elas()

const Spline * genie::INukeHadroData2018::XSecKpp_Elas ( void ) const

inline

Definition at line 100 of file INukeHadroData2018.h.

{ return fXSecKpp_Elas;     }

References fXSecKpp_Elas.

XSecNn_Elas()

const Spline * genie::INukeHadroData2018::XSecNn_Elas ( void ) const

inline

Definition at line 96 of file INukeHadroData2018.h.

{ return fXSecNn_Elas;      }

References fXSecNn_Elas.

XSecNn_Reac()

const Spline * genie::INukeHadroData2018::XSecNn_Reac ( void ) const

inline

Definition at line 97 of file INukeHadroData2018.h.

{ return fXSecNn_Reac;      }

References fXSecNn_Reac.

XSecNn_Tot()

const Spline * genie::INukeHadroData2018::XSecNn_Tot ( void ) const

inline

Definition at line 95 of file INukeHadroData2018.h.

{ return fXSecNn_Tot;       }

References fXSecNn_Tot.

XSecPi0n_CEx()

const Spline * genie::INukeHadroData2018::XSecPi0n_CEx ( void ) const

inline

Definition at line 81 of file INukeHadroData2018.h.

{ return fXSecPi0n_CEx;     }

References fXSecPi0n_CEx.

XSecPi0n_Elas()

const Spline * genie::INukeHadroData2018::XSecPi0n_Elas ( void ) const

inline

Definition at line 82 of file INukeHadroData2018.h.

{ return fXSecPi0n_Elas;    }

References fXSecPi0n_Elas.

XSecPi0n_Reac()

const Spline * genie::INukeHadroData2018::XSecPi0n_Reac ( void ) const

inline

Definition at line 83 of file INukeHadroData2018.h.

{ return fXSecPi0n_Reac;    }

References fXSecPi0n_Reac.

XSecPi0n_Tot()

const Spline * genie::INukeHadroData2018::XSecPi0n_Tot ( void ) const

inline

Definition at line 80 of file INukeHadroData2018.h.

{ return fXSecPi0n_Tot;     }

References fXSecPi0n_Tot.

XSecPi0p_Abs()

const Spline * genie::INukeHadroData2018::XSecPi0p_Abs ( void ) const

inline

Definition at line 88 of file INukeHadroData2018.h.

{ return fXSecPi0d_Abs;     }

References fXSecPi0d_Abs.

XSecPi0p_CEx()

const Spline * genie::INukeHadroData2018::XSecPi0p_CEx ( void ) const

inline

Definition at line 85 of file INukeHadroData2018.h.

{ return fXSecPi0p_CEx;     }

References fXSecPi0p_CEx.

XSecPi0p_Elas()

const Spline * genie::INukeHadroData2018::XSecPi0p_Elas ( void ) const

inline

Definition at line 86 of file INukeHadroData2018.h.

{ return fXSecPi0p_Elas;    }

References fXSecPi0p_Elas.

XSecPi0p_Reac()

const Spline * genie::INukeHadroData2018::XSecPi0p_Reac ( void ) const

inline

Definition at line 87 of file INukeHadroData2018.h.

{ return fXSecPi0p_Reac;    }

References fXSecPi0p_Reac.

XSecPi0p_Tot()

const Spline * genie::INukeHadroData2018::XSecPi0p_Tot ( void ) const

inline

Definition at line 84 of file INukeHadroData2018.h.

{ return fXSecPi0p_Tot;     }

References fXSecPi0p_Tot.

XSecPipn_CEx()

const Spline * genie::INukeHadroData2018::XSecPipn_CEx ( void ) const

inline

Definition at line 72 of file INukeHadroData2018.h.

{ return fXSecPipn_CEx;     }

References fXSecPipn_CEx.

XSecPipn_Elas()

const Spline * genie::INukeHadroData2018::XSecPipn_Elas ( void ) const

inline

Definition at line 73 of file INukeHadroData2018.h.

{ return fXSecPipn_Elas;    }

References fXSecPipn_Elas.

XSecPipn_Reac()

const Spline * genie::INukeHadroData2018::XSecPipn_Reac ( void ) const

inline

Definition at line 74 of file INukeHadroData2018.h.

{ return fXSecPipn_Reac;    }

References fXSecPipn_Reac.

XSecPipn_Tot()

const Spline * genie::INukeHadroData2018::XSecPipn_Tot ( void ) const

inline

Definition at line 71 of file INukeHadroData2018.h.

{ return fXSecPipn_Tot;     }

References fXSecPipn_Tot.

XSecPipp_Abs()

const Spline * genie::INukeHadroData2018::XSecPipp_Abs ( void ) const

inline

Definition at line 79 of file INukeHadroData2018.h.

{ return fXSecPipd_Abs;     }

References fXSecPipd_Abs.

XSecPipp_CEx()

const Spline * genie::INukeHadroData2018::XSecPipp_CEx ( void ) const

inline

Definition at line 76 of file INukeHadroData2018.h.

{ return fXSecPipp_CEx;     }

References fXSecPipp_CEx.

XSecPipp_Elas()

const Spline * genie::INukeHadroData2018::XSecPipp_Elas ( void ) const

inline

Definition at line 77 of file INukeHadroData2018.h.

{ return fXSecPipp_Elas;    }

References fXSecPipp_Elas.

XSecPipp_Reac()

const Spline * genie::INukeHadroData2018::XSecPipp_Reac ( void ) const

inline

Definition at line 78 of file INukeHadroData2018.h.

{ return fXSecPipp_Reac;    }

References fXSecPipp_Reac.

XSecPipp_Tot()

const Spline * genie::INukeHadroData2018::XSecPipp_Tot ( void ) const

inline

Definition at line 75 of file INukeHadroData2018.h.

{ return fXSecPipp_Tot;     }

References fXSecPipp_Tot.

XSecPn_Elas()

const Spline * genie::INukeHadroData2018::XSecPn_Elas ( void ) const

inline

Definition at line 93 of file INukeHadroData2018.h.

{ return fXSecPn_Elas;      }

References fXSecPn_Elas.

XSecPn_Reac()

const Spline * genie::INukeHadroData2018::XSecPn_Reac ( void ) const

inline

Definition at line 94 of file INukeHadroData2018.h.

{ return fXSecPn_Reac;      }

References fXSecPn_Reac.

XSecPn_Tot()

const Spline * genie::INukeHadroData2018::XSecPn_Tot ( void ) const

inline

Definition at line 92 of file INukeHadroData2018.h.

{ return fXSecPn_Tot;       }

References fXSecPn_Tot.

XSecPp_Elas()

const Spline * genie::INukeHadroData2018::XSecPp_Elas ( void ) const

inline

Definition at line 90 of file INukeHadroData2018.h.

{ return fXSecPp_Elas;      }

References fXSecPp_Elas.

XSecPp_Reac()

const Spline * genie::INukeHadroData2018::XSecPp_Reac ( void ) const

inline

Definition at line 91 of file INukeHadroData2018.h.

{ return fXSecPp_Reac;      }

References fXSecPp_Reac.

XSecPp_Tot()

const Spline * genie::INukeHadroData2018::XSecPp_Tot ( void ) const

inline

Definition at line 89 of file INukeHadroData2018.h.

{ return fXSecPp_Tot;       }

References fXSecPp_Tot.

Friends And Related Symbol Documentation

Cleaner

friend struct Cleaner

friend

Definition at line 253 of file INukeHadroData2018.h.

Member Data Documentation

fFracKA_Abs

Spline* genie::INukeHadroData2018::fFracKA_Abs

private

Definition at line 210 of file INukeHadroData2018.h.

Referenced by FracAIndep(), FracKA_Abs(), LoadCrossSections(), and ~INukeHadroData2018().

fFracKA_CEx

Spline* genie::INukeHadroData2018::fFracKA_CEx

private

Definition at line 208 of file INukeHadroData2018.h.

Referenced by FracKA_CEx(), LoadCrossSections(), and ~INukeHadroData2018().

fFracKA_Elas

Spline* genie::INukeHadroData2018::fFracKA_Elas

private

Definition at line 207 of file INukeHadroData2018.h.

Referenced by FracKA_Elas(), LoadCrossSections(), and ~INukeHadroData2018().

fFracKA_Inel

Spline* genie::INukeHadroData2018::fFracKA_Inel

private

Definition at line 209 of file INukeHadroData2018.h.

Referenced by FracAIndep(), FracKA_Inel(), LoadCrossSections(), and ~INukeHadroData2018().

fFracKA_Tot

Spline* genie::INukeHadroData2018::fFracKA_Tot

private

K+A x-section splines.

Definition at line 206 of file INukeHadroData2018.h.

Referenced by FracKA_Tot(), LoadCrossSections(), and ~INukeHadroData2018().

fFracNA_Abs

Spline* genie::INukeHadroData2018::fFracNA_Abs

private

Definition at line 204 of file INukeHadroData2018.h.

Referenced by FracAIndep(), FracNA_Abs(), LoadCrossSections(), and ~INukeHadroData2018().

fFracNA_CEx

Spline* genie::INukeHadroData2018::fFracNA_CEx

private

Definition at line 203 of file INukeHadroData2018.h.

Referenced by FracAIndep(), FracNA_CEx(), LoadCrossSections(), and ~INukeHadroData2018().

fFracNA_Cmp

Spline* genie::INukeHadroData2018::fFracNA_Cmp

private

Definition at line 218 of file INukeHadroData2018.h.

Referenced by FracAIndep(), FracNA_Cmp(), LoadCrossSections(), and ~INukeHadroData2018().

fFracNA_Elas

Spline* genie::INukeHadroData2018::fFracNA_Elas

private

Definition at line 201 of file INukeHadroData2018.h.

Referenced by FracNA_Elas().

fFracNA_Inel

Spline* genie::INukeHadroData2018::fFracNA_Inel

private

Definition at line 202 of file INukeHadroData2018.h.

Referenced by FracAIndep(), FracNA_Inel(), LoadCrossSections(), and ~INukeHadroData2018().

fFracNA_PiPro

Spline* genie::INukeHadroData2018::fFracNA_PiPro

private

Definition at line 205 of file INukeHadroData2018.h.

Referenced by FracAIndep(), FracNA_PiPro(), LoadCrossSections(), and ~INukeHadroData2018().

fFracNA_Tot

Spline* genie::INukeHadroData2018::fFracNA_Tot

private

Definition at line 200 of file INukeHadroData2018.h.

Referenced by FracNA_Tot(), LoadCrossSections(), and ~INukeHadroData2018().

fFracPA_Abs

Spline* genie::INukeHadroData2018::fFracPA_Abs

private

Definition at line 198 of file INukeHadroData2018.h.

Referenced by FracAIndep(), FracPA_Abs(), LoadCrossSections(), and ~INukeHadroData2018().

fFracPA_CEx

Spline* genie::INukeHadroData2018::fFracPA_CEx

private

Definition at line 197 of file INukeHadroData2018.h.

Referenced by FracAIndep(), FracPA_CEx(), LoadCrossSections(), and ~INukeHadroData2018().

fFracPA_Cmp

Spline* genie::INukeHadroData2018::fFracPA_Cmp

private

Definition at line 217 of file INukeHadroData2018.h.

Referenced by FracAIndep(), FracPA_Cmp(), LoadCrossSections(), and ~INukeHadroData2018().

fFracPA_Elas

Spline* genie::INukeHadroData2018::fFracPA_Elas

private

Definition at line 195 of file INukeHadroData2018.h.

Referenced by FracPA_Elas().

fFracPA_Inel

Spline* genie::INukeHadroData2018::fFracPA_Inel

private

Definition at line 196 of file INukeHadroData2018.h.

Referenced by FracAIndep(), FracPA_Inel(), LoadCrossSections(), and ~INukeHadroData2018().

fFracPA_PiPro

Spline* genie::INukeHadroData2018::fFracPA_PiPro

private

Definition at line 199 of file INukeHadroData2018.h.

Referenced by FracAIndep(), FracPA_PiPro(), LoadCrossSections(), and ~INukeHadroData2018().

fFracPA_Tot

Spline* genie::INukeHadroData2018::fFracPA_Tot

private

N+A x-section splines.

Definition at line 194 of file INukeHadroData2018.h.

Referenced by FracPA_Tot(), LoadCrossSections(), and ~INukeHadroData2018().

fhN2dXSecGamPi0N_Inelas

BLI2DNonUnifGrid* genie::INukeHadroData2018::fhN2dXSecGamPi0N_Inelas

private

Definition at line 239 of file INukeHadroData2018.h.

Referenced by hN2dXSecGamPi0N_Inelas(), LoadCrossSections(), XSec(), and ~INukeHadroData2018().

fhN2dXSecGamPi0P_Inelas

BLI2DNonUnifGrid* genie::INukeHadroData2018::fhN2dXSecGamPi0P_Inelas

private

Definition at line 238 of file INukeHadroData2018.h.

Referenced by hN2dXSecGamPi0P_Inelas(), LoadCrossSections(), XSec(), and ~INukeHadroData2018().

fhN2dXSecGamPimP_Inelas

BLI2DNonUnifGrid* genie::INukeHadroData2018::fhN2dXSecGamPimP_Inelas

private

Definition at line 241 of file INukeHadroData2018.h.

Referenced by hN2dXSecGamPimP_Inelas(), LoadCrossSections(), XSec(), and ~INukeHadroData2018().

fhN2dXSecGamPipN_Inelas

BLI2DNonUnifGrid* genie::INukeHadroData2018::fhN2dXSecGamPipN_Inelas

private

Definition at line 240 of file INukeHadroData2018.h.

Referenced by hN2dXSecGamPipN_Inelas(), LoadCrossSections(), XSec(), and ~INukeHadroData2018().

fhN2dXSecKpN_Abs

BLI2DNonUnifGrid* genie::INukeHadroData2018::fhN2dXSecKpN_Abs

private

Definition at line 235 of file INukeHadroData2018.h.

fhN2dXSecKpN_CEx

BLI2DNonUnifGrid* genie::INukeHadroData2018::fhN2dXSecKpN_CEx

private

Definition at line 234 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), and ~INukeHadroData2018().

fhN2dXSecKpN_Elas

BLI2DNonUnifGrid* genie::INukeHadroData2018::fhN2dXSecKpN_Elas

private

Definition at line 232 of file INukeHadroData2018.h.

Referenced by hN2dXSecKpN_Elas(), LoadCrossSections(), XSec(), and ~INukeHadroData2018().

fhN2dXSecKpP_Elas

BLI2DNonUnifGrid* genie::INukeHadroData2018::fhN2dXSecKpP_Elas

private

Definition at line 233 of file INukeHadroData2018.h.

Referenced by hN2dXSecKpP_Elas(), LoadCrossSections(), XSec(), and ~INukeHadroData2018().

fhN2dXSecNP_Elas

BLI2DNonUnifGrid* genie::INukeHadroData2018::fhN2dXSecNP_Elas

private

Definition at line 228 of file INukeHadroData2018.h.

Referenced by hN2dXSecNP_Elas(), LoadCrossSections(), XSec(), and ~INukeHadroData2018().

fhN2dXSecPi0N_Elas

BLI2DNonUnifGrid* genie::INukeHadroData2018::fhN2dXSecPi0N_Elas

private

Definition at line 230 of file INukeHadroData2018.h.

Referenced by hN2dXSecPi0N_Elas(), LoadCrossSections(), XSec(), and ~INukeHadroData2018().

fhN2dXSecPimN_Elas

BLI2DNonUnifGrid* genie::INukeHadroData2018::fhN2dXSecPimN_Elas

private

Definition at line 231 of file INukeHadroData2018.h.

Referenced by hN2dXSecPimN_Elas(), LoadCrossSections(), XSec(), and ~INukeHadroData2018().

fhN2dXSecPiN_Abs

BLI2DNonUnifGrid* genie::INukeHadroData2018::fhN2dXSecPiN_Abs

private

Definition at line 237 of file INukeHadroData2018.h.

Referenced by hN2dXSecPiN_Abs(), LoadCrossSections(), XSec(), and ~INukeHadroData2018().

fhN2dXSecPiN_CEx

BLI2DNonUnifGrid* genie::INukeHadroData2018::fhN2dXSecPiN_CEx

private

Definition at line 236 of file INukeHadroData2018.h.

Referenced by hN2dXSecPiN_CEx(), LoadCrossSections(), XSec(), and ~INukeHadroData2018().

fhN2dXSecPipN_Elas

BLI2DNonUnifGrid* genie::INukeHadroData2018::fhN2dXSecPipN_Elas

private

Definition at line 229 of file INukeHadroData2018.h.

Referenced by hN2dXSecPipN_Elas(), LoadCrossSections(), XSec(), and ~INukeHadroData2018().

fhN2dXSecPP_Elas

BLI2DNonUnifGrid* genie::INukeHadroData2018::fhN2dXSecPP_Elas

private

Definition at line 227 of file INukeHadroData2018.h.

Referenced by hN2dXSecPP_Elas(), LoadCrossSections(), XSec(), and ~INukeHadroData2018().

fInstance

INukeHadroData2018 * INukeHadroData2018::fInstance = 0

staticprivate

Definition at line 160 of file INukeHadroData2018.h.

Referenced by Instance(), INukeHadroData2018(), and genie::INukeHadroData2018::Cleaner::~Cleaner().

fMaxKinEnergyHA

double INukeHadroData2018::fMaxKinEnergyHA = 999.0

static

Definition at line 146 of file INukeHadroData2018.h.

Referenced by FracADep(), and FracAIndep().

fMaxKinEnergyHN

double INukeHadroData2018::fMaxKinEnergyHN = 1799.0

static

Definition at line 147 of file INukeHadroData2018.h.

Referenced by Frac(), genie::utils::intranuke2018::MeanFreePath(), and XSec().

fMinKinEnergy

double INukeHadroData2018::fMinKinEnergy = 1.0

static

Definition at line 145 of file INukeHadroData2018.h.

Referenced by Frac(), FracADep(), FracAIndep(), genie::utils::intranuke2018::MeanFreePath(), and XSec().

fXSecGamn_fs

Spline* genie::INukeHadroData2018::fXSecGamn_fs

private

Definition at line 212 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSecGamn_fs(), and ~INukeHadroData2018().

fXSecGamN_Tot

Spline* genie::INukeHadroData2018::fXSecGamN_Tot

private

Definition at line 213 of file INukeHadroData2018.h.

Referenced by Frac(), LoadCrossSections(), XSecGamN_Tot(), and ~INukeHadroData2018().

fXSecGamp_fs

Spline* genie::INukeHadroData2018::fXSecGamp_fs

private

gamma A x-section splines

Definition at line 211 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSecGamp_fs(), and ~INukeHadroData2018().

fXSecKpN_Abs

Spline* genie::INukeHadroData2018::fXSecKpN_Abs

private

Definition at line 192 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSecKpN_Abs(), and ~INukeHadroData2018().

fXSecKpn_CEx

Spline* genie::INukeHadroData2018::fXSecKpn_CEx

private

Definition at line 191 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), XSecKpn_CEx(), and ~INukeHadroData2018().

fXSecKpn_Elas

Spline* genie::INukeHadroData2018::fXSecKpn_Elas

private

K+N x-section splines.

Definition at line 189 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), XSecKpn_Elas(), and ~INukeHadroData2018().

fXSecKpN_Tot

Spline* genie::INukeHadroData2018::fXSecKpN_Tot

private

Definition at line 193 of file INukeHadroData2018.h.

Referenced by Frac(), LoadCrossSections(), XSecKpN_Tot(), and ~INukeHadroData2018().

fXSecKpp_Elas

Spline* genie::INukeHadroData2018::fXSecKpp_Elas

private

Definition at line 190 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSecKpp_Elas(), and ~INukeHadroData2018().

fXSecNn_Cmp

Spline* genie::INukeHadroData2018::fXSecNn_Cmp

private

Definition at line 216 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), and ~INukeHadroData2018().

fXSecNn_Elas

Spline* genie::INukeHadroData2018::fXSecNn_Elas

private

Definition at line 187 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), XSecNn_Elas(), and ~INukeHadroData2018().

fXSecNn_Reac

Spline* genie::INukeHadroData2018::fXSecNn_Reac

private

Definition at line 188 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), XSecNn_Reac(), and ~INukeHadroData2018().

fXSecNn_Tot

Spline* genie::INukeHadroData2018::fXSecNn_Tot

private

Definition at line 186 of file INukeHadroData2018.h.

Referenced by Frac(), LoadCrossSections(), XSecNn_Tot(), and ~INukeHadroData2018().

fXSecPi0d_Abs

Spline* genie::INukeHadroData2018::fXSecPi0d_Abs

private

Definition at line 179 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), XSecPi0p_Abs(), and ~INukeHadroData2018().

fXSecPi0n_CEx

Spline* genie::INukeHadroData2018::fXSecPi0n_CEx

private

Definition at line 172 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), XSecPi0n_CEx(), and ~INukeHadroData2018().

fXSecPi0n_Elas

Spline* genie::INukeHadroData2018::fXSecPi0n_Elas

private

Definition at line 173 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), XSecPi0n_Elas(), and ~INukeHadroData2018().

fXSecPi0n_Reac

Spline* genie::INukeHadroData2018::fXSecPi0n_Reac

private

Definition at line 174 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), XSecPi0n_Reac(), and ~INukeHadroData2018().

fXSecPi0n_Tot

Spline* genie::INukeHadroData2018::fXSecPi0n_Tot

private

pi0n hN x-section splines

Definition at line 171 of file INukeHadroData2018.h.

Referenced by Frac(), LoadCrossSections(), XSecPi0n_Tot(), and ~INukeHadroData2018().

fXSecPi0p_CEx

Spline* genie::INukeHadroData2018::fXSecPi0p_CEx

private

Definition at line 176 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), XSecPi0p_CEx(), and ~INukeHadroData2018().

fXSecPi0p_Elas

Spline* genie::INukeHadroData2018::fXSecPi0p_Elas

private

Definition at line 177 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), XSecPi0p_Elas(), and ~INukeHadroData2018().

fXSecPi0p_Reac

Spline* genie::INukeHadroData2018::fXSecPi0p_Reac

private

Definition at line 178 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), XSecPi0p_Reac(), and ~INukeHadroData2018().

fXSecPi0p_Tot

Spline* genie::INukeHadroData2018::fXSecPi0p_Tot

private

pi0p hN x-section splines

Definition at line 175 of file INukeHadroData2018.h.

Referenced by Frac(), LoadCrossSections(), XSecPi0p_Tot(), and ~INukeHadroData2018().

fXSecPipd_Abs

Spline* genie::INukeHadroData2018::fXSecPipd_Abs

private

Definition at line 170 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), XSecPipp_Abs(), and ~INukeHadroData2018().

fXSecPipn_CEx

Spline* genie::INukeHadroData2018::fXSecPipn_CEx

private

Definition at line 163 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), XSecPipn_CEx(), and ~INukeHadroData2018().

fXSecPipn_Elas

Spline* genie::INukeHadroData2018::fXSecPipn_Elas

private

Definition at line 164 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), XSecPipn_Elas(), and ~INukeHadroData2018().

fXSecPipn_Reac

Spline* genie::INukeHadroData2018::fXSecPipn_Reac

private

Definition at line 165 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), XSecPipn_Reac(), and ~INukeHadroData2018().

fXSecPipn_Tot

Spline* genie::INukeHadroData2018::fXSecPipn_Tot

private

pi+n hN x-section splines

Definition at line 162 of file INukeHadroData2018.h.

Referenced by Frac(), LoadCrossSections(), XSecPipn_Tot(), and ~INukeHadroData2018().

fXSecPipp_CEx

Spline* genie::INukeHadroData2018::fXSecPipp_CEx

private

Definition at line 167 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), XSecPipp_CEx(), and ~INukeHadroData2018().

fXSecPipp_Elas

Spline* genie::INukeHadroData2018::fXSecPipp_Elas

private

Definition at line 168 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), XSecPipp_Elas(), and ~INukeHadroData2018().

fXSecPipp_Reac

Spline* genie::INukeHadroData2018::fXSecPipp_Reac

private

Definition at line 169 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), XSecPipp_Reac(), and ~INukeHadroData2018().

fXSecPipp_Tot

Spline* genie::INukeHadroData2018::fXSecPipp_Tot

private

pi+p hN x-section splines

Definition at line 166 of file INukeHadroData2018.h.

Referenced by Frac(), LoadCrossSections(), XSecPipp_Tot(), and ~INukeHadroData2018().

fXSecPn_Cmp

Spline* genie::INukeHadroData2018::fXSecPn_Cmp

private

Definition at line 215 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), and ~INukeHadroData2018().

fXSecPn_Elas

Spline* genie::INukeHadroData2018::fXSecPn_Elas

private

Definition at line 184 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), XSecPn_Elas(), and ~INukeHadroData2018().

fXSecPn_Reac

Spline* genie::INukeHadroData2018::fXSecPn_Reac

private

Definition at line 185 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), XSecPn_Reac(), and ~INukeHadroData2018().

fXSecPn_Tot

Spline* genie::INukeHadroData2018::fXSecPn_Tot

private

Definition at line 183 of file INukeHadroData2018.h.

Referenced by Frac(), LoadCrossSections(), XSecPn_Tot(), and ~INukeHadroData2018().

fXSecPp_Cmp

Spline* genie::INukeHadroData2018::fXSecPp_Cmp

private

NN cmp (compound nucleus) fate.

Definition at line 214 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), and ~INukeHadroData2018().

fXSecPp_Elas

Spline* genie::INukeHadroData2018::fXSecPp_Elas

private

Definition at line 181 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), XSecPp_Elas(), and ~INukeHadroData2018().

fXSecPp_Reac

Spline* genie::INukeHadroData2018::fXSecPp_Reac

private

Definition at line 182 of file INukeHadroData2018.h.

Referenced by LoadCrossSections(), XSec(), XSecPp_Reac(), and ~INukeHadroData2018().

fXSecPp_Tot

Spline* genie::INukeHadroData2018::fXSecPp_Tot

private

p/nN x-section splines

Definition at line 180 of file INukeHadroData2018.h.

Referenced by Frac(), LoadCrossSections(), XSecPp_Tot(), and ~INukeHadroData2018().

TfracPipA_Abs

TGraph2D* genie::INukeHadroData2018::TfracPipA_Abs

private

Definition at line 224 of file INukeHadroData2018.h.

Referenced by FracADep(), LoadCrossSections(), and ~INukeHadroData2018().

TfracPipA_CEx

TGraph2D* genie::INukeHadroData2018::TfracPipA_CEx

private

Definition at line 221 of file INukeHadroData2018.h.

Referenced by FracADep(), LoadCrossSections(), and ~INukeHadroData2018().

TfracPipA_Inelas

TGraph2D* genie::INukeHadroData2018::TfracPipA_Inelas

private

Definition at line 223 of file INukeHadroData2018.h.

Referenced by FracADep(), LoadCrossSections(), and ~INukeHadroData2018().

TfracPipA_PiPro

TGraph2D* genie::INukeHadroData2018::TfracPipA_PiPro

private

Definition at line 225 of file INukeHadroData2018.h.

Referenced by FracADep(), LoadCrossSections(), and ~INukeHadroData2018().

---

The documentation for this class was generated from the following files:

• INukeHadroData2018.h
• INukeHadroData2018.cxx