genie::SmithMonizUtils Class Reference

Contains auxiliary functions for Smith-Moniz model.
. More...

#include <SmithMonizUtils.h>

Inheritance diagram for genie::SmithMonizUtils:

Classes
class	Func1D

Public Member Functions
	SmithMonizUtils ()
	SmithMonizUtils (string config)
virtual	~SmithMonizUtils ()
void	SetInteraction (const Interaction *i)
double	GetBindingEnergy (void) const
double	GetFermiMomentum (void) const
double	GetTheta_k (double v, double qv) const
double	GetTheta_p (double pv, double v, double qv, double &E_p) const
double	E_nu_thr_SM (void) const
Range1D_t	Q2QES_SM_lim (void) const
Range1D_t	vQES_SM_lim (double Q2) const
Range1D_t	kFQES_SM_lim (double nu, double Q2) const
double	PhaseSpaceVolume (KinePhaseSpace_t ps) const
void	Configure (const Registry &config)
void	Configure (string config)
Public Member Functions inherited from genie::Algorithm
virtual	~Algorithm ()
virtual void	FindConfig (void)
virtual const Registry &	GetConfig (void) const
Registry *	GetOwnedConfig (void)
virtual const AlgId &	Id (void) const
	Get algorithm ID. More...
virtual AlgStatus_t	GetStatus (void) const
	Get algorithm status. More...
virtual bool	AllowReconfig (void) const
virtual AlgCmp_t	Compare (const Algorithm *alg) const
	Compare with input algorithm. More...
virtual void	SetId (const AlgId &id)
	Set algorithm ID. More...
virtual void	SetId (string name, string config)
const Algorithm *	SubAlg (const RgKey &registry_key) const
void	AdoptConfig (void)
void	AdoptSubstructure (void)
virtual void	Print (ostream &stream) const
	Print algorithm info. More...

Static Public Member Functions
static double	rho (double P_Fermi, double T_Fermi, double p)
Static Public Member Functions inherited from genie::Algorithm
static string	BuildParamVectKey (const std::string &comm_name, unsigned int i)
static string	BuildParamVectSizeKey (const std::string &comm_name)

Private Member Functions
void	LoadConfig (void)
double	QEL_EnuMin_SM (double E_nu) const
double	Q2lim1_SM (double Q2) const
double	Q2lim2_SM (double Q2) const
double	LambdaFUNCTION (double a, double b, double c) const
void	DMINFC (Func1D< SmithMonizUtils > &F, double A, double B, double EPS, double DELTA, double &X, double &Y, bool &LLM) const
double	vQES_SM_low_bound (double Q2) const
double	vQES_SM_upper_bound (double Q2) const

Private Attributes
map< int, double >	fNucRmvE
string	fKFTable
bool	fUseParametrization
const Interaction *	fInteraction
double	E_nu
	Neutrino energy (GeV) More...
double	m_lep
	Mass of final charged lepton (GeV) More...
double	mm_lep
	Squared mass of final charged lepton (GeV) More...
double	m_ini
	Mass of initial hadron or hadron system (GeV) More...
double	mm_ini
	Sqared mass of initial hadron or hadron system (GeV) More...
double	m_fin
	Mass of final hadron or hadron system (GeV) More...
double	mm_fin
	Squared mass of final hadron or hadron system (GeV) More...
double	m_tar
	Mass of target nucleus (GeV) More...
double	mm_tar
	Squared mass of target nucleus (GeV) More...
double	m_rnu
	Mass of residual nucleus (GeV) More...
double	mm_rnu
	Squared mass of residual nucleus (GeV) More...
double	P_Fermi
	Maximum value of Fermi momentum of target nucleon (GeV) More...
double	E_BIN
	Binding energy (GeV) More...
double	Enu_in
	Running neutrino energy (GeV) More...

Additional Inherited Members
Protected Member Functions inherited from genie::Algorithm
	Algorithm ()
	Algorithm (string name)
	Algorithm (string name, string config)
void	Initialize (void)
void	DeleteConfig (void)
void	DeleteSubstructure (void)
Registry *	ExtractLocalConfig (const Registry &in) const
Registry *	ExtractLowerConfig (const Registry &in, const string &alg_key) const
	Split an incoming configuration Registry into a block valid for the sub-algo identified by alg_key. More...
template<class T >
bool	GetParam (const RgKey &name, T &p, bool is_top_call=true) const
template<class T >
bool	GetParamDef (const RgKey &name, T &p, const T &def) const
template<class T >
int	GetParamVect (const std::string &comm_name, std::vector< T > &v, bool is_top_call=true) const
	Handle to load vectors of parameters. More...
int	GetParamVectKeys (const std::string &comm_name, std::vector< RgKey > &k, bool is_top_call=true) const
int	AddTopRegistry (Registry *rp, bool owns=true)
	add registry with top priority, also update ownership More...
int	AddLowRegistry (Registry *rp, bool owns=true)
	add registry with lowest priority, also update ownership More...
int	MergeTopRegistry (const Registry &r)
int	AddTopRegisties (const vector< Registry * > &rs, bool owns=false)
	Add registries with top priority, also udated Ownerships. More...
Protected Attributes inherited from genie::Algorithm
bool	fAllowReconfig
bool	fOwnsSubstruc
	true if it owns its substructure (sub-algs,...) More...
AlgId	fID
	algorithm name and configuration set More...
vector< Registry * >	fConfVect
vector< bool >	fOwnerships
	ownership for every registry in fConfVect More...
AlgStatus_t	fStatus
	algorithm execution status More...
AlgMap *	fOwnedSubAlgMp
	local pool for owned sub-algs (taken out of the factory pool) More...

Detailed Description

Contains auxiliary functions for Smith-Moniz model.
.

[1] R.A.Smith and E.J.Moniz, Nuclear Physics B43, (1972) 605-622
[2] K.S. Kuzmin, V.V. Lyubushkin, V.A.Naumov Eur. Phys. J. C54, (2008) 517-538

Author

Igor Kakorin kakor.nosp@m.in@j.nosp@m.inr.r.nosp@m.u Joint Institute for Nuclear Research
adapted from fortran code provided by:
Konstantin Kuzmin kkuzm.nosp@m.in@t.nosp@m.heor..nosp@m.jinr.nosp@m..ru Joint Institute for Nuclear Research
Vladimir Lyubushkin Joint Institute for Nuclear Research
Vadim Naumov vnaum.nosp@m.ov@t.nosp@m.heor..nosp@m.jinr.nosp@m..ru Joint Institute for Nuclear Research
based on code of:
Costas Andreopoulos const.nosp@m.anti.nosp@m.nos.a.nosp@m.ndre.nosp@m.opoul.nosp@m.os@c.nosp@m.ern.c.nosp@m.h University of Liverpool & STFC Rutherford Appleton Laboratory

May 05, 2017

Copyright (c) 2003-2020, The GENIE Collaboration For the full text of the license visit http://copyright.genie-mc.org

Definition at line 51 of file SmithMonizUtils.h.

Constructor & Destructor Documentation

SmithMonizUtils::SmithMonizUtils

(

)

Definition at line 49 of file SmithMonizUtils.cxx.

                                 :
Algorithm("genie::SmithMonizUtils")
{

}

SmithMonizUtils::SmithMonizUtils

(

string

config

)

Definition at line 55 of file SmithMonizUtils.cxx.

                                              :
Algorithm("genie::SmithMonizUtils", config)
{

}

SmithMonizUtils::~SmithMonizUtils ( )

virtual

Definition at line 61 of file SmithMonizUtils.cxx.

{

}

Member Function Documentation

void SmithMonizUtils::Configure ( const Registry &  config )

virtual

methods overloading the Algorithm() interface implementation to build the fragmentation function from configuration data

Reimplemented from genie::Algorithm.

Definition at line 66 of file SmithMonizUtils.cxx.

{
  Algorithm::Configure(config);
  this->LoadConfig();
}

void SmithMonizUtils::Configure ( string  config )

virtual

Configure the algorithm from the AlgoConfigPool based on param_set string given in input An algorithm contains a vector of registries coming from different xml configuration files, which are loaded according a very precise prioriy This methods will load a number registries in order of priority: 1) "Tunable" parameter set from CommonParametes. This is loaded with the highest prioriry and it is designed to be used for tuning procedure Usage not expected from the user. 2) For every string defined in "CommonParame" the corresponding parameter set will be loaded from CommonParameter.xml 3) parameter set specified by the config string and defined in the xml file of the algorithm 4) if config is not "Default" also the Default parameter set from the same xml file will be loaded Effectively this avoids the repetion of a parameter when it is not changed in the requested configuration

Reimplemented from genie::Algorithm.

Definition at line 72 of file SmithMonizUtils.cxx.

{
  Algorithm::Configure(config);
  this->LoadConfig();
}

void SmithMonizUtils::DMINFC ( Func1D< SmithMonizUtils > &  F, double  A, double  B, double  EPS, double  DELTA, double &  X, double &  Y, bool &  LLM  ) const

private

Definition at line 418 of file SmithMonizUtils.cxx.

{
  const double W5 = TMath::Sqrt(5);
  const double HV = (3-W5)/2;
  const double HW = (W5-1)/2;
  const double R1 = 1.0;
  const double HF = R1/2;

  int N = -1;
  if(A!=B) N = TMath::Nint(2.08*TMath::Log(TMath::Abs((A-B)/EPS)));
  double C = A;
  double D = B;
  if(A>B)
  {
     C = B;
     D = A;
  }
  bool LLT = true;
  bool LGE = true;
  double V, FV, W, FW, H;
  while(true)
  {
     H = D-C;
     if(N<0)
     {
       X = HF*(C+D);
       Y = F(X);
       LLM = TMath::Abs(X-A)>DELTA && TMath::Abs(X-B)>DELTA;
       return;
     }
     if(LLT)
     {
       V = C+HV*H;
       FV = F(V);
     }
     if(LGE)
     {
        W = C+HW*H;
        FW = F(W);
     }
     if(FV<FW)
     {
        LLT = true;
        LGE = false;
        D = W;
        W = V;
        FW = FV;
     }
     else
     {
        LLT = false;
        LGE = true;
        C = V;
        V = W;
        FV = FW;
     }
     N = N-1;
  }
}

double SmithMonizUtils::E_nu_thr_SM

(

void

)

const

Definition at line 209 of file SmithMonizUtils.cxx.

{

  Func1D<SmithMonizUtils> QEL_EnuMin_SM_(*this, &SmithMonizUtils::QEL_EnuMin_SM);


  const int MFC = 10000;      //  Maximum of function call
  const double EPSABS = 0;
  const double EPSREL = 1.0e-08;
  double E_min = ((m_lep + m_rnu + m_fin)*(m_lep + m_rnu + m_fin) - mm_tar)/(2*m_tar);
  double Enu_2 = 5.0e+00;
  double Enu_rf;
  if (QEL_EnuMin_SM(E_min) > 0)
  {
    // C++ analog of fortran function Enu_rf= DZEROX(E_min,Enu_2,EPS,MFC,QEL_EnuMin_SM,1)
    ROOT::Math::RootFinder rfgb(ROOT::Math::RootFinder::kGSL_BRENT);
    rfgb.Solve(QEL_EnuMin_SM_, E_min, Enu_2, MFC, EPSABS, EPSREL); //convergence is reached using  tolerance = 2 *( epsrel * abs(x) + epsabs)
    Enu_rf = rfgb.Root();
  }
  else
  {
    Enu_rf = -1.0e+01;
  }
  E_min = TMath::Max(E_min,Enu_rf);
  if(E_min < 0)
  {
    E_min = 0;
    LOG("SmithMoniz", pDEBUG) << "E_min = " << E_min;
  }
  return E_min;

}

double SmithMonizUtils::GetBindingEnergy

(

void

)

const

Definition at line 493 of file SmithMonizUtils.cxx.

{
  return E_BIN;
}

double SmithMonizUtils::GetFermiMomentum

(

void

)

const

Definition at line 498 of file SmithMonizUtils.cxx.

{
  return P_Fermi;
}

double SmithMonizUtils::GetTheta_k	(	double	v,
		double	qv
	)		const

Definition at line 503 of file SmithMonizUtils.cxx.

{
  return TMath::ACos((v + (mm_lep-v*v+qv*qv)/2/E_nu)/qv);
}

double SmithMonizUtils::GetTheta_p	(	double	pv,
		double	v,
		double	qv,
		double &	E_p
	)		const

Definition at line 508 of file SmithMonizUtils.cxx.

{
  E_p = TMath::Sqrt(mm_ini+pv*pv)-E_BIN;
  if (pv != 0)
    return TMath::ACos(((v-E_BIN)*(2*E_p+v+E_BIN)-qv*qv+mm_ini-mm_fin)/(2*pv*qv));
  else
    return 0;
}

Range1D_t SmithMonizUtils::kFQES_SM_lim	(	double	nu,
		double	Q2
	)		const

Definition at line 395 of file SmithMonizUtils.cxx.

{
  double qv = TMath::Sqrt(nu*nu+Q2);
  double c_f = (nu-E_BIN)/qv;
  double d_f = (E_BIN*E_BIN-2*nu*E_BIN-Q2+mm_ini-mm_fin)/(2*qv*m_ini);
  double Ef_min= m_ini*(c_f*d_f+TMath::Sqrt(1.0-c_f*c_f+d_f*d_f))/(1.0-c_f*c_f);
  double kF_min= TMath::Sqrt(TMath::Max(Ef_min*Ef_min-mm_ini,0.0));
  double kF_max= P_Fermi;
  Range1D_t R;
  if (kF_min<=kF_max)
    R = Range1D_t(kF_min,kF_max);
  else
    R = Range1D_t(0.5*(kF_min+kF_max),0.5*(kF_min+kF_max));
  return R;

}

double SmithMonizUtils::LambdaFUNCTION ( double  a, double  b, double  c  ) const

private

Definition at line 412 of file SmithMonizUtils.cxx.

{
  return a*a+b*b+c*c-2*(a*b+b*c+a*c);
}

void SmithMonizUtils::LoadConfig ( void  )

private

Definition at line 78 of file SmithMonizUtils.cxx.

{


  GetParam( "FermiMomentumTable", fKFTable);
  GetParam( "RFG-UseParametrization", fUseParametrization);


  // Load removal energy for specific nuclei from either the algorithm's
  // configuration file or the UserPhysicsOptions file.
  // If none is used use Wapstra's semi-empirical formula.
  const std::string keyStart = "RFG-NucRemovalE@Pdg=";

  RgIMap entries = GetConfig().GetItemMap();

  for(RgIMap::const_iterator it = entries.begin(); it != entries.end(); ++it)
  {
    const std::string& key = it->first;
    int pdg = 0;
    int Z = 0;
    if (0 == key.compare(0, keyStart.size(), keyStart.c_str()))
    {
      pdg = atoi(key.c_str() + keyStart.size());
      Z = pdg::IonPdgCodeToZ(pdg);
    }
    if (0 != pdg && 0 != Z)
    {
      ostringstream key_ss ;
      key_ss << keyStart << pdg;
      RgKey rgkey   = key_ss.str();
      double eb;
      GetParam( rgkey, eb ) ;
      eb = TMath::Max(eb, 0.);
      fNucRmvE.insert(map<int,double>::value_type(Z,eb));
    }
  }


}

double SmithMonizUtils::PhaseSpaceVolume

(

KinePhaseSpace_t

ps

)

const

Definition at line 527 of file SmithMonizUtils.cxx.

{
   double vol = 0.0;
   if (ps == kPSQ2fE)
   {
     Range1D_t rQ2 = Q2QES_SM_lim();
     vol = rQ2.max - rQ2.min;
     return vol;
   }
   else if (ps == kPSQ2vfE)
   {
     Range1D_t rQ2 = Q2QES_SM_lim();
     const int kNQ2 = 101;
     double dQ2 = (rQ2.max-rQ2.min)/(kNQ2-1);
     for(int iq2=0; iq2<kNQ2; iq2++)
     {
       double Q2 = (rQ2.min + iq2*dQ2);
       Range1D_t rvlims  = vQES_SM_lim(Q2);
       double dv = (rvlims.max-rvlims.min);
       vol += (dQ2*dv);
     }
     return vol;
   }
   else
     return 0;
}

double SmithMonizUtils::Q2lim1_SM ( double  Q2 ) const

private

Definition at line 265 of file SmithMonizUtils.cxx.

{
  double s = 2*Enu_in*m_tar+mm_tar;
  double nu_max = (s-mm_tar-mm_lep*(s-mm_tar)/(Q2+mm_lep)-mm_tar*(Q2+mm_lep)/(s-mm_tar))/(2*m_tar);
  double E = sqrt(P_Fermi*P_Fermi+mm_ini);
  double b = (E-E_BIN)*(E-E_BIN)-P_Fermi*P_Fermi;
  double a = 0.5*(Q2+mm_fin-b);
  double nu_1 = (a*(E-E_BIN)-P_Fermi*TMath::Sqrt(a*a+Q2*b))/b;
  return nu_1-nu_max;

}

double SmithMonizUtils::Q2lim2_SM ( double  Q2 ) const

private

Definition at line 278 of file SmithMonizUtils.cxx.

{
  double nu_min = ((m_rnu+m_fin)*(m_rnu+m_fin)+Q2-mm_tar)/(2*m_tar);
  double E = sqrt(P_Fermi*P_Fermi+mm_ini);
  double b = (E-E_BIN)*(E-E_BIN)-P_Fermi*P_Fermi;
  double a = (Q2+mm_fin-b)*0.5;
  double nu_2  = (a*(E-E_BIN)+P_Fermi*TMath::Sqrt(a*a+Q2*b))/b;
  return nu_min-nu_2;

}

Range1D_t SmithMonizUtils::Q2QES_SM_lim

(

void

)

const

Definition at line 290 of file SmithMonizUtils.cxx.

{


  const int MFC = 1000;      //  Maximum of function call
  const double EPS  = 1.0e-08;
  const double Delta= 1.0e-14;
  const double EPSABS = 0;
  const double EPSREL = 1.0e-08;

  Enu_in = E_nu;
  double s = 2*E_nu*m_tar+mm_tar;
  double W2 = (m_rnu+m_fin)*(m_rnu+m_fin);
  double c = 0.5*(W2+mm_lep-mm_tar*(W2-mm_lep)/s);
  double sqrtD = TMath::Sqrt(LambdaFUNCTION(1.0, mm_lep/s, W2/s));
  double tmp = 0.5*(s-mm_tar);
  double Q2_min = tmp*(1.0-sqrtD)-c;
  double Q2_max = tmp*(1.0+sqrtD)-c;
  // if the nucleus is hydrogen or deuterium then there is no need in further calculation
  if (E_BIN == 0 && P_Fermi == 0)
  {
    Q2_min= TMath::Max(Q2_min,0.0);
    Range1D_t R(Q2_min,Q2_max);
    return R;
  }
  double F_MIN, Q2_0;
  bool LLM;
  // C++ analog of fortran function DMINFC(Q2lim1_SM,Q2_min,Q2_max,EPS,Delta,Q2_0,F_MIN,LLM)
  Func1D<SmithMonizUtils> Q2lim1_SM_(*this, &SmithMonizUtils::Q2lim1_SM);
  DMINFC(Q2lim1_SM_,Q2_min,Q2_max,EPS,Delta,Q2_0,F_MIN,LLM);
  if (F_MIN>0)
  {
     LOG("SmithMoniz", pFATAL)
     << "No overlapped area for energy " << E_nu << "\n" <<
     "Q2_min=" << Q2_min << " Q2_max=" << Q2_max << "\n" <<
     "Q2_0=" << Q2_0 << " F_MIN=" << F_MIN;
     exit(1);
  }
  if (Q2lim1_SM(Q2_min)>0)
  {
    //C++ analog of fortran function Q2_RF=DZEROX(Q2_min,Q2_0,EPS,MFC,Q2lim1_SM,1)
    ROOT::Math::RootFinder rfgb(ROOT::Math::RootFinder::kGSL_BRENT);
    rfgb.Solve(Q2lim1_SM_, Q2_min, Q2_0, MFC, EPSABS, EPSREL); //convergence is reached using  tolerance = 2 *( epsrel * abs(x) + epsabs)
    double Q2_RF = rfgb.Root();
    Q2_min= TMath::Max(Q2_min,Q2_RF);
  }
  if(Q2lim1_SM(Q2_max)>0)
  {
     // C++ analog of fortran function Q2_RF=DZEROX(Q2_0,Q2_max,Eps,MFC,Q2lim1_SM,1)
     ROOT::Math::RootFinder rfgb(ROOT::Math::RootFinder::kGSL_BRENT);
     rfgb.Solve(Q2lim1_SM_, Q2_0, Q2_max, MFC, EPSABS, EPSREL); //convergence is reached using  tolerance = 2 *( epsrel * abs(x) + epsabs)
     double Q2_RF = rfgb.Root();
     Q2_max= TMath::Min(Q2_max,Q2_RF);
  }
  Func1D<SmithMonizUtils> Q2lim2_SM_(*this, &SmithMonizUtils::Q2lim2_SM);
  if (Q2lim2_SM(Q2_min)>0)
  {
     if(Q2lim2_SM(Q2_max)>0)
     {
           LOG("SmithMoniz", pWARN) << "The RFG model is not applicable! The cross section is set zero!";
             Q2_min = Q2_max;
     }
     else
     {
       // C++ analog of fortran function Q2_RF = DZEROX(Q2_min,Q2_max,Eps,MFC,Q2lim2_SM,1)
       ROOT::Math::RootFinder rfgb(ROOT::Math::RootFinder::kGSL_BRENT);
       rfgb.Solve(Q2lim2_SM_, Q2_min,Q2_max, MFC, EPSABS, EPSREL); //convergence is reached using  tolerance = 2 *( epsrel * abs(x) + epsabs)
       double Q2_RF = rfgb.Root();
       Q2_min= TMath::Max(Q2_min,Q2_RF);
     }
  }

  Q2_min= TMath::Max(Q2_min,0.0);
  Range1D_t R(Q2_min,Q2_max);
  return R;

}

double SmithMonizUtils::QEL_EnuMin_SM ( double  E_nu ) const

private

Definition at line 243 of file SmithMonizUtils.cxx.

{
  const double EPS  = 1.0e-06;
  const double Delta= 1.0e-14;
  const double Precision = std::numeric_limits<double>::epsilon();
  Enu_in = Enu;
  double s = 2*Enu*m_tar+mm_tar;
  double W2 = (m_rnu+m_fin)*(m_rnu+m_fin);
  double c  = 0.5*(W2+mm_lep-mm_tar*(W2-mm_lep)/s);
  double sqrtD  = TMath::Sqrt(TMath::Max(Precision,LambdaFUNCTION(1.0, mm_lep/s, W2/s)));
  double tmp    = 0.5*(s-mm_tar);
  double Q2_lim1= tmp*(1.0-sqrtD)-c;
  double Q2_lim2= tmp*(1.0+sqrtD)-c;
  // C++ analog of fortran function DMINFC(Q2lim1_SM,Q2_lim1,Q2_lim2,EPS,Delta,Q2_0,F_MIN,LLM)
  Func1D<SmithMonizUtils> Q2lim1_SM_(*this, &SmithMonizUtils::Q2lim1_SM);
  double Q2_0,F_MIN;
  bool LLM;
  DMINFC(Q2lim1_SM_,Q2_lim1,Q2_lim2,EPS,Delta,Q2_0,F_MIN,LLM);
  return F_MIN;
}

double SmithMonizUtils::rho ( double  P_Fermi, double  T_Fermi, double  p  )

static

Definition at line 479 of file SmithMonizUtils.cxx.

{

  if (T_Fermi==0)                             //Pure Fermi gaz with T_Fermi=0
    if(p<=P_Fermi)
      return 1.0;
    else
      return 0.0;
  else                                        //Fermi-Dirac distribution
      return 1.0/(1.0 + TMath::Exp(-(P_Fermi-p)/T_Fermi));


}

void SmithMonizUtils::SetInteraction

(

const Interaction *

i

)

Definition at line 119 of file SmithMonizUtils.cxx.

{

  fInteraction = interaction;
  // Get kinematics & init-state parameters
  // unused // const Kinematics &   kinematics = interaction -> Kine();
  const InitialState & init_state = interaction -> InitState();
  const Target & target = init_state.Tgt();
  PDGLibrary * pdglib = PDGLibrary::Instance();

  E_nu = interaction->InitState().ProbeE(kRfLab);         //  Neutrino energy (GeV)

  assert(target.HitNucIsSet());
  // get lepton&nuclear masses (init & final state nucleus)
  m_lep = interaction->FSPrimLepton()->Mass();          //  Mass of final charged lepton (GeV)
  mm_lep     = TMath::Power(m_lep,    2);
  int nucl_pdg_ini = target.HitNucPdg();
  m_ini  = target.HitNucMass();
  mm_ini = TMath::Power(m_ini,    2);
  int nucl_pdg_fin = genie::pdg::SwitchProtonNeutron(nucl_pdg_ini);
  TParticlePDG * nucl_fin = pdglib->Find( nucl_pdg_fin );
  m_fin  = nucl_fin -> Mass();                         //  Mass of final hadron or hadron system (GeV)
  mm_fin = TMath::Power(m_fin,    2);
  m_tar = target.Mass();                             //  Mass of target nucleus (GeV)
  mm_tar = TMath::Power(m_tar,    2);

  // For hydrogen and deuterium RFG is not applied
  if (target.A()<3)
  {
    E_BIN = P_Fermi = m_rnu = mm_rnu = 0;
    return;
  }

  bool is_p = pdg::IsProton(nucl_pdg_ini);
  int Zi = target.Z();
  int Ai = target.A();
  int Zf = (is_p) ? Zi-1 : Zi;
  int Af = Ai-1;
  TParticlePDG * nucl_f = pdglib->Find( pdg::IonPdgCode(Af, Zf) );
  if(!nucl_f)
  {
    LOG("SmithMoniz", pFATAL)
    << "Unknwown nuclear target! No target with code: "
    << pdg::IonPdgCode(Af, Zf) << " in PDGLibrary!";
    exit(1);
  }

  m_rnu = nucl_f -> Mass();                          //  Mass of residual nucleus (GeV)
  mm_rnu = TMath::Power(m_rnu, 2);

  int Z = target.Z();
  int A = target.A();
  int N = A-Z;


  // Maximum value of Fermi momentum of target nucleon (GeV)
  if (A < 6 || !fUseParametrization)
  {
     //-- look up the Fermi momentum for this Target
     FermiMomentumTablePool * kftp = FermiMomentumTablePool::Instance();
     const FermiMomentumTable * kft = kftp->GetTable(fKFTable);
     P_Fermi = kft->FindClosestKF(pdg::IonPdgCode(A, Z), nucl_pdg_ini);
  }
  else
  {
    //-- define the Fermi momentum for this Target
    //
    P_Fermi = utils::nuclear::FermiMomentumForIsoscalarNucleonParametrization(target);
    //-- correct the Fermi momentum for the struck nucleon
    if(is_p) P_Fermi *= TMath::Power( 2.*Z/A, 1./3);
    else
           P_Fermi *= TMath::Power( 2.*N/A, 1./3);
  }

  // Neutrino binding energy (GeV)
  if (target.A() < 6 || !fUseParametrization)
  {
     map<int,double>::const_iterator it = fNucRmvE.find(Z);
     if(it != fNucRmvE.end()) E_BIN = it->second;
       else E_BIN = utils::nuclear::BindEnergyPerNucleon(target);
  }
  else
    E_BIN = utils::nuclear::BindEnergyPerNucleonParametrization(target);




}

Range1D_t SmithMonizUtils::vQES_SM_lim

(

double

Q2

)

const

Definition at line 369 of file SmithMonizUtils.cxx.

{
  double s = 2*E_nu*m_tar+mm_tar;
  double nu_min= ((m_rnu+m_fin)*(m_rnu+m_fin)+Q2-mm_tar)/(2*m_tar);
  double nu_max= (s-mm_tar-mm_lep*(s-mm_tar)/(Q2+mm_lep)-mm_tar*(Q2+mm_lep)/(s-mm_tar))/(2*m_tar);
  double E = TMath::Sqrt(P_Fermi*P_Fermi+mm_ini);
  double b = (E-E_BIN)*(E-E_BIN)-P_Fermi*P_Fermi;
  double a = (Q2+mm_fin-b)*0.5;
  double tmp1 = a*(E-E_BIN);
  double tmp2  = P_Fermi*TMath::Sqrt(a*a+Q2*b);
  double nu_1 = (tmp1-tmp2)/b;
  double nu_2 = (tmp1+tmp2)/b;
  nu_min= TMath::Max(nu_min,nu_1);
  nu_max= TMath::Min(nu_max,nu_2);
  Range1D_t R;
  if (E_BIN == 0 && P_Fermi == 0)
    nu_max=nu_min;
  if (nu_min<=nu_max)
    R = Range1D_t(nu_min,nu_max);
  else
    R = Range1D_t(0.5*(nu_min+nu_max),0.5*(nu_min+nu_max));
  return R;

}

double SmithMonizUtils::vQES_SM_low_bound ( double  Q2 ) const

private

Definition at line 517 of file SmithMonizUtils.cxx.

{
        return vQES_SM_lim(Q2).min;
}

double SmithMonizUtils::vQES_SM_upper_bound ( double  Q2 ) const

private

Definition at line 522 of file SmithMonizUtils.cxx.

{
  return -1.0*vQES_SM_lim(Q2).max;
}

Member Data Documentation

double genie::SmithMonizUtils::E_BIN

private

Binding energy (GeV)

Definition at line 117 of file SmithMonizUtils.h.

double genie::SmithMonizUtils::E_nu

private

Neutrino energy (GeV)

Definition at line 105 of file SmithMonizUtils.h.

double genie::SmithMonizUtils::Enu_in

mutableprivate

Running neutrino energy (GeV)

Definition at line 118 of file SmithMonizUtils.h.

const Interaction* genie::SmithMonizUtils::fInteraction

private

Definition at line 100 of file SmithMonizUtils.h.

string genie::SmithMonizUtils::fKFTable

private

Definition at line 97 of file SmithMonizUtils.h.

map<int, double> genie::SmithMonizUtils::fNucRmvE

private

Definition at line 96 of file SmithMonizUtils.h.

bool genie::SmithMonizUtils::fUseParametrization

private

Definition at line 98 of file SmithMonizUtils.h.

double genie::SmithMonizUtils::m_fin

private

Mass of final hadron or hadron system (GeV)

Definition at line 110 of file SmithMonizUtils.h.

double genie::SmithMonizUtils::m_ini

private

Mass of initial hadron or hadron system (GeV)

Definition at line 108 of file SmithMonizUtils.h.

double genie::SmithMonizUtils::m_lep

private

Mass of final charged lepton (GeV)

Definition at line 106 of file SmithMonizUtils.h.

double genie::SmithMonizUtils::m_rnu

private

Mass of residual nucleus (GeV)

Definition at line 114 of file SmithMonizUtils.h.

double genie::SmithMonizUtils::m_tar

private

Mass of target nucleus (GeV)

Definition at line 112 of file SmithMonizUtils.h.

double genie::SmithMonizUtils::mm_fin

private

Squared mass of final hadron or hadron system (GeV)

Definition at line 111 of file SmithMonizUtils.h.

double genie::SmithMonizUtils::mm_ini

private

Sqared mass of initial hadron or hadron system (GeV)

Definition at line 109 of file SmithMonizUtils.h.

double genie::SmithMonizUtils::mm_lep

private

Squared mass of final charged lepton (GeV)

Definition at line 107 of file SmithMonizUtils.h.

double genie::SmithMonizUtils::mm_rnu

private

Squared mass of residual nucleus (GeV)

Definition at line 115 of file SmithMonizUtils.h.

double genie::SmithMonizUtils::mm_tar

private

Squared mass of target nucleus (GeV)

Definition at line 113 of file SmithMonizUtils.h.

double genie::SmithMonizUtils::P_Fermi

private

Maximum value of Fermi momentum of target nucleon (GeV)

Definition at line 116 of file SmithMonizUtils.h.

---

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

• Generator/src/Physics/QuasiElastic/XSection/SmithMonizUtils.h
• Generator/src/Physics/QuasiElastic/XSection/SmithMonizUtils.cxx