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.

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

add Mori2020 model #160

Open
JostMigenda opened this issue Dec 7, 2021 · 0 comments
Open

add Mori2020 model #160

JostMigenda opened this issue Dec 7, 2021 · 0 comments
Assignees
Labels
SupernovaModel Implementing/correcting supernova model

Comments

@JostMigenda
Copy link
Member

Mori-san recently contacted me about adding the Mori2020 model to SNEWPY. It uses the same format as the Supernova Neutrino Database. Previously, we’ve converted those files to a different format assuming a pinched spectrum (@nuberoi, is the script to do this available publicly?); but I would suggest implementing the spectral interpolation described by Nakazato for improved accuracy.[1] I’ve already implemented that for sntools, so this shouldn’t be too much extra work; but don’t know yet when I’ll find the time.

[1] Comparison of original interpolation with pinched spectrum:
Running sntools <input file/format> --detector HyperK --distance 50 --starttime 100 --endtime 300 for both nakazato-shen-z0.02-t_rev100ms-s20.0.fits (SNEWPY’s fit to a pinched spectrum) and intp2001.data (original file, sntools’ implementation of the spectrum interpolation described by Nakazato), I get the following number of expected events:

channel orig. pinched
es-e 7.84 7.88
es-eb 3.31 3.33
es-x 2.01 2.02
es-xb 1.73 1.73
ibd 223.49 223.21
o16e 2.59 0.84
o16eb 3.24 2.35

So the difference is negligible for IBD and elastic scattering, but up to a factor of 3× for oxygen-16 CC interactions; presumably because the difference between spectra is only significant in the high-energy tail, which O16 interactions are particularly sensitive too.

@JostMigenda JostMigenda added the SupernovaModel Implementing/correcting supernova model label Dec 7, 2021
@JostMigenda JostMigenda self-assigned this Dec 7, 2021
@JostMigenda JostMigenda added this to the v1.4 milestone Sep 15, 2022
@JostMigenda JostMigenda removed this from the v1.4 milestone Apr 2, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
SupernovaModel Implementing/correcting supernova model
Projects
None yet
Development

No branches or pull requests

1 participant