Open idkakorin opened 1 year ago
There is a small logic mistake.
I would also ask to add a markdown file in
data/evgen/dcc/
to describe the content of the file: units, what it is, link to papers/authors of the file etc. Whatever information we can about it.
All links and authors, can be found in the appropriate header file. The table with amplitudes are used "as is", provided by author. I am not author of this theory, I think all information can be found in the papers, cited in header file.
The
DCCSPPPXSec
andDCCSPPXSec
algorithm are implemented. They are working for neutrinos and charged leptons.Now, the following commands are valid:
gmkspl -p 11L -t 1000000010 -n 30 -e 100 -o out.xml --tune DCC22_00a_00_000
--- calculates cross section for massless left polarized electron on neutron;gmkspl -p 11R -t 1000000010 -n 30 -e 100 -o out.xml --tune DCC22_00a_00_000
--- calculates cross section for massless right polarized electron on neutron;gmkspl -p 11 -t 1000000010 -n 30 -e 100 -o out.xml --tune DCC22_00a_00_000
--- calculates cross section for massive electron on neutron;gmkspl -p 12 -t 1000000010 -n 30 -e 100 -o out.xml --tune DCC22_00a_00_000
--- calculates cross section for electron neutrino on neutron;The helicity of massless lepton is added to the class
genie::InitialState
.DCCSPPPXSec:XSec(const Interaction * i, KinePhaseSpace_t k)
calculates: $\frac{d^2\sigma}{dWdQ^2}$ when k = kPSWQ2fE; $\frac{d^3\sigma}{dWdQ^2d\cos\theta\pi}$ when k = kPSWQ2ctpfE; $\frac{d^4\sigma}{dWdQ^2d\cos\theta\pi d\phi\pi}$ when k = kPSWQ2ctpphipfE, where $\theta\pi$ and $\phi_\pi$ are polar and azimuthal angles of final pion in the resonance rest frame.The single pion production reaction is defined by enumeration in
SppChannel.h
which was extended for charged lepton case.