import os,sys
import hippo
from numarray import *

app=hippo.HDApp()
canvas=app.canvas()

from hippo import NTupleController
ntc=NTupleController.instance()			#Instantiate an object
						
HistDataPath='/afs/slac/u/ki/lances/GEANT4/lsstgeant4/hist.dat'
HistDataNTup=ntc.createNTuple ( HistDataPath )  #Create an NTuple
HistDataArray=ntc.createDataArray(HistDataPath) #Create an Array
NumColumns=HistDataArray.columns
NumRows=HistDataArray.rows

x0=HistDataArray[0]
y0=HistDataArray[1]

#x0=HistDataNTup.getColumn('x0')
y0=HistDataNTup.getColumn('y0')
z0=HistDataNTup.getColumn('z0')
x1=HistDataNTup.getColumn('x1')
y1=HistDataNTup.getColumn('y1')
z1=HistDataNTup.getColumn('z1')
E=HistDataNTup.getColumn('Energy')
P=HistDataNTup.getColumn('Momentum')
M=HistDataNTup.getColumn('Mass')
EventID=HistDataNTup.getColumn('EventID')
TrackID=HistDataNTup.getColumn('TrackID')
ParentID=HistDataNTup.getColumn('ParentID')
PDGCode=HistDataNTup.getColumn('PDGCode')

NumColumns=len(x0)

canvas.setPlotMatrix(2,3)
from hippo import Display
hist1 = Display ( "Histogram", HistDataNTup, ("x0", ) )
canvas.addDisplay(hist1)
hist1.addDataRep ('Histogram', HistDataNTup, ['y0'])