Skip to content

Add algorithms to filter charged or brem calo clusters

Lucas Meyer Garcia requested to merge lmg_veto_matched_photons into 2024-patches

Use existing CaloDigit-Track matching and Brem recovery information to filter calo clusters. First attempt at addressing #484. With the current implementation, a veto for track-matched clusters is applied when constructing the NeutralBasicParticle objects, while a veto for Bremsstrahlung-matched clusters is applied in the TwoCaloClusters selection algorithm (which is the underlying algorithm of Hlt1Pi02GammaGamma).

This MR also includes a Parallel Reduction algorithm (device/utils/reduction) to perform parallel sum of arrays. This is implementation is based on these examples, and it is used to replace the sequential sum of ECAL cell energies performed in TotalECALEnergy.cu.

Some preliminary results based on the Hlt1Pi02GammaGamma monitoring line can be found below. The line is run with no prescaling and without the nPV cut on the upgrade_mc_minbias_scifi_v5_retinacluster_000_v1_newLHCbID.mdf sample (10k events). Mass fits are performed on the diphoton mass distrubution using a Gaussian to describe the \pi^0 signal and a power law to describe background.

With the default requirements, 16652 cadidate diphotons are selected, out of which 980 +- 68 are signal (see fit below).

diphoton

Removing diclusters for which at least one seed cell has been matched to a track, the number of selected candidates falls to 13140, out of which 945 +- 114 are signal (see fit below).

diphoton_tm

Removing diclusters for which at least one seed cell has been matched to Bremsstrahlung, the number of selected candidates is 14445, out of which 910 +- 134 are signal (see fit below).

diphoton_bm

Finally, applying both vetoes we get 11351 candiates and 922 +- 59 signal events (see fit below).

diphoton_tbm

Below you can alse find the distribution of vetoed diphotons fitted with background only.

diphoton_bg

In summary,

Candidates Variation S/B
Default 16652 - (5.9 +- 0.4)%
TM-only 13140 -21.1% (7.2 +- 0.9)%
BM-only 14445 -13.3% (6.3 +- 0.9)%
BM+TM 11351 -31.2% (8.1 +- 0.5)%

Edit: Updating comparison after modifications.

New comparison done with the same HLT1 configuration: running Hlt1Pi02GammaGamma on upgrade_mc_minbias_scifi_v5_retinacluster_000_v1_newLHCbID.mdf with no prescaling and no n_pvs cuts.

With the current reconstruction and selection of 2024-patches, 21162 cadidate diphotons are selected, out of which 1014 +- 73 are signal (see fit below). diphoton_2024_patches

With this MR, 14238 cadidate diphotons are selected, out of which 910 +- 63 are signal (see fit below). diphoton_mr

In summary,

Candidates Variation S/B
2024-patches 21162 - 5.0%
This MR 14238 -32.7% 6.8%
Edited by Rosen Matev

Merge request reports

Loading