EffectiveSF.h

Go to the documentation of this file.

//____________________________________________________________________________
/*!
 
\class    genie::EffectiveSF
 
\brief    An effective spectral function to match psi' superscaling.
          Implements the NuclearModelI interface.
 
\ref      http://arxiv.org/abs/1405.0583
 
\author   Brian Coopersmith, University of Rochester
 
\created  April 02, 2014
 
\cpright  Copyright (c) 2003-2025, The GENIE Collaboration
          For the full text of the license visit http://copyright.genie-mc.org
          
*/
//____________________________________________________________________________
 
#ifndef _EFFECTIVE_SF_H_
#define _EFFECTIVE_SF_H_
 
#include <map>
 
#include <TH1D.h>
 
#include "Physics/NuclearState/NuclearModelI.h"
 
using std::map;
 
namespace genie {
 
class EffectiveSF : public NuclearModelI {
 
public:
  EffectiveSF();
  EffectiveSF(string config);
  virtual ~EffectiveSF();
 
  using NuclearModelI::GenerateNucleon;  // inherit versions not overridden here
  using NuclearModelI::Prob;
 
  //-- implement the NuclearModelI interface
  bool           GenerateNucleon (const Target & t) const;
  double         Prob            (double mom, double w, const Target & t) const;
  NuclearModel_t ModelType       (const Target &) const
  {
    return kNucmEffSpectralFunc;
  }
 
  //-- override the Algorithm::Configure methods to load configuration
  //   data to private data members
  void Configure (const Registry & config);
  void Configure (string param_set);
 
private:
  TH1D * ProbDistro (const Target & t) const;
 
  TH1D * MakeEffectiveSF(const Target & target) const;
 
  TH1D * MakeEffectiveSF(double bs, double bp, double alpha, double beta,
                         double c1, double c2, double c3,
                         const Target & target) const;
 
  double ReturnBindingEnergy(const Target & target) const;
  double GetTransEnh1p1hMod(const Target& target) const;
 
  double Returnf1p1h(const Target & target) const;
  void   LoadConfig (void);
 
  mutable map<string, TH1D *> fProbDistroMap;
  double fPMax;
  double fPCutOff;
  bool   fEjectSecondNucleon2p2h;
  bool   fUseElFFTransEnh;
 
  // Map from PDG code to spectral function parameters
  map<int, double> fNucRmvE;
  map<int, double> f1p1hMap;
  map<int, std::vector<double> > fProbDistParams;
  map<int, double> fTransEnh1p1hMods;
 
  // Map from range of A (pair<lowA, highA> inclusive> to spectral
  // function parameters.
  map<pair<int, int>, double> fRangeNucRmvE;
  map<pair<int, int>, double> fRange1p1hMap;
  map<pair<int, int>, std::vector<double> > fRangeProbDistParams;
  map<pair<int, int>, double> fRangeTransEnh1p1hMods;
};
 
}         // genie namespace
#endif    // _EFFECTIVE_SF_H_