Skip to content

Unified QEL model with spectral function #485

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.

Sign up for GitHub

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

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

Unified QEL model with spectral function #485

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
e9b9fcd
Unified QEL model with spectral function
LiangLiu212 May 3, 2026
e6e40c3
remove the test function
LiangLiu212 May 3, 2026
815ba9a
Updated the Unified QE cross section
LiangLiu212 May 7, 2026
b183a3b
Remove Eigen, create a TensorUtil to do Matrix
LiangLiu212 May 7, 2026
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
89 changes: 89 additions & 0 deletions config/KellyELFormFactorsModel.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,89 @@
<?xml version="1.0" encoding="ISO-8859-1"?>

<alg_conf>

<!--
Configuration for the Kelly elastic form factors model

Form factors based on J.J. Kelly, Phys.Rev.C 70 (2004) 068202

G(Q^2) \propto \frac{\sum_{k=0}^{n} a_k \tau^k}{1 + \sum_{k=1}^{n+2} b_k \tau^k}

Both numerator and denominator are polynomials in tau = Q^2 / 4mp^2
a and b are parameters for numerator and denominator.
This is just an ancient baseline for spectral function, do not use it for production.

Configurable Parameters:
.........................................................................................
Name Type Optional Comment Default
.........................................................................................
Kelly-Gep-a0 double No coefficient -
Kelly-Gep-a1 double No coefficient -
Kelly-Gep-a2 double No coefficient -
Kelly-Gep-b1 double No coefficient -
Kelly-Gep-b2 double No coefficient -
Kelly-Gep-b3 double No coefficient -
Kelly-Gep-b4 double No coefficient -
Kelly-Gmp-a0 double No coefficient -
Kelly-Gmp-a1 double No coefficient -
Kelly-Gmp-a2 double No coefficient -
Kelly-Gmp-b1 double No coefficient -
Kelly-Gmp-b2 double No coefficient -
Kelly-Gmp-b3 double No coefficient -
Kelly-Gmp-b4 double No coefficient -
Kelly-Gen-a0 double No coefficient -
Kelly-Gen-a1 double No coefficient -
Kelly-Gen-a2 double No coefficient -
Kelly-Gen-b1 double No coefficient -
Kelly-Gen-b2 double No coefficient -
Kelly-Gen-b3 double No coefficient -
Kelly-Gen-b4 double No coefficient -
Kelly-Gmn-a0 double No coefficient -
Kelly-Gmn-a1 double No coefficient -
Kelly-Gmn-a2 double No coefficient -
Kelly-Gmn-b1 double No coefficient -
Kelly-Gmn-b2 double No coefficient -
Kelly-Gmn-b3 double No coefficient -
Kelly-Gmn-b4 double No coefficient -
-->

<param_set name="Default">

<param type="string" name="CommonParam"> MagnMoments,QuasiElastic </param>
<!--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Kelly elastic form factor parameters
-->
<param type="double" name="Kelly-Gep-a0"> 1. </param>
<param type="double" name="Kelly-Gep-a1"> -0.24 </param>
<param type="double" name="Kelly-Gep-a2"> 0.0 </param>
<param type="double" name="Kelly-Gep-b1"> 10.98 </param>
<param type="double" name="Kelly-Gep-b2"> 12.82 </param>
<param type="double" name="Kelly-Gep-b3"> 21.97 </param>
<param type="double" name="Kelly-Gep-b4"> 0.0 </param>
<param type="double" name="Kelly-Gmp-a0"> 1. </param>
<param type="double" name="Kelly-Gmp-a1"> 0.12 </param> <!-- -->
<param type="double" name="Kelly-Gmp-a2"> 0.0 </param> <!-- -->
<param type="double" name="Kelly-Gmp-b1"> 10.97 </param> <!-- -->
<param type="double" name="Kelly-Gmp-b2"> 18.86 </param> <!-- -->
<param type="double" name="Kelly-Gmp-b3"> 6.55 </param> <!-- -->
<param type="double" name="Kelly-Gmp-b4"> 0.0 </param> <!-- -->
<param type="double" name="Kelly-Gen-a0"> 0.0 </param>
<param type="double" name="Kelly-Gen-a1"> 1.7 </param>
<param type="double" name="Kelly-Gen-a2"> 0.0 </param>
<param type="double" name="Kelly-Gen-b1"> 3.3 </param> <!-- -->
<param type="double" name="Kelly-Gen-b2"> 0.0 </param> <!-- -->
<param type="double" name="Kelly-Gen-b3"> 0.0 </param> <!-- -->
<param type="double" name="Kelly-Gen-b4"> 0.0 </param> <!-- -->
<param type="double" name="Kelly-Gmn-a0"> 1. </param>
<param type="double" name="Kelly-Gmn-a1"> 2.33 </param> <!-- -->
<param type="double" name="Kelly-Gmn-a2"> 0.0 </param> <!-- -->
<param type="double" name="Kelly-Gmn-b1"> 14.72 </param> <!-- -->
<param type="double" name="Kelly-Gmn-b2"> 24.20 </param> <!-- -->
<param type="double" name="Kelly-Gmn-b3"> 84.1 </param> <!-- -->
<param type="double" name="Kelly-Gmn-b4"> 0.0 </param> <!-- -->

</param_set>

</alg_conf>

33 changes: 33 additions & 0 deletions config/LwlynSmithFFEM.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,33 @@
<?xml version="1.0" encoding="ISO-8859-1"?>

<alg_conf>

<!--
Configuration sets for Llewellyn-Smith QEL EM form factors

Configurable Parameters:
..........................................................................................................
Name Type Optional Comment Default
..........................................................................................................
ElasticFormFactorsModel alg No Elastic form factors model CommonParam[ElasticFF]
UseElFFTransverseEnhancement bool No CommonParam[ElasticFF]
TransverseEnhancement alg No TransverseEnhancement model CommonParam[ElasticFF]
AxialFormFactorModel alg No Axial form factor model
AnomMagnMoment-P double No CommonParam[MagnMoments]
AnomMagnMoment-N double No CommonParam[MagnMoments]
WeinbergAngle double No Weinberg angle CommonParam[WeakInt]
SU3-D double No CommonParam[StrongInt]
SU3-F double No CommonParam[StrongInt]
-->

<param_set name="Default">

<param type="string" name="CommonParam"> WeakInt,MagnMoments,StrongInt,ElasticFF </param>

<!-- Dummy axial form factor model required since this class inherits from LwlynSmithFF.
All axial form factors are set to zero by LwlynSmithFFEM. -->
<param type="alg" name="AxialFormFactorModel"> genie::DipoleAxialFormFactorModel/Default </param>

</param_set>

</alg_conf>
3 changes: 3 additions & 0 deletions config/Messenger.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -112,6 +112,7 @@
<priority msgstream="GXSecFunc"> WARN </priority>
<priority msgstream="HadTransp"> WARN </priority>
<priority msgstream="HadronicVtx"> WARN </priority>
<priority msgstream="HadronTensors"> WARN </priority>
<priority msgstream="HybridHad"> WARN </priority>
<priority msgstream="IBD"> NOTICE </priority>
<priority msgstream="IMDXSec"> NOTICE </priority>
Expand Down Expand Up @@ -146,6 +147,7 @@
<priority msgstream="MArunAxialFormFactorModel"> NOTICE </priority>
<priority msgstream="KPhaseSpace"> NOTICE </priority>
<priority msgstream="LeptonicVertex"> WARN </priority>
<priority msgstream="LeptonTensor"> WARN </priority>
<priority msgstream="LocalFGM"> WARN </priority>
<priority msgstream="LwlynSmith"> NOTICE </priority>
<priority msgstream="LHAPDF6"> NOTICE </priority>
Expand Down Expand Up @@ -178,6 +180,7 @@
<priority msgstream="PathL"> NOTICE </priority>
<priority msgstream="PhysUtil"> WARN </priority>
<priority msgstream="QELXSec"> NOTICE </priority>
<priority msgstream="UnifiedQE"> WARN </priority>
<priority msgstream="QELCharmXSec"> NOTICE </priority>
<priority msgstream="QELEvent"> NOTICE </priority>
<priority msgstream="QELHadronicVtx"> WARN </priority>
Expand Down
3 changes: 3 additions & 0 deletions config/Messenger_inuke_verbose.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -91,6 +91,7 @@
<priority msgstream="GLRES"> WARN </priority>
<priority msgstream="HadTransp"> INFO </priority>
<priority msgstream="HadronicVtx"> NOTICE </priority>
<priority msgstream="HadronTensors"> WARN </priority>
<priority msgstream="HybridHad"> INFO </priority>
<priority msgstream="IBD"> NOTICE </priority>
<priority msgstream="IMDXSec"> NOTICE </priority>
Expand Down Expand Up @@ -156,6 +157,7 @@
<priority msgstream="QELEvent"> INFO </priority>
<priority msgstream="QELHadronicVtx"> INFO </priority>
<priority msgstream="QELKinematics"> INFO </priority>
<priority msgstream="UnifiedQE"> WARN </priority>
<priority msgstream="SPPEventGen"> INFO </priority>
<priority msgstream="RESXSec"> INFO </priority>
<priority msgstream="RESPXSec"> WARN </priority>
Expand Down Expand Up @@ -203,6 +205,7 @@
<priority msgstream="HEDISKinematics"> WARN </priority>
<priority msgstream="HEDISStrucFunc"> WARN </priority>
<priority msgstream="LeptoHad"> ERROR </priority>
<priority msgstream="LeptonTensor"> WARN </priority>

<priority msgstream="gevgen"> NOTICE </priority>
<priority msgstream="gevgen_ndcy"> NOTICE </priority>
Expand Down
3 changes: 3 additions & 0 deletions config/Messenger_laconic.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -107,6 +107,7 @@
<priority msgstream="GIBUU"> WARN </priority>
<priority msgstream="HadTransp"> WARN </priority>
<priority msgstream="HadronicVtx"> WARN </priority>
<priority msgstream="HadronTensors"> WARN </priority>
<priority msgstream="HybridHad"> WARN </priority>
<priority msgstream="IBD"> WARN </priority>
<priority msgstream="IMDXSec"> WARN </priority>
Expand Down Expand Up @@ -177,6 +178,7 @@
<priority msgstream="QELEvent"> WARN </priority>
<priority msgstream="QELHadronicVtx"> WARN </priority>
<priority msgstream="QELKinematics"> WARN </priority>
<priority msgstream="UnifiedQE"> WARN </priority>
<priority msgstream="SPPEventGen"> WARN </priority>
<priority msgstream="RESXSec"> WARN </priority>
<priority msgstream="RESPXSec"> WARN </priority>
Expand Down Expand Up @@ -224,6 +226,7 @@
<priority msgstream="HEDISKinematics"> WARN </priority>
<priority msgstream="HEDISStrucFunc"> WARN </priority>
<priority msgstream="LeptoHad"> WARN </priority>
<priority msgstream="LeptonTensor"> WARN </priority>

<priority msgstream="gevgen"> WARN </priority>
<priority msgstream="gevgen_dm"> WARN </priority>
Expand Down
3 changes: 3 additions & 0 deletions config/Messenger_rambling.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -107,6 +107,7 @@
<priority msgstream="GLRES"> WARN </priority>
<priority msgstream="HadTransp"> INFO </priority>
<priority msgstream="HadronicVtx"> NOTICE </priority>
<priority msgstream="HadronTensors"> WARN </priority>
<priority msgstream="HybridHad"> INFO </priority>
<priority msgstream="IBD"> NOTICE </priority>
<priority msgstream="IMDXSec"> NOTICE </priority>
Expand Down Expand Up @@ -175,6 +176,7 @@
<priority msgstream="QELEvent"> INFO </priority>
<priority msgstream="QELHadronicVtx"> INFO </priority>
<priority msgstream="QELKinematics"> INFO </priority>
<priority msgstream="UnifiedQE"> WARN </priority>
<priority msgstream="SPPEventGen"> INFO </priority>
<priority msgstream="RESXSec"> INFO </priority>
<priority msgstream="RESPXSec"> WARN </priority>
Expand Down Expand Up @@ -220,6 +222,7 @@
<priority msgstream="HEDISKinematics"> WARN </priority>
<priority msgstream="HEDISStrucFunc"> WARN </priority>
<priority msgstream="LeptoHad"> ERROR </priority>
<priority msgstream="LeptonTensor"> WARN </priority>

<priority msgstream="gevgen"> NOTICE </priority>
<priority msgstream="gevgen_dm"> NOTICE </priority>
Expand Down
3 changes: 3 additions & 0 deletions config/Messenger_whisper.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -107,6 +107,7 @@
<priority msgstream="GLRES"> FATAL </priority>
<priority msgstream="HadTransp"> FATAL </priority>
<priority msgstream="HadronicVtx"> FATAL </priority>
<priority msgstream="HadronTensors"> FATAL </priority>
<priority msgstream="HybridHad"> FATAL </priority>
<priority msgstream="IBD"> FATAL </priority>
<priority msgstream="IMDXSec"> FATAL </priority>
Expand Down Expand Up @@ -174,6 +175,7 @@
<priority msgstream="QELEvent"> FATAL </priority>
<priority msgstream="QELHadronicVtx"> FATAL </priority>
<priority msgstream="QELKinematics"> FATAL </priority>
<priority msgstream="UnifiedQE"> FATAL </priority>
<priority msgstream="SPPEventGen"> FATAL </priority>
<priority msgstream="RESXSec"> FATAL </priority>
<priority msgstream="RESPXSec"> FATAL </priority>
Expand Down Expand Up @@ -219,6 +221,7 @@
<priority msgstream="HEDISKinematics"> FATAL </priority>
<priority msgstream="HEDISStrucFunc"> FATAL </priority>
<priority msgstream="LeptoHad"> FATAL </priority>
<priority msgstream="LeptonTensor"> FATAL </priority>

<priority msgstream="gevgen"> FATAL </priority>
<priority msgstream="gevgen_dm"> FATAL </priority>
Expand Down
22 changes: 19 additions & 3 deletions config/SpectralFunc.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3,7 +3,7 @@
<alg_conf>

<!--
Configuration for the SpectralFunc NuclerModelI implementation
Configuration for the SpectralFunc NuclearModelI implementation

Configurable Parameters:
....................................................................................................................................
Expand All @@ -14,8 +14,24 @@ FermiMomentumTable string No Table of Fermi momentum (kF) constan
-->

<param_set name="Default">
<param type="string" name="CommonParam"> FermiGas </param>

<param type="string" name="CommonParam"> FermiGas </param>

<!-- The data file path may be specified relative to the GENIE root folder or
as an absolute path -->
<param type="string" name="DataPathType"> relative </param>
<param type="string" name="DataPath"> data/evgen/nucl/spectral_functions </param>

<!-- <param type="string" name="SpectFuncTable@Pdg=1000060120"> benhar-sf-12c.data </param> -->
<param type="string" name="SpectFuncTable@Pdg=1000060120_2212"> pke12_tot.data </param>
Comment thread
nusense marked this conversation as resolved.
<param type="string" name="SpectFuncTable@Pdg=1000060120_2112"> pke12_tot.data </param>
<param type="string" name="SpectFuncTable@Pdg=1000080160_2212"> pke16_tot.data </param>
<param type="string" name="SpectFuncTable@Pdg=1000080160_2112"> pke16_tot.data </param>
<param type="string" name="SpectFuncTable@Pdg=1000260560_2212"> benhar-sf-56fe.data </param>
<param type="string" name="SpectFuncTable@Pdg=1000260560_2112"> benhar-sf-56fe.data </param>
<param type="string" name="SpectFuncTable@Pdg=1000180400_2212"> pke40p_tot.data </param>
<param type="string" name="SpectFuncTable@Pdg=1000180400_2112"> pke40n_tot.data </param>

</param_set>

</alg_conf>

69 changes: 69 additions & 0 deletions config/UnifiedQELPXSec.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,69 @@
<?xml version="1.0" encoding="ISO-8859-1"?>

<alg_conf>

<!--
Configuration for the CBF spectral function QEL differential xsec algorithm.

Configurable Parameters:
.....................................................................................................
Name Type Optional Comment Default
.....................................................................................................
QEL-CC-XSecScale double Yes Scaling factor GPL value
CCFormFactorsAlg alg No QEL-CC form factors algorithm
NCFormFactorsAlg alg No QEL-NC form factors algorithm
EMFormFactorsAlg alg No QEL-EM form factors algorithm
XSec-Integrator alg No
IntegralNuclearModel alg
IntegralNucleonBindingMode string
TensorModel string No Type of hadron tensor Noemi-hadron-tensor
Options:
- Noemi-hadron-tensor
- Noemi-hadron-tensor-cc
PauliBlockerAlg alg No PauliBlocking algorithm
DoPauliBlocking bool Yes Do Pauli blocking true
DoRotate_qAlong_z bool Yes Rotates system so q is along z


.....................................................................................................
Parameters needed when Integrating with this model to generate splines:
...........................................................................................................................
Name Type Optional Comment Default
...........................................................................................................................
IntegralNuclearModel alg No
IntegralNucleonBindingMode alg No
-->

<param_set name="Default">

<param type="string" name="CommonParam"> CKM </param>

<param type="double" name="QEL-CC-XSecScale"> 1.000 </param>
<param type="double" name="IntegralNuclearInfluenceCutoffEnergy"> 100.0 </param>

<param type="alg" name="XSec-Integrator"> genie::NewQELXSec/Default </param>
<param type="alg" name="IntegralNuclearModel"> genie::SpectralFunc/Default </param>
<param type="string" name="IntegralNucleonBindingMode"> UseNuclearModel </param>

<param type="string" name="TensorModel"> Noemi-hadron-tensor </param>

@nusense nusense May 4, 2026

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

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

Looking at the code below that uses the "TensorModel" string, it looks like this can only currently take on one of two values. It would be good to document that above in the Configurable Parameters comment.

@LiangLiu212 LiangLiu212 May 7, 2026

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

I update the comments


<param type="alg" name="PauliBlockerAlg"> genie::PauliBlocker/Default </param>
<param type="bool" name="DoPauliBlocking"> true </param>

<param type="bool" name="DoRotate_qAlong_z"> false </param>

<param type="alg" name="NCFormFactorsAlg"> genie::LwlynSmithFFNC/Default </param>
<param type="alg" name="EMFormFactorsAlg"> genie::LwlynSmithFFEM/Default </param>

</param_set>

<param_set name="Dipole">
<param type="alg" name="CCFormFactorsAlg"> genie::LwlynSmithFFCC/Dipole </param>
</param_set>


<param_set name="ZExp">
<param type="alg" name="CCFormFactorsAlg"> genie::LwlynSmithFFCC/ZExp </param>
</param_set>

</alg_conf>
Loading