BergerSehgalCOHPiPXSec2015.h

Go to the documentation of this file.

//____________________________________________________________________________
/*!
 
\class   genie::BergerSehgalCOHPiPXSec2015
 
\brief   Computes the double differential cross section for CC & NC coherent
         pion production according to the \b Berger-Sehgal model.
         v(vbar)A->v(vbar)Api0, vA->l-Api+, vbarA->l+Api-
 
         The t-dependence of the triple differential cross (d^3xsec/dxdydt)
         is integrated out.
 
         Is a concrete implementation of the XSecAlgorithmI interface.
 
\ref     PRD 79, 053003 (2009) by Berger and Sehgal
 
\author  G. Perdue, H. Gallagher, D. Cherdack
 
\created 2014
 
\cpright Copyright (c) 2003-2025, The GENIE Collaboration
         For the full text of the license visit http://copyright.genie-mc.org
*/
//____________________________________________________________________________
 
#ifndef _BERGER_SEHGAL_COHPI_PXSEC_2015_H_
#define _BERGER_SEHGAL_COHPI_PXSEC_2015_H_
 
#include "Framework/EventGen/XSecAlgorithmI.h"
 
namespace genie {
 
  class XSecIntegratorI;
 
  class BergerSehgalCOHPiPXSec2015 : public XSecAlgorithmI {
 
    public:
      BergerSehgalCOHPiPXSec2015();
      BergerSehgalCOHPiPXSec2015(string config);
      virtual ~BergerSehgalCOHPiPXSec2015();
 
      //-- XSecAlgorithmI interface implementation
      double XSec            (const Interaction * i, KinePhaseSpace_t k) const;
      double Integral        (const Interaction * i) const;
      bool   ValidProcess    (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);
 
      double ExactKinematicTerm(const Interaction * i) const;
      double PionCOMAbsMomentum(const Interaction * i) const;
 
      //-- private data members loaded from config Registry or set to defaults
      double fMa;          ///< axial mass
      double fRo;          ///< nuclear size scale parameter
      double fCos8c2;      ///< cos^2(Cabibbo angle)
      bool fRSPionXSec;    ///< Use Rein-Sehgal "style" pion-nucleon xsecs
 
      const XSecIntegratorI * fXSecIntegrator;
  };
 
}       // genie namespace
 
#endif  // _BERGER_SEHGAL_COHPI_PXSEC_2015_H_