while checking pt of decay products of sister particle, allow multiple leptons and tracks

Author: mgratti

GeneratorInterface/GenFilters/plugins/PythiaFilterMotherSister.cc

@@ -68,20 +68,19 @@ bool PythiaFilterMotherSister::filter(edm::StreamID, edm::Event& iEvent, const e
                ++dau) {
             // find the daugther you're interested in
             if (abs((*dau)->pdg_id()) == abs(sisterID)) {
-              bool passTrackPt = false;
-              bool passLeptonPt = false;
+              int failTrackPt = 0;
+              int failLeptonPt = 0;
               // check pt of the nephews
               for (HepMC::GenVertex::particle_iterator nephew = (*dau)->end_vertex()->particles_begin(HepMC::children);
                    nephew != (*dau)->end_vertex()->particles_end(HepMC::children);
                    ++nephew) {
                 int nephew_pdgId = abs((*nephew)->pdg_id());
-                // implicit requirement that only one newphew is a lepton
-                if (nephew_pdgId == 11 or nephew_pdgId == 13 or nephew_pdgId == 15)
-                  passLeptonPt = ((*nephew)->momentum().perp() > minLeptonPt);
-                if (nephew_pdgId == 211)
-                  passTrackPt = ((*nephew)->momentum().perp() > minTrackPt);
+                if (minLeptonPt > 0. and (nephew_pdgId == 11 or nephew_pdgId == 13 or nephew_pdgId == 15))
+                  failLeptonPt += ((*nephew)->momentum().perp() < minLeptonPt);
+                if (minTrackPt > 0. and nephew_pdgId == 211)
+                  failTrackPt += ((*nephew)->momentum().perp() < minTrackPt);
               }
-              if (not passLeptonPt or not passTrackPt)
+              if (failLeptonPt > 0 or failTrackPt > 0)
                 return false;
               // calculate displacement of the sister particle, from production to decay
               HepMC::GenVertex* v1 = (*dau)->production_vertex();