Add Truth plugins to fuse

[HenningSE]

This PR will add truth output to fuse. The logic differs from the one used in WFSim so one has to be a bit careful with the interpretation of the output.

The code in the PR works as follows:

• The cluster_id from the microphysics clustering will be propagated up to the propagated_photons level. For each photon, it is now possible to tell which cluster it came from.
• The PMTResponseAndDAQ puts these photons into waveforms and splits the waveforms into raw_records. It will now save an additional output called contributing_clusters. This output contains the information which cluster contributed with how many photons to a given raw_record.
• In the next step a new PeakTruth plugin evaluates which clusters contribute to a peak and calculates truth data for each peak. One of the provided columns is the observable_energy_truth (we can rename columns if needed), indicating how large the energy was that went into a certain peak. Some examples:
  • Let's assume there would be a cluster of 100 keV creating 100 S2 photons. If just 95 photons make it into a peak, the observable energy for this peak will be given as 95 keV.
  • If strax would split the S2 originating to this cluster in two peaks each with 50 photons, the observable energy for each peak would be listed as 50 keV.

[HenningSE]

The columns provided by PeakTruth are now:

• s1_photon_number_truth: Number of S1 photons in the peak
• s2_photon_number_truth: Number of S2 photons in the peak
• ap_photon_number_truth: Number of PMT afterpulse 'photons' in the peak
• raw_area_truth: Sum of the photon_gains for the peak decided by the channels PMT gain
• observable_energy_truth: Energy associated to the peak. (see description and example above)
• number_of_contributing_clusters: Number of clusters contributing to the peak
• average_x_of_contributing_clusters: Weighted average of the contributing clusters x position. The weights are the number of photons in the peak from a certain cluster divided by the total number of photons in the peak.
• average_y_of_contributing_clusters: See average_x_of_contributing_clusters, but for y
• average_z_of_contributing_clusters: See average_x_of_contributing_clusters, but for z

If there are more columns needed as output for the truth plugin, please comment below.

[HenningSE]

These are the things that I want to add next:

• Plugin that returns the clusters that do not make it into a peak
• Example notebook explaining how to use the truth plugins

[HenningSE]

Right now the weighted x and y position in the PeakTruth plugin are computed from the "real" position in the TPC. We should add the observed positions as well to allow a better comparison with our PosRec output.

[HenningSE]

I will now bump the plugin versions & update the documentation. When merging this PR we should also bump the overall fuse version to 1.1.0.

[mhliu0001]

@HenningSE In tagged clusters, can you also include how many S1/S2 photons created by a certain cluster are detected? Just by counting the propagated_photons, should be easy. Maybe call them photons_s1 and photons_s2?

[HenningSE]

@HenningSE In tagged clusters, can you also include how many S1/S2 photons created by a certain cluster are detected? Just by counting the propagated_photons, should be easy. Maybe call them photons_s1 and photons_s2?

We can use n_s1_photon_hits of the S1PhotonHits plugin and sum_s2_photons of the SecondaryScintillation plugin for these numbers or? They can be loaded together with the clusters.