genie::KineGeneratorWithCache Class Reference

Abstract class. Provides a data caching mechanism for for concrete implementations of the EventRecordVisitorI interface, generating kinematics and wishing to cache maximum differential xsecs. More...

#include <KineGeneratorWithCache.h>

Inheritance diagram for genie::KineGeneratorWithCache:

Collaboration diagram for genie::KineGeneratorWithCache:

Protected Member Functions
	KineGeneratorWithCache ()
	KineGeneratorWithCache (string name)
	KineGeneratorWithCache (string name, string config)
	~KineGeneratorWithCache ()
virtual double	ComputeMaxXSec (const Interaction *in) const =0
virtual double	ComputeMaxXSec (const Interaction *in, const int nkey) const
virtual double	MaxXSec (GHepRecord *evrec, const int nkey=0) const
virtual double	FindMaxXSec (const Interaction *in, const int nkey=0) const
virtual void	CacheMaxXSec (const Interaction *in, double xsec, const int nkey=0) const
virtual double	Energy (const Interaction *in) const
virtual CacheBranchFx *	AccessCacheBranch (const Interaction *in, const int nkey=0) const
virtual void	AssertXSecLimits (const Interaction *in, double xsec, double xsec_max) const
Protected Member Functions inherited from genie::EventRecordVisitorI
	EventRecordVisitorI ()
	EventRecordVisitorI (string name)
	EventRecordVisitorI (string name, string config)
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.
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.
int	GetParamVectKeys (const std::string &comm_name, std::vector< RgKey > &k, bool is_top_call=true) const
template<class T>
int	GetParamMat (const std::string &comm_name, TMatrixT< T > &mat, bool is_top_call=true) const
	Handle to load matrix of parameters.
template<class T>
int	GetParamMatSym (const std::string &comm_name, TMatrixTSym< T > &mat, bool is_top_call=true) const
int	GetParamMatKeys (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
int	AddLowRegistry (Registry *rp, bool owns=true)
	add registry with lowest priority, also update ownership
int	MergeTopRegistry (const Registry &r)
int	AddTopRegisties (const vector< Registry * > &rs, bool owns=false)
	Add registries with top priority, also udated Ownerships.

Protected Attributes
const XSecAlgorithmI *	fXSecModel
double	fSafetyFactor
	ComputeMaxXSec -> ComputeMaxXSec * fSafetyFactor.
std::vector< double >	vSafetyFactors
	MaxXSec -> MaxXSec * fSafetyFactors[nkey].
int	fNumOfSafetyFactors
	Number of given safety factors.
std::vector< string >	vInterpolatorTypes
	Type of interpolator for each key in a branch.
int	fNumOfInterpolatorTypes
	Number of given interpolators types.
double	fMaxXSecDiffTolerance
	max{100*(xsec-maxxsec)/.5*(xsec+maxxsec)} if xsec>maxxsec
double	fEMin
	min E for which maxxsec is cached - forcing explicit calc.
bool	fGenerateUniformly
	uniform over allowed phase space + event weight?
Protected Attributes inherited from genie::Algorithm
bool	fAllowReconfig
bool	fOwnsSubstruc
	true if it owns its substructure (sub-algs,...)
AlgId	fID
	algorithm name and configuration set
vector< Registry * >	fConfVect
vector< bool >	fOwnerships
	ownership for every registry in fConfVect
AlgStatus_t	fStatus
	algorithm execution status
AlgMap *	fOwnedSubAlgMp
	local pool for owned sub-algs (taken out of the factory pool)

Additional Inherited Members
Public Member Functions inherited from genie::EventRecordVisitorI
virtual	~EventRecordVisitorI ()
virtual void	ProcessEventRecord (GHepRecord *event_rec) const =0
Public Member Functions inherited from genie::Algorithm
virtual	~Algorithm ()
virtual void	Configure (const Registry &config)
virtual void	Configure (string config)
virtual void	FindConfig (void)
virtual const Registry &	GetConfig (void) const
Registry *	GetOwnedConfig (void)
virtual const AlgId &	Id (void) const
	Get algorithm ID.
virtual AlgStatus_t	GetStatus (void) const
	Get algorithm status.
virtual bool	AllowReconfig (void) const
virtual AlgCmp_t	Compare (const Algorithm *alg) const
	Compare with input algorithm.
virtual void	SetId (const AlgId &id)
	Set algorithm ID.
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.
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)
static string	BuildParamMatKey (const std::string &comm_name, unsigned int i, unsigned int j)
static string	BuildParamMatRowSizeKey (const std::string &comm_name)
static string	BuildParamMatColSizeKey (const std::string &comm_name)

Detailed Description

Abstract class. Provides a data caching mechanism for for concrete implementations of the EventRecordVisitorI interface, generating kinematics and wishing to cache maximum differential xsecs.

This class provides some common implementation for handling (retrieving, creating, searching, adding to) the cache. The various super-classes should implement the ComputeMaxXSec(...) method for computing the maximum xsec in case it has not already being pushed into the cache at a previous iteration.
Update May 15, 2022 IK: It makes possible to cache several values having different keys. The example of using this opportunity see in the class QELEventGeneratorSM.

Author

Costas Andreopoulos <c.andreopoulos \at cern.ch> University of Liverpool
Igor Kakorin kakor.nosp@m.in@j.nosp@m.inr.r.nosp@m.u Joint Institute for Nuclear Research

Created:\n December 15, 2004

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 49 of file KineGeneratorWithCache.h.

Constructor & Destructor Documentation

KineGeneratorWithCache() [1/3]

KineGeneratorWithCache::KineGeneratorWithCache ( )

protected

Definition at line 34 of file KineGeneratorWithCache.cxx.

                                               : 
EventRecordVisitorI(), fSafetyFactor(1.), fNumOfSafetyFactors(-1), fNumOfInterpolatorTypes(-1)
{
 
}

References genie::EventRecordVisitorI::EventRecordVisitorI(), fNumOfInterpolatorTypes, fNumOfSafetyFactors, and fSafetyFactor.

Referenced by genie::COHKinematicsGenerator::COHKinematicsGenerator(), genie::COHKinematicsGenerator::COHKinematicsGenerator(), genie::DFRKinematicsGenerator::DFRKinematicsGenerator(), genie::DFRKinematicsGenerator::DFRKinematicsGenerator(), genie::DISKinematicsGenerator::DISKinematicsGenerator(), genie::DISKinematicsGenerator::DISKinematicsGenerator(), genie::DMDISKinematicsGenerator::DMDISKinematicsGenerator(), genie::DMDISKinematicsGenerator::DMDISKinematicsGenerator(), genie::DMEKinematicsGenerator::DMEKinematicsGenerator(), genie::DMEKinematicsGenerator::DMEKinematicsGenerator(), genie::DMELEventGenerator::DMELEventGenerator(), genie::DMELEventGenerator::DMELEventGenerator(), genie::DMELKinematicsGenerator::DMELKinematicsGenerator(), genie::DMELKinematicsGenerator::DMELKinematicsGenerator(), genie::HEDISKinematicsGenerator::HEDISKinematicsGenerator(), genie::HEDISKinematicsGenerator::HEDISKinematicsGenerator(), genie::HELeptonKinematicsGenerator::HELeptonKinematicsGenerator(), genie::HELeptonKinematicsGenerator::HELeptonKinematicsGenerator(), genie::IBDKinematicsGenerator::IBDKinematicsGenerator(), genie::IBDKinematicsGenerator::IBDKinematicsGenerator(), genie::NuEKinematicsGenerator::NuEKinematicsGenerator(), genie::NuEKinematicsGenerator::NuEKinematicsGenerator(), genie::QELEventGenerator::QELEventGenerator(), genie::QELEventGenerator::QELEventGenerator(), genie::QELEventGeneratorSM::QELEventGeneratorSM(), genie::QELEventGeneratorSM::QELEventGeneratorSM(), genie::QELEventGeneratorSuSA::QELEventGeneratorSuSA(), genie::QELEventGeneratorSuSA::QELEventGeneratorSuSA(), genie::QELKinematicsGenerator::QELKinematicsGenerator(), genie::QELKinematicsGenerator::QELKinematicsGenerator(), genie::RESKinematicsGenerator::RESKinematicsGenerator(), genie::RESKinematicsGenerator::RESKinematicsGenerator(), genie::SKKinematicsGenerator::SKKinematicsGenerator(), genie::SKKinematicsGenerator::SKKinematicsGenerator(), genie::SPPEventGenerator::SPPEventGenerator(), and genie::SPPEventGenerator::SPPEventGenerator().

KineGeneratorWithCache() [2/3]

KineGeneratorWithCache::KineGeneratorWithCache ( string name )

protected

Definition at line 40 of file KineGeneratorWithCache.cxx.

                                                          : 
EventRecordVisitorI(name), fSafetyFactor(1.), fNumOfSafetyFactors(-1), fNumOfInterpolatorTypes(-1)
{
 
}

References genie::EventRecordVisitorI::EventRecordVisitorI(), fNumOfInterpolatorTypes, fNumOfSafetyFactors, and fSafetyFactor.

KineGeneratorWithCache() [3/3]

KineGeneratorWithCache::KineGeneratorWithCache ( string name, string config )

protected

Definition at line 46 of file KineGeneratorWithCache.cxx.

                                                                         : 
EventRecordVisitorI(name, config), fSafetyFactor(1.), fNumOfSafetyFactors(-1), fNumOfInterpolatorTypes(-1)
{
 
}

References genie::EventRecordVisitorI::EventRecordVisitorI(), fNumOfInterpolatorTypes, fNumOfSafetyFactors, and fSafetyFactor.

~KineGeneratorWithCache()

KineGeneratorWithCache::~KineGeneratorWithCache ( )

protected

Definition at line 52 of file KineGeneratorWithCache.cxx.

{
 
}

Member Function Documentation

AccessCacheBranch()

CacheBranchFx * KineGeneratorWithCache::AccessCacheBranch ( const Interaction * in, const int nkey = 0 ) const

protectedvirtual

Definition at line 224 of file KineGeneratorWithCache.cxx.

{
// Returns the cache branch for this algorithm and this interaction. If no
// branch is found then one is created.
 
  Cache * cache = Cache::Instance();
 
  // build the cache branch key as: namespace::algorithm/config/interaction/nkey
  string algkey = this->Id().Key();
  string intkey = interaction->AsString();
  string key    = cache->CacheBranchKey(algkey, intkey, std::to_string(nkey));
 
  CacheBranchFx * cache_branch =
              dynamic_cast<CacheBranchFx *> (cache->FindCacheBranch(key));
  if(!cache_branch) {
    //-- create the cache branch at the first pass
    LOG("Kinematics", pINFO) << "No cache branch found";
    LOG("Kinematics", pINFO) << "Creating cache branch - key = " << key;
 
    cache_branch = new CacheBranchFx("Max over phase space");
    cache->AddCacheBranch(key, cache_branch);
  }
  assert(cache_branch);
 
  return cache_branch;
}

References genie::Cache::AddCacheBranch(), genie::Interaction::AsString(), genie::Cache::CacheBranchKey(), genie::Cache::FindCacheBranch(), genie::Algorithm::Id(), genie::Cache::Instance(), genie::AlgId::Key(), LOG, and pINFO.

Referenced by CacheMaxXSec(), and FindMaxXSec().

AssertXSecLimits()

void KineGeneratorWithCache::AssertXSecLimits ( const Interaction * in, double xsec, double xsec_max ) const

protectedvirtual

Definition at line 252 of file KineGeneratorWithCache.cxx.

{
  // check the computed cross section for the current kinematics against the
  // maximum cross section used in the rejection MC method for the current
  // interaction at the current energy.
  if(xsec>xsec_max) {
    double f = 200*(xsec-xsec_max)/(xsec_max+xsec);
    if(f>fMaxXSecDiffTolerance) {
       LOG("Kinematics", pFATAL)
          << "xsec: (curr) = " << xsec
                 << " > (max) = " << xsec_max << "\n for " << *interaction;
       LOG("Kinematics", pFATAL)
          << "*** Exceeding estimated maximum differential cross section";
       std::terminate();
    } else {
       LOG("Kinematics", pWARN)
          << "xsec: (curr) = " << xsec
                 << " > (max) = " << xsec_max << "\n for " << *interaction;
       LOG("Kinematics", pWARN)
            << "*** The fractional deviation of " << f << " % was allowed";
    }
  }
 
  // this should never happen - print an error mesg just in case...
  if(xsec<0) {
    LOG("Kinematics", pERROR)
     << "Negative cross section for current kinematics!! \n" << *interaction;
  }
}

References fMaxXSecDiffTolerance, LOG, pERROR, pFATAL, and pWARN.

Referenced by genie::COHKinematicsGenerator::CalculateKin_AlvarezRuso(), genie::SKKinematicsGenerator::CalculateKin_AtharSingleKaon(), genie::COHKinematicsGenerator::CalculateKin_ReinSehgal(), genie::DFRKinematicsGenerator::ProcessEventRecord(), genie::DISKinematicsGenerator::ProcessEventRecord(), genie::DMDISKinematicsGenerator::ProcessEventRecord(), genie::DMEKinematicsGenerator::ProcessEventRecord(), genie::DMELEventGenerator::ProcessEventRecord(), genie::DMELKinematicsGenerator::ProcessEventRecord(), genie::HEDISKinematicsGenerator::ProcessEventRecord(), genie::HELeptonKinematicsGenerator::ProcessEventRecord(), genie::IBDKinematicsGenerator::ProcessEventRecord(), genie::NuEKinematicsGenerator::ProcessEventRecord(), genie::QELEventGenerator::ProcessEventRecord(), genie::QELEventGeneratorSM::ProcessEventRecord(), genie::QELKinematicsGenerator::ProcessEventRecord(), genie::RESKinematicsGenerator::ProcessEventRecord(), genie::SPPEventGenerator::ProcessEventRecord(), genie::QELEventGeneratorSuSA::SelectLeptonKinematics(), genie::DMELKinematicsGenerator::SpectralFuncExperimentalCode(), and genie::QELKinematicsGenerator::SpectralFuncExperimentalCode().

CacheMaxXSec()

void KineGeneratorWithCache::CacheMaxXSec ( const Interaction * in, double xsec, const int nkey = 0 ) const

protectedvirtual

Definition at line 190 of file KineGeneratorWithCache.cxx.

{
  LOG("Kinematics", pINFO)
                       << "Adding the computed max value to cache";
  CacheBranchFx * cb = this->AccessCacheBranch(interaction, nkey);
 
  double E = this->Energy(interaction);
  if (E<fEMin) return;
  if(max_xsec>0) cb->AddValues(E,max_xsec);
 
  if(! cb->Spl() ) {
    if( cb->Map().size() > 40 ) 
      cb->CreateSpline(nkey<=fNumOfInterpolatorTypes-1?vInterpolatorTypes[nkey]:"");
  }
 
  if( cb->Spl() ) {
     if( E < cb->Spl()->XMin() || E > cb->Spl()->XMax() ) {
        cb->CreateSpline(nkey<=fNumOfInterpolatorTypes-1?vInterpolatorTypes[nkey]:"");
     }
  }
}

References AccessCacheBranch(), genie::CacheBranchFx::AddValues(), genie::CacheBranchFx::CreateSpline(), Energy(), fEMin, fNumOfInterpolatorTypes, LOG, genie::CacheBranchFx::Map(), pINFO, genie::CacheBranchFx::Spl(), vInterpolatorTypes, and genie::Spline::XMax().

Referenced by MaxXSec().

ComputeMaxXSec() [1/2]

virtual double genie::KineGeneratorWithCache::ComputeMaxXSec ( const Interaction * in ) const

protectedpure virtual

Implemented in genie::COHKinematicsGenerator, genie::DFRKinematicsGenerator, genie::DISKinematicsGenerator, genie::DMDISKinematicsGenerator, genie::DMEKinematicsGenerator, genie::DMELEventGenerator, genie::DMELKinematicsGenerator, genie::HEDISKinematicsGenerator, genie::HELeptonKinematicsGenerator, genie::IBDKinematicsGenerator, genie::NuEKinematicsGenerator, genie::QELEventGenerator, genie::QELEventGeneratorSM, genie::QELEventGeneratorSuSA, genie::QELKinematicsGenerator, genie::RESKinematicsGenerator, genie::SKKinematicsGenerator, and genie::SPPEventGenerator.

Referenced by ComputeMaxXSec(), and MaxXSec().

ComputeMaxXSec() [2/2]

double KineGeneratorWithCache::ComputeMaxXSec ( const Interaction * in, const int nkey ) const

protectedvirtual

Reimplemented in genie::QELEventGeneratorSM.

Definition at line 283 of file KineGeneratorWithCache.cxx.

{
  if (nkey == 0)
  {
     return this->ComputeMaxXSec(in);
  }
  else
  {
    return -1;
  }
}

References ComputeMaxXSec().

Energy()

double KineGeneratorWithCache::Energy ( const Interaction * in ) const

protectedvirtual

Reimplemented in genie::COHKinematicsGenerator, genie::DMEKinematicsGenerator, genie::HELeptonKinematicsGenerator, genie::NuEKinematicsGenerator, and genie::SKKinematicsGenerator.

Definition at line 213 of file KineGeneratorWithCache.cxx.

{
// Returns the neutrino energy at the struck nucleon rest frame. Kinematic
// generators should override this method if they need to cache the max-xsec
// values for another energy value (eg kinematic generators for IMD or COH)
 
  const InitialState & init_state = interaction->InitState();
  double E = init_state.ProbeE(kRfHitNucRest);
  return E;
}

References genie::Interaction::InitState(), genie::kRfHitNucRest, and genie::InitialState::ProbeE().

Referenced by CacheMaxXSec(), and FindMaxXSec().

FindMaxXSec()

double KineGeneratorWithCache::FindMaxXSec ( const Interaction * in, const int nkey = 0 ) const

protectedvirtual

Definition at line 105 of file KineGeneratorWithCache.cxx.

{
// Find a cached max xsec for the specified xsec algorithm & interaction and
// close to the specified energy
 
  // get neutrino energy
  double E = this->Energy(interaction);
  LOG("Kinematics", pINFO) << "E = " << E;
 
  if(E < fEMin) {
     LOG("Kinematics", pINFO)
         << "Below minimum energy - Forcing explicit calculation";
     return -1.;
  }
 
  // access the the cache branch
  CacheBranchFx * cb = this->AccessCacheBranch(interaction, nkey);
 
  // if there are enough points stored in the cache buffer to build a
  // spline, then intepolate
  if( cb->Spl() ) {
     if( E >= cb->Spl()->XMin() && E <= cb->Spl()->XMax()) {
       double spl_max_xsec = cb->Spl()->Evaluate(E);
       LOG("Kinematics", pINFO)
          << "\nInterpolated: max (E=" << E << ") = " << spl_max_xsec;
       return spl_max_xsec;
     }
     LOG("Kinematics", pINFO)
          << "Outside spline boundaries - Forcing explicit calculation";
     return -1.;
  }
 
  // if there are not enough points at the cache buffer to have a spline,
  // look whether there is another point that is sufficiently close
  double dE = TMath::Min(0.25, 0.05*E);
  const map<double,double> & fmap = cb->Map();
  map<double,double>::const_iterator iter = fmap.lower_bound(E);
  if(iter != fmap.end()) {
     if(TMath::Abs(E - iter->first) < dE) return iter->second;
  }
 
  return -1;
 
/*
  // build the search rule
  double dE = TMath::Min(0.25, 0.05*E);
  ostringstream search;
  search << "(x-" << E << " < " << dE << ") && (x>=" << E << ")";
 
  // query for all the entries at a window around the current energy
  TSQLResult * result = cb->Ntuple()->Query("x:y", search.str().c_str());
  int nrows = result->GetRowCount();
  LOG("Kinematics", pDEBUG)
            << "Found " << nrows << " rows with " << search.str();
  if(nrows <= 0) {
     delete result;
     return -1;
  }
 
  // and now select the entry with the closest energy
  double max_xsec = -1.0;
  double Ep       = 0;
  double dEmin    = 999;
  TSQLRow * row   = 0;
  while( (row = result->Next()) ) {
     double cE    = atof( row->GetField(0) );
     double cxsec = atof( row->GetField(1) );
     double dE    = TMath::Abs(E-cE);
     if(dE < dEmin) {
        max_xsec = cxsec;
        Ep       = cE;
        dEmin    = TMath::Min(dE,dEmin);
     }
     delete row;
  }
  delete result;
 
  LOG("Kinematics", pINFO)
     << "\nRetrieved: max xsec = " << max_xsec << " cached at E = " << Ep;
 
  return max_xsec;
*/
}

References AccessCacheBranch(), Energy(), genie::Spline::Evaluate(), fEMin, LOG, genie::CacheBranchFx::Map(), pINFO, genie::CacheBranchFx::Spl(), and genie::Spline::XMin().

Referenced by MaxXSec().

MaxXSec()

double KineGeneratorWithCache::MaxXSec ( GHepRecord * evrec, const int nkey = 0 ) const

protectedvirtual

Definition at line 57 of file KineGeneratorWithCache.cxx.

{
  LOG("Kinematics", pINFO)
                << "Getting max. for the rejection method";
 
  double xsec_max = -1;
  Interaction * interaction = event_rec->Summary();
 
  LOG("Kinematics", pINFO)
                  << "Attempting to find a cached max value";
  xsec_max = this->FindMaxXSec(interaction, nkey);
  if(xsec_max>0) return nkey<=fNumOfSafetyFactors-1?vSafetyFactors[nkey]*xsec_max:xsec_max;
 
  LOG("Kinematics", pINFO)
                  << "Attempting to compute the max value";
  if (nkey == 0)
  {
    xsec_max = this->ComputeMaxXSec(interaction);
  }
  else
  {
    xsec_max = this->ComputeMaxXSec(interaction, nkey);
  }
  
  if(xsec_max>0) {
     LOG("Kinematics", pINFO) << "max = " << xsec_max;
     this->CacheMaxXSec(interaction, xsec_max, nkey);
     return nkey<=fNumOfSafetyFactors-1?vSafetyFactors[nkey]*xsec_max:xsec_max;
  }
 
  LOG("Kinematics", pNOTICE)
            << "Can not generate event kinematics max_xsec<=0)";
  // xsec for selected kinematics = 0
  event_rec->SetDiffXSec(0,kPSNull);
  // switch on error flag
  event_rec->EventFlags()->SetBitNumber(kKineGenErr, true);
  // reset 'trust' bits
  interaction->ResetBit(kISkipProcessChk);
  interaction->ResetBit(kISkipKinematicChk);
  // throw exception
  genie::exceptions::EVGThreadException exception;
  exception.SetReason("kinematics generation: max_xsec<=0");
  exception.SwitchOnFastForward();
  throw exception;
 
  return 0;
}

References CacheMaxXSec(), ComputeMaxXSec(), genie::GHepRecord::EventFlags(), FindMaxXSec(), fNumOfSafetyFactors, genie::kISkipKinematicChk, genie::kISkipProcessChk, genie::kKineGenErr, genie::kPSNull, LOG, pINFO, pNOTICE, genie::GHepRecord::SetDiffXSec(), genie::exceptions::EVGThreadException::SetReason(), genie::GHepRecord::Summary(), genie::exceptions::EVGThreadException::SwitchOnFastForward(), and vSafetyFactors.

Referenced by genie::COHKinematicsGenerator::CalculateKin_AlvarezRuso(), genie::SKKinematicsGenerator::CalculateKin_AtharSingleKaon(), genie::COHKinematicsGenerator::CalculateKin_BergerSehgal(), genie::COHKinematicsGenerator::CalculateKin_BergerSehgalFM(), genie::COHKinematicsGenerator::CalculateKin_ReinSehgal(), genie::DFRKinematicsGenerator::ProcessEventRecord(), genie::DISKinematicsGenerator::ProcessEventRecord(), genie::DMDISKinematicsGenerator::ProcessEventRecord(), genie::DMEKinematicsGenerator::ProcessEventRecord(), genie::DMELEventGenerator::ProcessEventRecord(), genie::DMELKinematicsGenerator::ProcessEventRecord(), genie::HELeptonKinematicsGenerator::ProcessEventRecord(), genie::IBDKinematicsGenerator::ProcessEventRecord(), genie::NuEKinematicsGenerator::ProcessEventRecord(), genie::QELEventGenerator::ProcessEventRecord(), genie::QELEventGeneratorSM::ProcessEventRecord(), genie::QELKinematicsGenerator::ProcessEventRecord(), genie::RESKinematicsGenerator::ProcessEventRecord(), genie::SPPEventGenerator::ProcessEventRecord(), genie::QELEventGeneratorSuSA::SelectLeptonKinematics(), genie::DMELKinematicsGenerator::SpectralFuncExperimentalCode(), and genie::QELKinematicsGenerator::SpectralFuncExperimentalCode().

Member Data Documentation

fEMin

double genie::KineGeneratorWithCache::fEMin

protected

min E for which maxxsec is cached - forcing explicit calc.

Definition at line 76 of file KineGeneratorWithCache.h.

Referenced by CacheMaxXSec(), FindMaxXSec(), genie::COHKinematicsGenerator::LoadConfig(), genie::DFRKinematicsGenerator::LoadConfig(), genie::DISKinematicsGenerator::LoadConfig(), genie::DMDISKinematicsGenerator::LoadConfig(), genie::DMEKinematicsGenerator::LoadConfig(), genie::DMELEventGenerator::LoadConfig(), genie::DMELKinematicsGenerator::LoadConfig(), genie::IBDKinematicsGenerator::LoadConfig(), genie::NuEKinematicsGenerator::LoadConfig(), genie::QELEventGenerator::LoadConfig(), genie::QELEventGeneratorSM::LoadConfig(), genie::QELEventGeneratorSuSA::LoadConfig(), genie::QELKinematicsGenerator::LoadConfig(), genie::RESKinematicsGenerator::LoadConfig(), genie::SKKinematicsGenerator::LoadConfig(), and genie::SPPEventGenerator::LoadConfig().

fGenerateUniformly

bool genie::KineGeneratorWithCache::fGenerateUniformly

protected

uniform over allowed phase space + event weight?

Definition at line 77 of file KineGeneratorWithCache.h.

Referenced by genie::COHKinematicsGenerator::CalculateKin_AlvarezRuso(), genie::SKKinematicsGenerator::CalculateKin_AtharSingleKaon(), genie::COHKinematicsGenerator::CalculateKin_ReinSehgal(), genie::COHKinematicsGenerator::LoadConfig(), genie::DFRKinematicsGenerator::LoadConfig(), genie::DISKinematicsGenerator::LoadConfig(), genie::DMDISKinematicsGenerator::LoadConfig(), genie::DMEKinematicsGenerator::LoadConfig(), genie::DMELEventGenerator::LoadConfig(), genie::DMELKinematicsGenerator::LoadConfig(), genie::IBDKinematicsGenerator::LoadConfig(), genie::NuEKinematicsGenerator::LoadConfig(), genie::QELEventGenerator::LoadConfig(), genie::QELEventGeneratorSM::LoadConfig(), genie::QELEventGeneratorSuSA::LoadConfig(), genie::QELKinematicsGenerator::LoadConfig(), genie::RESKinematicsGenerator::LoadConfig(), genie::SKKinematicsGenerator::LoadConfig(), genie::SPPEventGenerator::LoadConfig(), genie::COHKinematicsGenerator::ProcessEventRecord(), genie::DFRKinematicsGenerator::ProcessEventRecord(), genie::DISKinematicsGenerator::ProcessEventRecord(), genie::DMDISKinematicsGenerator::ProcessEventRecord(), genie::DMEKinematicsGenerator::ProcessEventRecord(), genie::DMELEventGenerator::ProcessEventRecord(), genie::DMELKinematicsGenerator::ProcessEventRecord(), genie::HELeptonKinematicsGenerator::ProcessEventRecord(), genie::IBDKinematicsGenerator::ProcessEventRecord(), genie::NuEKinematicsGenerator::ProcessEventRecord(), genie::QELEventGenerator::ProcessEventRecord(), genie::QELEventGeneratorSM::ProcessEventRecord(), genie::QELKinematicsGenerator::ProcessEventRecord(), genie::RESKinematicsGenerator::ProcessEventRecord(), genie::SKKinematicsGenerator::ProcessEventRecord(), and genie::SPPEventGenerator::ProcessEventRecord().

fMaxXSecDiffTolerance

double genie::KineGeneratorWithCache::fMaxXSecDiffTolerance

protected

max{100*(xsec-maxxsec)/.5*(xsec+maxxsec)} if xsec>maxxsec

Definition at line 75 of file KineGeneratorWithCache.h.

Referenced by AssertXSecLimits(), genie::COHKinematicsGenerator::LoadConfig(), genie::DFRKinematicsGenerator::LoadConfig(), genie::DISKinematicsGenerator::LoadConfig(), genie::DMDISKinematicsGenerator::LoadConfig(), genie::DMEKinematicsGenerator::LoadConfig(), genie::DMELEventGenerator::LoadConfig(), genie::DMELKinematicsGenerator::LoadConfig(), genie::HEDISKinematicsGenerator::LoadConfig(), genie::HELeptonKinematicsGenerator::LoadConfig(), genie::IBDKinematicsGenerator::LoadConfig(), genie::NuEKinematicsGenerator::LoadConfig(), genie::QELEventGenerator::LoadConfig(), genie::QELEventGeneratorSM::LoadConfig(), genie::QELEventGeneratorSuSA::LoadConfig(), genie::QELKinematicsGenerator::LoadConfig(), genie::RESKinematicsGenerator::LoadConfig(), genie::SKKinematicsGenerator::LoadConfig(), and genie::SPPEventGenerator::LoadConfig().

fNumOfInterpolatorTypes

int genie::KineGeneratorWithCache::fNumOfInterpolatorTypes

protected

Number of given interpolators types.

Definition at line 74 of file KineGeneratorWithCache.h.

Referenced by CacheMaxXSec(), KineGeneratorWithCache(), KineGeneratorWithCache(), KineGeneratorWithCache(), and genie::QELEventGeneratorSM::LoadConfig().

fNumOfSafetyFactors

int genie::KineGeneratorWithCache::fNumOfSafetyFactors

protected

Number of given safety factors.

Definition at line 72 of file KineGeneratorWithCache.h.

Referenced by KineGeneratorWithCache(), KineGeneratorWithCache(), KineGeneratorWithCache(), genie::QELEventGeneratorSM::LoadConfig(), and MaxXSec().

fSafetyFactor

double genie::KineGeneratorWithCache::fSafetyFactor

protected

ComputeMaxXSec -> ComputeMaxXSec * fSafetyFactor.

Definition at line 70 of file KineGeneratorWithCache.h.

Referenced by genie::COHKinematicsGenerator::ComputeMaxXSec(), genie::DMEKinematicsGenerator::ComputeMaxXSec(), genie::DMELEventGenerator::ComputeMaxXSec(), genie::DMELKinematicsGenerator::ComputeMaxXSec(), genie::HELeptonKinematicsGenerator::ComputeMaxXSec(), genie::IBDKinematicsGenerator::ComputeMaxXSec(), genie::NuEKinematicsGenerator::ComputeMaxXSec(), genie::QELEventGenerator::ComputeMaxXSec(), genie::QELEventGeneratorSuSA::ComputeMaxXSec(), genie::QELKinematicsGenerator::ComputeMaxXSec(), genie::RESKinematicsGenerator::ComputeMaxXSec(), genie::SKKinematicsGenerator::ComputeMaxXSec(), genie::SPPEventGenerator::ComputeMaxXSec(), KineGeneratorWithCache(), KineGeneratorWithCache(), KineGeneratorWithCache(), genie::COHKinematicsGenerator::LoadConfig(), genie::DFRKinematicsGenerator::LoadConfig(), genie::DISKinematicsGenerator::LoadConfig(), genie::DMDISKinematicsGenerator::LoadConfig(), genie::DMEKinematicsGenerator::LoadConfig(), genie::DMELEventGenerator::LoadConfig(), genie::DMELKinematicsGenerator::LoadConfig(), genie::HEDISKinematicsGenerator::LoadConfig(), genie::HELeptonKinematicsGenerator::LoadConfig(), genie::IBDKinematicsGenerator::LoadConfig(), genie::NuEKinematicsGenerator::LoadConfig(), genie::QELEventGenerator::LoadConfig(), genie::QELEventGeneratorSuSA::LoadConfig(), genie::QELKinematicsGenerator::LoadConfig(), genie::RESKinematicsGenerator::LoadConfig(), genie::SKKinematicsGenerator::LoadConfig(), genie::SPPEventGenerator::LoadConfig(), genie::COHKinematicsGenerator::MaxXSec_AlvarezRuso(), and genie::HEDISKinematicsGenerator::ProcessEventRecord().

fXSecModel

const XSecAlgorithmI* genie::KineGeneratorWithCache::fXSecModel

mutableprotected

Definition at line 68 of file KineGeneratorWithCache.h.

Referenced by genie::COHKinematicsGenerator::CalculateKin_AlvarezRuso(), genie::SKKinematicsGenerator::CalculateKin_AtharSingleKaon(), genie::COHKinematicsGenerator::CalculateKin_BergerSehgal(), genie::COHKinematicsGenerator::CalculateKin_BergerSehgalFM(), genie::COHKinematicsGenerator::CalculateKin_ReinSehgal(), genie::COHKinematicsGenerator::ComputeMaxXSec(), genie::DFRKinematicsGenerator::ComputeMaxXSec(), genie::DISKinematicsGenerator::ComputeMaxXSec(), genie::DMDISKinematicsGenerator::ComputeMaxXSec(), genie::DMEKinematicsGenerator::ComputeMaxXSec(), genie::DMELEventGenerator::ComputeMaxXSec(), genie::DMELKinematicsGenerator::ComputeMaxXSec(), genie::HELeptonKinematicsGenerator::ComputeMaxXSec(), genie::IBDKinematicsGenerator::ComputeMaxXSec(), genie::NuEKinematicsGenerator::ComputeMaxXSec(), genie::QELEventGenerator::ComputeMaxXSec(), genie::QELEventGeneratorSM::ComputeMaxXSec(), genie::QELEventGeneratorSM::ComputeMaxXSec(), genie::QELEventGeneratorSuSA::ComputeMaxXSec(), genie::QELKinematicsGenerator::ComputeMaxXSec(), genie::RESKinematicsGenerator::ComputeMaxXSec(), genie::SKKinematicsGenerator::ComputeMaxXSec(), genie::SPPEventGenerator::ComputeMaxXSec(), genie::COHKinematicsGenerator::MaxXSec_AlvarezRuso(), genie::COHKinematicsGenerator::MaxXSec_BergerSehgal(), genie::COHKinematicsGenerator::MaxXSec_BergerSehgalFM(), genie::COHKinematicsGenerator::MaxXSec_ReinSehgal(), genie::COHKinematicsGenerator::ProcessEventRecord(), genie::DFRKinematicsGenerator::ProcessEventRecord(), genie::DISKinematicsGenerator::ProcessEventRecord(), genie::DMDISKinematicsGenerator::ProcessEventRecord(), genie::DMEKinematicsGenerator::ProcessEventRecord(), genie::DMELEventGenerator::ProcessEventRecord(), genie::DMELKinematicsGenerator::ProcessEventRecord(), genie::HEDISKinematicsGenerator::ProcessEventRecord(), genie::HELeptonKinematicsGenerator::ProcessEventRecord(), genie::IBDKinematicsGenerator::ProcessEventRecord(), genie::NuEKinematicsGenerator::ProcessEventRecord(), genie::QELEventGenerator::ProcessEventRecord(), genie::QELEventGeneratorSM::ProcessEventRecord(), genie::QELEventGeneratorSuSA::ProcessEventRecord(), genie::QELKinematicsGenerator::ProcessEventRecord(), genie::RESKinematicsGenerator::ProcessEventRecord(), genie::SKKinematicsGenerator::ProcessEventRecord(), genie::SPPEventGenerator::ProcessEventRecord(), genie::HEDISKinematicsGenerator::Scan(), genie::QELEventGeneratorSuSA::SelectLeptonKinematics(), genie::DMELKinematicsGenerator::SpectralFuncExperimentalCode(), and genie::QELKinematicsGenerator::SpectralFuncExperimentalCode().

vInterpolatorTypes

std::vector<string> genie::KineGeneratorWithCache::vInterpolatorTypes

protected

Type of interpolator for each key in a branch.

Definition at line 73 of file KineGeneratorWithCache.h.

Referenced by CacheMaxXSec(), and genie::QELEventGeneratorSM::LoadConfig().

vSafetyFactors

std::vector<double> genie::KineGeneratorWithCache::vSafetyFactors

protected

MaxXSec -> MaxXSec * fSafetyFactors[nkey].

Definition at line 71 of file KineGeneratorWithCache.h.

Referenced by genie::QELEventGeneratorSM::LoadConfig(), and MaxXSec().

---

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

• KineGeneratorWithCache.h
• KineGeneratorWithCache.cxx