Photon-recoil lag correlation ($\ell =t_{\gamma }-t_e$)

A photon-recoil lag correlation asks how often a photon record occurs at a specified time displacement from an Electron record. Let $r_e,r_{\gamma }$ be observed separate-channel rates, $S$ the rate of detected genuine pairs before a timing cut, and $f(\ell )$ their normalized timing probability density function. For an original dilute, stationary model with observation time $T$, neglecting end effects and detector dead time, the expected pair density isThis counts all pairs of records, rather than at most one detector trigger per Electron. It uses independent unrelated detections for the flat term and an added same-event contribution. A normal distribution for timing spread is a chosen approximation, not a claim about the historical instrument. Intrinsic emission delay, counter response and readout uncertainty can all affect the observed probability density function.

Define the accepted lag interval by $|\ell |<w/2$: $w$ is its full width. Integrating the flat baseline gives the existing accidental coincidence expectation $T{r}_e{r}_{\gamma }w$. The additional area is $TS{\int }_{-w/2}^{w/2}f(\ell )d\ell$. Moving a window of the same full width away from the peak tests the baseline under the model's stationarity assumption. It does not remove correlated instrumental backgrounds or changing rates automatically.

For the lag model, choose $T=1000\mathrm{s}$, $r_e=2{\mathrm{s}}^{-1}$, $r_{\gamma }=0.5{\mathrm{s}}^{-1}$, $S=0.1{\mathrm{s}}^{-1}$ and timing standard deviation $\sigma =1\mathrm{ms}$. A full width $w=2\mathrm{ms}$ gives $2$ expected accidental coincidences and $100\mathrm{erf}(1/\sqrt{2})\simeq 68.27$ accepted genuine pairs, using the error function. The total is approximately $70.27$, whereas the displaced interval from $4$ to $6\mathrm{ms}$ has approximately $2.003$ expected pairs. A finite sample fluctuates around these expectations; these values are not measured counts or a significance calculation.

The right panel isolates the information lost by retaining singles totals alone. In a separate idealized model, stipulate $M=10^6$ independent $1\mathrm{ms}$ bins, with Bernoulli distribution indicators $X,Y$ for Electron and photon detection in a bin. Set $p_e=\mathrm{Pr}(X=1)=0.002$ and $p_{\gamma }=\mathrm{Pr}(Y=1)=0.0005$. Independence predicts $\mathrm{Pr}(X=1,Y=1)=p_e{p}_{\gamma }=10^{-6}$, hence one expected double detection. Alternatively chooseThese nonnegative probabilities sum to one and preserve both marginal probabilities, yet predict $101$ double detections in expectation. Their covariance is $0.000101-p_e{p}_{\gamma }=0.000100$. This is an exact joint probability distribution illustrating pairing; it does not model the timing spread in the left panel or reconstruct either historical hypothesis quantitatively. Its lesson is why two singles totals cannot decide whether the observed events share a cause.