Add existing file

#!/usr/bin/env python3

# -*- coding: utf-8 -*-

"""

Created on Wed Jul 24 12:21:39 2019

@author: christian

"""

import gcn

import healpy as hp

from tables import *

import numpy as np

import ligo.skymap

import os

class Event(IsDescription):

    GraceID=StringCol(30)

    AlertType=StringCol(30)

    Instruments = StringCol(30)

    FAR = StringCol(30)

    skymap=StringCol(30)

    Group=StringCol(30)

    BNS=StringCol(30)

    NSBH=StringCol(30)

    BBH=StringCol(30)

    Terrestrial = StringCol(30)

    HasNS=StringCol(30)

    HasRemnant=StringCol(30)

# Function to call every time a GCN is received.

# Run only for notices of type

# LVC_PRELIMINARY, LVC_INITIAL, or LVC_UPDATE.

@gcn.handlers.include_notice_types(

    gcn.notice_types.LVC_PRELIMINARY,

    gcn.notice_types.LVC_INITIAL,

    gcn.notice_types.LVC_UPDATE,

    gcn.notice_types.LVC_RETRACTION)

def process_gcn(payload, root):

    # Respond only to 'test' events.

    # VERY IMPORTANT! Replace with the following code

    # to respond to only real 'observation' events.

    # if root.attrib['role'] != 'observation':

    #    return

    if root.attrib['role'] != 'test':

        return

    h5file = open_file("Event Database",mode="r+",title="eventinfo")

    # Read all of the VOEvent parameters from the "What" section.

    params = {elem.attrib['name']:

              elem.attrib['value']

              for elem in root.iterfind('.//Param')}

    # Respond only to 'CBC' events. Change 'CBC' to "Burst'

    # to respond to only unmodeled burst events.

    if params['Group'] != 'CBC':

        return

    filepath = params['skymap_fits']

    skymap = params['GraceDB']+'.png'

    os.sys("curl -0 "+filepath)

    os.sys("ligo-skymap-plot "+os.path.basename(filepath)+" -o "+skymap+" --annotate --contour 50 90 --geo")

    #group name: events

    table = h5file.create_table(h5file.root.events,params['GraceID'],Event,'CBC event')

    det_event = table.row

    # Print all parameters.

    for key, value in params.items():

        print(key, '=', value)

        try:

            det_event[key] = value

        except:

            continue

    det_event['skymap'] = skymap

 #   if 'skymap_fits' in params:

  #      # Read the HEALPix sky map and the FITS header.

    #    skymap, header = hp.read_map(params['skymap_fits'],

   #                                  h=True, verbose=False)

  #      header = dict(header)

#

        # Print some values from the FITS header.

    #    print('Distance =', header['DISTMEAN'], '+/-', header['DISTSTD'])

    table.flush()

    h5file.close()

gcn.listen(handler=process_gcn)

 No newline at end of file