WIP: Add RD strange four body HLT2 lines
Implements HLT2 lines for the following 4-body strange decay lines:
- KS0 -> 2(mu+ mu-)
- KS0 -> 2(e+ e-)
- KS0 -> pi+ pi- e+ e-
- KS0 -> pi+ pi- mu+ mu-
- KS0 -> mu+ mu- e+ e-
running over upgrade_minbias_hlt1_filtered
gives (after some updates and fixing the muon PID bug)
----------------------------------------------------------------------------------------------------
INFO: Starting /scratch03/titus.mombacher/HLT2Dev/stack/MooreAnalysis/HltEfficiencyChecker/scripts/hlt_calculate_rates.py...
----------------------------------------------------------------------------------------------------
INFO: No lines specified. Defaulting to all...
HLT rates:
----------------------------------------------------------------------------------------------------
Line: Hlt2RD_KS0To4E_LineDecision Incl: 0.13348461589801777 +/- 0.09 kHz, Excl: 0.13348461589801777 +/- 0.09 kHz
Line: Hlt2RD_KS0To4Mu_LineDecision Incl: 0.06674230794900889 +/- 0.06 kHz, Excl: 0.0 +/- 0.0 kHz
Line: Hlt2RD_KS0ToMuMuEE_LineDecision Incl: 0.06674230794900889 +/- 0.06 kHz, Excl: 0.0 +/- 0.0 kHz
Line: Hlt2RD_KS0ToPiPiEE_LineDecision Incl: 0.13348461589801777 +/- 0.09 kHz, Excl: 0.13348461589801777 +/- 0.09 kHz
Line: Hlt2RD_KS0ToPiPiMuMu_LineDecision Incl: 0.06674230794900889 +/- 0.06 kHz, Excl: 0.06674230794900889 +/- 0.06 kHz
----------------------------------------------------------------------------------------------------
Hlt2 Total: Rate: 0 +/- 0.2 kHz
----------------------------------------------------------------------------------------------------
TOS Efficiency from running on signal MC (outdated, but the numbers shouldn't have changed significantly):
- KS24mu:
Integrated HLT efficiencies for the lines with denominator: AllEvents
------------------------------------------------------------------------------------
Line: Hlt2RD_KS0To4E_LineDecision Efficiency: 0.000 +/- 0.000
Line: Hlt2RD_KS0To4Mu_LineDecision Efficiency: 4.4e-4
Line: Hlt2RD_KS0ToMuMuEE_LineDecision Efficiency: 0.000 +/- 0.000
Line: Hlt2RD_KS0ToPiPiEE_LineDecision Efficiency: 0.000 +/- 0.000
Line: Hlt2RD_KS0ToPiPiMuMu_LineDecision Efficiency: 0.002 +/- 0.000
Line: KS_Hlt2RD_KS0To4E_LineDecisionTrueSim Efficiency: 0.000 +/- 0.000
Line: KS_Hlt2RD_KS0To4Mu_LineDecisionTrueSim Efficiency: 2.2e-4
Line: KS_Hlt2RD_KS0ToMuMuEE_LineDecisionTrueSim Efficiency: 0.000 +/- 0.000
Line: KS_Hlt2RD_KS0ToPiPiEE_LineDecisionTrueSim Efficiency: 0.000 +/- 0.000
Line: KS_Hlt2RD_KS0ToPiPiMuMu_LineDecisionTrueSim Efficiency: 0.002 +/- 0.000
------------------------------------------------------------------------------------
Integrated HLT efficiencies for the lines with denominator: CanRecoChildren
------------------------------------------------------------------------------------
Line: Hlt2RD_KS0To4E_LineDecision Efficiency: 0.000 +/- 0.000
Line: Hlt2RD_KS0To4Mu_LineDecision Efficiency: 0.081 +/- 0.013
Line: Hlt2RD_KS0ToMuMuEE_LineDecision Efficiency: 0.009 +/- 0.004
Line: Hlt2RD_KS0ToPiPiEE_LineDecision Efficiency: 0.005 +/- 0.003
Line: Hlt2RD_KS0ToPiPiMuMu_LineDecision Efficiency: 0.041 +/- 0.009
Line: KS_Hlt2RD_KS0To4E_LineDecisionTrueSim Efficiency: 0.000 +/- 0.000
Line: KS_Hlt2RD_KS0To4Mu_LineDecisionTrueSim Efficiency: 0.074 +/- 0.012
Line: KS_Hlt2RD_KS0ToMuMuEE_LineDecisionTrueSim Efficiency: 0.000 +/- 0.000
Line: KS_Hlt2RD_KS0ToPiPiEE_LineDecisionTrueSim Efficiency: 0.002 +/- 0.002
Line: KS_Hlt2RD_KS0ToPiPiMuMu_LineDecisionTrueSim Efficiency: 0.036 +/- 0.009
------------------------------------------------------------------------------------
- KS22pi2e:
Integrated HLT efficiencies for the lines with denominator: AllEvents
------------------------------------------------------------------------------------
Line: Hlt2RD_KS0To4E_LineDecision Efficiency: 0.000 +/- 0.000
Line: Hlt2RD_KS0To4Mu_LineDecision Efficiency: 0.000 +/- 0.000
Line: Hlt2RD_KS0ToMuMuEE_LineDecision Efficiency: 0.000 +/- 0.000
Line: Hlt2RD_KS0ToPiPiEE_LineDecision Efficiency: 7.5e-4
Line: Hlt2RD_KS0ToPiPiMuMu_LineDecision Efficiency: 0.000 +/- 0.000
Line: KS_Hlt2RD_KS0To4E_LineDecisionTrueSim Efficiency: 0.000 +/- 0.000
Line: KS_Hlt2RD_KS0To4Mu_LineDecisionTrueSim Efficiency: 0.000 +/- 0.000
Line: KS_Hlt2RD_KS0ToMuMuEE_LineDecisionTrueSim Efficiency: 0.000 +/- 0.000
Line: KS_Hlt2RD_KS0ToPiPiEE_LineDecisionTrueSim Efficiency: 4.2e-4
Line: KS_Hlt2RD_KS0ToPiPiMuMu_LineDecisionTrueSim Efficiency: 0.000 +/- 0.000
------------------------------------------------------------------------------------
------------------------------------------------------------------------------------
Integrated HLT efficiencies for the lines with denominator: CanRecoChildren
------------------------------------------------------------------------------------
Line: Hlt2RD_KS0To4E_LineDecision Efficiency: 0.004 +/- 0.002
Line: Hlt2RD_KS0To4Mu_LineDecision Efficiency: 0.000 +/- 0.000
Line: Hlt2RD_KS0ToMuMuEE_LineDecision Efficiency: 0.004 +/- 0.002
Line: Hlt2RD_KS0ToPiPiEE_LineDecision Efficiency: 0.040 +/- 0.007
Line: Hlt2RD_KS0ToPiPiMuMu_LineDecision Efficiency: 0.001 +/- 0.001
Line: KS_Hlt2RD_KS0To4E_LineDecisionTrueSim Efficiency: 0.000 +/- 0.000
Line: KS_Hlt2RD_KS0To4Mu_LineDecisionTrueSim Efficiency: 0.000 +/- 0.000
Line: KS_Hlt2RD_KS0ToMuMuEE_LineDecisionTrueSim Efficiency: 0.003 +/- 0.002
Line: KS_Hlt2RD_KS0ToPiPiEE_LineDecisionTrueSim Efficiency: 0.037 +/- 0.007
Line: KS_Hlt2RD_KS0ToPiPiMuMu_LineDecisionTrueSim Efficiency: 0.001 +/- 0.001
------------------------------------------------------------------------------------
Merge request reports
Activity
assigned to @jubermud
added RTA label
assigned to @sargueda
So I don't know what your plans are @msaur (if you just write up your plans here, we can restructure also in this MR), but there are a few things we should do here:
- physics-wise: one could play a bit with the PID cuts. The PIDK cut on the pions is probably too tight. Requiring isMuon is also quite strict. Replacing this by PIDmu>0 could be possible as well. The PIDe cuts can be tightened a bit to increase efficiencies.
- organisatorial:
- probably the builders and line registration should be split in two files
- the lines have to be added to the RD-level
__init__.py
- the ghostprob cuts are supposed to be removed for now. They are organised centrally
Is there a reason why you don't include downstream lines as well? For KS0 decays you should be able to double your yield this way, and the backgrounds will also be lower so you can probably relax some of your cuts. Obviously we won't have downstream tracks for the first part of R3 but perhaps it would be good to have these in as placeholders, even if they are switched off for now?
I agree with Titus that requiring isMuon on all four muons is likely to be rather harsh. One option might be to add a second line that takes three muons satisfying isMuon with a fourth track that fails isMuon (it would need to explicitly fail the cut otherwise you will end up duplicating candidates in the 4mu line). Something like this was done for in the stripping for B0->4mu in R1/R2. In that case signal was swamped by hadronic background from decays like B0->Jpsi(mm)pipi but in your case I imagine it would be less of an issue since KS0->pipimumu is right at the kinematic limit.
right. I added some downstream placeholders, but they have trouble atm (not sure which functor is causing it). So in favour of getting these lines ready now-ish I would prefer to not fiddle too much with the downstream stuff. Indeed probably the PID is not needed for the pipimumu and 4mu because of the kinematic limit
Basically the rate and efficiencies are zero, from a sort-of quick check I didn't see error messages pointing me to what is buggy and atm I don't have time to check every individual functor. So I thought it's best to get the long track versions merged first (I definitely want at least one version of these lines in the system, so it's well prepared for the early data), as we won't have downstream tracking for a while. Sorry to not be able to give you more details atm.
Thanks @tmombach. It looks good. Happy to approve the merge request on the assumption that some further tweaks may be needed once the muon PID issue is resolved.
After all I noticed that the muon ID for the KS24mu ones can be relaxed to PIDmu>-1 (because of the kinematic edge). In that case the extra line with even looser muon requirements doesn't really add anything. So I'm leaving it out again. Will update the signal efficiencies on simulation in the coming days and then hopefully this can go.
mentioned in merge request !1353 (closed)
added 260 commits
-
527e069d...4f907403 - 257 commits from branch
master
- 87e90cbb - Merge branch 'master' into jubermud_RD_4body_lines
- 719aecb3 - Merge branch 'master' into jubermud_RD_4body_lines
- 58cf4c1d - modify a bit to improve rates and include rd prefilters
Toggle commit list-
527e069d...4f907403 - 257 commits from branch
mentioned in merge request !1399 (merged)
added 118 commits
-
a6861765...b76cc1cc - 116 commits from branch
master
- ced17dd6 - Merge branch 'master' into jubermud_RD_4body_lines
- 3fd35b30 - add 3mupi line
-
a6861765...b76cc1cc - 116 commits from branch
added 219 commits
-
3fd35b30...af3a0dd6 - 217 commits from branch
master
- e1563f91 - fix merge conflict
- 2fd50bea - fix merge conflict
-
3fd35b30...af3a0dd6 - 217 commits from branch