# How are the HEP experiments' invariant mass plots generated?

I am experimenting and playing around with some data, and I'm having trouble seeing how to generate invariant mass plots.

The data I have has a bunch of events, and variables such as $P,P_T,\eta,\phi$ etc.. but no energy. There is a histogram which generates the invariant mass, but no macro provided for this. For example, I know there are tracks for electrons and protons (TPC), and I want to see if there are any pairs as a result of $\gamma\rightarrow e^+ e^-$. Conceptually I need to loop over all positive negative pairs, and plot the histogram. I should get a peak around zero for $\gamma$ pair production. The invariant mass formula is $$m_\gamma^2 = (E_{e^+}+E_{e^-})^2 - (\vec{p}_{e^+}+\vec{p}_{e^-})^2$$

But there is no data for the energy, how do the high energy physicists do it? Do they guess some value for the energy and make plots with that guess?

• I have thought about it, and I think experimenters must have some idea of the identity of the individual tracks. Otherwise, any combinatorial algorithm is meaningless. – yayu Aug 11 '11 at 22:06

• If you have some kind of PID reason to believe that the tracks might be electron you just assume that they are and compute the energy from the momentum and $m_e$.
• In you have a calorimeter in the detector stack you measure the energy, and project the maximum likelihood energy loss back to the vertex. Note, however that the measured signal may have error much larger than $m_e$, so it may be better to use the device for PID and fall back on the previous suggestion.