#include <MuonCaloCleanerAllCrossed.h>

Inheritance diagram for MuonCaloCleanerAllCrossed:

Public Member Functions
	MuonCaloCleanerAllCrossed (const edm::ParameterSet &)

	~MuonCaloCleanerAllCrossed ()

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
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 (edm::Event &, const edm::EventSetup &, const reco::Candidate *, detIdToFloatMap &)

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

Private Attributes
edm::InputTag	srcESrecHits_

edm::InputTag	srcSelectedMuons_

TrackDetectorAssociator	trackAssociator_

TrackAssociatorParameters	trackAssociatorParameters_

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 collections of calorimeter recHits from which energy deposits of muons are subtracted.

This module attributes to the muon all energy deposited in calorimeter cells crossed by the muon.

WARNING: this module is to be used only for calorimeter types for which the assumption that no other particle deposited energy close to the muon track is a good approximation (e.g. HO)

Author: Tomasz Maciej Frueboes; Christian Veelken, LLR

Version

Revision:: 1.2

Id:: MuonCaloCleanerAllCrossed.h,v 1.2 2012/11/25 15:43:12 veelken Exp

Definition at line 37 of file MuonCaloCleanerAllCrossed.h.

Member Typedef Documentation

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

private

Definition at line 46 of file MuonCaloCleanerAllCrossed.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 18 of file MuonCaloCleanerAllCrossed.cc.

References edm::EDConsumerBase::consumesCollector(), edm::ParameterSet::getParameter(), TrackAssociatorParameters::loadParameters(), trackAssociator_, trackAssociatorParameters_, and TrackDetectorAssociator::useDefaultPropagator().

   : srcSelectedMuons_(cfg.getParameter<edm::InputTag>("selectedMuons")),
     srcESrecHits_(cfg.getParameter<edm::InputTag>("esRecHits"))
 {
   edm::ParameterSet cfgTrackAssociator = cfg.getParameter<edm::ParameterSet>("trackAssociator");
   edm::ConsumesCollector iC = consumesCollector();
   trackAssociatorParameters_.loadParameters(cfgTrackAssociator, iC);
   trackAssociator_.useDefaultPropagator();
 
   // maps of detId to energy deposit attributed to muon
   produces<detIdToFloatMap>("energyDepositsMuPlus");
   produces<detIdToFloatMap>("energyDepositsMuMinus");
 }

MuonCaloCleanerAllCrossed::~MuonCaloCleanerAllCrossed ( )

Definition at line 32 of file MuonCaloCleanerAllCrossed.cc.

 {
 // nothing to be done yet...
 }

Member Function Documentation

void MuonCaloCleanerAllCrossed::fillEnergyDepositMap	(	edm::Event &	evt,
		const edm::EventSetup &	es,
		const reco::Candidate *	muon,
		detIdToFloatMap &	energyDepositMap
	)

private

Definition at line 53 of file MuonCaloCleanerAllCrossed.cc.

References TrackDetMatchInfo::crossedEcalRecHits, TrackDetMatchInfo::crossedHcalRecHits, TrackDetMatchInfo::crossedHORecHits, TrackDetMatchInfo::crossedPreshowerIds, edm::Event::getByLabel(), getTrackDetMatchInfo(), srcESrecHits_, trackAssociator_, and trackAssociatorParameters_.

Referenced by produce().

 {
   TrackDetMatchInfo trackDetMatchInfo = getTrackDetMatchInfo(evt, es, trackAssociator_, trackAssociatorParameters_, muon);
   
   for ( std::vector<const EcalRecHit*>::const_iterator rh = trackDetMatchInfo.crossedEcalRecHits.begin();
     rh != trackDetMatchInfo.crossedEcalRecHits.end(); ++rh ) {
     energyDepositMap[(*rh)->detid().rawId()] += (*rh)->energy();
   }
   
   for ( std::vector<const HBHERecHit*>::const_iterator rh = trackDetMatchInfo.crossedHcalRecHits.begin();
     rh != trackDetMatchInfo.crossedHcalRecHits.end(); ++rh ) {
     energyDepositMap[(*rh)->detid().rawId()] += (*rh)->energy();
   }
   
   for ( std::vector<const HORecHit*>::const_iterator rh = trackDetMatchInfo.crossedHORecHits.begin();
     rh != trackDetMatchInfo.crossedHORecHits.end(); ++rh ) {
     energyDepositMap[(*rh)->detid().rawId()] += (*rh)->energy();
   }
   
   // TF: there exists no better way
   edm::Handle<ESRecHitCollection> esRecHits;
   evt.getByLabel(srcESrecHits_, esRecHits);
   for ( std::vector<DetId>::const_iterator detId = trackDetMatchInfo.crossedPreshowerIds.begin();
     detId != trackDetMatchInfo.crossedPreshowerIds.end(); ++detId ) {
     ESRecHitCollection::const_iterator rh = esRecHits->find(*detId);
     if ( rh != esRecHits->end() ) {
       energyDepositMap[rh->detid().rawId()] += rh->energy();
     }
   }
 }

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

privatevirtual

Implements edm::EDProducer.

Definition at line 37 of file MuonCaloCleanerAllCrossed.cc.

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

 {
   std::auto_ptr<detIdToFloatMap> energyDepositsMuPlus(new detIdToFloatMap());
   std::auto_ptr<detIdToFloatMap> energyDepositsMuMinus(new detIdToFloatMap());
   
   std::vector<reco::CandidateBaseRef> selMuons = getSelMuons(evt, srcSelectedMuons_);
   const reco::CandidateBaseRef muPlus  = getTheMuPlus(selMuons);
   const reco::CandidateBaseRef muMinus = getTheMuMinus(selMuons);
 
   if ( muPlus.isNonnull()  ) fillEnergyDepositMap(evt, es, &(*muPlus), *energyDepositsMuPlus);
   if ( muMinus.isNonnull() ) fillEnergyDepositMap(evt, es, &(*muMinus), *energyDepositsMuMinus);
 
   evt.put(energyDepositsMuPlus, "energyDepositsMuPlus");
   evt.put(energyDepositsMuMinus, "energyDepositsMuMinus");
 }