MuonDetCleaner Class Reference

#include <MuonDetCleaner.h>

Inheritance diagram for MuonDetCleaner:

Public Member Functions
	MuonDetCleaner (const edm::ParameterSet &)
	~MuonDetCleaner ()
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
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) 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, int >	detIdToIntMap

Private Member Functions
void	fillHitMap (edm::Event &, const edm::EventSetup &, const reco::Candidate *, detIdToIntMap &)
virtual void	produce (edm::Event &, const edm::EventSetup &)

Private Attributes
edm::InputTag	srcSelectedMuons_
TrackDetectorAssociator	trackAssociator_
TrackAssociatorParameters	trackAssociatorParameters_
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 collections of recHits in muon detectors from which the hits of the two muons from the Z -> mu+ mu- decay are subtracted (the remaining hits may be due to pile-up/heavy quark decays in Z+jets events)

Author

Christian Veelken, LLR

Version

Revision:

1.2

Id:

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

Definition at line 31 of file MuonDetCleaner.h.

Member Typedef Documentation

typedef std::map<uint32_t, int> MuonDetCleaner::detIdToIntMap

private

Definition at line 40 of file MuonDetCleaner.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 13 of file MuonDetCleaner.cc.

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

  : srcSelectedMuons_(cfg.getParameter<edm::InputTag>("selectedMuons"))
{
  edm::ParameterSet cfgTrackAssociator = cfg.getParameter<edm::ParameterSet>("trackAssociator");
  edm::ConsumesCollector iC = consumesCollector();
  trackAssociatorParameters_.loadParameters(cfgTrackAssociator, iC);
  trackAssociator_.useDefaultPropagator();

  // maps of detId to number of hits attributed to muon
  produces<detIdToIntMap>("hitsMuPlus");
  produces<detIdToIntMap>("hitsMuMinus");

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

MuonDetCleaner::~MuonDetCleaner

(

)

Definition at line 29 of file MuonDetCleaner.cc.

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

Member Function Documentation

void MuonDetCleaner::fillHitMap ( edm::Event &  evt, const edm::EventSetup &  es, const reco::Candidate *  muon, detIdToIntMap &  hitMap  )

private

Definition at line 76 of file MuonDetCleaner.cc.

References TrackDetectorAssociator::Any, TrackDetectorAssociator::associate(), TrackDetMatchInfo::chambers, reco::Candidate::charge(), gather_cfg::cout, reco::Candidate::eta(), edm::Ref< C, T, F >::isNonnull(), trackingTruthProducerFastSim_cfi::muon, reco::Muon::outerTrack(), reco::Candidate::phi(), reco::Candidate::pt(), reco::Candidate::px(), reco::Candidate::py(), reco::Candidate::pz(), TrackingRecHit::rawId(), AlCaHLTBitMon_QueryRunRegistry::string, trackAssociator_, trackAssociatorParameters_, verbosity_, and reco::Candidate::vertex().

Referenced by produce().

{
  int numHits = 0;
  
  const reco::Muon* recoMuon = dynamic_cast<const reco::Muon*>(muon);
  if ( recoMuon && recoMuon->outerTrack().isNonnull() ) {
    const reco::TrackRef muonOuterTrack = recoMuon->outerTrack();   
    TrackDetMatchInfo trackDetMatchInfo = trackAssociator_.associate(evt, es, *muonOuterTrack, trackAssociatorParameters_, TrackDetectorAssociator::Any);
    for ( std::vector<TAMuonChamberMatch>::const_iterator rh = trackDetMatchInfo.chambers.begin();
      rh != trackDetMatchInfo.chambers.end(); ++rh ) {
      ++hitMap[rh->id.rawId()];
      ++numHits;
    }
    for ( trackingRecHit_iterator rh = muonOuterTrack->recHitsBegin();
      rh != muonOuterTrack->recHitsEnd(); ++rh ) {
      fillHitMapRH(**rh, hitMap, numHits);
    }
  } else {
    GlobalVector muonP3(muon->px(), muon->py(), muon->pz()); 
    GlobalPoint muonVtx(muon->vertex().x(), muon->vertex().y(), muon->vertex().z());
    TrackDetMatchInfo trackDetMatchInfo = trackAssociator_.associate(evt, es, muonP3, muonVtx, muon->charge(), trackAssociatorParameters_);
    for ( std::vector<TAMuonChamberMatch>::const_iterator rh = trackDetMatchInfo.chambers.begin();
      rh != trackDetMatchInfo.chambers.end(); ++rh ) {
      ++hitMap[rh->id.rawId()];
      ++numHits;
    }
  }
  
  if ( verbosity_ ) {
    std::string muonCharge_string = "";
    if      ( muon->charge() > +0.5 ) muonCharge_string = "+";
    else if ( muon->charge() < -0.5 ) muonCharge_string = "-";
    std::cout << "Mu" << muonCharge_string << ": Pt = " << muon->pt() << ", eta = " << muon->eta() << ", phi = " << muon->phi() 
          << " --> #Hits = " << numHits << std::endl;
    if ( verbosity_ >= 2 ) printHitMapRH(es, hitMap);    
  }
}

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

privatevirtual

Implements edm::EDProducer.

Definition at line 34 of file MuonDetCleaner.cc.

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

{
  std::auto_ptr<detIdToIntMap> hitsMuPlus(new detIdToIntMap());
  std::auto_ptr<detIdToIntMap> hitsMuMinus(new detIdToIntMap());
  
  std::vector<reco::CandidateBaseRef> selMuons = getSelMuons(evt, srcSelectedMuons_);
  const reco::CandidateBaseRef muPlus  = getTheMuPlus(selMuons);
  const reco::CandidateBaseRef muMinus = getTheMuMinus(selMuons);

  if ( muPlus.isNonnull()  ) fillHitMap(evt, es, &(*muPlus), *hitsMuPlus);
  if ( muMinus.isNonnull() ) fillHitMap(evt, es, &(*muMinus), *hitsMuMinus);

  evt.put(hitsMuPlus, "hitsMuPlus");
  evt.put(hitsMuMinus, "hitsMuMinus");
}

Member Data Documentation

edm::InputTag MuonDetCleaner::srcSelectedMuons_

private

Definition at line 43 of file MuonDetCleaner.h.

Referenced by produce().

TrackDetectorAssociator MuonDetCleaner::trackAssociator_

private

Definition at line 45 of file MuonDetCleaner.h.

Referenced by fillHitMap(), and MuonDetCleaner().

TrackAssociatorParameters MuonDetCleaner::trackAssociatorParameters_

private

Definition at line 46 of file MuonDetCleaner.h.

Referenced by fillHitMap(), and MuonDetCleaner().

int MuonDetCleaner::verbosity_

private

Definition at line 48 of file MuonDetCleaner.h.

Referenced by fillHitMap(), and MuonDetCleaner().