DTSegmentSelector Class Reference

#include <DTSegmentSelector.h>

Public Member Functions
	DTSegmentSelector (edm::ParameterSet const &pset, edm::ConsumesCollector &iC)
bool	operator() (DTRecSegment4D const &, edm::Event const &, edm::EventSetup const &)
	~DTSegmentSelector ()

Private Member Functions
bool	checkNoisySegment (edm::ESHandle< DTStatusFlag > const &, std::vector< DTRecHit1D > const &)

Private Attributes
bool	checkNoisyChannels_
double	maxAnglePhi_
double	maxAngleZ_
double	maxChi2_
int	minHitsPhi_
int	minHitsZ_
edm::InputTag	muonTags_
edm::EDGetTokenT< reco::MuonCollection >	muonToken_

Detailed Description

Definition at line 24 of file DTSegmentSelector.h.

Constructor & Destructor Documentation

DTSegmentSelector::DTSegmentSelector	(	edm::ParameterSet const &	pset,
		edm::ConsumesCollector &	iC
	)

Definition at line 16 of file DTSegmentSelector.cc.

References edm::ConsumesCollector::consumes(), muonTags_, and muonToken_.

    : muonTags_(pset.getParameter<edm::InputTag>("Muons")),
      checkNoisyChannels_(pset.getParameter<bool>("checkNoisyChannels")),
      minHitsPhi_(pset.getParameter<int>("minHitsPhi")),
      minHitsZ_(pset.getParameter<int>("minHitsZ")),
      maxChi2_(pset.getParameter<double>("maxChi2")),
      maxAnglePhi_(pset.getParameter<double>("maxAnglePhi")),
      maxAngleZ_(pset.getParameter<double>("maxAngleZ")) {
  muonToken_ = iC.consumes<reco::MuonCollection>(muonTags_);
}

DTSegmentSelector::~DTSegmentSelector ( )

inline

Definition at line 27 of file DTSegmentSelector.h.

References checkNoisySegment(), and operator()().

{}

Member Function Documentation

bool DTSegmentSelector::checkNoisySegment ( edm::ESHandle< DTStatusFlag > const &  statusMap, std::vector< DTRecHit1D > const &  dtHits  )

private

Definition at line 125 of file DTSegmentSelector.cc.

References DTStatusFlag::cellStatus(), and LogTrace.

Referenced by operator()(), and ~DTSegmentSelector().

                                                                               {
  bool segmentNoisy = false;

  std::vector<DTRecHit1D>::const_iterator dtHit = dtHits.begin();
  std::vector<DTRecHit1D>::const_iterator dtHits_end = dtHits.end();
  for (; dtHit != dtHits_end; ++dtHit) {
    DTWireId wireId = dtHit->wireId();
    // Check for noisy channels to skip them
    bool isNoisy = false, isFEMasked = false, isTDCMasked = false, isTrigMask = false, isDead = false, isNohv = false;
    statusMap->cellStatus(wireId, isNoisy, isFEMasked, isTDCMasked, isTrigMask, isDead, isNohv);
    if (isNoisy) {
      LogTrace("Calibration") << "Wire: " << wireId << " is noisy, skipping!";
      segmentNoisy = true;
      break;
    }
  }
  return segmentNoisy;
}

bool DTSegmentSelector::operator()	(	DTRecSegment4D const &	segment,
		edm::Event const &	event,
		edm::EventSetup const &	setup
	)

Definition at line 27 of file DTSegmentSelector.cc.

References checkNoisyChannels_, checkNoisySegment(), DTRecSegment4D::chi2(), simKBmtfDigis_cfi::chiSquare, DTRecSegment4D::degreesOfFreedom(), flavorHistoryFilter_cfi::dr, PVValHelper::dx, PVValHelper::dy, edm::EventSetup::get(), DTRecSegment4D::hasPhi(), DTRecSegment4D::hasZed(), edm::InputTag::label(), DTRecSegment4D::localDirection(), DTRecSegment4D::localPosition(), muonTagProbeFilters_cff::matched, maxAnglePhi_, maxAngleZ_, maxChi2_, minHitsPhi_, minHitsZ_, PDWG_BPHSkim_cff::muons, muonTags_, muonToken_, DTRecSegment4D::phiSegment(), Geom::pi(), edm::Handle< T >::product(), TrackingRecHit::rawId(), mps_fire::result, DTRecSegment2D::specificRecHits(), mathSSE::sqrt(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::z(), and DTRecSegment4D::zSegment().

Referenced by ~DTSegmentSelector().

                                                               {
  bool result = true;

  /* get the muon collection if one is specified
    (not specifying a muon collection switches off muon matching */
  if (!muonTags_.label().empty()) {
    edm::Handle<reco::MuonCollection> muonH;
    event.getByToken(muonToken_, muonH);
    const std::vector<reco::Muon>* muons = muonH.product();
    //std::cout << " Muon collection size: " << muons->size() << std::endl;
    if (muons->empty())
      return false;

    DTChamberId segId(segment.rawId());

    bool matched = false;
    for (auto& imuon : *muons)
      for (const auto& ch : imuon.matches()) {
        DetId chId(ch.id.rawId());
        if (chId != segId)
          continue;
        if (imuon.pt() < 15 || !imuon.isGlobalMuon())
          continue;

        int nsegs = ch.segmentMatches.size();
        if (!nsegs)
          continue;

        LocalPoint posHit = segment.localPosition();
        float dx = (posHit.x() ? posHit.x() - ch.x : 0);
        float dy = (posHit.y() ? posHit.y() - ch.y : 0);
        float dr = sqrt(dx * dx + dy * dy);
        if (dr < 5)
          matched = true;
      }

    if (!matched)
      result = false;
  }

  edm::ESHandle<DTStatusFlag> statusMap;
  if (checkNoisyChannels_)
    setup.get<DTStatusFlagRcd>().get(statusMap);

  // Get the Phi 2D segment
  int nPhiHits = -1;
  bool segmentNoisyPhi = false;
  if (segment.hasPhi()) {
    const DTChamberRecSegment2D* phiSeg = segment.phiSegment();  // phiSeg lives in the chamber RF
    std::vector<DTRecHit1D> phiRecHits = phiSeg->specificRecHits();
    nPhiHits = phiRecHits.size();
    if (checkNoisyChannels_)
      segmentNoisyPhi = checkNoisySegment(statusMap, phiRecHits);
    //} else result = false;
  }
  // Get the Theta 2D segment
  int nZHits = -1;
  bool segmentNoisyZ = false;
  if (segment.hasZed()) {
    const DTSLRecSegment2D* zSeg = segment.zSegment();  // zSeg lives in the SL RF
    std::vector<DTRecHit1D> zRecHits = zSeg->specificRecHits();
    nZHits = zRecHits.size();
    if (checkNoisyChannels_)
      segmentNoisyZ = checkNoisySegment(statusMap, zRecHits);
  }

  // Segment selection
  // Discard segment if it has a noisy cell
  if (segmentNoisyPhi || segmentNoisyZ)
    result = false;

  // 2D-segment number of hits
  if (segment.hasPhi() && nPhiHits < minHitsPhi_)
    result = false;

  if (segment.hasZed() && nZHits < minHitsZ_)
    result = false;

  // Segment chi2
  double chiSquare = segment.chi2() / segment.degreesOfFreedom();
  if (chiSquare > maxChi2_)
    result = false;

  // Segment angle
  LocalVector segment4DLocalDir = segment.localDirection();
  double angleZ = fabs(atan(segment4DLocalDir.y() / segment4DLocalDir.z()) * 180. / Geom::pi());
  if (angleZ > maxAngleZ_)
    result = false;

  double anglePhi = fabs(atan(segment4DLocalDir.x() / segment4DLocalDir.z()) * 180. / Geom::pi());
  if (anglePhi > maxAnglePhi_)
    result = false;

  return result;
}

Member Data Documentation

bool DTSegmentSelector::checkNoisyChannels_

private

Definition at line 35 of file DTSegmentSelector.h.

Referenced by operator()().

double DTSegmentSelector::maxAnglePhi_

private

Definition at line 39 of file DTSegmentSelector.h.

Referenced by operator()().

double DTSegmentSelector::maxAngleZ_

private

Definition at line 40 of file DTSegmentSelector.h.

Referenced by operator()().

double DTSegmentSelector::maxChi2_

private

Definition at line 38 of file DTSegmentSelector.h.

Referenced by operator()().

int DTSegmentSelector::minHitsPhi_

private

Definition at line 36 of file DTSegmentSelector.h.

Referenced by operator()().

int DTSegmentSelector::minHitsZ_

private

Definition at line 37 of file DTSegmentSelector.h.

Referenced by operator()().

edm::InputTag DTSegmentSelector::muonTags_

private

Definition at line 33 of file DTSegmentSelector.h.

Referenced by DTSegmentSelector(), and operator()().

edm::EDGetTokenT<reco::MuonCollection> DTSegmentSelector::muonToken_

private

Definition at line 34 of file DTSegmentSelector.h.

Referenced by DTSegmentSelector(), and operator()().