ldf2fits.py

#!/usr/local/bin/python
#
#                               Copyright 2002
#                                     by
#                        The Board of Trustees of the
#                     Leland Stanford Junior University.
#                            All rights reserved.
#

__facility__ = "Online"
__abstract__ = "Archive to FITS converter"
__author__   = "S. Tuvi <stuvi@SLAC.Stanford.edu> SLAC - GLAST LAT I&T/Online"
__date__     = "02/27/2003"
__version__  = "$Revision: 2.2 $"
__credits__  = "SLAC"

import LATTE.copyright_SLAC

import sys
import getopt

def __usage():
  print "\n%s Usage:" % (sys.argv[0])
  print "\t-h or --help for Help"
  print "\t-a or --archive to specify an input archive file (*.ldf)"
  print "\t-f or --fits to specify the output FITS file (append .gz to the filename for compression)"
  print "\t-r or --runreport to specify an input run report filename (ie. rcReport.out)"
  print "\t-v or --verbose to specify progress messages during the conversion"


if __name__ == '__main__':
  archiveFile=''
  fitsFile=''
  verbose=0
  runReportFile=''
  try:
    opts, args = getopt.getopt(sys.argv[1:], "ha:f:vr:",
                               ["help", "archive=", "fits=", "verbose", "runreport="])
  except getopt.GetoptError:
    # print help information and exit:
    __usage()
    sys.exit(2)
  for o, a in opts:
    if o in ("-h", "--help"):
      __usage()
      sys.exit()
    if o in ("-a", "--archive"):
      archiveFile = a
    if o in ("-f", "--fits"):
      fitsFile = a
    if o in ("-v", "--verbose"):
      verbose = 1
    if o in ("-r", "--runreport"):
      runReportFile = a

  if archiveFile == '':
    print "\nERROR: Missing archive file parameter"
    __usage()
    sys.exit(2)

  if fitsFile == '':
    print "\nERROR: Missing FITS file parameter"
    __usage()
    sys.exit(2)

  if runReportFile == '':
    print "\nERROR: Missing Run Report file parameter"
    __usage()
    sys.exit(2)

  import rcArchiver
  arch = rcArchiver.rcArchiver(fileName=archiveFile, mode=rcArchiver.MODE_READONLY)

  import rcReportGen
  reportGen = rcReportGen.rcReportGen(fileName=runReportFile)
  reportGen.initialize()
  reportGen.readReportXML()

  import rcFitsWriter
  fits = rcFitsWriter.rcFitsWriter(fileName=fitsFile, format=rcFitsWriter.FORMAT_BYTE, reportGen=reportGen)

  cnt=0
  while(True):
    dat = arch.read()
    if dat == None: break
    cnt+=1
    if verbose and (cnt % 100) == 0: print "Progress line:", cnt
    fits.write(dat)
  if verbose and (cnt % 100) != 0: print "Progress line:", cnt

  arch.close()
  fits.close()

  if verbose: print "Conversion completed successfully!"

---