35 TNtuple * nt =
new TNtuple(
"nt",
"",
"Z:A:nucl:Q2:Rdef:Rm10p:Rp10p");
46 const double logQ2min = TMath::Log(Q2min);
47 const double logQ2max = TMath::Log(Q2max);
48 const double dlogQ2 = (logQ2max-logQ2min)/(N-1);
55 for(
int itgt=0; itgt<ntgt; itgt++) {
56 for(
int inuc=0; inuc<nnuc; inuc++) {
64 <<
"nuclear target = " << target[itgt] <<
", nucleon = " << nucleon[inuc];
66 for(
int i=0; i<
N; i++) {
67 double Q2 = TMath::Exp(logQ2min + i*dlogQ2);
73 <<
"Q2 = " << Q2 <<
" GeV, Rdef = " << Rdef
74 <<
", R(kF->0.9kF) = " << Rm10p <<
", R(kF->1.1*kF) = " << Rp10p;
76 nt->Fill(Z,A,nucleon[inuc],Q2,Rdef,Rm10p,Rp10p);
82 TFile
f(
"./fg_pauli_suppression_factors.root",
"recreate");
static const double kNucleonMass
double Q2(const Interaction *const i)
int IonPdgCodeToA(int pdgc)
A table of Fermi momentum constants.
double RQEFG_generic(double q2, double Mn, double kFi, double kFf, double pmax)
Singleton class to load & serve tables of Fermi momentum constants.
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
const FermiMomentumTable * GetTable(string name)
static constexpr double GeV
static constexpr double GeV2
int IonPdgCodeToZ(int pdgc)
double FindClosestKF(int target_pdgc, int nucleon_pdgc) const
const int kPdgTgtDeuterium