Add additional capability (filtering by decay R/Z) to MCMultiParticle… #33033
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
…Filter; fill up vectors with defaults if necessary for MCSmartSingleParticleFilter
PR description:
Additional capability (filtering according to the decay region, the same functionality which has been already implemented in the MCSmartSingleParticleFilter) is added to the MCMultiParticleFilter. A small typo (bug, but no real effect) for loading motherID_ in the original filter has been fixed.
Fill up vectors with defaults if necessary, fix both MCMultiParticleFilter and MCSmartSingleParticleFilter filters.
PR validation:
Managed to filter the B(s,d) -> hh -> mumu events by adding such a filter (requiring two muon decaying from K/pi, and within the detector volume) to the process:
mugenfilter = cms.EDFilter(
"MCMultiParticleFilter",
ParticleID = cms.vint32(-13, 13, -13, 13),
MotherID = cms.untracked.vint32(321, -321, 211, -211),
MaxDecayRadius = cms.untracked.vdouble(2000.0, 2000.0, 2000.0, 2000.0),
MaxDecayZ = cms.untracked.vdouble(4000.0, 4000.0, 4000.0, 4000.0),
MinDecayZ = cms.untracked.vdouble(-4000.0, -4000.0, -4000.0, -4000.0),
EtaMax = cms.vdouble(2.5, 2.5, 2.5, 2.5),
PtMin = cms.vdouble(3.5, 3.5, 3.5, 3.5),
Status = cms.vint32(1,1,1,1),
NumRequired = cms.int32(2),
AcceptMore = cms.bool(True)
)
if this PR is a backport please specify the original PR and why you need to backport that PR:
Backport from PR #32999, in order to enable the generation for Run-2 UL