demo2.py

Go to the documentation of this file.

# Create an oscillatable spectrum

import cafana

import ROOT

fname = 'prod_sumdecaf_R17-03-01-prod3reco.k_fd_genie_nonswap_fhc_nova_v08_full_nue_or_numu_or_nus_contain_v1'

loader = cafana.SpectrumLoader(fname)

bins = cafana.Binning.Simple(100, 0, 1000)

# In fact kNHit was already defined within CAFAna. It's better to use that

nhit = cafana.Spectrum('Number of hits in slice', bins, loader, cafana.kNHit, cafana.kIsNumuCC)

# With true energy axis too
sosc = cafana.OscillatableSpectrum('Number of hits in slice', bins, loader, kNHit, cafana.kIsNumuCC)


loader.Go()

pot = 18e20

nhit.ToTH1(pot).Draw('hist')

# And 2D histograms
canv = ROOT.TCanvas()
sosc.ToTH2(pot).Draw('colz')