demo1.py

Go to the documentation of this file.

# Make two simple spectrum plots

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)

kTrackLen = cafana.CVar('''
if(sr.trk.kalman.ntracks == 0) return 0.0f;
return float(sr.trk.kalman.tracks[0].len);
''')

# For the simplest cases, can just provide the variable name to return
kNHit = cafana.SimpleVar('slc.nhit')

# Spectrum to be filled from the loader
length = cafana.Spectrum('Track length (cm)', bins, loader, kTrackLen, cafana.kIsNumuCC)

# And a second one...
nhit = cafana.Spectrum("Number of hits in slice", bins, loader, kNHit, cafana.kIsNumuCC)

# Do it!
loader.Go()

# How to scale histograms
pot = 18e20

# We have histograms
length.ToTH1(pot).Draw('hist')
canv = ROOT.TCanvas()
nhit.ToTH1(pot).Draw('hist')