Skip to content

MEC dial updates#51

Open
sjgardiner wants to merge 27 commits into
masterfrom
feature/mec-dial-updates
Open

MEC dial updates#51
sjgardiner wants to merge 27 commits into
masterfrom
feature/mec-dial-updates

Conversation

@sjgardiner

Copy link
Copy Markdown
Member

Replacement for PR #49 that has been rebased onto master and cleaned up. Contributions originally developed by Lars Bathe-Peters (@larsb-p).

Updated list of 2p2h/MEC uncertainties handled by Reweight:

  • NormCCMEC
  • DecayAngMEC, DecayAng2MEC
  • DecayAngMECLegendre, DecayAngMECLegendre2, DecayAngMECLegendre3, DecayAngMECLegendre4, DecayAngMECLegendre5, DecayAngMECLegendre6
  • FracPN_CCMEC
  • FracDelta_CCMEC
  • XSecShape_CCMEC
  • XSecShape_CCMEC_Empirical
  • XSecShape_CCMEC_Martini
  • EnergyDependence_CCMEC

Input CCMEC models that are compatible with the updated dials include Valencia, SuSAv2, and Martini-Ericson-Chanfray-Marteau.

larsb-p added 20 commits July 15, 2026 10:57
…enabled actual reweighting to the 2p-2h Martini model by inputting MartiniEricsonChanfrayMarteauMECPXSec2024 into placeholder for XSecShape_CCMEC_Martini dial
…onsidering the 2p2h Martini model and added normalsation for DecayAngMEC parameter with two tweak dials
…cia model and the energy transfer cut-off in the Martini model
…ecShape_CCMEC and XSecShape_CCMEC_Martini dials
…Martini-Ericson_Chanfray-Marteau model is input as default
…ied instead of divided (input default model gets up- or down-reweighted according to renormalised cross section ratio being greater or smaller to 1.)
…ng hard-coded xsec ratio tables to it reading the xsecs directly from the user's input cross-section spline xml-file. The user can now input multiple xml-files via the cross-sections option of grwght1p.
…RwghtNCorrelatedParams.cxx (matches original implementation)
@sjgardiner
sjgardiner force-pushed the feature/mec-dial-updates branch from 0de7241 to 855182e Compare July 15, 2026 15:57
…s may

be supplied via the --cross-sections command-line option
special "Reweight" configuration of the empirical MEC xsec algorithm
@sjgardiner

Copy link
Copy Markdown
Member Author

This is mostly in good shape but still needs some attention. I note that the reweighting of Valencia CCMEC events that involve an internal Δ line vs. those that do not appears to be broken by these changes. I will have a look at restoring that functionality and double-checking everything else.

@sjgardiner
sjgardiner force-pushed the feature/mec-dial-updates branch from 5e0bbfa to 768e6ed Compare July 16, 2026 16:12
@sjgardiner
sjgardiner force-pushed the feature/mec-dial-updates branch from 768e6ed to d50764c Compare July 16, 2026 22:23
…ns to

the CC MEC cross section under the Valencia model. This functionality was
removed when implementing pn-fraction reweighting for the SuSAv2 and
Martini MEC models. It now lives in a separate member function of
GReWeightXSecMEC. For models other than Valencia, the new function simply
returns a unit weight.
@sjgardiner

Copy link
Copy Markdown
Member Author

The "Δ / no Δ" reweighting for the Valencia model is now fixed on some new commits. I moved the relevant code into a dedicated member function to make the implementation cleaner.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants