Skip to content

Cannot filter photons to use as input of ParticleFlow

FYI @mveghel @dovombru @mvesteri @ausachov @mramospe

We are trying to reduce the BW of some jets lines, after the changes reported in !4426 (comment 9377752). Since https://gitlab.cern.ch/lhcb/Moore/-/blob/master/Hlt/RecoConf/python/RecoConf/calorimeter_reconstruction.py?ref_type=heads#L237 is thought to be not configurable, we are trying to place some filters on photons right before ParticleFlow, in order to get back to the previous configuration. We are doing something like this in standard_particles.py (then we can create a make_photons_PF function)

@configurable
def make_photons(
    make_neutral_protoparticles=_make_neutral_protoparticles,
    pv_maker=_make_pvs,
    **kwargs,
):
    """creates photon LHCb::Particles from LHCb::ProtoParticles (PVs are optional)"""
    pvs = pv_maker()
    code = F.require_all(F.CLUSTERMATCH_CHI2>10)
    particles = PhotonMaker(
        InputProtoParticles=make_neutral_protoparticles(),
        InputPrimaryVertices=pvs,
        **kwargs,
    ).Particles
    particles = ParticleFilter(particles, F.FILTER(code))
    return particles

This seems not to work, and the message we get is

DeRich2RadiatorGas                     INFO Current GasParameters : P=605490536.714180 T=293.978897 Scale=0.996228
DeRich1RadiatorGas                     INFO Current GasParameters : P=606041349.940539 T=293.411403 Scale=1.038453
DependencyHandler                      INFO Inserted 3 [3] conditions to pool-iov: run(0):[0-9223372036854775807]   [0.00001 seconds]
DependencyHandler                      INFO Inserted 26 [26] conditions to pool-iov: run(0):[291593-9223372036854775807]   [0.00002 seconds]
DependencyHandler                      INFO Inserted 82 [82] conditions to pool-iov: run(0):[307586-307586]   [0.00004 seconds]
DetectorDataService                    INFO +  Created/Accessed a total of 6633 conditions (S:  3261,L:  3261,C:   111,M:0)  Load:4.36206s/4132KB Compute:4.20402s/196608KB
FTRawBankDecoder                    WARNING Suppressing message: 'No end fragment found. Ignoring the cluster.'
HLTControlFlowMgr                     FATAL Event failed in Node ParticleRangeFilter/SimpleJets_296edc46
HLTControlFlowMgr                      INFO Application return code set to 3
HLTControlFlowMgr                     FATAL Event 2 on slot 2 failed!
HLTControlFlowMgr                     FATAL *** Too many consecutive failures 1, stopping now ***
HLTControlFlowMgr                     FATAL Event failed in Node ParticleRangeFilter/SimpleJets_296edc46
HLTControlFlowMgr                      INFO Application return code set to 3
HLTControlFlowMgr                     FATAL Event 6 on slot 6 failed!
HLTControlFlowMgr                     FATAL Event failed in Node ParticleRangeFilter/SimpleJets_296edc46
HLTControlFlowMgr                      INFO Application return code set to 3
HLTControlFlowMgr                     FATAL Event 1 on slot 1 failed!
HLTControlFlowMgr                     FATAL Event failed in Node ParticleRangeFilter/SimpleJets_296edc46
HLTControlFlowMgr                      INFO Application return code set to 3
HLTControlFlowMgr                     FATAL Event 7 on slot 7 failed!
HLTControlFlowMgr                     FATAL Event failed in Node ParticleRangeFilter/SimpleJets_296edc46
HLTControlFlowMgr                      INFO Application return code set to 3
HLTControlFlowMgr                     FATAL Event 0 on slot 0 failed!
HLTControlFlowMgr                     FATAL Event failed in Node ParticleRangeFilter/SimpleJets_296edc46
HLTControlFlowMgr                      INFO Application return code set to 3
HLTControlFlowMgr                     FATAL Event 8 on slot 2 failed!
HLTControlFlowMgr                     FATAL Event failed in Node ParticleRangeFilter/SimpleJets_296edc46
HLTControlFlowMgr                      INFO Application return code set to 3
HLTControlFlowMgr                     FATAL Event 4 on slot 4 failed!
HLTControlFlowMgr                     FATAL Event failed in Node ParticleRangeFilter/SimpleJets_296edc46
HLTControlFlowMgr                      INFO Application return code set to 3
HLTControlFlowMgr                     FATAL Event 3 on slot 3 failed!
Gaudi__Histograming__Sink__Root        INFO Completed update of ROOT histograms in: hlt2_qee_bw_histos.root

One thing to notice is that, we trying to (re)build jets at DaVinci level, this trick seems to work.

Any help is appreciated!

Edited by Davide Zuliani