MuonCaloDistanceProducer Class Reference

#include <MuonCaloDistanceProducer.h>

Inheritance diagram for MuonCaloDistanceProducer:

Public Member Functions
	MuonCaloDistanceProducer (const edm::ParameterSet &)
	~MuonCaloDistanceProducer ()
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 (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &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	fillDistanceMap (edm::Event &, const edm::EventSetup &, const reco::Candidate *, detIdToFloatMap &, detIdToFloatMap &)
virtual void	produce (edm::Event &, const edm::EventSetup &)

Private Attributes
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
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

Compute distance traversed by muon through calorimeter cells.

NOTE: The output of this class is used as input for MuonCaloCleanerByDistance.

Author

Tomasz Maciej Frueboes; Christian Veelken, LLR

Version

Revision:

1.3

Id:

MuonCaloDistanceProducer.h,v 1.3 2013/01/30 10:10:49 aburgmei Exp

Definition at line 30 of file MuonCaloDistanceProducer.h.

Member Typedef Documentation

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

private

Definition at line 39 of file MuonCaloDistanceProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 24 of file MuonCaloDistanceProducer.cc.

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

  : srcSelectedMuons_(cfg.getParameter<edm::InputTag>("selectedMuons"))
{
  // maps of detId to distance traversed by muon through detector volume
  produces<reco::CandidateCollection>("muons");
  produces<detIdToFloatMap>("distancesMuPlus");
  produces<detIdToFloatMap>("distancesMuMinus");
  produces<detIdToFloatMap>("depositsMuPlus");
  produces<detIdToFloatMap>("depositsMuMinus");
  
  edm::ParameterSet cfgTrackAssociator = cfg.getParameter<edm::ParameterSet>("trackAssociator");
  edm::ConsumesCollector iC = consumesCollector();
  trackAssociatorParameters_.loadParameters(cfgTrackAssociator, iC);
  trackAssociator_.useDefaultPropagator();
}

MuonCaloDistanceProducer::~MuonCaloDistanceProducer

(

)

Definition at line 40 of file MuonCaloDistanceProducer.cc.

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

Member Function Documentation

void MuonCaloDistanceProducer::fillDistanceMap ( edm::Event &  evt, const edm::EventSetup &  es, const reco::Candidate *  muon, detIdToFloatMap &  distanceMap, detIdToFloatMap &  depositMap  )

private

Definition at line 75 of file MuonCaloDistanceProducer.cc.

References assert(), TrackDetMatchInfo::crossedEcalIds, TrackDetMatchInfo::crossedEcalRecHits, TrackDetMatchInfo::crossedHcalIds, TrackDetMatchInfo::crossedHcalRecHits, TrackDetMatchInfo::crossedHOIds, TrackDetMatchInfo::crossedHORecHits, TrackDetMatchInfo::crossedPreshowerIds, CaloRecHit::detid(), EcalRecHit::detid(), CaloRecHit::energy(), EcalRecHit::energy(), edm::EventSetup::get(), TrackDetectorAssociator::getCachedTrajector(), CachedTrajectory::getEcalTrajectory(), CaloSubdetectorGeometry::getGeometry(), CachedTrajectory::getHcalTrajectory(), CachedTrajectory::getHOTrajectory(), CachedTrajectory::getPreshowerTrajectory(), getTrackDetMatchInfo(), CaloCellGeometry::inside(), point, DetId::rawId(), mathSSE::sqrt(), trackAssociator_, trackAssociatorParameters_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by produce().

{
  TrackDetMatchInfo trackDetMatchInfo = getTrackDetMatchInfo(evt, es, trackAssociator_, trackAssociatorParameters_, muon);

  BOOST_FOREACH(const EcalRecHit * rh, trackDetMatchInfo.crossedEcalRecHits)
    depositMap[rh->detid().rawId()]+=rh->energy();
  
  BOOST_FOREACH(const HBHERecHit * rh, trackDetMatchInfo.crossedHcalRecHits)
    depositMap[rh->detid().rawId()]+=rh->energy();

  BOOST_FOREACH(const HORecHit * rh, trackDetMatchInfo.crossedHORecHits)
    depositMap[rh->detid().rawId()]+=rh->energy();

  typedef std::map<std::string, const std::vector<DetId>*> CaloToDetIdMap;
  CaloToDetIdMap caloToDetIdMap;
  caloToDetIdMap["ecal"] = &(trackDetMatchInfo.crossedEcalIds);
  caloToDetIdMap["hcal"] = &(trackDetMatchInfo.crossedHcalIds);
  caloToDetIdMap["ho"]   = &(trackDetMatchInfo.crossedHOIds);
  caloToDetIdMap["es"]   = &(trackDetMatchInfo.crossedPreshowerIds);

  edm::ESHandle<CaloGeometry> caloGeo;
  es.get<CaloGeometryRecord>().get(caloGeo); 

  for ( CaloToDetIdMap::const_iterator caloToDetIdEntry = caloToDetIdMap.begin();
    caloToDetIdEntry != caloToDetIdMap.end(); ++caloToDetIdEntry ) {
    std::vector<SteppingHelixStateInfo>::const_iterator itHelixState_first, itHelixState_last;
    if ( caloToDetIdEntry->first == "ecal" ) {
      itHelixState_first = trackAssociator_.getCachedTrajector().getEcalTrajectory().begin();
      itHelixState_last  = trackAssociator_.getCachedTrajector().getEcalTrajectory().end();
    } else if ( caloToDetIdEntry->first == "hcal" ) {
      itHelixState_first = trackAssociator_.getCachedTrajector().getHcalTrajectory().begin();
      itHelixState_last  = trackAssociator_.getCachedTrajector().getHcalTrajectory().end();
    } else if ( caloToDetIdEntry->first == "ho" ) {
      itHelixState_first = trackAssociator_.getCachedTrajector().getHOTrajectory().begin();
      itHelixState_last  = trackAssociator_.getCachedTrajector().getHOTrajectory().end();
    } else if ( caloToDetIdEntry->first == "es" ) {
      itHelixState_first = trackAssociator_.getCachedTrajector().getPreshowerTrajectory().begin();
      itHelixState_last  = trackAssociator_.getCachedTrajector().getPreshowerTrajectory().end();
    } else assert(0);
    
    // copy trajectory points
    std::vector<GlobalPoint> trajectory;
    for ( std::vector<SteppingHelixStateInfo>::const_iterator helixState = itHelixState_first;
      helixState != itHelixState_last; ++helixState ) {
      trajectory.push_back(helixState->position());
    }
    
    // iterate over crossed detIds
    for ( std::vector<DetId>::const_iterator detId = caloToDetIdEntry->second->begin();
      detId != caloToDetIdEntry->second->end(); ++detId ) {
      if ( detId->rawId() == 0 ) continue;
      
      const CaloSubdetectorGeometry* subDetGeo = caloGeo->getSubdetectorGeometry(*detId);
      const CaloCellGeometry* caloCellGeo = subDetGeo->getGeometry(*detId);
      GlobalPoint previousPoint;
      bool previousPoint_initialized;
      float distanceWithinDetId = 0;
      for ( std::vector<GlobalPoint>::const_iterator point = trajectory.begin();
        point != trajectory.end(); ++point ) {
    if ( previousPoint_initialized ) {
      float dx = point->x() - previousPoint.x();
      float dy = point->y() - previousPoint.y();
      float dz = point->z() - previousPoint.z();
      float distanceBetweenPoints = sqrt(dx*dx + dy*dy + dz*dz);
      int numSteps = 100;
      int numStepsWithinDetId = 0;
      for ( int iStep = 0; iStep <= numSteps; ++iStep ){
        float stepX = previousPoint.x() + iStep*dx/numSteps;
        float stepY = previousPoint.y() + iStep*dy/numSteps;
        float stepZ = previousPoint.z() + iStep*dz/numSteps;
        GlobalPoint stepPoint(stepX, stepY, stepZ);     
        bool isWithinDetId = caloCellGeo->inside(stepPoint);
        if ( isWithinDetId ) ++numStepsWithinDetId;
      }
      distanceWithinDetId += (numStepsWithinDetId/float(numSteps + 1))*distanceBetweenPoints;
    }
    previousPoint = (*point);
    previousPoint_initialized = true;
      } 
      distanceMap[detId->rawId()] = distanceWithinDetId;
    }
  }
}

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

privatevirtual

Implements edm::EDProducer.

Definition at line 45 of file MuonCaloDistanceProducer.cc.

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

{
  std::auto_ptr<detIdToFloatMap> distanceMuPlus(new detIdToFloatMap());
  std::auto_ptr<detIdToFloatMap> distanceMuMinus(new detIdToFloatMap());
  std::auto_ptr<detIdToFloatMap> depositMuPlus(new detIdToFloatMap());
  std::auto_ptr<detIdToFloatMap> depositMuMinus(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()  ) fillDistanceMap(evt, es, &(*muPlus), *distanceMuPlus, *depositMuPlus);
  if ( muMinus.isNonnull() ) fillDistanceMap(evt, es, &(*muMinus), *distanceMuMinus, *depositMuMinus);

  std::auto_ptr<reco::CandidateCollection> muons(new reco::CandidateCollection);
  if ( muPlus.isNonnull()  ) muons->push_back(new reco::ShallowCloneCandidate(muPlus));
  if ( muMinus.isNonnull() ) muons->push_back(new reco::ShallowCloneCandidate(muMinus));

  // References to the muons themselves
  evt.put(muons, "muons");

  // maps of detId to distance traversed by muon through calorimeter cell
  evt.put(distanceMuPlus, "distancesMuPlus");
  evt.put(distanceMuMinus, "distancesMuMinus");

  // maps of detId to energy deposited in calorimeter cell
  evt.put(depositMuPlus, "depositsMuPlus");
  evt.put(depositMuMinus, "depositsMuMinus");
}

Member Data Documentation

edm::InputTag MuonCaloDistanceProducer::srcSelectedMuons_

private

Definition at line 42 of file MuonCaloDistanceProducer.h.

Referenced by produce().

TrackDetectorAssociator MuonCaloDistanceProducer::trackAssociator_

private

Definition at line 44 of file MuonCaloDistanceProducer.h.

Referenced by fillDistanceMap(), and MuonCaloDistanceProducer().

TrackAssociatorParameters MuonCaloDistanceProducer::trackAssociatorParameters_

private

Definition at line 45 of file MuonCaloDistanceProducer.h.

Referenced by fillDistanceMap(), and MuonCaloDistanceProducer().