Skip to content

ICARUS NuMI QE-like cross section (arXiv:2604.24925)#112

Open
jedori0228 wants to merge 6 commits into
NUISANCEMC:mainfrom
jedori0228:feature/ICARUSNuMIXsec
Open

ICARUS NuMI QE-like cross section (arXiv:2604.24925)#112
jedori0228 wants to merge 6 commits into
NUISANCEMC:mainfrom
jedori0228:feature/ICARUSNuMIXsec

Conversation

@jedori0228
Copy link
Copy Markdown
Contributor

@jedori0228 jedori0228 commented May 13, 2026

Adding ICARUS NuMI QE-like (1muNp0pi) cross section measurement: https://arxiv.org/abs/2604.24925

Unique features of this measurement are:

  1. Signal definition contains both numu and numubar
    1. When a sample is generated with numu+numubar mixed flux, there is no issue. But in most cases, user produce them separately. To resolve this issue, I am saving flux/evt/Nevent for each input files separately in InputHandler, and re-calculate numu+numubar averaged values inside the measurement constructor. In the card, user must give N-number of numu files first in the input list, and then numubar. And also give N in the card with numubar_start_index key.
  2. Correlated fit: Multiple (two) 1D distributions fitted simultaneously, which results in correlated 1D cross sections.
    1. Overriding FillHistograms when correlated result is requested.
  3. Log-normal distribution of the cross section
    1. Support filling histograms and calculating chi2s in log(xsec)

I have tested that we can reproduce the data and MC results from the original paper with nuiscomp. Testing with nuismin will follow.

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.

1 participant