@@ -27,12 +27,11 @@ PythiaFilterMotherSister::PythiaFilterMotherSister(const edm::ParameterSet& iCon
2727 sisterID(iConfig.getUntrackedParameter(" SisterID" , 0 )),
2828 maxSisDisplacement(iConfig.getUntrackedParameter(" MaxSisterDisplacement" , -1 .)),
2929 nephewIDs(iConfig.getUntrackedParameter(" NephewIDs" , std::vector<int >{0 })),
30- minNephewPts(iConfig.getUntrackedParameter(" MinNephewPts" , std::vector<double >{0 .}))
31- {
32- if (nephewIDs.size () != minNephewPts.size ()) {
33- throw cms::Exception (" BadConfig" ) << " PythiaFilterMotherSister: "
34- << " 'nephewIDs' and 'minNephewPts' need same length." ;
35- }
30+ minNephewPts(iConfig.getUntrackedParameter(" MinNephewPts" , std::vector<double >{0 .})) {
31+ if (nephewIDs.size () != minNephewPts.size ()) {
32+ throw cms::Exception (" BadConfig" ) << " PythiaFilterMotherSister: "
33+ << " 'nephewIDs' and 'minNephewPts' need same length." ;
34+ }
3635}
3736
3837PythiaFilterMotherSister::~PythiaFilterMotherSister () {
@@ -79,11 +78,12 @@ bool PythiaFilterMotherSister::filter(edm::StreamID, edm::Event& iEvent, const e
7978 ++nephew) {
8079 int nephew_pdgId = abs ((*nephew)->pdg_id ());
8180 for (unsigned int i = 0 ; i < nephewIDs.size (); i++) {
82- if (nephew_pdgId == abs (nephewIDs.at (i)))
83- failNephewPt += ((*nephew)->momentum ().perp () < minNephewPts.at (i));
81+ if (nephew_pdgId == abs (nephewIDs.at (i)))
82+ failNephewPt += ((*nephew)->momentum ().perp () < minNephewPts.at (i));
8483 }
8584 }
86- if (failNephewPt > 0 ) return false ;
85+ if (failNephewPt > 0 )
86+ return false ;
8787 // calculate displacement of the sister particle, from production to decay
8888 HepMC::GenVertex* v1 = (*dau)->production_vertex ();
8989 HepMC::GenVertex* v2 = (*dau)->end_vertex ();
0 commit comments