Skip to content

Conversation

@idkakorin
Copy link
Contributor

The DCCSPPPXSec and DCCSPPXSec 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.

Copy link
Member

@mroda88 mroda88 left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

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.

@idkakorin
Copy link
Contributor Author

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.

@idkakorin idkakorin requested a review from mroda88 June 15, 2023 12:45
Copy link
Contributor Author

@idkakorin idkakorin left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Address all requests

Copy link
Contributor Author

@idkakorin idkakorin left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Address all requests

@idkakorin idkakorin requested review from mroda88 and removed request for mroda88 February 14, 2024 20:30
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants