def from_file_hessio(filename, item): """ reads the Telescope data out of the open hessio file Parameters ---------- filename: string name of the hessio file (must be a hessio file!) """ tel_id = h.get_telescope_ids() tel_num = h.get_num_telescope() tel_posX = [-1] * u.m tel_posY = [-1] * u.m tel_posZ = [-1] * u.m return (tel_id, tel_num, tel_posX, tel_posY, tel_posZ)
def from_hessio_to_astropytable(filename): h.file_open(filename) next(h.move_to_next_event()) tel_id = h.get_telescope_ids() tel_table = Table() tel_table["TelID"] = tel_id tel_posX = [h.get_telescope_position(i)[0] for i in tel_id] * u.m tel_posY = [h.get_telescope_position(i)[1] for i in tel_id] * u.m tel_posZ = [h.get_telescope_position(i)[2] for i in tel_id] * u.m tel_table["TelX"] = tel_posX tel_table["TelX"].unit = u.m tel_table["TelY"] = tel_posY tel_table["TelY"].unit = u.m tel_table["TelZ"] = tel_posZ tel_table["TelZ"].unit = u.m for t in range(len(tel_id)): table = Table() pix_posX = h.get_pixel_position(tel_id[t])[0] * u.m pix_posY = h.get_pixel_position(tel_id[t])[1] * u.m pix_id = np.arange(len(pix_posX)) pix_area = h.get_pixel_area(tel_id[t]) * u.m ** 2 table["TelID"] = [tel_id[t] for i in range(len(pix_posX))] table["PixelID"] = pix_id table["PixX"] = pix_posX table["PixX"].unit = u.m table["PixY"] = pix_posY table["PixY"].unit = u.m table["PixArea"] = pix_area.value table["PixArea"].unit = pix_area.unit if t == 0: cam_table = table else: cam_table = join(cam_table, table, join_type="outer") return [tel_table, cam_table]
def from_hessio_to_astropytable(filename): h.file_open(filename) next(h.move_to_next_event()) tel_id = h.get_telescope_ids() tel_table = Table() tel_table['TelID'] = tel_id tel_posX = [h.get_telescope_position(i)[0] for i in tel_id] * u.m tel_posY = [h.get_telescope_position(i)[1] for i in tel_id] * u.m tel_posZ = [h.get_telescope_position(i)[2] for i in tel_id] * u.m tel_table['TelX'] = tel_posX tel_table['TelX'].unit = u.m tel_table['TelY'] = tel_posY tel_table['TelY'].unit = u.m tel_table['TelZ'] = tel_posZ tel_table['TelZ'].unit = u.m for t in range(len(tel_id)): table = Table() pix_posX = h.get_pixel_position(tel_id[t])[0] * u.m pix_posY = h.get_pixel_position(tel_id[t])[1] * u.m pix_id = np.arange(len(pix_posX)) pix_area = h.get_pixel_area(tel_id[t]) * u.m**2 table['TelID'] = [tel_id[t] for i in range(len(pix_posX))] table['PixelID'] = pix_id table['PixX'] = pix_posX table['PixX'].unit = u.m table['PixY'] = pix_posY table['PixY'].unit = u.m table['PixArea'] = pix_area.value table['PixArea'].unit = pix_area.unit if t == 0: cam_table = table else: cam_table = join(cam_table, table, join_type='outer') return [tel_table, cam_table]
def load_hessio(filename): """ Function to open and load a hessio file Parameters ---------- filename: string name of the file """ print("Hessio file will be opened.") h.file_open(filename) next(h.move_to_next_event()) #version = h.get... version = 'Feb2016' #Creating 3 dictionaries where the instrument configuration will be stored #The dictionaries themselves contain astropy.table.Table objects telescope = {} camera = {} optics = {} #-------------------------------------------------------------------------- #Telescope configuration tel_table_prime = Table() tel_table_prime.meta = {'VERSION': version} try: tel_id = h.get_telescope_ids() tel_table_prime['TelID']= tel_id except: pass try: tel_posX = [h.get_telescope_position(i)[0] for i in tel_id] tel_table_prime['TelX'] = tel_posX tel_table_prime['TelX'].unit = u.m tel_table_prime.meta['TelX_description'] =\ 'x-position of the telescope measured by...' except: pass try: tel_posY = [h.get_telescope_position(i)[1] for i in tel_id] tel_table_prime['TelY'] = tel_posY tel_table_prime['TelY'].unit = u.m except: pass try: tel_posZ = [h.get_telescope_position(i)[2] for i in tel_id] tel_table_prime['TelZ'] = tel_posZ tel_table_prime['TelZ'].unit = u.m except: pass try: tel_table['CameraClass'] = [h.get_camera_class(i) for i in tel_id] except: pass try: tel_table_prime['MirA'] = [h.get_mirror_area(i) for i in tel_id] tel_table_prime['MirA'].unit = u.m**2 except: pass try: tel_table_prime['MirN'] = [h.get_mirror_number(i) for i in tel_id] except: pass try: tel_table_prime['FL'] = [h.get_optical_foclen(i) for i in tel_id] tel_table_prime['FL'].unit = u.m except: pass try: tel_table_prime.meta['TelNum'] = len(tel_posX) except: pass #Beside other tables containimng telescope configuration data, the main #telescope table is written into the telescope dictionary. telescope['TelescopeTable_Version%s' % version] = tel_table_prime #-------------------------------------------------------------------------- #Camera and Optics configuration try: for t in range(len(tel_id)): cam_table_prime = Table() cam_table_prime.meta = {'TELID': tel_id[t], 'VERSION': version} opt_table_prime = Table() opt_table_prime.meta = {'TELID': tel_id[t], 'VERSION': version} try: pix_posX = h.get_pixel_position(tel_id[t])[0] pix_id = np.arange(len(pix_posX)) cam_table_prime['PixID'] = pix_id cam_table_prime['PixX'] = pix_posX cam_table_prime['PixX'].unit = u.m cam_table_prime.meta['PixXDescription'] =\ 'x-position of the pixel measured by...' except: pass try: pix_posY = h.get_pixel_position(tel_id[t])[1] cam_table_prime['PixY'] = pix_posY cam_table_prime['PixY'].unit = u.m except: pass try: camera_class = CD.guess_camera_geometry(pix_posX*u.m,pix_posY*u.m) pix_area_prime = camera_class.pix_area pix_type_prime = camera_class.pix_type pix_neighbors_prime = camera_class.pix_neighbors except: pass try: pix_area = h.get_pixel_area(tel_id[t]) cam_table_prime['PixA'] = pix_area cam_table_prime['PixA'].unit = u.mm**2 except: try: cam_table_prime['PixA'] = pix_area_prime cam_table_prime['PixA'].unit = u.mm**2 except: pass try: pix_type = h.get_pixel_type(tel_id[t]) except: try: pix_type = pix_type_prime except: pix_type = 'unknown' cam_table_prime.meta['PixType'] = pix_type try: pix_neighbors = h.get_pixel_neighbor(tel_id[t]) cam_table_prime['PixNeig'] = pix_neighbors except: try: cam_table_prime['PixNeig'] = pix_neighbors_prime except: pass #as long as no mirror IDs are given, use the following: opt_table_prime['MirrID'] = [1,2] try: opt_table_prime.meta['MirNum'] = h.get_mirror_number(tel_id[t]) except: pass try: opt_table_prime['MirArea'] = h.get_mirror_area(tel_id[t]) opt_table_prime['MirArea'].unit = u.m**2 opt_table_prime.meta['MirAreaDescription'] =\ 'Area of all mirrors' except: pass try: opt_table_prime['OptFocLen'] = h.get_optical_foclen(tel_id[t]) opt_table_prime['OptFocLen'].unit = u.m except: pass #Beside other tables containing camera and optics configuration #data, the main tables are written into the camera and optics #dictionary. camera['CameraTable_Version%s_TelID%i' % (version,tel_id[t])] \ = cam_table_prime optics['OpticsTable_Version%s_TelID%i' % (version,tel_id[t])] \ = opt_table_prime except: pass print('Astropy tables have been created.') h.close_file() print("Hessio file has been closed.") return(telescope,camera,optics)
def load_hessio(filename): """ Function to open and load a hessio file Parameters ---------- filename: string name of the file """ print("Hessio file will be opened.") h.file_open(filename) next(h.move_to_next_event()) #version = h.get... version = 'Feb2016' #Creating 3 dictionaries where the instrument configuration will be stored #The dictionaries themselves contain astropy.table.Table objects telescope = {} camera = {} optics = {} #-------------------------------------------------------------------------- #Telescope configuration tel_table_prime = Table() tel_table_prime.meta = {'VERSION': version} try: tel_id = h.get_telescope_ids() tel_table_prime['TelID'] = tel_id except: pass try: tel_posX = [h.get_telescope_position(i)[0] for i in tel_id] tel_table_prime['TelX'] = tel_posX tel_table_prime['TelX'].unit = u.m tel_table_prime.meta['TelX_description'] =\ 'x-position of the telescope measured by...' except: pass try: tel_posY = [h.get_telescope_position(i)[1] for i in tel_id] tel_table_prime['TelY'] = tel_posY tel_table_prime['TelY'].unit = u.m except: pass try: tel_posZ = [h.get_telescope_position(i)[2] for i in tel_id] tel_table_prime['TelZ'] = tel_posZ tel_table_prime['TelZ'].unit = u.m except: pass try: tel_table['CameraClass'] = [h.get_camera_class(i) for i in tel_id] except: pass try: tel_table_prime['MirA'] = [h.get_mirror_area(i) for i in tel_id] tel_table_prime['MirA'].unit = u.m**2 except: pass try: tel_table_prime['MirN'] = [h.get_mirror_number(i) for i in tel_id] except: pass try: tel_table_prime['FL'] = [h.get_optical_foclen(i) for i in tel_id] tel_table_prime['FL'].unit = u.m except: pass try: tel_table_prime.meta['TelNum'] = len(tel_posX) except: pass #Beside other tables containimng telescope configuration data, the main #telescope table is written into the telescope dictionary. telescope['TelescopeTable_Version%s' % version] = tel_table_prime #-------------------------------------------------------------------------- #Camera and Optics configuration try: for t in range(len(tel_id)): cam_table_prime = Table() cam_table_prime.meta = {'TELID': tel_id[t], 'VERSION': version} opt_table_prime = Table() opt_table_prime.meta = {'TELID': tel_id[t], 'VERSION': version} try: pix_posX = h.get_pixel_position(tel_id[t])[0] pix_id = np.arange(len(pix_posX)) cam_table_prime['PixID'] = pix_id cam_table_prime['PixX'] = pix_posX cam_table_prime['PixX'].unit = u.m cam_table_prime.meta['PixXDescription'] =\ 'x-position of the pixel measured by...' except: pass try: pix_posY = h.get_pixel_position(tel_id[t])[1] cam_table_prime['PixY'] = pix_posY cam_table_prime['PixY'].unit = u.m except: pass try: camera_class = CD.guess_camera_geometry( pix_posX * u.m, pix_posY * u.m) pix_area_prime = camera_class.pix_area pix_type_prime = camera_class.pix_type pix_neighbors_prime = camera_class.pix_neighbors except: pass try: pix_area = h.get_pixel_area(tel_id[t]) cam_table_prime['PixA'] = pix_area cam_table_prime['PixA'].unit = u.mm**2 except: try: cam_table_prime['PixA'] = pix_area_prime cam_table_prime['PixA'].unit = u.mm**2 except: pass try: pix_type = h.get_pixel_type(tel_id[t]) except: try: pix_type = pix_type_prime except: pix_type = 'unknown' cam_table_prime.meta['PixType'] = pix_type try: pix_neighbors = h.get_pixel_neighbor(tel_id[t]) cam_table_prime['PixNeig'] = pix_neighbors except: try: cam_table_prime['PixNeig'] = pix_neighbors_prime except: pass #as long as no mirror IDs are given, use the following: opt_table_prime['MirrID'] = [1, 2] try: opt_table_prime.meta['MirNum'] = h.get_mirror_number(tel_id[t]) except: pass try: opt_table_prime['MirArea'] = h.get_mirror_area(tel_id[t]) opt_table_prime['MirArea'].unit = u.m**2 opt_table_prime.meta['MirAreaDescription'] =\ 'Area of all mirrors' except: pass try: opt_table_prime['OptFocLen'] = h.get_optical_foclen(tel_id[t]) opt_table_prime['OptFocLen'].unit = u.m except: pass #Beside other tables containing camera and optics configuration #data, the main tables are written into the camera and optics #dictionary. camera['CameraTable_Version%s_TelID%i' % (version,tel_id[t])] \ = cam_table_prime optics['OpticsTable_Version%s_TelID%i' % (version,tel_id[t])] \ = opt_table_prime except: pass print('Astropy tables have been created.') h.close_file() print("Hessio file has been closed.") return (telescope, camera, optics)
def load_hessio(filename): """ Function to open and load a hessio file Parameters ---------- filename: string name of the file """ h.file_open(filename) print("Hessio file %s has been opened" % filename) next(h.move_to_next_event()) tel_id = h.get_telescope_ids() tel_table = Table() tel_table['TelID']= tel_id tel_posX = [h.get_telescope_position(i)[0] for i in tel_id] tel_posY = [h.get_telescope_position(i)[1] for i in tel_id] tel_posZ = [h.get_telescope_position(i)[2] for i in tel_id] tel_table['TelX'] = tel_posX tel_table['TelX'].unit = u.m tel_table['TelY'] = tel_posY tel_table['TelY'].unit = u.m tel_table['TelZ'] = tel_posZ tel_table['TelZ'].unit = u.m #tel_table['CameraClass'] = [h.get_camera_class(i) for i in tel_id] tel_table['MirrorArea'] = [h.get_mirror_area(i) for i in tel_id] tel_table['MirrorArea'].unit = u.m**2 tel_table['NMirrors'] = [h.get_mirror_number(i) for i in tel_id] tel_table['FL'] = [h.get_optical_foclen(i) for i in tel_id] tel_table['FL'].unit = u.m for t in range(len(tel_id)): table = Table() pix_posX = h.get_pixel_position(tel_id[t])[0] pix_posY = h.get_pixel_position(tel_id[t])[1] pix_id = np.arange(len(pix_posX)) pix_area = h.get_pixel_area(tel_id[t]) table['TelID'] = [tel_id[t] for i in range(len(pix_posX))] table['PixelID'] = pix_id table['PixX'] = pix_posX table['PixX'].unit = u.m table['PixY'] = pix_posY table['PixY'].unit = u.m table['PixArea'] = pix_area table['PixArea'].unit = u.mm**2 if t == 0: cam_table = table else: cam_table = join(cam_table,table,join_type='outer') #for t in range(len(tel_id)): # table = Table() # if t == 0: # opt_table = table # else: # opt_table = join(opt_table,table,join_type='outer') opt_table = Table() print('Astropy tables have been created.') #to use this, one has to go through every event of the run... #n_channel = h.get_num_channel(tel_id) #ld.channel_num = n_channel #for chan in range(n_channel): # ld.adc_samples.append(h.get_adc_sample(tel_id,chan).tolist()) h.close_file() print("Hessio file has been closed.") return [tel_table,cam_table,opt_table]