genie::FermiMomentumTable Class Reference

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 <constantinos.andreopoulos cern.ch> University of Liverpool & STFC Rutherford Appleton Laboratory

August 18, 2005

Copyright (c) 2003-2020, 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.

{
  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.

{
  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.

---

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

• Generator/src/Physics/NuclearState/FermiMomentumTable.h
• Generator/src/Physics/NuclearState/FermiMomentumTable.cxx