Reconstruction To Do List

Generation
1. Use genr8 for a2 production
2. Skim off (hd_filter) the wanted events by cutting on the total neutral energy (> 6 GeV) but record the total number of events before skimming. The resulting MC data set will be stored for general use.
3. Use Pythia for background generation after removing a₂ → η π⁰ from the MC data stream as well as events with total neutral energy <= 6GeV.
4. Use a2 charge exchange production xsection vs pythia total xsection to estimate number of signal and noise events
Analysis
1. After pruning the Pythia events, check the generated vs. reconstructed energy in a plot.
2. Look at the a2 mass peak - adjust cuts on the eta and pi0 masses then look at the effect on the different angular distributions (i.e. Gottfried-Jackson angles, etc.)
3. look at difference between reconstructed and generated z position for less than 65 cm
4. look at single photon reconstruction efficiencies
5. adjust the error matrix for the BCAL (look in PID/DPhoton_factory.cc)
Gather all relevant DocDB documents and cross section publications at one spot on the Wiki.

Documents to Review

There are 5000 thrown events (γp -> ηπ0p -> 4γp) therefore containing 5000 protons with an average of 4-5 showers produced in the GlueX calorimeters (shown below).

Looking at the difference between the azimuthal angle at the BCAL of the swum proton and the reconstructed showers (DPhoton) there is a definite peak at Δφ = 0.

Making a cut on Δφ < 0.1 radian (~5.7 degrees) there are 2860 surviving showers (~57%).

Looking at the difference between the final proton Z position and the possible showers, there is a clear peak around zero (shown in Red below) after making the cut on Δφ < 0.1 radians.