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.

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

Jump to bottom

Noah's Fortran interface to Noemi's QE spectral function model #338

Draft
wants to merge 3 commits into
base: master
Choose a base branch
from
Draft

Noah's Fortran interface to Noemi's QE spectral function model #338

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
3de2ccd
Merge Noah's Fortran interface to Noemi's spectral function treatment of
sjgardiner Oct 18, 2023
7a5dd5f
Added abstract interface class HadronTensorInterfaceI and concrete im…
nastein Oct 26, 2023
7689a6f
Changed name of config parameter which selects hadron tensor models
nastein Oct 31, 2023
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
6 changes: 6 additions & 0 deletions Makefile
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -97,6 +97,7 @@ physics-neutrino-scattering-modes: FORCE
cd ${GENIE}/src/Physics/DeepInelastic/EventGen && $(MAKE) && \
cd ${GENIE}/src/Physics/Diffractive/XSection && $(MAKE) && \
cd ${GENIE}/src/Physics/Diffractive/EventGen && $(MAKE) && \
cd ${GENIE}/src/Physics/Fortran/XSection && $(MAKE) && \
cd ${GENIE}/src/Physics/HELepton/XSection && $(MAKE) && \
cd ${GENIE}/src/Physics/HELepton/EventGen && $(MAKE) && \
cd ${GENIE}/src/Physics/InverseBetaDecay/XSection && $(MAKE) && \
Expand Down Expand Up @@ -393,6 +394,7 @@ endif
mkdir ${GENIE_INC_INSTALLATION_PATH}/Physics/Diffractive
mkdir ${GENIE_INC_INSTALLATION_PATH}/Physics/Diffractive/XSection
mkdir ${GENIE_INC_INSTALLATION_PATH}/Physics/Diffractive/EventGen
mkdir ${GENIE_INC_INSTALLATION_PATH}/Physics/Fortran/XSection
mkdir ${GENIE_INC_INSTALLATION_PATH}/Physics/HELepton
mkdir ${GENIE_INC_INSTALLATION_PATH}/Physics/HELepton/XSection
mkdir ${GENIE_INC_INSTALLATION_PATH}/Physics/HELepton/EventGen
Expand Down Expand Up @@ -465,6 +467,7 @@ copy-install-files: FORCE
cd ${GENIE}/src/Physics/DeepInelastic/EventGen && $(MAKE) install && \
cd ${GENIE}/src/Physics/Diffractive/XSection && $(MAKE) install && \
cd ${GENIE}/src/Physics/Diffractive/EventGen && $(MAKE) install && \
cd ${GENIE}/src/Physics/Fortran/XSection && $(MAKE) install && \
cd ${GENIE}/src/Physics/HELepton/XSection && $(MAKE) install && \
cd ${GENIE}/src/Physics/HELepton/EventGen && $(MAKE) install && \
cd ${GENIE}/src/Physics/Hadronization && $(MAKE) install && \
Expand Down Expand Up @@ -527,6 +530,7 @@ purge: FORCE
cd ${GENIE}/src/Physics/DeepInelastic/EventGen && $(MAKE) purge && \
cd ${GENIE}/src/Physics/Diffractive/XSection && $(MAKE) purge && \
cd ${GENIE}/src/Physics/Diffractive/EventGen && $(MAKE) purge && \
cd ${GENIE}/src/Physics/Fortran/XSection && $(MAKE) purge && \
cd ${GENIE}/src/Physics/HELepton/XSection && $(MAKE) purge && \
cd ${GENIE}/src/Physics/HELepton/EventGen && $(MAKE) purge && \
cd ${GENIE}/src/Physics/Hadronization && $(MAKE) purge && \
Expand Down Expand Up @@ -590,6 +594,7 @@ clean-files: FORCE
cd ${GENIE}/src/Physics/DeepInelastic/EventGen && $(MAKE) clean && \
cd ${GENIE}/src/Physics/Diffractive/XSection && $(MAKE) clean && \
cd ${GENIE}/src/Physics/Diffractive/EventGen && $(MAKE) clean && \
cd ${GENIE}/src/Physics/Fortran/XSection && $(MAKE) clean && \
cd ${GENIE}/src/Physics/HELepton/XSection && $(MAKE) clean && \
cd ${GENIE}/src/Physics/HELepton/EventGen && $(MAKE) clean && \
cd ${GENIE}/src/Physics/Hadronization && $(MAKE) clean && \
Expand Down Expand Up @@ -678,6 +683,7 @@ endif
cd ${GENIE}/src/Physics/DeepInelastic/EventGen && $(MAKE) distclean && \
cd ${GENIE}/src/Physics/Diffractive/XSection && $(MAKE) distclean && \
cd ${GENIE}/src/Physics/Diffractive/EventGen && $(MAKE) distclean && \
cd ${GENIE}/src/Physics/Fortran/XSection && $(MAKE) distclean && \
cd ${GENIE}/src/Physics/HELepton/XSection && $(MAKE) distclean && \
cd ${GENIE}/src/Physics/HELepton/EventGen && $(MAKE) distclean && \
cd ${GENIE}/src/Physics/Hadronization && $(MAKE) distclean && \
Expand Down
1 change: 1 addition & 0 deletions config/CommonParam.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -236,6 +236,7 @@ Or changing the name of this parameter set
Elastic form factors used for QEL CC cross section calculation.
Options are:
- genie::DipoleELFormFactorsModel
- genie::KellyELFormFactorsModel, J.J. Kelly, Phys.Rev.C 70 (2004) 068202
- genie::BBA03ELFormFactorsModel, H.Budd, NuINT-02 proceedings
- genie::BBA05ELFormFactorsModel, R.Bradford, A.Bodek, H.Budd and J.Arrington, hep-ex/0602017
- genie::BBA07ELFormFactorsModel, R.Bradford, A.Bodek, H.Budd and S.Avvakumov, Euro.Phys.J.C53 (2008);[arXiv:0708.1946 [hep-ex]]
Expand Down
332 changes: 332 additions & 0 deletions config/GTEST23_40a/CommonParam.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,332 @@
<?xml version="1.0" encoding="ISO-8859-1"?>

<!--
***************************************************************************************************
This file defines parameters sets for parameters which are common between algorithms and therefore
has to be the same for those parameters to be meaningful.
This is file is read from the AlgConfigPool and the parameters set are expected to be used by the
algorithms to configure these common parameters.

Author:
Costas Andreopoulos <costas.andreopoulos \at stfc.ac.uk>
STFC, Rutherford Appleton Laboratory

and
Marco Roda <marco.roda \at liverpool.ac.uk>
University of Liverpool
***************************************************************************************************
-->



<common_Param_list>


<param_set name="Tunable">
<!-- Reserved register for tuning -->

<!-- <param type="double" name="DIS-XSecScale"> 1.032 </param> -->
<!-- <param type="double" name="MEC-CC-XSecScale"> 1.000 </param> -->


</param_set>


<param_set name="WeakInt">
<param type="double" name="WeinbergAngle"> 0.501716712132 </param> <!-- 0.501568 -->

</param_set>

<param_set name="StrongInt">
<param type="double" name="SU3-D"> 0.804 </param>
<param type="double" name="SU3-F"> 0.463 </param>

</param_set>


<param_set name="CKM">
<!-- ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
CKM quark mixing parameters
Chinese Physics C Vol.40, No. 10 (2016) 100001
Review of Particle Physics
-->
<param type="double" name="CKM-Vud"> 0.97417 </param>
<param type="double" name="CKM-Vus"> 0.2248 </param>
<param type="double" name="CKM-Vcd"> 0.220 </param> <!-- 0.2239 -->
<param type="double" name="CKM-Vcs"> 0.995 </param>

<param type="double" name="CabibboAngle"> 0.227780466682 </param> <!-- 0.226893 -->

</param_set>

<param_set name="NUCL">

<!-- ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
NUCL-Ro (in fm) is a scale parameter driving the effective nuclear sizes (Ro in R=Ro*A^1/3)
-->
<param type="double" name="NUCL-R0"> 1.4 </param>
<param type="double" name="NUCL-NR"> 3.0 </param>

</param_set>

<param_set name="Masses">
<param type="double" name="Charm-Mass"> 1.430 </param> <!-- GeV -->
</param_set>


<param_set name="NonResBackground">
<!--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
RES/DIS joining scheme
At this point there is a single scheme, the one originally developed in neugen3.
For W > Wcut : RES -> 0, + DIS -> full
For W <= Wcut : RES -> full + `DIS' (non-RES bkg) -> modified by DIS-HMultWgt-* params
-
- ...

-->
<param type="bool" name="UseDRJoinScheme"> true </param>
<param type="double" name="Wcut"> 1.7 </param>

<!--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
NEUGEN parameters applied to DIS hadronic multiplicity distributions
to avoid 2-ble counting with RES
-->
<param type="double" name="DIS-HMultWgt-vp-CC-m2"> 0.100 </param>
<param type="double" name="DIS-HMultWgt-vp-CC-m3"> 1.000 </param>
<param type="double" name="DIS-HMultWgt-vp-NC-m2"> 0.100 </param>
<param type="double" name="DIS-HMultWgt-vp-NC-m3"> 1.000 </param>
<param type="double" name="DIS-HMultWgt-vn-CC-m2"> 0.300 </param>
<param type="double" name="DIS-HMultWgt-vn-CC-m3"> 1.000 </param>
<param type="double" name="DIS-HMultWgt-vn-NC-m2"> 0.300 </param>
<param type="double" name="DIS-HMultWgt-vn-NC-m3"> 1.000 </param>
<param type="double" name="DIS-HMultWgt-vbp-CC-m2"> 0.300 </param>
<param type="double" name="DIS-HMultWgt-vbp-CC-m3"> 1.000 </param>
<param type="double" name="DIS-HMultWgt-vbp-NC-m2"> 0.300 </param>
<param type="double" name="DIS-HMultWgt-vbp-NC-m3"> 1.000 </param>
<param type="double" name="DIS-HMultWgt-vbn-CC-m2"> 0.100 </param>
<param type="double" name="DIS-HMultWgt-vbn-CC-m3"> 1.000 </param>
<param type="double" name="DIS-HMultWgt-vbn-NC-m2"> 0.100 </param>
<param type="double" name="DIS-HMultWgt-vbn-NC-m3"> 1.000 </param>

</param_set>


<param_set name="FermiGas">
<!--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Parameters related with GENIE's Fermi Gas (RFG) model implementation
- RFG-Momentum-CutOff is a momentum cut-off for the NN correlation tail
- RFG-NucRemovalE@Pdg=10LZZZAAAI is the removal energy for the nucleus with the specified pdg
code. If none is used then the average binding energy will be computed from Wapstra's semi-
empirical formula.
Currently, if you explicitly specify a binding energy for a nucleus then the same value will
be used for all isotopes.
-->
<param type="double" name="RFG-MomentumCutOff"> 0.5000 </param>
<param type="double" name="RFG-NucRemovalE@Pdg=1000030060"> 0.0170 </param> <!-- Li6 -->
<param type="double" name="RFG-NucRemovalE@Pdg=1000060120"> 0.0250 </param> <!-- C12 -->
<param type="double" name="RFG-NucRemovalE@Pdg=1000080160"> 0.0270 </param> <!-- O16 -->
<param type="double" name="RFG-NucRemovalE@Pdg=1000120240"> 0.0320 </param> <!-- Mg24 -->
<param type="double" name="RFG-NucRemovalE@Pdg=1000180400"> 0.0295 </param> <!-- Ar40 -->
<param type="double" name="RFG-NucRemovalE@Pdg=1000200400"> 0.0280 </param> <!-- Ca40 -->
<param type="double" name="RFG-NucRemovalE@Pdg=1000220480"> 0.0300 </param> <!-- Ti48 -->
<param type="double" name="RFG-NucRemovalE@Pdg=1000260560"> 0.0360 </param> <!-- Fe56 -->
<param type="double" name="RFG-NucRemovalE@Pdg=1000280580"> 0.0360 </param> <!-- Ni58 -->
<param type="double" name="RFG-NucRemovalE@Pdg=1000822080"> 0.0440 </param> <!-- Pb208 -->
<param type="bool" name="RFG-UseParametrization"> false </param>

<!--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Table of Fermi momentum (kF) constants for various nuclei
The tables can be found in $GENIE/config/FermiMomentumTables.xml
-->
<param type="string" name="FermiMomentumTable"> Default </param>

</param_set>


<!--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
INTRANUKE-specific parameters:
- Mode options are: hA, hN
- NucRemovalE is the binding E to subtract from cascade nucleons (in GeV)
- typical values for pion, nucleon DelR are 0.5, 0.7 (hA) and 0.2, 0.2 (hN)

Probably these should all be moved in the local configurations
-->

<param_set name="INUKE">
<param type="double" name="INUKE-NucRemovalE"> 0.00 </param>
<param type="double" name="INUKE-HadStep"> 0.05 </param>
<param type="double" name="INUKE-NucAbsFac"> 1.0 </param>
<param type="double" name="INUKE-NucQEFac"> 1.0 </param>
<param type="double" name="INUKE-NucCEXFac"> 1.0 </param>
<param type="double" name="INUKE-Energy_Pre_Eq"> 0.041 </param>
<param type="double" name="INUKE-FermiFac"> 1.0 </param>
<param type="double" name="INUKE-FermiMomentum"> 0.250 </param>
<param type="bool" name="INUKE-DoFermi"> true </param>
<param type="bool" name="INUKE-DoCompoundNucleus"> true </param>
<param type="bool" name="INUKE-XsecNNCorr"> true </param>
</param_set>

<param_set name="HAINUKE">
<param type="double" name="HAINUKE-DelRPion"> 0.5 </param>
<param type="double" name="HAINUKE-DelRNucleon"> 0.7 </param>
</param_set>

<param_set name="HNINUKE">
<param type="double" name="HNINUKE-DelRPion"> 0.0 </param>
<param type="double" name="HNINUKE-DelRNucleon"> 0.0 </param>
<param type="bool" name="HNINUKE-UseOset"> true </param>
</param_set>

<param_set name="Coherent">

<param type="double" name="COH-Ro"> 1.000 </param>

<!--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Minimum and maximum considered Q^2 for Berger-Sehgal coherent reactions when estimating
the max cross section.
Units in GeV^2.
-->
<param type="double" name="COH-Q2-min"> 0.000 </param>
<param type="double" name="COH-Q2-max"> 1.000 </param>

<!--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Maximum considered t for Berger-Sehgal coherent reactions when estimating the
max cross section.
Units in GeV^2.
-->
<param type="double" name="COH-t-max"> 0.250 </param>

</param_set>


<param_set name="Diffractive">

<!--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Slope parameter beta for diffractive scattering (GeV^-2) (b ~ 0.333 * nucleon_size^2)
-->
<param type="double" name="DFR-Beta"> 7.0 </param>

<!--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Maximum considered t for diffractive scattering reactions when estimating the max cross section.
Units in GeV^2.

This value is read in KpHase space!!! It's mental, but we have to be careful before removing this
Or changing the name of this parameter set
-->
<param type="double" name="DFR-t-max"> 0.35 </param>

</param_set>

<param_set name="QuasiElastic">
<!--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Value of axial form factor at Q2=0
-->
<param type="double" name="QEL-FA0"> -1.2694 </param>
<!--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Axial and vector masses for quasi-elastic scattering
-->
<param type="double" name="QEL-Ma"> 1.0 </param>
<param type="double" name="QEL-Mv"> 0.840 </param>

</param_set>

<param_set name="MagnMoments">
<!--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Proton and neutron anomalous magnetic moments
-->
<param type="double" name="AnomMagnMoment-P"> 2.79278 </param> <!-- 2.7928473 -->
<param type="double" name="AnomMagnMoment-N"> -1.91315 </param>

</param_set>

<param_set name="MultiNucleons">
<!-- Q3 max for 2p2h model -->
<param type="double" name="NSV-Q3Max"> 1.2 </param>
</param_set>

<param_set name="ElasticFF">
<!--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Elastic form factors used for QEL CC cross section calculation.
Options are:
- genie::DipoleELFormFactorsModel
- genie::KellyELFormFactorsModel, J.J. Kelly, Phys.Rev.C 70 (2004) 068202
- genie::BBA03ELFormFactorsModel, H.Budd, NuINT-02 proceedings
- genie::BBA05ELFormFactorsModel, R.Bradford, A.Bodek, H.Budd and J.Arrington, hep-ex/0602017
- genie::BBA07ELFormFactorsModel, R.Bradford, A.Bodek, H.Budd and S.Avvakumov, Euro.Phys.J.C53 (2008);[arXiv:0708.1946 [hep-ex]]
-->
<param type="alg" name="ElasticFormFactorsModel"> genie::KellyELFormFactorsModel/Default </param>

<!-- Option for turning on Transverse Enhancement by Elastic Form Factor adjustment.
See http://arxiv.org/abs/1405.0583 and http://arxiv.org/abs/1106.0340
-->
<param type="bool" name="UseElFFTransverseEnhancement"> false </param>
<param type="alg" name="TransverseEnhancement"> genie::TransverseEnhancementFFModel/Default </param>

</param_set>

<param_set name="Resonances">
<!--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Default list of baryon resonances included in cross section models and generation threads
-->
<param type="string" name="ResonanceNameList">
P33(1232),S11(1535),D13(1520),S11(1650),D13(1700),D15(1675),
S31(1620),D33(1700),P11(1440),P33(1600),P13(1720),F15(1680),
P31(1910),P33(1920),F35(1905),F37(1950),P11(1710),F17(1970)
</param>

</param_set>

<param_set name="Lepton">
<param type="bool" name="ApplyCoulombCorrection"> false </param>

</param_set>


<param_set name="BoostedDarkMatter">
<!--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Set the coupling of dark matter to nucleons
-->
<param type="double" name="ZpCoupling"> 0.1 </param>

</param_set>


<param_set name="KNO2Pythia">
<!--

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Invariant mass window for the transition from KNO model to PYTHIA
-->

<param type="double" name="KNO2PYTHIA-Wmin"> 2.30 </param>
<param type="double" name="KNO2PYTHIA-Wmax"> 3.00 </param>
</param_set>



<param_set name="Validation">
<!--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
GENIE validity range
-->
<param type="double" name="GVLD-Emin"> 0.010 </param>
<param type="double" name="GVLD-Emax"> 1000.000 </param>

</param_set>

</common_Param_list>
Loading