Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Dynamic coupled channels model by Toru Sato #292

Open
wants to merge 38 commits into
base: master
Choose a base branch
from

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.

src/Framework/Interaction/Interaction.cxx Outdated Show resolved Hide resolved
src/Framework/Interaction/Interaction.cxx Outdated Show resolved Hide resolved
@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