Skip to content

Harvard-Neutrino/atmospheric-neutrino-experiment-mc

Repository files navigation

Monte Carlo Simulations

We provide the MC simlation files of atmospheric neutrinos for the following experiments in HDF5 format. All of them can be downloaded from the Harvard Dataverse

  • SuperK: No neutron tagging covering the first 3 phases of the experiment.
  • SuperK with neutron tagging on Hydrogen: For phases 4 and 5 of SuperK and also used as HyperK's simulation file.
  • SuperK with neutron tagging on Gd: Assumes the SuperK detector is loaded with Gd at the goal concentration of 0.1 % .
  • IceCube - Upgrade: Copied from publibly available simulation from the IceCube Collaboration converted into .hdf5 format for completeness.
  • ORCA: Projected simulation based on IceCube - Upgrade's upgrade Monte Carlo applying reported detector response.

Plotting some variables

We provide a simple python code to plot variables from these Monte Carlo files with different option to get you started. The code uses pretty standard requirements which can be intalled via pip3:

  • h5py
  • matplotlib
  • numpy
pip3 install -r requirements.txt

Follow the instruction below to produce your first plots:

usage:

plot_sim.py [-h] --experiment [{SK,SK-Htag,SK-Gd,ORCA,ICUp,HK}] --fname [FNAME]
				[--variables VARIABLES [VARIABLES ...]] [--flavor [{e,mu,e+mu,tau}]] [--CP [{nu,antinu,both}]]
               [--interaction [{CC,NC,ALL,False}]] [--samples [SAMPLES]]

required arguments:

  --experiment [{SK,SK-Htag,SK-Gd,ORCA,ICUp,HK}] Experiment you would like to use.
  --fname [FNAME] Path to simulation file of the experiment you would like to use (SK, SK-Htag, SK-Gd, ORCA, ICUp or HK).

options:

  -h, --help show this help message and exit
  --variables VARIABLES [VARIABLES ...] Set of variables you want to plot. Deault variables are the true neutrino energy, the reconstructed cosine of zenith angle and the reconstructed neutrino energy.
  --flavor [{e,mu,e+mu,tau}] Flavor cut and breakdown. Default is break down of $\nu_{e}$ and $\nu_{\mu}$.
  --CP [{nu,antinu,both}] Flavor cut and breakdown. Default is break down of both.
  --interaction [{CC,NC,ALL,False}] Interaction mode(s) cut and breakdown. Dafult is no breakdown or cut.
  --samples [SAMPLES] Comma separated set of event samples you want to plot (for IC: tracks or cascades; for ORCA: tracks, intermediate or cascades; for SK and HK numerical indeces (displayed when calling these detectors)). Default is plotting all samples.

SuperK and HyperK files

Variable name(s) Description Name in file
neutrino,nu Neutrino flavor ipnu
Enu,E Neutrino energy (GeV) pnu
dirnu_x,nudir_x Direction of neutrino, x dirnuX
dirnu_y,nudir_y Direction of neutrino, y dirnuY
dirnu_z,nudir_z,cos_zen,cos_zentih Neutrino cosine zenith dirnuZ
azi,azimuth Neutrino azimuth angle azi
Plep,plepton,lepton_mom True lepton momentum (GeV/c) plep
dirlep_x,lepton_dir_x True lepton direction, x dirlepX
dirlep_y,lepton_dir_y True lepton direction, y dirlepY
dirlep_z,lepton_dir_z True lepton direction, z dirlepZ
interaction,int_mode,mode Interaction mode mode
inv_mass Invariant mass of pi0 imass
mom_reco_mer Reco. momentum of most energetic ring (GeV/c) pmax
reco_energy Reconstructed energy (GeV) evis
recodir_x,reco_dir_x Reconstructed direction, x recodirX
recodir_y,reco_dir_y Reconstructed direction, y recodirY
recodir_z,reco_dir_z,reco_coszen Reconstructed direction, z (zenith) recodirZ
ring_ip,reco_ring_ip,ip Reco. ring ID ip
number_of_rings,nring Number of reconstructed rings nring
muedk,reco_decay_e Number of tagged decay electrons muedk
number_neutrons,neutrons Number of tagged * neutrons neutron
itype,event_type Event ID (sample index) itype
------------------------------------------- --------------------------------------------- ------------
  • The number of tagged neutrons is only available in SK-Htag, SK-Gd and HK.

Note

There are more variables in these files that may be useful for the user, please look at their usage in kamioka.py.

IceCube-Upgrade files

Variable name(s) Description Name in file
neutrino,nu Neutrino flavor pdg
Enu,E Neutrino energy (GeV) true_energy
dirnu_z,nudir_z,cos_zen,cos_zentih Neutrino cosine zenith (cos_zenith derived from) true_zenith
azi,azimuth Neutrino azimuth angle true_azimuth
interaction,int_mode,mode Interaction mode: QE (0), RES (1), DIS (2), other (3) interaction_type
current,current_type Neutral (0) or charged (1) current interaction interaction_type
reco_energy Reconstructed energy (GeV) reco_energy
reco_azi,reco_azimuth Reconstructed azimuth reco_azimuth
recodir_z,reco_dir_z,reco_coszen Reconstructed direction, z (cos_zenith derived from) reco_zenith
itype,event_type,ip Event ID (sample index) pid
Q2,mom_transf Momentum transfer Q^2 (GeV^2/c^2) Q2
W,inv_mass_had Hadronic invariant mass (GeV/c^2) W
x_bjorken,x Bjorken x variable x
y,inelasticity Inelasticity of interaction y
cross_section,xsec Cross section xsec
diff_cross_section,dxsec Differential cross section dxsec
------------------------------------------- ---------------------------------------------------- ----------------

Note

There are more variables in these files that may be useful for the user, please look at their usage in south_pole.py.

ORCA file

The simulation file for ORCA follows the same notation and contains the same variables as the IceCube Upgrade one.

Other

For further questions and details, please check our joint analysis paper or reach out to the authors.

About

No description, website, or topics provided.

Resources

License

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages