def initialize_fits(filename,tel_id, file_closed):
hdulist = file_closed
if file_closed == 1:
hdulist = load_fits(filename)
teles = hdulist[1].data
ld.telescope = teles["TelID"].tolist()
else:
pass
index = ld.telescope_id.index(tel_id)
ld.mirror_area.append(hdulist[1].data[index]["MirrorArea"])
if file_closed == 1:
util_functions.close_fits(hdulist)
def initialize_telescope(filename, file_closed = True):
ld.clear_lists_telescope()
ld.clear_lists_camera()
ld.clear_lists_optics()
if 'simtel.gz' in filename:
file_closed = TD.initialize_hessio(filename,file_closed)
for tel in ld.telescope_id:
CD.initialize_hessio(filename,tel,file_closed)
OD.initialize_hessio(filename,tel,file_closed)
util_functions.close_hessio()
elif 'fits' in filename:
file_closed = TD.initialize_fits(filename,file_closed)
for tel in ld.telescope_id:
CD.initialize_fits(filename,tel,file_closed)
OD.initialize_fits(filename,tel,file_closed)
util_functions.close_fits(file_closed)
def initialize_telescope(filename, file_closed=True):
ld.clear_lists_telescope()
ld.clear_lists_camera()
ld.clear_lists_optics()
if 'simtel.gz' in filename:
file_closed = TD.initialize_hessio(filename, file_closed)
for tel in ld.telescope_id:
CD.initialize_hessio(filename, tel, file_closed)
OD.initialize_hessio(filename, tel, file_closed)
util_functions.close_hessio()
elif 'fits' in filename:
file_closed = TD.initialize_fits(filename, file_closed)
for tel in ld.telescope_id:
CD.initialize_fits(filename, tel, file_closed)
OD.initialize_fits(filename, tel, file_closed)
util_functions.close_fits(file_closed)
def initialize_fits(filename, tel_id, file_closed):
hdulist = file_closed
if file_closed == 1:
hdulist = util_functions.load_fits(filename)
teles = hdulist[1].data
ld.telescope = teles["TelID"].tolist()
else:
pass
index = ld.telescope_id.index(tel_id)
ld.camera_fov.append(hdulist[1].data[index]["FOV"])
pixel_id_cam = []
index2 = np.where(hdulist[2].data['L0ID'] == index)[0]
for i in index2:
pixel_id_cam.append(hdulist[2].data[i]['PixelID'])
ld.pixel_id.append(pixel_id_cam)
if file_closed == 1:
util_functions.close_fits(hdulist)
def initialize_fits(filename,tel_id, file_closed):
hdulist = file_closed
if file_closed == 1:
hdulist = util_functions.load_fits(filename)
teles = hdulist[1].data
ld.telescope = teles["TelID"].tolist()
else:
pass
index = ld.telescope_id.index(tel_id)
ld.camera_fov.append(hdulist[1].data[index]["FOV"])
pixel_id_cam = []
index2 = np.where(hdulist[2].data['L0ID'] == index)[0]
for i in index2:
pixel_id_cam.append(hdulist[2].data[i]['PixelID'])
ld.pixel_id.append(pixel_id_cam)
if file_closed == 1:
util_functions.close_fits(hdulist)