Allen merge requestshttps://gitlab.cern.ch/lhcb/Allen/-/merge_requests2024-03-28T19:50:49+01:00https://gitlab.cern.ch/lhcb/Allen/-/merge_requests/1414Add algorithms to filter charged or brem calo clusters2024-03-28T19:50:49+01:00Lucas Meyer GarciaAdd algorithms to filter charged or brem calo clustersUse existing CaloDigit-Track matching and Brem recovery information to filter calo clusters. First attempt at addressing Allen#484.
With the current implementation, a veto for track-matched clusters is applied when constructing the `Neut...Use existing CaloDigit-Track matching and Brem recovery information to filter calo clusters. First attempt at addressing Allen#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`).
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 \pm 68$ are signal (see fit below).
![diphoton](/uploads/7176ac388e3196e78a6476fabe696ab5/diphoton.png)
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 \pm 114$ are signal (see fit below).
![diphoton_tm](/uploads/577488a89052288d6e6a7e35b8289f96/diphoton_tm.png)
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 \pm 134$ are signal (see fit below).
![diphoton_bm](/uploads/1fd2f39ec610a554b095df94aac66b4b/diphoton_bm.png)
Finally, applying both vetoes we get $11351$ candiates and $922 \pm 59$ signal events (see fit below).
![diphoton_tbm](/uploads/c171047c43f596dc47dcbae96f375968/diphoton_tbm.png)
Below you can alse find the distribution of vetoed diphotons fitted with background only.
![diphoton_bg](/uploads/c8b0409918b3c57782754937c4bf13a3/diphoton_bg.png)
In summary,
| | Candidates | Variation | S/B |
| ------- | ---------- | --------- | ----------------- |
| Default | $16652$ | $-$ | $(5.9 \pm 0.4)\%$ |
| TM-only | $13140$ | $-21.1\%$ | $(7.2 \pm 0.9)\%$ |
| BM-only | $14445$ | $-13.3\%$ | $(6.3 \pm 0.9)\%$ |
| BM+TM | $11351$ | $-31.2\%$ | $(8.1 \pm 0.5)\%$ |Anfeng LiAnfeng Lihttps://gitlab.cern.ch/lhcb/Allen/-/merge_requests/1246Allen ecal llt decoding and ecal llt eTot lumi counter2024-03-18T08:58:42+01:00Shu XianAllen ecal llt decoding and ecal llt eTot lumi counter* merge @baudurie's [calo_llt_audurier](https://gitlab.cern.ch/lhcb/Allen/-/tree/calo_llt_audurier) branch
* some minor changes to make it runs on GPU
* Adds Ecal llt decoding algorithm: `calo_decode_llt`
* ET `m_et` is used for ...* merge @baudurie's [calo_llt_audurier](https://gitlab.cern.ch/lhcb/Allen/-/tree/calo_llt_audurier) branch
* some minor changes to make it runs on GPU
* Adds Ecal llt decoding algorithm: `calo_decode_llt`
* ET `m_et` is used for lumi counter
* Address, ET max and multiplicity is also decoded, but not used in this MR
* Adds sequence `calo_llt` for ECal llt decode testing
* Add `ECalLLTET` and `ECalLLTEtot` to the lumi ecal counters list
* The Calo LLT counters are computed in the FPGAs, they are treated as different Calo Lumi counters from the Allen-decoded ones.
* The very same LLT counters implemented in run1 and run2 worked well
* Both LLT and HLT1 Calo counters are supposed to be collected, in parallel
* Could be used as a cross-check, further optimization possible after validationDaniel Charles CraikAnfeng LiDaniel Charles Craik