A table of Fermi momentum constants. More...

#include <FermiMomentumTable.h>

Public Member Functions
	FermiMomentumTable ()
	FermiMomentumTable (const FermiMomentumTable &fmt)
virtual	~FermiMomentumTable ()
double	FindClosestKF (int target_pdgc, int nucleon_pdgc) const
void	AddTableEntry (int target_pdgc, KF_t kf)

Private Attributes
map< int, KF_t >	fKFSets

Detailed Description

A table of Fermi momentum constants.

Author: Costas Andreopoulos <c.andreopoulos \at cern.ch> University of Liverpool

Created:\n August 18, 2005

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

Definition at line 33 of file FermiMomentumTable.h.

Constructor & Destructor Documentation

◆ FermiMomentumTable() [1/2]

FermiMomentumTable::FermiMomentumTable ( )

Definition at line 26 of file FermiMomentumTable.cxx.

27{

28}

Referenced by FermiMomentumTable().

◆ FermiMomentumTable() [2/2]

FermiMomentumTable::FermiMomentumTable ( const FermiMomentumTable & fmt )

Definition at line 30 of file FermiMomentumTable.cxx.

{
 
}

References FermiMomentumTable().

◆ ~FermiMomentumTable()

FermiMomentumTable::~FermiMomentumTable ( )

virtual

Definition at line 35 of file FermiMomentumTable.cxx.

{
 
}

Member Function Documentation

◆ AddTableEntry()

void FermiMomentumTable::AddTableEntry	(	int	target_pdgc,
		KF_t	kf )

Definition at line 40 of file FermiMomentumTable.cxx.

{
  fKFSets.insert(map<int, KF_t>::value_type(tgt_pdgc, kf));
}

References fKFSets.

Referenced by genie::FermiMomentumTablePool::ParseXMLTables().

◆ FindClosestKF()

double FermiMomentumTable::FindClosestKF	(	int	target_pdgc,
		int	nucleon_pdgc ) const

Definition at line 45 of file FermiMomentumTable.cxx.

{
  LOG("FermiP", pINFO)
       << "Finding Fermi momenta table entry for (tgt = "
                      << tgt_pdgc << ", nucl = " << nucleon_pdgc << ")";
 
  if(fKFSets.size()==0) {
      LOG("FermiP", pWARN)
         << "The Fermi momenta table is empty! Returning kf(tgt = "
                  << tgt_pdgc << ", nucl = " << nucleon_pdgc << ") = 0";
      return 0;
  }
 
  double kf=0;
  bool isp = pdg::IsProton(nucleon_pdgc);
 
  if(fKFSets.count(tgt_pdgc) == 1) {
     LOG("FermiP", pDEBUG) << "Got exact match in Fermi momenta table";
     map<int, KF_t>::const_iterator table_iter = fKFSets.find(tgt_pdgc);
     if(isp) kf = table_iter->second.p;
     else    kf = table_iter->second.n;
     LOG("FermiP", pINFO) << "kF = " << kf;
     return kf;
  }
  LOG("FermiP", pINFO) << "Couldn't find exact match in Fermi momenta table";
 
  int  Z   = pdg::IonPdgCodeToZ(tgt_pdgc);
  int    Ac=9999, Zc=9999, dZmin=9999;
  map<int, KF_t>::const_iterator kfiter;
  for(kfiter=fKFSets.begin(); kfiter!=fKFSets.end(); ++kfiter) {
    int pdgc = kfiter->first;
    int Zt = pdg::IonPdgCodeToZ(pdgc);
    int dZ = TMath::Abs(Zt-Z);
    if(dZ<dZmin) {
      dZmin = dZ;
      Zc = Zt;
      Ac = pdg::IonPdgCodeToA(pdgc);
      KF_t kft  = kfiter->second;
      if(isp) kf=kft.p;
      else    kf=kft.n;
    }
  }
  LOG("FermiP", pINFO)
       << "The closest nucleus in table is pdgc = " << pdg::IonPdgCode(Ac,Zc);
  LOG("FermiP", pINFO) << "kF = " << kf;
  return kf;
}

References fKFSets, genie::pdg::IonPdgCode(), genie::pdg::IonPdgCodeToA(), genie::pdg::IonPdgCodeToZ(), genie::pdg::IsProton(), LOG, genie::EKF_t::n, genie::EKF_t::p, pDEBUG, pINFO, and pWARN.

Referenced by main(), genie::utils::nuclear::NuclQELXSecSuppression(), genie::SmithMonizUtils::SetInteraction(), genie::BSKLNBaseRESPXSec2014::XSec(), genie::MKSPPPXSec2020::XSec(), genie::ReinSehgalRESPXSec::XSec(), genie::SuSAv2MECPXSec::XSec(), and genie::SuSAv2QELPXSec::XSecScaling().

Member Data Documentation

◆ fKFSets

map<int, KF_t> genie::FermiMomentumTable::fKFSets

private

Definition at line 44 of file FermiMomentumTable.h.

Referenced by AddTableEntry(), and FindClosestKF().

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

Public Member Functions

Private Attributes