source: fact/tools/pyscripts/sandbox/dneise/plotsome.py

#!/usr/bin/python -tti
#
# Example
#    * looping over RawData class object
# 

from pyfact   import RawData
from plotters import Plotter
from drs_spikes import DRSSpikes
import sys

data_filename = '/home/dominik/py/test/20120304_018.fits.gz'
calib_filename = '/home/dominik/py/test/20120304_012.drs.fits.gz'

if len(sys.argv) < 2:
    print 'call: scriptname number-of-events-to-process'

run = RawData(data_filename, calib_filename, return_dict = True)

p = Plotter('test')

for event in run:
    index = event['event_id'].value
    print event['trigger_type'].value ,
    p(event['acal_data'][0])
    if index == int(sys.argv[1]):
        sys.exit(0)