def collection_procedure(filename, number_of_events=0, min_photons=np.NaN):
# File import -----------------------------------------------------------
importer = ImporterRoot()
importer.open_root_file(filename)
event_collection = importer.import_all_spad_events(number_of_events)
print("#### Opening file ####")
print(filename)
print(event_collection.qty_spad_triggered)
# Energy discrimination -------------------------------------------------
event_collection.remove_events_with_too_many_photons()
CEnergyDiscrimination.discriminate_by_energy(
event_collection, low_threshold_kev=0, high_threshold_kev=700)
# Filtering of unwanted photon types ------------------------------------
event_collection.remove_unwanted_photon_types(
remove_thermal_noise=False,
remove_after_pulsing=False,
remove_crosstalk=False,
remove_masked_photons=True)
event_collection.save_for_hardware_simulator()
# Sharing of TDCs --------------------------------------------------------
# event_collection.apply_tdc_sharing(pixels_per_tdc_x=1, pixels_per_tdc_y=1)
# First photon discriminator ---------------------------------------------
# DiscriminatorMultiWindow.DiscriminatorMultiWindow(event_collection)
DiscriminatorDualWindow.DiscriminatorDualWindow(event_collection,
min_photons)
#event_collection.remove_events_with_fewer_photons(100)
# Apply TDC - Must be applied after making the coincidences because the
# coincidence adds a random time offset to pairs of events
#tdc = CTdc(system_clock_period_ps=5000, tdc_bin_width_ps=1, tdc_jitter_std=1)
#tdc.get_sampled_timestamps(event_collection)
#tdc.get_sampled_timestamps(coincidence_collection.detector2)
# Making of coincidences -------------------------------------------------
coincidence_collection = CCoincidenceCollection(event_collection)
return event_collection, coincidence_collection
