zMCLeptonDaughters.cc File Reference

#include "PhysicsTools/CandUtils/interface/zMCLeptonDaughters.h"
#include "DataFormats/Candidate/interface/Candidate.h"
#include "FWCore/Utilities/interface/EDMException.h"
#include <cassert>

Go to the source code of this file.

Functions
std::pair< const Candidate *, const Candidate * >	zMCLeptonDaughters (const Candidate &z, int leptonPdgId)

Function Documentation

std::pair<const Candidate*, const Candidate*> zMCLeptonDaughters	(	const Candidate &	z,
		int	leptonPdgId
	)

Definition at line 8 of file zMCLeptonDaughters.cc.

References funct::abs(), reco::Candidate::daughter(), reco::Candidate::numberOfDaughters(), reco::Candidate::pdgId(), and reco::Candidate::status().

                                                                                                     {
  if(z.numberOfDaughters()<2)
    throw cms::Exception("RuntimeError") <<
      "calling helper function reco::zMCLeptonDaughters passing a Z candidate"
      "with less than 2 daughters (" << z.numberOfDaughters() << ").\n";
  const Candidate * dau0 = z.daughter(0);
  const Candidate * dau1 = z.daughter(1);
  for(size_t i0 = 0; i0 < dau0->numberOfDaughters(); ++i0) {
    const Candidate * ddau0 = dau0->daughter(i0);
    if(std::abs(ddau0->pdgId())==leptonPdgId && ddau0->status()==1) {
      dau0 = ddau0; break;
    }
  }
  for(size_t i1 = 0; i1 < dau1->numberOfDaughters(); ++i1) {
    const Candidate * ddau1 = dau1->daughter(i1);
    if(std::abs(ddau1->pdgId())==leptonPdgId && ddau1->status()==1) {
      dau1 = ddau1; break;
    }
  }
  assert(std::abs(dau0->pdgId())==leptonPdgId && dau0->status()==1);
  assert(std::abs(dau1->pdgId())==leptonPdgId && dau1->status()==1);
  return std::make_pair(dau0, dau1);
}