GENIEGenerator
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
NuElectronXSec.cxx
Go to the documentation of this file.
1 //____________________________________________________________________________
2 /*
3  Copyright (c) 2003-2024, The GENIE Collaboration
4  For the full text of the license visit http://copyright.genie-mc.org
5 
6  Costas Andreopoulos <c.andreopoulos \at cern.ch>
7  University of Liverpool
8 */
9 //____________________________________________________________________________
10 
11 #include <TMath.h>
12 #include <Math/IFunction.h>
13 #include <Math/Integrator.h>
14 
15 #include "Framework/Conventions/GBuild.h"
23 
24 using namespace genie;
25 using namespace genie::constants;
26 
27 //____________________________________________________________________________
29 XSecIntegratorI("genie::NuElectronXSec")
30 {
31 
32 }
33 //____________________________________________________________________________
35 XSecIntegratorI("genie::NuElectronXSec", config)
36 {
37 
38 }
39 //____________________________________________________________________________
41 {
42 
43 }
44 //____________________________________________________________________________
46  const XSecAlgorithmI * model, const Interaction * in) const
47 {
48  if(! model->ValidProcess(in) ) return 0.;
49 
50  const KPhaseSpace & kps = in->PhaseSpace();
51  if(!kps.IsAboveThreshold()) {
52  LOG("NuEXSec", pDEBUG) << "*** Below energy threshold";
53  return 0;
54  }
55  Range1D_t yl = kps.Limits(kKVy);
56 
57  LOG("NuEXSec", pDEBUG) << "y = (" << yl.min << ", " << yl.max << ")";
58 
59  Interaction * interaction = new Interaction(*in);
60  interaction->SetBit(kISkipProcessChk);
61  //interaction->SetBit(kISkipKinematicChk);
62 
63  ROOT::Math::IBaseFunctionOneDim * func =
64  new utils::gsl::dXSec_dy_E(model, interaction);
65  ROOT::Math::IntegrationOneDim::Type ig_type =
67  ROOT::Math::Integrator ig(*func,ig_type,1,fGSLRelTol,fGSLMaxEval);
68  double xsec = ig.Integral(yl.min, yl.max) * (1E-38 * units::cm2);
69 
70  //LOG("NuEXSec", pDEBUG) << "*** XSec[ve-] (E=" << E << ") = " << xsec;
71 
72  delete interaction;
73  delete func;
74  return xsec;
75 }
76 //____________________________________________________________________________
77 void NuElectronXSec::Configure(const Registry & config)
78 {
79  Algorithm::Configure(config);
80  this->LoadConfig();
81 }
82 //____________________________________________________________________________
83 void NuElectronXSec::Configure(string config)
84 {
85  Algorithm::Configure(config);
86  this->LoadConfig();
87 }
88 //____________________________________________________________________________
90 {
91  // Get GSL integration type & relative tolerance
92  GetParamDef( "gsl-integration-type", fGSLIntgType, string( "adaptive" ) ) ;
93  GetParamDef( "gsl-relative-tolerance", fGSLRelTol, 0.01 ) ;
94  int max;
95  GetParamDef( "gsl-max-eval", max, 100000 ) ;
96  fGSLMaxEval = (unsigned int) max ;
97 }
98 //____________________________________________________________________________
Cross Section Calculation Interface.
const KPhaseSpace & PhaseSpace(void) const
Definition: Interaction.h:73
ROOT::Math::IntegrationOneDim::Type Integration1DimTypeFromString(string type)
Definition: GSLUtils.cxx:23
string fGSLIntgType
name of GSL numerical integrator
Cross Section Integrator Interface.
A simple [min,max] interval for doubles.
Definition: Range1.h:42
Range1D_t Limits(KineVar_t kvar) const
Return the kinematical variable limits.
Summary information for an interaction.
Definition: Interaction.h:56
#define LOG(stream, priority)
A macro that returns the requested log4cpp::Category appending a string (using the FILE...
Definition: Messenger.h:96
static constexpr double cm2
Definition: Units.h:69
virtual void Configure(const Registry &config)
Definition: Algorithm.cxx:62
Kinematical phase space.
Definition: KPhaseSpace.h:33
double func(double x, double y)
virtual bool ValidProcess(const Interaction *i) const =0
Can this cross section algorithm handle the input process?
double max
Definition: Range1.h:53
int fGSLMaxEval
GSL max evaluations.
double Integrate(const XSecAlgorithmI *model, const Interaction *i) const
XSecIntegratorI interface implementation.
A registry. Provides the container for algorithm configuration parameters.
Definition: Registry.h:65
bool IsAboveThreshold(void) const
Checks whether the interaction is above the energy threshold.
void Configure(const Registry &config)
double min
Definition: Range1.h:52
bool GetParamDef(const RgKey &name, T &p, const T &def) const
const UInt_t kISkipProcessChk
if set, skip process validity checks
Definition: Interaction.h:47
#define pDEBUG
Definition: Messenger.h:63
double fGSLRelTol
required relative tolerance (error)