Closed jtenavidal closed 8 months ago
No FSI case for DIS has no neutrons:
with FSI
From gst file check we can conclude genie only makes 1pi- events from en interactions. No ep events observed in gst file
No true neutrons
I think it is the 1pi (2 hadron) events where the AGKY has big issue with xF asymmetry. The more hadrons the more it is smeared, I guess
Fig 2 in https://arxiv.org/pdf/0904.4043.pdf
An easy test is to go to https://github.com/GENIE-MC/Generator/blob/master/src/Physics/Hadronization/AGKYLowW2019.cxx and comment out lines 919-940 and replace them with accepted=true
This will remove the biasing to match the input distributions and would just accept an isotropic phase space decay. This is probably as far back as you can go in the opposite direction of what we tried to do, so it would bracket the possible effect of mismodelling
The nucleon gets this xF distribution and the pion is pushed forward in the hadronic CM. We knew this is perhaps too strong and the event ought to be more symmetric (as similar final states from resonances). Could it be that this SIS mismodelling appears as an excess at low pmiss? Don't know. But seems weird that the suspect appears to be SIS/DIS... There is not much space for large mismodelling there, unless you have nfortunate sensitivity to some of the weird features of AGKY
if it has an effect here, the modification (and what is configurable) needs to be a bit subtler
It is easily possible to make this so that it is a gradual change, and not a step change between current <-> isotropic
Can we add a new knob? Maybe also useful for TKI
Running sample with:
python eAScatteringGridSubmitter.py --git-location https://github.com/e4nu/Generator_Q08.git --git-branch jtena/hadphasespace --jobs-topdir /pnfs/genie/scratch/users/jtenavid/GENIE_e4nu_Generations/G18_10a/Q2_08/haddecaytest/ --total-xsec /pnfs/genie/persistent/users/jtenavid/e4nu_files/GENIE_Files/Splines/G18_10a_00_000/G18_10a_00_000_eNuclei_Q2min_08_total_xsec.xml --config-dir /pnfs/genie/persistent/users/jtenavid/e4nu_files/e4nu-GENIE-config/ --ebeam-energy 4.461 --e-tgt-list 1000060120 --gst-output --no-ghep-output --starting-point 4 --job-lifetime-generation 6 --e-ntotevents 20000000 --store-comitinfo --tune G18_10a_00_000
./plot_e4nuanalysis --mc_location /Users/juliatenavidal/Desktop/Postdoc/e4nu/PionAnalysis/1p1pi/mc_files/ --input_mc_files e4nuanalysis_1p1pimanalysis_G18_10a_Q2_04_e_on_1000060120_2261MeV_NoRad,e4nuanalysis_1p1pimanalysis_G18_10b_Q2_04_e_on_1000060120_2261MeV_NoRad,e4nuanalysis_1p1pimanalysis_G18_10c_Q2_04_e_on_1000060120_2261MeV_NoRad --analysis_id 1p1pim --observable_list MissingMomentum --nofsi_file e4nuanalysis_1p1pimanalysis_G18_10a_NoFSI_Q2_04_e_on_1000060120_2261MeV_NoRad --output_location /Users/juliatenavidal/Desktop/Postdoc/e4nu/PionAnalysis/1p1pi/ --output_name "test"
Fig1: NoFSI available
Costas suggestion does not improve the agreement
I think the problem is only the magnitude. No need to re-run samples