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 <constantinos.andreopoulos \at cern.ch>
          University of Liverpool & STFC Rutherford Appleton Laboratory

\created  July 12, 2019

\cpright  Copyright (c) 2003-2020, 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_