PattonCEvNSPXSec.h

Go to the documentation of this file.

//____________________________________________________________________________
/*!
 
\class    genie::PattonCEvNSPXSec
 
\brief    Differential cross section for v+As coherent elastic scattering.
          Is a concrete implementation of the XSecAlgorithmI interface.
 
\ref      K.Patton, J.Engel, G.McLaughlin, and N.Schunck, arXiv:1207.0693v1
 
\author   Costas Andreopoulos <c.andreopoulos \at cern.ch>
          University of Liverpool
 
\created  July 12, 2019
 
\cpright  Copyright (c) 2003-2025, The GENIE Collaboration
          For the full text of the license visit http://copyright.genie-mc.org
*/
//____________________________________________________________________________
 
#ifndef _PATTON_ET_AL_COHERENT_ELASTIC_PXSEC_H_
#define _PATTON_ET_AL_COHERENT_ELASTIC_PXSEC_H_
 
#include "Framework/EventGen/XSecAlgorithmI.h"
 
namespace genie {
 
class XSecIntegratorI;
 
class PattonCEvNSPXSec : public XSecAlgorithmI {
 
public:
  PattonCEvNSPXSec();
  PattonCEvNSPXSec(string config);
  virtual ~PattonCEvNSPXSec();
 
  // XSecAlgorithmI interface implementation
  double XSec            (const Interaction * i, KinePhaseSpace_t k) const;
  double Integral        (const Interaction * i) const;
  bool   ValidProcess    (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);
 
  // Calculate nuclear density moments
  double NuclearDensityMoment(int A, int k) const;
 
  const XSecIntegratorI * fXSecIntegrator;  ///< cross section integrator
  double fSin2thw;                          ///< sin^2(weinberg angle)
 
  // Parameters used for the numerical integration yielding nuclear density moments
  double fNuclDensMomentCalc_UpperIntegrationLimit; // upper integration limit (in units of nuclear radii R0*A^1/3)
  double fNuclDensMomentCalc_RelativeTolerance;     // relative tolerance for numerical integrator
  double fNuclDensMomentCalc_AbsoluteTolerance;     // absolute tolerance for numerical integrator
  int    fNuclDensMomentCalc_MaxNumOfEvaluations;   // maximum number of integran evaluations in numerical integration
 
};
 
}       // genie namespace
 
#endif  // _PATTON_ET_AL_COHERENT_ELASTIC_PXSEC_H_