MuonPFCandidateCleaner Class Reference

Inheritance diagram for MuonPFCandidateCleaner:

Public Member Functions
	MuonPFCandidateCleaner (const edm::ParameterSet &)
	~MuonPFCandidateCleaner ()
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
virtual void	produce (edm::Event &, const edm::EventSetup &)

Private Attributes
double	dRmatch_
int	maxWarnings_tooFew_
int	maxWarnings_tooMany_
int	numWarnings_tooFew_
int	numWarnings_tooMany_
bool	removeDuplicates_
edm::InputTag	srcPFCandidates_
edm::InputTag	srcSelectedMuons_
int	verbosity_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Produce collection of PFCandidates in Z –> mu+ mu- event from which the two muons are removed (later to be replaced by simulated tau decay products)

Authors

Tomasz Maciej Frueboes; Christian Veelken, LLR

Version

Revision:

1.1

Id:

MuonPFCandidateCleaner.cc,v 1.1 2013/02/05 19:59:05 veelken Exp

Definition at line 33 of file MuonPFCandidateCleaner.cc.

Constructor & Destructor Documentation

MuonPFCandidateCleaner::MuonPFCandidateCleaner ( const edm::ParameterSet &  cfg )

explicit

Definition at line 56 of file MuonPFCandidateCleaner.cc.

References edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), and verbosity_.

  : srcSelectedMuons_(cfg.getParameter<edm::InputTag>("selectedMuons")),
    srcPFCandidates_(cfg.getParameter<edm::InputTag>("pfCands")),
    dRmatch_(cfg.getParameter<double>("dRmatch")),
    removeDuplicates_(cfg.getParameter<bool>("removeDuplicates")),
    maxWarnings_tooMany_(100),
    numWarnings_tooMany_(0),
    maxWarnings_tooFew_(3),
    numWarnings_tooFew_(0)
{
  verbosity_ = ( cfg.exists("verbosity") ) ?
    cfg.getParameter<int>("verbosity") : 0;

  produces<reco::PFCandidateCollection>();
}

MuonPFCandidateCleaner::~MuonPFCandidateCleaner ( )

inline

Definition at line 37 of file MuonPFCandidateCleaner.cc.

{}

Member Function Documentation

void MuonPFCandidateCleaner::produce ( edm::Event &  evt, const edm::EventSetup &  es  )

privatevirtual

Implements edm::EDProducer.

Definition at line 118 of file MuonPFCandidateCleaner.cc.

References gather_cfg::cout, reco::deltaR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, dRmatch_, reco::Candidate::eta(), edm::Event::getByLabel(), getSelMuons(), getTheMuMinus(), getTheMuPlus(), customizeTrackingMonitorSeedNumber::idx, reco::isMuon(), edm::RefToBase< T >::isNonnull(), maxWarnings_tooFew_, maxWarnings_tooMany_, numWarnings_tooFew_, numWarnings_tooMany_, reco::Candidate::p4(), reco::tau::pfCandidates(), reco::Candidate::phi(), reco::Candidate::pt(), edm::Event::put(), removeDuplicates_, python.multivaluedict::sort(), srcPFCandidates_, srcSelectedMuons_, and verbosity_.

{
  if ( verbosity_ ) std::cout << "<MuonPFCandidateCleaner::produce>:" << std::endl;

  std::vector<reco::CandidateBaseRef> selMuons = getSelMuons(evt, srcSelectedMuons_);
  const reco::CandidateBaseRef muPlus  = getTheMuPlus(selMuons);
  const reco::CandidateBaseRef muMinus = getTheMuMinus(selMuons);

  std::vector<reco::Particle::LorentzVector> selMuonP4s;
  if ( muPlus.isNonnull()  ) {
    if ( verbosity_ ) std::cout << " muPlus: Pt = " << muPlus->pt() << ", eta = " << muPlus->eta() << ", phi = " << muPlus->phi() << std::endl;
    selMuonP4s.push_back(muPlus->p4());
  }
  if ( muMinus.isNonnull() ) {
    if ( verbosity_ ) std::cout << " muMinus: Pt = " << muMinus->pt() << ", eta = " << muMinus->eta() << ", phi = " << muMinus->phi() << std::endl;
    selMuonP4s.push_back(muMinus->p4());
  }

//--- produce collection of PFCandidate excluding muons
  edm::Handle<reco::PFCandidateCollection> pfCandidates;
  evt.getByLabel(srcPFCandidates_, pfCandidates);

  std::auto_ptr<reco::PFCandidateCollection> pfCandidates_woMuons(new reco::PFCandidateCollection());   
   
//--- iterate over list of reconstructed PFCandidates, 
//    add PFCandidate to output collection in case it does not correspond to any selected muon
  std::vector<muonToPFCandMatchInfoType> selMuonToPFCandMatches;
  for ( std::vector<reco::Particle::LorentzVector>::const_iterator selMuonP4 = selMuonP4s.begin();
    selMuonP4 != selMuonP4s.end(); ++selMuonP4 ) {
    std::vector<muonToPFCandMatchInfoType> tmpMatches;
    for ( reco::PFCandidateCollection::const_iterator pfCandidate = pfCandidates->begin();
      pfCandidate != pfCandidates->end(); ++pfCandidate ) {
      double dR = reco::deltaR(pfCandidate->p4(), *selMuonP4);
      if ( dR < dRmatch_ ) tmpMatches.push_back(muonToPFCandMatchInfoType(*selMuonP4, &(*pfCandidate), dR));
    }
    // rank muon-to-pfCandidate matches by quality
    std::sort(tmpMatches.begin(), tmpMatches.end(), SortMuonToPFCandMatchInfosDescendingMatchQuality());
    if ( tmpMatches.size() > 0 ) selMuonToPFCandMatches.push_back(tmpMatches.front());
    if ( removeDuplicates_ ) {
      // CV: remove all high Pt charged PFCandidates very close to muon direction
      //    (duplicate tracks arise in case muon track in SiStrip + Pixel detector is reconstructed as 2 disjoint segments)
      for ( std::vector<muonToPFCandMatchInfoType>::const_iterator tmpMatch =  tmpMatches.begin();
        tmpMatch !=  tmpMatches.end(); ++tmpMatch ) {
    if ( tmpMatch->dR_ < 1.e-3 && fabs(tmpMatch->pfCandidateCharge_) > 0.5 && tmpMatch->pfCandidatePt_ > (0.33*tmpMatch->muonPt_) ) selMuonToPFCandMatches.push_back(*tmpMatch);
      }
    }
  }

  std::vector<const reco::PFCandidate*> removedPFCandidates;
  for ( reco::PFCandidateCollection::const_iterator pfCandidate = pfCandidates->begin();
    pfCandidate != pfCandidates->end(); ++pfCandidate ) {
    bool isMuon = false;
    for ( std::vector<muonToPFCandMatchInfoType>::const_iterator muonMatchInfo = selMuonToPFCandMatches.begin();
      muonMatchInfo != selMuonToPFCandMatches.end(); ++muonMatchInfo ) {
      if ( muonMatchInfo->pfCandidate_ == &(*pfCandidate) ) isMuon = true;
    }
    if ( verbosity_ && pfCandidate->pt() > 10. && fabs(pfCandidate->charge()) > 0.5 ) {
      std::cout << "pfCandidate: Pt = " << pfCandidate->pt() << ", eta = " << pfCandidate->eta() << ", phi = " << pfCandidate->phi() << ", isMuon = " << isMuon << std::endl;
    }
    if ( isMuon ) removedPFCandidates.push_back(&(*pfCandidate)); // pfCandidate belongs to a selected muon, do not copy
    else pfCandidates_woMuons->push_back(*pfCandidate);
  }
  if ( (removedPFCandidates.size() > selMuons.size() && numWarnings_tooMany_ < maxWarnings_tooMany_) &&
       (removedPFCandidates.size() < selMuons.size() && numWarnings_tooFew_  < maxWarnings_tooFew_ ) ) {
    edm::LogWarning("MuonPFCandidateCleaner") 
      << " (" << runLumiEventNumbers_to_string(evt) << ")" << std::endl
      << " Removed " << removedPFCandidates.size() << " PF-candidates from event containing " << selMuons.size() << " muons !!" << std::endl;
    if ( muPlus.isNonnull() ) std::cout << " muPlus: Pt = " << muPlus->pt() << ", eta = " << muPlus->eta() << ", phi = " << muPlus->phi() << std::endl;
    if ( muMinus.isNonnull() ) std::cout << " muMinus: Pt = " << muMinus->pt() << ", eta = " << muMinus->eta() << ", phi = " << muMinus->phi() << std::endl;
    int idx = 0;
    for ( std::vector<const reco::PFCandidate*>::const_iterator removedPFCandidate = removedPFCandidates.begin();
      removedPFCandidate != removedPFCandidates.end(); ++removedPFCandidate ) {
      std::cout << "PF-candidate #" << idx << " (charge = " << (*removedPFCandidate)->charge() << "):" 
        << " Pt = " << (*removedPFCandidate)->pt() << ", eta = " << (*removedPFCandidate)->eta() << ", phi = " << (*removedPFCandidate)->phi() << std::endl;
      ++idx;
    }
    if ( removedPFCandidates.size() > selMuons.size() ) ++numWarnings_tooMany_;
    if ( removedPFCandidates.size() < selMuons.size() ) ++numWarnings_tooFew_;
  }

  evt.put(pfCandidates_woMuons);
}

Member Data Documentation

double MuonPFCandidateCleaner::dRmatch_

private

Definition at line 45 of file MuonPFCandidateCleaner.cc.

Referenced by produce().

int MuonPFCandidateCleaner::maxWarnings_tooFew_

private

Definition at line 50 of file MuonPFCandidateCleaner.cc.

Referenced by produce().

int MuonPFCandidateCleaner::maxWarnings_tooMany_

private

Definition at line 48 of file MuonPFCandidateCleaner.cc.

Referenced by produce().

int MuonPFCandidateCleaner::numWarnings_tooFew_

private

Definition at line 51 of file MuonPFCandidateCleaner.cc.

Referenced by produce().

int MuonPFCandidateCleaner::numWarnings_tooMany_

private

Definition at line 49 of file MuonPFCandidateCleaner.cc.

Referenced by produce().

bool MuonPFCandidateCleaner::removeDuplicates_

private

Definition at line 46 of file MuonPFCandidateCleaner.cc.

Referenced by produce().

edm::InputTag MuonPFCandidateCleaner::srcPFCandidates_

private

Definition at line 43 of file MuonPFCandidateCleaner.cc.

Referenced by produce().

edm::InputTag MuonPFCandidateCleaner::srcSelectedMuons_

private

Definition at line 42 of file MuonPFCandidateCleaner.cc.

Referenced by produce().

int MuonPFCandidateCleaner::verbosity_

private

Definition at line 53 of file MuonPFCandidateCleaner.cc.

Referenced by MuonPFCandidateCleaner(), and produce().