genie::FermiMomentumTable Class Reference

A table of Fermi momentum constants. More...

#include <FermiMomentumTable.h>

Collaboration diagram for genie::FermiMomentumTable:

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 cern.ch> University of Liverpool

Created:

August 18, 2005

License:

Copyright (c) 2003-2024, 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::FermiMomentumTable

(

)

Definition at line 26 of file FermiMomentumTable.cxx.

{
}

FermiMomentumTable::FermiMomentumTable

(

const FermiMomentumTable &

fmt

)

Definition at line 30 of file FermiMomentumTable.cxx.

{

}

FermiMomentumTable::~FermiMomentumTable ( )

virtual

Definition at line 35 of file FermiMomentumTable.cxx.

{

}

Member Function Documentation

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

Definition at line 40 of file FermiMomentumTable.cxx.

References fKFSets.

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

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

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

Definition at line 45 of file FermiMomentumTable.cxx.

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 genie::NievesQELCCPXSec::CNCTCLimUcalc(), genie::PauliBlocker::GetFermiMomentum(), main(), genie::utils::nuclear::NuclQELXSecSuppression(), genie::SmithMonizUtils::SetInteraction(), genie::SuSAv2MECPXSec::XSec(), genie::ReinSehgalRESPXSec::XSec(), genie::BSKLNBaseRESPXSec2014::XSec(), genie::MKSPPPXSec2020::XSec(), and genie::SuSAv2QELPXSec::XSecScaling().

{
  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;
}

Member Data Documentation

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:

• FermiMomentumTable.h
• FermiMomentumTable.cxx