PaisQELLambdaPXSec.h

Go to the documentation of this file.

//____________________________________________________________________________
/*!
 
\class    genie::PaisQELLambdaPXSec
 
\brief    Implementation of the quasi-elastic scattering formula for
          production of particles with different masses than the target.
 
\ref      Weak Interactions at High Energies
          A. Pais, Annals of Physics 63, 361-392 (1971)
          Implemented here is equation 2.34 of the Pais paper, but ignoring
          lepton mass terms.  This equation is given also as Equation 3.37 in
          the Llewellyn-Smith paper, though this paper uses slightly
          different notation than that used in the Pais paper, and introduces
          a small error in the kinematic coefficient of the w2 term.  The
          notation here by and large follows that of the Llewelyn-Smith paper.
 
\author   Hugh Gallagher
          Tufts University
 
\created  June 10, 2004
 
\cpright  Copyright (c) 2003-2025, The GENIE Collaboration
          For the full text of the license visit http://copyright.genie-mc.org        
*/
//____________________________________________________________________________
 
#ifndef _PAIS_QEL_LAMBDA_PARTIAL_XSEC_H_
#define _PAIS_QEL_LAMBDA_PARTIAL_XSEC_H_
 
#include "Framework/EventGen/XSecAlgorithmI.h"
#include "Physics/QuasiElastic/XSection/QELFormFactors.h"
 
namespace genie {
 
class QELFormFactorsModelI;
class XSecIntegratorI;
 
class PaisQELLambdaPXSec : public XSecAlgorithmI {
 
public:
  PaisQELLambdaPXSec();
  PaisQELLambdaPXSec(string config);
  virtual ~PaisQELLambdaPXSec();
 
  //-- XSecAlgorithmI interface implementation
  double XSec            (const Interaction * i, KinePhaseSpace_t k) const;
  double Integral        (const Interaction * i) const;
  bool   ValidProcess    (const Interaction * i) const;
  bool   ValidKinematics (const Interaction * i) const;
 
  //-- override the Algorithm::Configure methods to load configuration
  //   data to private data members
  void Configure (const Registry & config);
  void Configure (string param_set);
 
private:
  void  LoadConfig (void);
  double MHyperon(const Interaction * interaction) const;
 
  mutable QELFormFactors          fFormFactors;
  const   QELFormFactorsModelI *  fFormFactorsModel;
  const   XSecIntegratorI *       fXSecIntegrator;
  double                          fSin8c2;
 
};
 
} // genie namespace
#endif  // _PAIS_QEL_LAMBDA_PARTIAL_XSEC_H_