PFMuonCaloCleaner Class Reference

#include <PFMuonCaloCleaner.h>

Inheritance diagram for PFMuonCaloCleaner:

Public Member Functions
	PFMuonCaloCleaner (const edm::ParameterSet &)
	~PFMuonCaloCleaner ()
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	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
std::vector< ConsumesInfo >	consumesInfo () const
	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
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Types
typedef std::map< uint32_t, float >	detIdToFloatMap

Private Member Functions
void	fillEnergyDepositMap (const reco::Muon *, const edm::View< reco::PFCandidate > &, detIdToFloatMap &)
virtual void	produce (edm::Event &, const edm::EventSetup &)

Private Attributes
double	dRmatch_
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
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
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::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

Definition at line 17 of file PFMuonCaloCleaner.h.

Member Typedef Documentation

typedef std::map<uint32_t, float> PFMuonCaloCleaner::detIdToFloatMap

private

Definition at line 26 of file PFMuonCaloCleaner.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 23 of file PFMuonCaloCleaner.cc.

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

  : srcSelectedMuons_(cfg.getParameter<edm::InputTag>("selectedMuons")),
    srcPFCandidates_(cfg.getParameter<edm::InputTag>("pfCandidates")),
    dRmatch_(cfg.getParameter<double>("dRmatch"))
{
  // maps of detId to energy deposit attributed to muon
  produces<detIdToFloatMap>("energyDepositsMuPlus");
  produces<detIdToFloatMap>("energyDepositsMuMinus");

  verbosity_ = ( cfg.exists("verbosity") ) ?
    cfg.getParameter<int>("verbosity") : 0;
}

PFMuonCaloCleaner::~PFMuonCaloCleaner ( )

inline

Definition at line 21 of file PFMuonCaloCleaner.h.

{}

Member Function Documentation

void PFMuonCaloCleaner::fillEnergyDepositMap ( const reco::Muon *  muon, const edm::View< reco::PFCandidate > &  pfCandidates, detIdToFloatMap &  energyDepositMap  )

private

Definition at line 72 of file PFMuonCaloCleaner.cc.

References edm::OwnVector< T, P >::begin(), edm::View< T >::begin(), gather_cfg::cout, deltaR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, dRmatch_, edm::OwnVector< T, P >::end(), edm::View< T >::end(), reco::LeafCandidate::eta(), edm::hlt::Exception, customizeTrackingMonitorSeedNumber::idx, trackerHitRTTI::isMatched(), edm::InputTag::label(), reco::PFCandidate::mu, reco::LeafCandidate::p4(), reco::LeafCandidate::phi(), reco::LeafCandidate::pt(), srcSelectedMuons_, and verbosity_.

Referenced by produce().

{
  if ( !muon )
    throw cms::Exception("InvalidData") 
      << "Collection 'selectedMuons' = " << srcSelectedMuons_.label() << " is not of valid type !!\n";
  
  bool isMatched = false;
  for ( edm::View<reco::PFCandidate>::const_iterator pfCandidate = pfCandidates.begin();
    pfCandidate != pfCandidates.end(); ++pfCandidate ) {
    if ( pfCandidate->particleId() == reco::PFCandidate::mu ) {
      double dR = deltaR(pfCandidate->p4(), muon->p4());
      if ( dR < dRmatch_ ) {
    // NOTE: particleFlow sequence needs to be rerun for Zmumu event
    //       in order to recreate the links from PFCandidates->PFBlocks->PFClusters->PFRecHitFractions
    //      (not stored in AOD/RECO event content)
    const reco::PFCandidate::ElementsInBlocks& pfBlocks = pfCandidate->elementsInBlocks();
    for ( reco::PFCandidate::ElementsInBlocks::const_iterator pfBlock = pfBlocks.begin();
          pfBlock != pfBlocks.end(); ++pfBlock ) {
      const edm::OwnVector<reco::PFBlockElement>& pfBlockElements = pfBlock->first->elements();
      for ( edm::OwnVector<reco::PFBlockElement>::const_iterator pfBlockElement = pfBlockElements.begin();
        pfBlockElement != pfBlockElements.end(); ++pfBlockElement ) {
        if ( pfBlockElement->clusterRef().isNonnull() ) {
          reco::PFClusterRef pfCluster = pfBlockElement->clusterRef();
          const std::vector<reco::PFRecHitFraction>& pfRecHitFractions = pfCluster->recHitFractions();
          for ( std::vector<reco::PFRecHitFraction>::const_iterator pfRecHitFraction = pfRecHitFractions.begin();
            pfRecHitFraction != pfRecHitFractions.end(); ++pfRecHitFraction ) {
        const reco::PFRecHitRef& pfRecHit = pfRecHitFraction->recHitRef();
        energyDepositMap[pfRecHit->detId()] += pfRecHitFraction->fraction();
          }
        }
      }
    }
    isMatched = true;
      }
    }
  }
  if ( !isMatched ) {
    edm::LogWarning("PFMuonCaloCleaner") 
      << "Failed to match Muon to PFCandidate: Pt = " << muon->pt() << ", eta = " << muon->eta() << ", phi = " << muon->phi() << " !!" << std::endl;
    if ( verbosity_ ) {
      int idx = 0;
      for ( edm::View<reco::PFCandidate>::const_iterator pfCandidate = pfCandidates.begin();
        pfCandidate != pfCandidates.end(); ++pfCandidate ) {
    std::cout << "PFCandidate #" << idx << " (" << getPFCandidateType(pfCandidate->particleId()) << "):"
          << " Pt = " << pfCandidate->pt() << ", eta = " << pfCandidate->eta() << ", phi = " << pfCandidate->phi() << std::endl;
    ++idx;
      }
    }
  }
}

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

privatevirtual

Implements edm::EDProducer.

Definition at line 36 of file PFMuonCaloCleaner.cc.

References fillEnergyDepositMap(), edm::Event::getByLabel(), getSelMuons(), getTheMuMinus(), getTheMuPlus(), edm::RefToBase< T >::isNonnull(), reco::tau::pfCandidates(), edm::Event::put(), srcPFCandidates_, and srcSelectedMuons_.

{
  std::auto_ptr<detIdToFloatMap> energyDepositMuPlus(new detIdToFloatMap());
  std::auto_ptr<detIdToFloatMap> energyDepositMuMinus(new detIdToFloatMap());

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

  typedef edm::View<reco::PFCandidate> PFCandidateView;
  edm::Handle<PFCandidateView> pfCandidates;
  evt.getByLabel(srcPFCandidates_, pfCandidates);

  if ( muPlus.isNonnull()  ) fillEnergyDepositMap(dynamic_cast<const reco::Muon*>(&*muPlus), *pfCandidates, *energyDepositMuPlus);
  if ( muMinus.isNonnull() ) fillEnergyDepositMap(dynamic_cast<const reco::Muon*>(&*muMinus), *pfCandidates, *energyDepositMuMinus);

  evt.put(energyDepositMuPlus, "energyDepositsMuPlus");
  evt.put(energyDepositMuMinus, "energyDepositsMuMinus");
}

Member Data Documentation

double PFMuonCaloCleaner::dRmatch_

private

Definition at line 31 of file PFMuonCaloCleaner.h.

Referenced by fillEnergyDepositMap().

edm::InputTag PFMuonCaloCleaner::srcPFCandidates_

private

Definition at line 30 of file PFMuonCaloCleaner.h.

Referenced by produce().

edm::InputTag PFMuonCaloCleaner::srcSelectedMuons_

private

Definition at line 29 of file PFMuonCaloCleaner.h.

Referenced by fillEnergyDepositMap(), and produce().

int PFMuonCaloCleaner::verbosity_

private

Definition at line 33 of file PFMuonCaloCleaner.h.

Referenced by fillEnergyDepositMap(), and PFMuonCaloCleaner().