I analyzed 10 files of run 5242 which was a pi0 production run with the carbon target. I selected events with 2 clusters only and reconstructed the pi0 from the 2 clusters. Next, I flagged the veto status for each cluster. Comparisons are shown below with a few veto/cluster combinations.
Here is a legend of the different colors:
Black - No veto required.
Red - Veto says both clusters are photons.
Green - Veto says both clusters are charged.
Blue - Veto says 1 cluster is charged, and the other cluster is a photon.
Invariant Mass of 2 photons
Invariant Mass of 2 photons (better view of pi0 peak)
In this section, I analyzed the same files in the previous section, but this time I selected events with only 3 clusters. I made calculated the invariant mass for different 2 cluster combinations and plotted the mass with the same veto status cuts listed above.
Number of clusters per event
Invariant Mass of 2 photons
Invariant Mass of 2 photons (better view of pi0 peak)
I analyzed pi0 production data for events with clusters from 2 to 10. The plots below are from run 5242 (Carbon target) and 5267 (Lead target). The statistics in this example are low (1 run); however, the majority of the contributions to the pi0 peak are from 2, 3, and 4 clusters.
Carbon target
Lead target
In this figure, I made ratios of the mass spectra with veto cuts by the mass spectrum without veto cuts. The number of clusters selected was 2 and 3 per event.
This time I analyzed all of the Pb target data and plotted the 2 photon invariant mass. With the full statistics, the veto cuts do a better job of reducing the background.
Lead target (all files)
Below, I fitted the mass spectrum for the 2 cluster events when there is no veto cut(left) and when the veto cut is set for both clusters to be photons(right). The mass of the peak is 0.131 GeV with a FWHM = 5.8 MeV.
Inv Mass w/ 2 cluster events and no veto
Inv. Mass w/ 2 cluster events and both neutrals from veto
I focused on the pi0 peak of the right plot above to see what the fit was missing. The peak has a high mass shoulder and is not exactly gaussian.