IntTimeLst = [806.,242.,73.,22.,7.]	# In "GALPROP" units (1500 ???).

IntTimeLst = [539.]

ExpectToPoissTupleList = [
('ModelDataPaddedIntTime7.0.060529.fits',   'PoisDatonPaddedIntTime7.0.060604.fits'  ),
('ModelDataPaddedIntTime22.0.060529.fits',  'PoisDatonPaddedIntTime22.0.060604.fits' ),
('ModelDataPaddedIntTime73.0.060529.fits',  'PoisDatonPaddedIntTime73.0.060604.fits' ),
('ModelDataPaddedIntTime242.0.060529.fits', 'PoisDatonPaddedIntTime242.0.060604.fits'),
('ModelDataPaddedIntTime806.0.060529.fits', 'PoisDatonPaddedIntTime806.0.060604.fits')
]

ExpectToPoissTupleList = [('ModelDataPaddedIntTime539.0.060529.fits','PoisDatonPaddedIntTime539.X.060612.fits')]

ExpectToPoissTupleList = [
    ('counts_g1-9_g004gt1GeVpadded64x64b.fits','poiss_g1-9_g004gt1GeVpadded64x64b.fits')
    ]

ExpectToPoissTupleList = [
    ('Nu3CompDataPaddedIntTime133.0.070905.fits','Nu3CompPoisDatonPaddedIntTime133.0.070905.fits')
    ]

#from numarray import * ## Changed to numpy Aug 2007
from numpy import *
import pyfits
#import random
import numpy
from run_simpler_skymap_datons import *
from simpler_skymap_datons import *

for this_expecteddata_tuple in ExpectToPoissTupleList:
    this_modeldata_fitsfile, this_OutPoisDatonsFitsFile = this_expecteddata_tuple
    this_modeldata_HDU = pyfits.open(this_modeldata_fitsfile)
#   "Touch" header and data to make sure they are read in:
    print 'Input expected data cnts at 012:', this_modeldata_HDU[0].data[0][1][2]
    print 'Input expected data cnts header:', this_modeldata_HDU[0].header
#
# 4/ Get some Poisson counts! what fun!!
#
    this_poisdatons_HDU = PoisDatonsHDU(this_modeldata_HDU[0],this_OutPoisDatonsFitsFile)