code: def event_table(self):
          
          # Create event list table for each source
          tables = []
          for source in self.source_table:
              lon = source['GLON'] * np.ones(source['COUNTS'])
              lat = source['GLAT'] * np.ones(source['COUNTS'])
              coord = SkyCoord(lon, lat, unit='deg', frame='galactic')
  
              # TODO: scatter positions assuming Gaussian PSF on the sky
              # using SkyOffsetFrame.
  
              table = Table()
              table['GLON'] = lon
              table['GLAT'] = lat
              table['SOURCE_IDX'] = source.index
  
              tables.append(table)
  
          # Stack all tables together
          table = table_vstack(tables)
  
          return table
docs: Event table (`~astropy.table.Table`).
  
          Columns: GLON, GLAT, SOURCE_IDX