]
mag_output_sex = []
mag_output_wayback = []
mag_output_error_sex = []
mag_output_error_wayback = []
flux_input = []
flux_output = []
flux_output_error = []
flux_output_max = []
flux_output_max_error = []
number_lost_objects = []
# Call Source Extractor
catalog_name = sex_caller("{}.fits".format(PICTURE), "{}_{}".format(PICTURE, FILTER))
# Read catalog
catag = CatalogReader(catalog_name)
catag.read()
# Create object for simulation
simulation = ObjectCreator(catag.fcat)
simulation.general_histograms(catag.fcat)
print "\nMasking . . . {}\n"
# For-loop for 1 to 30 mag
for mag in mag_input:
simulation.matrix_data = []
simulation.masking_matrix("{}.fits".format(PICTURE))
