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
Open
Changes from 1 commit
Commits
Show all changes
38 commits
Select commit Hold shift + click to select a range
a89ca26
begining of work on DCC
Jun 27, 2022
7c179b7
rename files
Jun 27, 2022
c24bceb
continue change spp
Jun 27, 2022
547b6de
change spp channel definitions and methods
Jun 27, 2022
d765d44
revert sign of Clebsch–Gordan coefficients to wrong ones just for con…
Jun 28, 2022
ab106ef
add single pion scattering type
Aug 3, 2022
5c0e3f2
minute changes
Aug 4, 2022
4d015c5
save current work
Aug 15, 2022
73455c6
save progress
Aug 18, 2022
d07539a
save current work
Aug 24, 2022
03b7436
save current work
Aug 27, 2022
f2c668e
savr current work
Aug 28, 2022
baafe0a
save current work
Aug 29, 2022
89a5120
save current work
Aug 29, 2022
c979983
save current work
Aug 30, 2022
c1cd010
save current work
Aug 31, 2022
a836775
save current work
Sep 1, 2022
64bec6e
save current work
Sep 4, 2022
0decfc3
save current work
Sep 4, 2022
f599359
fix very small energy case
Oct 20, 2022
f87308d
use literal for imaginary unit
Oct 21, 2022
03aa8d1
change integrator parameters
Oct 21, 2022
c6572fb
change integration accuracy
Oct 22, 2022
b24e1f9
change formatting
Oct 23, 2022
f45a484
small changes
Oct 23, 2022
100b2bd
take into account the non-isoscalar nature of nucleons, but this is w…
Apr 9, 2023
9155faf
implement XSec DCC class, everything is compiled
Apr 22, 2023
d51c893
add some rearranges
Apr 23, 2023
87521fe
fix SPP kinematics
Apr 26, 2023
bfad77e
implement DXsec algorithm for DCC
Jun 8, 2023
4c9afa0
add dcc integrator
Jun 9, 2023
726f8de
add comment
Jun 10, 2023
045cb32
fix factor for massless lepton case
Jun 12, 2023
e65e972
fix mistake
Jun 12, 2023
3c9242f
Merge branch 'master' into DCC_Sato
idkakorin Jun 15, 2023
acfcb39
Merge branch 'master' into DCC_Sato
idkakorin Jul 7, 2023
6f830fc
Merge branch 'master' into DCC_Sato
idkakorin Feb 14, 2024
1fedf8f
sync with master
idkakorin Feb 14, 2024
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
Next Next commit
small changes
Igor Kakorin committed Oct 23, 2022
commit f45a4847a8f270fe7ad5bc1cb8515f3af3c4b1f3
2 changes: 1 addition & 1 deletion src/Physics/Resonance/XSection/DCCSPPPXSec.cxx
Original file line number Diff line number Diff line change
@@ -146,7 +146,7 @@ double DCCSPPPXSec::XSec(const Interaction * interaction, KinePhaseSpace_t kps)
double Gamma = kAem*qL_gamma*Elf_L/2./kPi2/Q2/Eli_L/(1. - epsilon);


std::complex<double> F1, F2,F3, F4, F7, F8;
std::complex<double> F1, F2, F3, F4, F7, F8;
for (unsigned int L = 0; L<5; L++)
{
VAmpl vampl = Amplitudes(W, Q2, L, spp_channel);
6 changes: 1 addition & 5 deletions src/Physics/Resonance/XSection/DCCSPPXSecWithCache.cxx
Original file line number Diff line number Diff line change
@@ -263,12 +263,8 @@ genie::utils::gsl::d3XSecSPP_dWQ2CosTheta_E::d3XSecSPP_dWQ2CosTheta_E(

Wl = kps->WLim_SPP();
// model restrictions
if (Wl.min < 1.08 && Wl.max < 1.08)
{
isZero = true;
return;
}
Wl.min = TMath::Max (Wl.min, 1.08);
Wl.max = TMath::Max (Wl.max, 1.08);
Wl.max = TMath::Min (Wl.max, 2.00);

if (fWcut >= Wl.min)