5 #if !defined(__CINT__) || defined(__MAKECINT__) 11 #include "EVGCore/InteractionList.h" 12 #include "VLE/StrumiaVissaniIBDPXSec.h" 13 #include "PDG/PDGCodeList.h" 14 #include "PDG/PDGUtils.h" 19 #include "PDG/PDGCodes.h" 26 2.01, 2.25, 2.51, 2.80,
27 3.12, 3.48, 3.89, 4.33, 4.84,
28 5.40, 6.02, 6.72, 7.49, 8.36,
29 8.83, 9.85, 11.0, 12.3, 13.7,
30 15.3, 17.0, 19.0, 21.2, 23.6,
31 26.4, 29.4, 32.8, 36.6, 40.9,
32 43.2, 48.2, 53.7, 59.9, 66.9,
33 74.6, 83.2, 92.9, 104, 116,
34 129, 144, 160, 179, 200
37 0.00351, 0.00735, 0.0127, 0.0202,
38 0.0304, 0.0440, 0.06190, 0.0854, 0.116,
39 0.155, 0.205, 0.269, 0.349, 0.451,
40 0.511, 0.654, 0.832, 1.05, 1.33,
41 1.67, 2.09, 2.61, 3.24, 4.01,
42 4.95, 6.08, 7.44, 9.08, 11.0,
43 12.1, 14.7, 17.6, 21.0, 25.0,
44 29.6, 34.8, 40.7, 47.3, 54.6,
45 62.7, 71.5, 81.0, 91.3, 102.
49 using namespace genie;
54 void testXsec(
const Char_t* outfn=
"VLExsecNT.root") {
55 gROOT->Macro(
"loadlibs.C");
68 pool->
FindRegistry(
"genie::StrumiaVissaniIBDPXSec/Default");
75 outf = TFile::Open(outfn,
"recreate");
76 nt =
new TNtuple(
"nt",
"nt",
"Ev:xsec:xsnun:xspaper");
88 for (i = 0; i <
nknots; i++) {
89 TLorentzVector
p4(0,0,
E[i]*1
e-3,
E[i]*1
e-3);
102 for (i = 0; i <
nknots; i++) {
103 TLorentzVector n4(0,0,
E[i]*1
e-3,
E[i]*1
e-3);
108 for (i = 0; i <
nknots; i++) {
109 Printf(
"xsec(E=%g MeV) = %g \t xsec_n = %g",
E[i],
Cross Section Calculation Interface.
void SetProbeP4(const TLorentzVector &P4)
THE MAIN GENIE PROJECT NAMESPACE
A singleton class holding all configuration registries built while parsing all loaded XML configurati...
An implementation of the neutrino - (free) nucleon [inverse beta decay] cross section, valid from the threshold energy (1.806MeV) up to hundreds of MeV. Currently cut off at 1/2 nucleon mass. Based on the Strumia/Vissani paper Phys.Lett.B564:42-54,2003.
Summary information for an interaction.
void testXsec(const Char_t *outfn="VLExsecNT.root")
A class encapsulating an enumeration of interaction types (EM, Weak-CC, Weak-NC) and scattering types...
virtual double Integral(const Interaction *i) const =0
A Neutrino Interaction Target. Is a transparent encapsulation of quite different physical systems suc...
virtual void Configure(const Registry &config)
A registry. Provides the container for algorithm configuration parameters.
void SetHitNucPdg(int pdgc)
Target * TgtPtr(void) const
InitialState * InitStatePtr(void) const
Registry * FindRegistry(string key) const
static AlgConfigPool * Instance()
Initial State information.