Semi-Merged Diphoton Systems in Higgs Decays

• Semi-merged diphoton systems are defined by an intermediate photon separation, combining a resolved and a merged photon cluster in exotic decay chains.
• Advanced ECAL clustering and shower-shape analyses, complemented by machine learning mass regression, enable accurate reconstruction of these complex photon signatures.
• These systems provide new search avenues at the LHC for light boson decays and extended resonance models while offering precise efficiency and background modeling insights.

A semi-merged diphoton system is a composite photon-like object arising from the decay of a light boson—commonly in exotic decay chains such as $H \to \mathcal{AA} \to \gamma\gamma\gamma\gamma$—where one $\mathcal{A} \to \gamma\gamma$ decay produces two resolved photon candidates, while the other yields a highly collimated photon pair reconstructed as a single merged photon cluster in the electromagnetic calorimeter (ECAL). This regime is distinguished by an opening angle, $\Delta R_{\gamma\gamma}$, intermediate between the limits set by ECAL granularity and cluster-separation criteria, specifically defined for $1 < m_{\mathcal{A}} < 15$ GeV in the context of LHC searches for new light resonances (Collaboration, 1 Jan 2026, Collaboration, 2024, Bi et al., 2015).

1. Kinematic Regimes and Definition

For a scalar $\mathcal{A}$ produced in Higgs decays, the photon pair from $\mathcal{A} \to \gamma\gamma$ acquires an opening angle $$\Delta R_{\gamma\gamma} \simeq 2 m_{\mathcal{A}} / p_T^{\mathcal{A}}$$ in the lab frame. Detailed studies by CMS show:

• $m_{\mathcal{A}} \lesssim 1$ GeV: $\Delta R_{\gamma\gamma} \lesssim 0.02$–$0.035$ (typical ECAL Molière radius); photons form a nearly indistinguishable shower overlap ("fully merged").
• $\mathcal{A} \to \gamma\gamma$0 GeV: $\mathcal{A} \to \gamma\gamma$1–$\mathcal{A} \to \gamma\gamma$2; partially resolved with strong overlap.
• $\mathcal{A} \to \gamma\gamma$3 GeV: $\mathcal{A} \to \gamma\gamma$4–$\mathcal{A} \to \gamma\gamma$5; photons are fully resolved as separate clusters.

The “semi-merged” diphoton regime is operationally defined for $\mathcal{A} \to \gamma\gamma$6, corresponding to intermediate photon separations. In this interval, one $\mathcal{A} \to \gamma\gamma$7 decay is reconstructed as two resolved photons ("resolved leg") and the other as a merged photon-like object ("merged leg"). For high-mass parents ($\mathcal{A} \to \gamma\gamma$8 in $\mathcal{A} \to \gamma\gamma$9), semi-merged diphotons occur for $\Delta R_{\gamma\gamma}$0–$\Delta R_{\gamma\gamma}$1 (Collaboration, 2024, Bi et al., 2015).

2. ECAL Clustering and Photon Identification

Photon clustering in CMS ECAL follows a seed-based algorithm:

• Seed crystals are required to have $\Delta R_{\gamma\gamma}$2 MeV.
• Basic clusters aggregate energy from adjacent crystals within local $\Delta R_{\gamma\gamma}$3 for superclusters (for $\Delta R_{\gamma\gamma}$4 GeV).
• For merged configurations, two collimated photon showers coalesce into an extended energy distribution, which standard PF (particle-flow) algorithms may not split.

Discrimination between single photons, merged diphotons, and hadronic backgrounds employs several shower-shape variables:

• $\Delta R_{\gamma\gamma}$5
• $\Delta R_{\gamma\gamma}$6
• $\Delta R_{\gamma\gamma}$7 (ratio of hadronic to electromagnetic energy)
• Charged-hadron isolation, $\Delta R_{\gamma\gamma}$8, and electron veto

Merged photon-like PF candidates that satisfy loose photon ID (tight $\Delta R_{\gamma\gamma}$9, $1 < m_{\mathcal{A}} < 15$0, $1 < m_{\mathcal{A}} < 15$1) but are not split into two PF photons are classified as the merged leg in semi-merged event selections (Collaboration, 1 Jan 2026).

3. Event Selection and Categorization

Selection of semi-merged diphoton topologies in dedicated searches (e.g., $1 < m_{\mathcal{A}} < 15$2) is performed as follows:

• Trigger requires diphoton events with $1 < m_{\mathcal{A}} < 15$3 GeV, leading $1 < m_{\mathcal{A}} < 15$4 GeV, subleading $1 < m_{\mathcal{A}} < 15$5 GeV; all photons must satisfy $1 < m_{\mathcal{A}} < 15$6 (restricted to ECAL barrel).
• Offline, exactly three photon-like PF candidates must pass preselection ($1 < m_{\mathcal{A}} < 15$7, $1 < m_{\mathcal{A}} < 15$8, $1 < m_{\mathcal{A}} < 15$9, $\mathcal{A}$0 GeV, no pixel seed).
• Triphoton invariant mass $\mathcal{A}$1 GeV.
• Resolved $\mathcal{A}$2 candidate comprises the closest photon pair; merged $\mathcal{A}$3 leg's mass is reconstructed via machine learning regression.

In extended resonance searches ($\mathcal{A}$4), semi-merged objects are classified by CNNs as either diphoton, single photon, or hadron, using a normalized $\mathcal{A}$5 ECAL energy image. Mass regression CNNs, taking crystal energy images and candidate $\mathcal{A}$6, predict $\mathcal{A}$7 to determine the cluster mass (Collaboration, 2024).

4. Mass Reconstruction and Machine Learning Techniques

Mass reconstruction for merged diphoton objects employs dedicated machine learning models:

• In $\mathcal{A}$8 searches, a graph neural network (GNN) regresses $\mathcal{A}$9 of the merged leg, producing linear response ($\mathcal{A} \to \gamma\gamma$0) across $\mathcal{A} \to \gamma\gamma$1–$\mathcal{A} \to \gamma\gamma$2 GeV with typical resolution $\mathcal{A} \to \gamma\gamma$3–$\mathcal{A} \to \gamma\gamma$4 GeV and scale uncertainty $\mathcal{A} \to \gamma\gamma$5–$\mathcal{A} \to \gamma\gamma$6.
• For $\mathcal{A} \to \gamma\gamma$7, the mass regression CNN outputs $\mathcal{A} \to \gamma\gamma$8; cluster mass is $\mathcal{A} \to \gamma\gamma$9. Resolution matches simulation within $$\Delta R_{\gamma\gamma} \simeq 2 m_{\mathcal{A}} / p_T^{\mathcal{A}}$$0; energy-scale uncertainty per cluster is $$\Delta R_{\gamma\gamma} \simeq 2 m_{\mathcal{A}} / p_T^{\mathcal{A}}$$1.

Selection efficiency for semi-merged topologies rises with $$\Delta R_{\gamma\gamma} \simeq 2 m_{\mathcal{A}} / p_T^{\mathcal{A}}$$2—$$\Delta R_{\gamma\gamma} \simeq 2 m_{\mathcal{A}} / p_T^{\mathcal{A}}$$3–$$\Delta R_{\gamma\gamma} \simeq 2 m_{\mathcal{A}} / p_T^{\mathcal{A}}$$4 for $$\Delta R_{\gamma\gamma} \simeq 2 m_{\mathcal{A}} / p_T^{\mathcal{A}}$$5 GeV—and the CNN classifier achieves $$\Delta R_{\gamma\gamma} \simeq 2 m_{\mathcal{A}} / p_T^{\mathcal{A}}$$655% efficiency for true merged $$\Delta R_{\gamma\gamma} \simeq 2 m_{\mathcal{A}} / p_T^{\mathcal{A}}$$7 (Collaboration, 1 Jan 2026, Collaboration, 2024).

5. Backgrounds and Statistical Modeling

Dominant backgrounds include:

• QCD multijet events with jets faking photons (“$$\Delta R_{\gamma\gamma} \simeq 2 m_{\mathcal{A}} / p_T^{\mathcal{A}}$$8+jets”, “jet+jet”)
• Prompt $$\Delta R_{\gamma\gamma} \simeq 2 m_{\mathcal{A}} / p_T^{\mathcal{A}}$$9+jet production

Background shapes are extracted from multiple sideband regions in the two-dimensional plane of merged vs. resolved masses, $m_{\mathcal{A}} \lesssim 1$0. Empirical functions (dijet-like, modified dijet, diphoton, power×exp, four-parameter) are fit to the invariant mass spectra in binned $m_{\mathcal{A}} \lesssim 1$1 categories, with discrete profiling and floating nuisance parameters (Collaboration, 2024).

Validation compares predicted 2D backgrounds to data in sideband regions, with residuals fit by Chebyshev polynomials to assign shape uncertainties. All systematic sources (luminosity, trigger, ID efficiency, energy scale/resolution, ML calibration, background function choice) are treated as nuisance parameters in final profile-likelihood fits (Collaboration, 1 Jan 2026, Collaboration, 2024).

6. Theoretical Interpretations and Model Significance

The semi-merged diphoton signature naturally arises in models with new light scalars coupling to the Higgs or heavy sector. One illustrative example is the $m_{\mathcal{A}} \lesssim 1$2 extension of the Standard Model (Bi et al., 2015), with $m_{\mathcal{A}} \lesssim 1$3 and $m_{\mathcal{A}} \lesssim 1$4 GeV. Highly boosted $m_{\mathcal{A}} \lesssim 1$5 yields photon pairs with $m_{\mathcal{A}} \lesssim 1$6; for $m_{\mathcal{A}} \lesssim 1$7 GeV, the pair is fully merged in the ECAL, satisfying $m_{\mathcal{A}} \lesssim 1$8.

Experimental constraints—dijet cross section bounds, photon-jet searches, ECAL granularity—are respected by restricting the relevant parameter space; prediction efficiency factors $m_{\mathcal{A}} \lesssim 1$9–$\Delta R_{\gamma\gamma} \lesssim 0.02$0 and production rates $\Delta R_{\gamma\gamma} \lesssim 0.02$1–$\Delta R_{\gamma\gamma} \lesssim 0.02$2 fb are demonstrated to explain observed anomalies without contradicting negative searches in the broader parameter space.

7. Experimental Limits and Outlook

Dedicated analyses at CMS have set stringent limits on the cross section times branching ratio: $\Delta R_{\gamma\gamma} \lesssim 0.02$3 is constrained to $\Delta R_{\gamma\gamma} \lesssim 0.02$4–$\Delta R_{\gamma\gamma} \lesssim 0.02$5 pb at $\Delta R_{\gamma\gamma} \lesssim 0.02$6 CL for $\Delta R_{\gamma\gamma} \lesssim 0.02$7–$\Delta R_{\gamma\gamma} \lesssim 0.02$8 GeV—the strongest bounds to date in the $\Delta R_{\gamma\gamma} \lesssim 0.02$9–$0.035$0 GeV regime (Collaboration, 1 Jan 2026). Analogous searches for $0.035$1 set 95% CL bounds from $0.035$2 to $0.035$3 fb for $0.035$4–$0.035$5 GeV and $0.035$6–$0.035$7 (Collaboration, 2024). No significant excess has been observed; future searches will benefit from enhanced ECAL granularity, improved machine learning mass regression, and further event topology exploitation, which will extend sensitivity into new regimes of collimated photon emission.