Feature/incl tune h potential - #13
Open
LiangLiu212 wants to merge 108 commits into
Open
Conversation
… 2p2h - 3p3h separation yet. No A-scaling
… make it consist with FIS
… make it consist with FIS
Implement of NucleusGenerator to replace the old FermiMover and VertexGenerator. NucleusGenerator is inherited from EventRecordVisitorI. This interface facilitates the use of various implementations of nuclear models. Now I create two implementations: NucleusGenINCL and NucleusGenTraditional. NucleusGenINCL nuclear model is used for both initial state and FSI model to have a consistent physics results. NucleusGenTraditional combines the FermiMover and VertexGenerator.
…ematics. Add in optional bias towards forward nucleon momentum.
modify the INCL nucleus to use the separation energy (M - E_separation) On branch inclxx Your branch is up to date with 'origin/inclxx'. Changes to be committed: new file: config/GINCL00_00a/EventGeneratorListAssembler.xml new file: config/GINCL00_00a/QELEventGenerator.xml modified: config/GINCL00_00a/TuneGeneratorList.xml modified: src/Apps/gNtpConv.cxx modified: src/Physics/NuclearState/INCLNucleus.cxx modified: src/Physics/NuclearState/NucleusGenINCL.cxx
it is the first implement the INCL in QEL-CC INCL nuclear model is different from nuclear model in GENIE and it is not work well with GENIE QEL-CC model. In this commit, I just use the vertex of hit nucleon from INCL. I use nuclear model in GENIE to simulate the fermi momentum for the hit nucleon.
The basic framework of QEL, RES, MEC and DIS has been completed. I need to adjust the detail of INCL nuclear model to setup the correct configuration.
Now I can run QEL, RES, MEC, and DIS.
channel with INCL. I think it will be better to create the new class to handle the CCQEL and CCMEC channel, so, we won't mess up with current models.
Refactor the QEL and MEC channel to using my new inferface of nuclear model. This interface will simulate the fermi motion and vertex simultaneously. So, I create a duplications of QEL/MEC EventGen to use the new interface without interference with old nuclear interface.
Add new QEL XSec integral and QEL Event Generator to using my new interface NucleusGenI. I use NucleusGenTraditional to test the refactor. The new interface works with old GENIE nuclear model and vertex generator.
Apply all of the physical consideration of hadronic binary collision in INCL nuclear model for neutrino vertices 1. Out of well: mark pions and kaons that have been created outside their well 2. Enforce energy conservation 3. Check the mass of outgoing delta resonances above the threshold 4. Pauli blocking 5. CDPP blocking 6. Bring particle inside the nucleus 7. Biasing of the final state 8. Make target spectator and decrement cascading; make participant and increment cascading
When using hybrid nuclear model with GENIE vertex and GENIE nuclear model, INCL will initialize the nuclear twice by mistake.
…/Multinucleon/XSection/MartiniEricsonChanfrayMarteauMECPXSec2024.cxx src/Physics/QuasiElastic/XSection/MartiniQELPXSec.cxx
Previous method: Generate a neutron from GENIE vertex generator, find the closest one in INCL and replace it. New method: Generate a neutron from GENIE vertex generator, randomly drop a neutron in INCL nucleus.
Nieves CCQE model uses pre-calculated hadron tensor in the rest frame of initial nucleon and q should be in z direction. Now we use coordinate system in the rest frame of nucleus. We need to do boost and rotation to do the transfermation.
Here I converted Noah's fortran code (calculated hadron tensor) to c++ version. Also add Kelly FF, and unified qel xsec model.
Fix a bug in Nieves CCQE hadron tensor
Setup the tune to turn off the delta transportation
C++ version of QEL model from Noah's Fortran code
Sync martini2p2h with upstream master before pull request.
1. remove test configurations 2. setup conditional compilation for INCL 3. setup configuration for alternative MEC model
Martini2p2h
1. Add EM MEC configuration in LFG and SF test tunes 2. Add configuration options for INCL cluster production
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
5 participants
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
No description provided.