ReinSehgalRESXSecFast.h

Go to the documentation of this file.

//____________________________________________________________________________
/*!

\class    genie::ReinSehgalRESXSecFast

\brief    Computes the cross section for an exclusive 1pi reaction through
          resonance neutrinoproduction according to the Rein-Sehgal model.

          This algorithm produces in principle what you could also get from
          the genie::RESXSec algorithm (RES cross section integrator) by
          specifying the genie::ReinSehgalRESPXSec as the differential
          cross section model. However, ReinSehgalRESXSecFast offers a faster
          alternative. Before computing any RES cross section this algorithm
          computes and caches splines for resonance neutrino-production cross
          sections. This improves the speed of the GENIE spline construction
          phase if splines for multiple nuclear targets are to be computed.
          Also this class integrates cross sections faster, than
          ReinSehgalRESXSec because of integration area transformation.

          Is a concrete implementation of the XSecAlgorithmI interface.\n

\ref      D.Rein and L.M.Sehgal, Neutrino Excitation of Baryon Resonances
          and Single Pion Production, Ann.Phys.133, 79 (1981)

\author   Igor Kakorin <kakorin@jinr.ru>
          Joint Institute for Nuclear Research

          based on code of
          Costas Andreopoulos <c.andreopoulos \at cern.ch>
          University of Liverpool

\created  March 01, 2017

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

*/
//____________________________________________________________________________

#ifndef _REIN_SEHGAL_RES_XSEC_FAST_H_
#define _REIN_SEHGAL_RES_XSEC_FAST_H_

#include "Physics/Resonance/XSection/ReinSehgalRESXSecWithCacheFast.h"

namespace genie {

class ReinSehgalRESXSecFast : public ReinSehgalRESXSecWithCacheFast {

public:
  ReinSehgalRESXSecFast();
  ReinSehgalRESXSecFast(string param_set);
  virtual ~ReinSehgalRESXSecFast();

  // XSecIntegratorI interface implementation
  double Integrate(const XSecAlgorithmI * model, const Interaction * i) const;

  // Overload the Algorithm::Configure() methods to load private data
  // members from configuration options
  void Configure(const Registry & config);
  void Configure(string config);

private:
  void LoadConfig(void);

  bool fUsePauliBlocking;      ///< account for Pauli blocking?
};

}       // genie namespace
#endif  // _REIN_SEHGAL_RES_XSEC_H_