#!/usr/bin/env python

"""
  options:

  -o -output_file :  
        The output file, root file
        default=auto_pwd 
  -E : energy           default=1e6
  -D : distance         default=100
  -C : cosine angle     default=0.8
  -P : pmt cosine angle default=-1.0
"""

print "plotPdf starting..."

import ROOT,os,sys

try :
    import aa
except ImportError:
    print "Failed to import aa; Make sure AADIR is in your PYHTONPATH!"
    sys.exit()
import rec

options = aa.Options( __doc__, sys.argv[1:] )
#options.read_env()
print options

from ROOT import vector, Hit, Det, Evt, Trk, Vec, EventFile, Timer
from ROOT import MestShowerPdf, DigitalShowerPdf, ShowerFit
from ROOT import JppShowerPdf, JppShowerFit
from ROOT import PlotPdf

#-----------------------------------------------------
# Fancy-smancy setting of members
#-----------------------------------------------------

#for T in Trk, JppShowerPdf, DigitalShowerPdf: T.set = aa.__set

# input & output file 
#-----------------------------------------

pdffile    = "./pdf_data/genhen_JUNE_2___.hist.root"
jshowerfile    = ["J13p.dat", "J14p.dat"]
estart     = 5e5
print " output file =", options.o

# plotting
#-----------------------------------------
aapdf = DigitalShowerPdf()
aapdf.init()

myplotter = PlotPdf();
myplotter.set_all(options.E, options.D, options.C, options.P);
myplotter.plot_cd(options.o+"_cd.root");
myplotter.plot_D(options.o+"_D.root");
myplotter.plot_pmta(options.o+"_pmta.root");

print "plotPdf done."