Correction to free NN xsec in nuclear matter. More...

#include <INukeNucleonCorr.h>

Collaboration diagram for INukeNucleonCorr:

Public Member Functions
double	getAvgCorrection (const double rho, const double A, const double Ek)
	get the correction for given four-momentum and density
void	OutputFiles (int A, int Z)
double	AvgCorrection (const double rho, const int A, const int Z, const int pdg, const double Ek)
	generate kinematics fRepeat times to calculate average correction

Static Public Member Functions
static INukeNucleonCorr *	getInstance ()
	get single instance of INukeNucleonCorr; create if necessary

Private Member Functions
	INukeNucleonCorr ()
	private constructor (called only by getInstance())
	INukeNucleonCorr (const INukeNucleonCorr &)
	block copy constructor
INukeNucleonCorr &	operator= (const INukeNucleonCorr &)
	block assignment operator
double	beta (const double rho)
	potential component (Eq. 2.18)
double	lambda (const double rho)
	potential component (Eq. 2.19)
void	setFermiLevel (const double rho, const int A, const int Z)
double	fermiMomentum (const int pdg)
	return proper Fermi momentum based on nucleon PDG
double	mstar (const double rho, const double k2)
	m* calculated based on Eqs. 2.6 and 2.16
double	localFermiMom (const double rho, const int A, const int Z, const int pdg)
	calculate local Fermi momentum
TLorentzVector	generateTargetNucleon (const double mass, const double fermiMomentum)
	generate target nucleon
double	getCorrection (const double mass, const double rho, const TVector3 &k1, const TVector3 &k2, const TVector3 &k3, const TVector3 &k4)
	calculate xsec correction

Private Attributes
double	fFermiMomProton
double	fFermiMomNeutron

Static Private Attributes
static INukeNucleonCorr *	fInstance = NULL
	single instance of INukeNucleonCorr
static const unsigned int	fRepeat = 1000
	number of repetition to get average correction
static const double	fRho0 = 0.16
	equilibrium density
static const double	fAlpha1
	alpha coefficient as defined by Eq. 2.17
static const double	fAlpha2
	alpha coefficient as defined by Eq. 2.17
static const double	fBeta1 = -116.00 / fRho0 / 1000.0
	beta coefficient as defined by Eq. 2.18
static const double	fLambda0 = 3.29 / (units::fermi)
	lambda coefficient as defined by Eq. 2.19
static const double	fLambda1 = -0.373 / (units::fermi) / fRho0
	lambda coefficient as defined by Eq. 2.19
static const int	fNDensityBins
	cache binning for density
static const double	fDensityStep
	within this density step correction is assumed to be constant
static const int	fNEnergyBins1
	cache binning for energy
static const double	fMaxEnergy1
	above this energy correction is assumed to be constant
static const double	fEnergyStep1
	within this energy step correction is assumed to be constant
static const int	fNEnergyBins2
	cache binning for energy
static const double	fMaxEnergy2
	above this energy correction is assumed to be constant
static const double	fEnergyStep2
	within this energy step correction is assumed to be constant
static const int	fNEnergyBins3
	cache binning for energy
static const double	fMaxEnergy3
	above this energy correction is assumed to be constant
static const double	fEnergyStep3
	within this energy step correction is assumed to be constant

Detailed Description

Correction to free NN xsec in nuclear matter.

Author: Kyle Bachinski, Tomasz Golan

Date: 2015

Remarks: V.R. Pandharipande and S. C. Pieper, Phys. Rev. C45 (1992) 791

Definition at line 18 of file INukeNucleonCorr.h.

Constructor & Destructor Documentation

◆ INukeNucleonCorr() [1/2]

INukeNucleonCorr::INukeNucleonCorr ( )

inlineprivate

private constructor (called only by getInstance())

Definition at line 78 of file INukeNucleonCorr.h.

Referenced by getInstance(), INukeNucleonCorr(), and operator=().

◆ INukeNucleonCorr() [2/2]

INukeNucleonCorr::INukeNucleonCorr ( const INukeNucleonCorr & )

private

block copy constructor

References INukeNucleonCorr().

Member Function Documentation

◆ AvgCorrection()

double INukeNucleonCorr::AvgCorrection	(	const double	rho,
		const int	A,
		const int	Z,
		const int	pdg,
		const double	Ek )

generate kinematics fRepeat times to calculate average correction

Definition at line 147 of file INukeNucleonCorr.cxx.

{
 
  RandomGen * rnd = RandomGen::Instance();
 
  setFermiLevel (rho, A, Z); // set Fermi momenta for protons and neutrons
 
  const double mass   = PDGLibrary::Instance()->Find(pdg)->Mass(); // mass of incoming nucleon
  const double energy = Ek + mass;
 
  TLorentzVector p (0.0, 0.0, sqrt (energy * energy - mass * mass), energy); // incoming particle 4-momentum
  GHepParticle incomingParticle (pdg, kIStUndefined, -1,-1,-1,-1, p, TLorentzVector ()); // for IntBounce
 
  double corrPauliBlocking = 0.0; // correction coming from Pauli blocking
  double     corrPotential = 0.0; // correction coming from potential
 
  for (unsigned int i = 0; i < fRepeat; i++) // generate kinematics fRepeat times to get avg corrections
  {
    // get proton vs neutron randomly based on Z/A
    const int targetPdg = rnd->RndGen().Rndm() < (double) Z / A ? kPdgProton : kPdgNeutron;
 
    const double targetMass = PDGLibrary::Instance()->Find(targetPdg)->Mass(); // set nucleon mass
 
    const TLorentzVector target = generateTargetNucleon (targetMass, fermiMomentum (targetPdg)); // generate target nucl
 
    TLorentzVector outNucl1, outNucl2, RemnP4; // final 4-momenta
 
    // random scattering angle
    double C3CM = INukeHadroData2018::Instance()->IntBounce (&incomingParticle, targetPdg, pdg, kIHNFtElas);
 
    // generate kinematics
    utils::intranuke2018::TwoBodyKinematics (mass, targetMass, p, target, outNucl1, outNucl2, C3CM, RemnP4);
 
    // update Pauli blocking correction
    corrPauliBlocking += (outNucl1.Vect().Mag() > fermiMomentum (pdg) and outNucl2.Vect().Mag() > fermiMomentum (targetPdg));
 
    // update potential-based correction
    corrPotential += getCorrection (mass, rho, p.Vect(), target.Vect(), outNucl1.Vect(), outNucl2.Vect());
  }
 
  corrPauliBlocking /= fRepeat;
      corrPotential /= fRepeat;
 
      return corrPotential * corrPauliBlocking;
 
}

References fermiMomentum(), genie::PDGLibrary::Find(), fRepeat, generateTargetNucleon(), getCorrection(), genie::INukeHadroData2018::Instance(), genie::PDGLibrary::Instance(), genie::RandomGen::Instance(), genie::INukeHadroData2018::IntBounce(), genie::kIHNFtElas, genie::kIStUndefined, genie::kPdgNeutron, genie::kPdgProton, genie::RandomGen::RndGen(), setFermiLevel(), and genie::utils::intranuke2018::TwoBodyKinematics().

Referenced by OutputFiles().

◆ beta()

double INukeNucleonCorr::beta ( const double rho )

inlineprivate

potential component (Eq. 2.18)

Definition at line 84 of file INukeNucleonCorr.h.

References fBeta1.

Referenced by mstar().

◆ fermiMomentum()

double INukeNucleonCorr::fermiMomentum ( const int pdg )

inlineprivate

return proper Fermi momentum based on nucleon PDG

Definition at line 94 of file INukeNucleonCorr.h.

94{return pdg == genie::kPdgProton ? fFermiMomProton : fFermiMomNeutron;}

INukeNucleonCorr::fFermiMomProton

double fFermiMomProton

Definition INukeNucleonCorr.h:73

INukeNucleonCorr::fFermiMomNeutron

double fFermiMomNeutron

Definition INukeNucleonCorr.h:74

References fFermiMomNeutron, fFermiMomProton, and genie::kPdgProton.

Referenced by AvgCorrection().

◆ generateTargetNucleon()

TLorentzVector INukeNucleonCorr::generateTargetNucleon	(	const double	mass,
		const double	fermiMomentum )

private

generate target nucleon

generate random momentum direction and return 4-momentum of target nucleon

Definition at line 117 of file INukeNucleonCorr.cxx.

{
  RandomGen * rnd = RandomGen::Instance();
 
  // get random momentum direction
  const double costheta = 2.0 * rnd->RndGen().Rndm() - 1.0;  // random cos (theta)
  const double sintheta = sqrt (1.0 - costheta * costheta);  // sin (theta)
  const double      phi = 2.0 * kPi * rnd->RndGen().Rndm(); // random phi
 
  // set nucleon 4-momentum
  const double p = rnd->RndGen().Rndm() * fermiMom; // random nucleon momentum up to Fermi level
 
  const TVector3   p3 = TVector3 (p * sintheta * cos (phi), p * sintheta * sin (phi), p * costheta); // 3-momentum
  const double energy = sqrt (p3.Mag2() + mass * mass); // energy
 
  return TLorentzVector (p3, energy);
}

References genie::RandomGen::Instance(), genie::constants::kPi, and genie::RandomGen::RndGen().

Referenced by AvgCorrection().

◆ getAvgCorrection()

double INukeNucleonCorr::getAvgCorrection	(	const double	rho,
		const double	A,
		const double	Ek )

get the correction for given four-momentum and density

Definition at line 237 of file INukeNucleonCorr.cxx.

{
  //Read in energy and density to determine the row and column of the correction table - adjust for variable binning - throws away some of the accuracy
   int Column = round(rho*100);
   if(rho<.01) Column = 1;
   if (Column>=NColumns) Column = NColumns-1;
   int Row = 0;
   if(ke<=.002) Row = 1;
   if(ke>.002&&ke<=.1) Row = round(ke*1000.);
   if(ke>.1&&ke<=.5) Row = round(.1*1000.+(ke-.1)*200);
   if(ke>.5&&ke<=1) Row = round(.1*1000.+(.5-.1)*200+(ke-.5)*40);
   if(ke>1) Row = NRows-1;
   //LOG ("INukeNucleonCorr",pNOTICE)
   //  << "row, column = " << Row << "   " << Column;
  //If the table of correction values has already been created
  // return a value. Else, interpolate the needed correction table//
//  static double cache[NRows][NColumns] = {{-1}};
  static bool ReadFile = false;
  if( ReadFile == true ) {
   int Npoints = 6;
   TGraph * Interp = new TGraph(Npoints);
   //LOG("INukeNucleonCorr",pNOTICE)
   //  << HeliumValues[Row][Column];
   Interp->SetPoint(0,4,HeliumValues[Row][Column]);
   Interp->SetPoint(1,12,CarbonValues[Row][Column]);
   Interp->SetPoint(2,40,CalciumValues[Row][Column]);
   Interp->SetPoint(3,56,IronValues[Row][Column]);
   Interp->SetPoint(4,120,TinValues[Row][Column]);
   Interp->SetPoint(5,238,UraniumValues[Row][Column]);
 
   //   Interpolated[e][r] = Interp->Eval(A);
        //      LOG("INukeNucleonCorr",pNOTICE)
        //        << "e,r,value= " << e << "   " << r << "   " << Interpolated[e][r];
   double returnval = Interp->Eval(A);
   delete Interp;
   LOG("INukeNucleonCorr",pINFO)
      << "Nucleon Corr interpolated correction factor = "
      << returnval  //cache[Row][Column]
      << " for rho, KE, A= "<<  rho << "  " << ke << "   " << A;
    //    return cache[Row][Column];
   return returnval;
  } else {
    //Reading in correction files//
    //    string dir = genie_dir + string("/data/evgen/nncorr/");
 
    read_file(dir+"NNCorrection_2_4.txt");
    HeliumValues = infile_values;
    infile_values = clear;
 
    read_file(dir+"NNCorrection_6_12.txt");
    CarbonValues = infile_values;
    infile_values = clear;
 
    read_file(dir+"NNCorrection_20_40.txt");
    CalciumValues = infile_values;
    infile_values = clear;
 
    read_file(dir+"NNCorrection_26_56.txt");
    IronValues = infile_values;
    infile_values = clear;
 
    read_file(dir+"NNCorrection_50_120.txt");
    TinValues = infile_values;
    infile_values = clear;
 
    read_file(dir+"NNCorrection_92_238.txt");
    UraniumValues = infile_values;
    infile_values = clear;
 
    LOG("INukeNucleonCorr",pNOTICE)
      << "Nucleon Corr interpolation files read in successfully";
    //get interpolated value for first event.
    int Npoints = 6;
    TGraph * Interp = new TGraph(Npoints);
    //LOG("INukeNucleonCorr",pNOTICE)
    //  << HeliumValues[Row][Column];
    Interp->SetPoint(0,4,HeliumValues[Row][Column]);
    Interp->SetPoint(1,12,CarbonValues[Row][Column]);
    Interp->SetPoint(2,40,CalciumValues[Row][Column]);
    Interp->SetPoint(3,56,IronValues[Row][Column]);
    Interp->SetPoint(4,120,TinValues[Row][Column]);
    Interp->SetPoint(5,238,UraniumValues[Row][Column]);
 
    //  LOG("INukeNucleonCorr",pNOTICE)
    //    << "Row,Column,value= " << Row << "   " << Column << "   " << Interp->Eval(A);
    double returnval = Interp->Eval(A);
    delete Interp;
    ReadFile = true;
    LOG("INukeNucleonCorr",pINFO)
      << "Nucleon Corr interpolated correction factor = "
      << returnval
      << " for rho, KE, A= "<<  rho << "  " << ke << "   " << A;
    return returnval;
  }
}

References CalciumValues, CarbonValues, clear, dir, HeliumValues, infile_values, IronValues, LOG, NColumns, NRows, pINFO, pNOTICE, read_file(), TinValues, and UraniumValues.

◆ getCorrection()

double INukeNucleonCorr::getCorrection	(	const double	mass,
		const double	rho,
		const TVector3 &	k1,
		const TVector3 &	k2,
		const TVector3 &	k3,
		const TVector3 &	k4 )

private

calculate xsec correction

calculate correction given by eq. 2.9

Definition at line 136 of file INukeNucleonCorr.cxx.

{
  const double num = (k1 - k2).Mag() * mstar (rho, (k3.Mag2() + k4.Mag2()) / 2.0) / mass / mass;
  const double den = (k1 * (1.0 / mstar (rho, k1.Mag2())) - k2 * (1.0 / mstar (rho, k2.Mag2()))).Mag();
 
  return num / den;
}

References mstar().

Referenced by AvgCorrection().

◆ getInstance()

INukeNucleonCorr * INukeNucleonCorr::getInstance ( )

inlinestatic

get single instance of INukeNucleonCorr; create if necessary

Definition at line 23 of file INukeNucleonCorr.h.

23{return fInstance ? fInstance : (fInstance = new INukeNucleonCorr);}

INukeNucleonCorr::INukeNucleonCorr

INukeNucleonCorr()

private constructor (called only by getInstance())

Definition INukeNucleonCorr.h:78

INukeNucleonCorr::fInstance

static INukeNucleonCorr * fInstance

single instance of INukeNucleonCorr

Definition INukeNucleonCorr.h:33

References fInstance, and INukeNucleonCorr().

Referenced by genie::utils::intranuke2018::MeanFreePath(), genie::utils::intranuke2025::MeanFreePath(), and OutputFiles().

◆ lambda()

double INukeNucleonCorr::lambda ( const double rho )

inlineprivate

potential component (Eq. 2.19)

Definition at line 85 of file INukeNucleonCorr.h.

References fLambda0, and fLambda1.

Referenced by mstar().

◆ localFermiMom()

double INukeNucleonCorr::localFermiMom	(	const double	rho,
		const int	A,
		const int	Z,
		const int	pdg )

private

calculate local Fermi momentum

$k_F = (\frac{3}{2}\pi^2\rho)^{1/3}$

Definition at line 101 of file INukeNucleonCorr.cxx.

{
  static double factor = 3.0 * kPi * kPi / 2.0;
 
  return pdg == kPdgProton ? pow (factor * rho * Z / A, 1.0 / 3.0) / (units::fermi) :
                             pow (factor * rho * (A - Z) / A, 1.0 / 3.0) / (units::fermi);
}

References genie::units::fermi, genie::kPdgProton, and genie::constants::kPi.

Referenced by setFermiLevel().

◆ mstar()

double INukeNucleonCorr::mstar	(	const double	rho,
		const double	k2 )

private

m* calculated based on Eqs. 2.6 and 2.16

$m^* (k,\rho) = m \frac{(\Lambda^2 + k^2)^2}{\Lambda^2 + k^2)^2 - 2\Lambda^2\beta m}$

Definition at line 83 of file INukeNucleonCorr.cxx.

{
  // density [fm^-3], momentum square [GeV^2]
 
  static const double m = (PDGLibrary::Instance()->Find(kPdgProton)->Mass() +
                           PDGLibrary::Instance()->Find(kPdgNeutron)->Mass()) / 2.0;
 
  const double L = lambda (rho); // potential coefficient lambda
  const double B =   beta (rho); // potential coefficient beta
 
  const double L2 = L * L; // lambda^2 used twice in equation
 
  const double num = (L2 + k2) * (L2 + k2); // numerator
 
  return m * num / (num - 2.0 * L2 * B * m);
}

References beta(), genie::PDGLibrary::Find(), genie::PDGLibrary::Instance(), genie::kPdgNeutron, genie::kPdgProton, and lambda().

Referenced by getCorrection().

◆ operator=()

INukeNucleonCorr & INukeNucleonCorr::operator= ( const INukeNucleonCorr & )

private

block assignment operator

References INukeNucleonCorr().

◆ OutputFiles()

void INukeNucleonCorr::OutputFiles	(	int	A,
		int	Z )

Definition at line 334 of file INukeNucleonCorr.cxx.

{
  string outputdir = genie_dir + string("/data/evgen/nncorr/");
  double pdgc;
  string file;
  string header;
 
 
    file = outputdir+"NNCorrection.txt";
    header = "##Correction values for protons (density(horizontal) and energy(vertical))";
    pdgc = 2212;
 
 
  double output[1002][18];
  //label densities (columns)//
  for(int r = 1; r < 18; r++)
    { double rho = (r-1)*0.01;
      output[0][r] = rho;}
 
  //label energies (rows) //
  for(int e = 1; e < 1002; e++)
    { double energy = (e-1)*0.001;
      output[e][0] = energy;}
 
  //loop over each energy and density to get corrections and build the correction table//
  for(int e = 1; e < 1002; e++){
    for(int r = 1; r < 18; r++){
      double energy  = (e-1)*0.001;
      double density = (r-1)*0.01;
      double correction = INukeNucleonCorr::getInstance()-> AvgCorrection (density, A, Z, pdgc, energy);
      output[e][r] = correction;
    }
  }
  //output the new correction table //
  ofstream outfile;
  outfile.open((char*)file.c_str(), ios::trunc);
  if (outfile.is_open()) {
 
    outfile << header << endl;
    outfile << left << setw(12) << "##";
    for (int e = 0; e < 1002; e++) {
      for (int r = 0; r < 18; r++) {
        if((e == 0) && (r == 0)) {}
        else{outfile << left << setw(12) << output[e][r];}
      }
      outfile << endl;
 
    }
  }
  outfile.close();
}

References AvgCorrection(), e, genie_dir(), and getInstance().

◆ setFermiLevel()

void INukeNucleonCorr::setFermiLevel	(	const double	rho,
		const int	A,
		const int	Z )

inlineprivate

Parameters

Z	set up Fermi momenta

Definition at line 87 of file INukeNucleonCorr.h.

    {
      fFermiMomProton  = localFermiMom (rho, A, Z, genie::kPdgProton);  // local Fermi momentum for protons
      fFermiMomNeutron = localFermiMom (rho, A, Z, genie::kPdgNeutron); // local Fermi momentum for neutrons
    }

References fFermiMomNeutron, fFermiMomProton, genie::kPdgNeutron, genie::kPdgProton, and localFermiMom().

Referenced by AvgCorrection().

Member Data Documentation

◆ fAlpha1

const double INukeNucleonCorr::fAlpha1

staticprivate

alpha coefficient as defined by Eq. 2.17

Definition at line 43 of file INukeNucleonCorr.h.

◆ fAlpha2

const double INukeNucleonCorr::fAlpha2

staticprivate

alpha coefficient as defined by Eq. 2.17

Definition at line 44 of file INukeNucleonCorr.h.

◆ fBeta1

const double INukeNucleonCorr::fBeta1 = -116.00 / fRho0 / 1000.0

staticprivate

beta coefficient as defined by Eq. 2.18

Definition at line 45 of file INukeNucleonCorr.h.

Referenced by beta().

◆ fDensityStep

const double INukeNucleonCorr::fDensityStep

staticprivate

within this density step correction is assumed to be constant

Definition at line 53 of file INukeNucleonCorr.h.

◆ fEnergyStep1

const double INukeNucleonCorr::fEnergyStep1

staticprivate

within this energy step correction is assumed to be constant

Definition at line 61 of file INukeNucleonCorr.h.

◆ fEnergyStep2

const double INukeNucleonCorr::fEnergyStep2

staticprivate

within this energy step correction is assumed to be constant

Definition at line 64 of file INukeNucleonCorr.h.

◆ fEnergyStep3

const double INukeNucleonCorr::fEnergyStep3

staticprivate

within this energy step correction is assumed to be constant

Definition at line 67 of file INukeNucleonCorr.h.

◆ fFermiMomNeutron

double INukeNucleonCorr::fFermiMomNeutron

private

Definition at line 74 of file INukeNucleonCorr.h.

Referenced by fermiMomentum(), and setFermiLevel().

◆ fFermiMomProton

double INukeNucleonCorr::fFermiMomProton

private

Definition at line 73 of file INukeNucleonCorr.h.

Referenced by fermiMomentum(), and setFermiLevel().

◆ fInstance

INukeNucleonCorr * INukeNucleonCorr::fInstance = NULL

staticprivate

single instance of INukeNucleonCorr

Definition at line 33 of file INukeNucleonCorr.h.

Referenced by getInstance().

◆ fLambda0

const double INukeNucleonCorr::fLambda0 = 3.29 / (units::fermi)

staticprivate

lambda coefficient as defined by Eq. 2.19

Definition at line 46 of file INukeNucleonCorr.h.

Referenced by lambda().

◆ fLambda1

const double INukeNucleonCorr::fLambda1 = -0.373 / (units::fermi) / fRho0

staticprivate

lambda coefficient as defined by Eq. 2.19

Definition at line 47 of file INukeNucleonCorr.h.

Referenced by lambda().

◆ fMaxEnergy1

const double INukeNucleonCorr::fMaxEnergy1

staticprivate

above this energy correction is assumed to be constant

Definition at line 60 of file INukeNucleonCorr.h.

◆ fMaxEnergy2

const double INukeNucleonCorr::fMaxEnergy2

staticprivate

above this energy correction is assumed to be constant

Definition at line 63 of file INukeNucleonCorr.h.

◆ fMaxEnergy3

const double INukeNucleonCorr::fMaxEnergy3

staticprivate

above this energy correction is assumed to be constant

Definition at line 66 of file INukeNucleonCorr.h.

◆ fNDensityBins

const int INukeNucleonCorr::fNDensityBins

staticprivate

cache binning for density

Definition at line 49 of file INukeNucleonCorr.h.

◆ fNEnergyBins1

const int INukeNucleonCorr::fNEnergyBins1

staticprivate

cache binning for energy

Definition at line 59 of file INukeNucleonCorr.h.

◆ fNEnergyBins2

const int INukeNucleonCorr::fNEnergyBins2

staticprivate

cache binning for energy

Definition at line 62 of file INukeNucleonCorr.h.

◆ fNEnergyBins3

const int INukeNucleonCorr::fNEnergyBins3

staticprivate

cache binning for energy

Definition at line 65 of file INukeNucleonCorr.h.

◆ fRepeat

const unsigned int INukeNucleonCorr::fRepeat = 1000

staticprivate

number of repetition to get average correction

Definition at line 37 of file INukeNucleonCorr.h.

Referenced by AvgCorrection().

◆ fRho0

const double INukeNucleonCorr::fRho0 = 0.16

staticprivate

equilibrium density

Definition at line 41 of file INukeNucleonCorr.h.

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

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ INukeNucleonCorr() [1/2]

◆ INukeNucleonCorr() [2/2]

Member Function Documentation

◆ AvgCorrection()

◆ beta()

◆ fermiMomentum()

◆ generateTargetNucleon()

◆ getAvgCorrection()

◆ getCorrection()

◆ getInstance()

◆ lambda()

◆ localFermiMom()

◆ mstar()

◆ operator=()

◆ OutputFiles()

◆ setFermiLevel()

Member Data Documentation

◆ fAlpha1

◆ fAlpha2

◆ fBeta1

◆ fDensityStep

◆ fEnergyStep1

◆ fEnergyStep2

◆ fEnergyStep3

◆ fFermiMomNeutron

◆ fFermiMomProton

◆ fInstance

◆ fLambda0

◆ fLambda1

◆ fMaxEnergy1

◆ fMaxEnergy2

◆ fMaxEnergy3

◆ fNDensityBins

◆ fNEnergyBins1

◆ fNEnergyBins2

◆ fNEnergyBins3

◆ fRepeat

◆ fRho0