CheckHitPattern Class Reference

#include <CheckHitPattern.h>

Classes
struct	Result

Public Types
typedef std::pair< uint32_t, uint32_t >	DetInfo

Public Member Functions
void	init (const edm::EventSetup &iSetup)
Result	operator() (const reco::Track &track, const VertexState &vert) const

Static Public Member Functions
static bool	barrel (uint32_t subDet)
static DetInfo	interpretDetId (DetId detId, const TrackerTopology *tTopo)
static void	print (const reco::Track &track)
static void	print (const reco::HitPattern::HitCategory category, const reco::HitPattern &hp)

Private Types
typedef std::map< DetInfo, std::pair< double, double > >	RZrangeMap

Private Attributes
bool	geomInitDone_ = false
RZrangeMap	rangeRorZ_
edm::ESHandle< TransientTrackBuilder >	trkTool_

Detailed Description

Definition at line 25 of file CheckHitPattern.h.

Member Typedef Documentation

DetInfo

typedef std::pair<uint32_t, uint32_t> CheckHitPattern::DetInfo

Definition at line 46 of file CheckHitPattern.h.

RZrangeMap

typedef std::map<DetInfo, std::pair<double, double> > CheckHitPattern::RZrangeMap

private

Definition at line 60 of file CheckHitPattern.h.

Member Function Documentation

barrel()

bool CheckHitPattern::barrel ( uint32_t  subDet )

static

Definition at line 88 of file CheckHitPattern.cc.

References PixelSubdetector::PixelBarrel, StripSubdetector::TIB, and StripSubdetector::TOB.

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

                                            {
  // Determines if given sub-detector is in the barrel.
  return (subDet == StripSubdetector::TIB || subDet == StripSubdetector::TOB ||
          subDet == PixelSubdetector::PixelBarrel);
}

init()

void CheckHitPattern::init

(

const edm::EventSetup &

iSetup

)

Definition at line 21 of file CheckHitPattern.cc.

References barrel(), gather_cfg::cout, ztail::d, TrackerGeometry::dets(), geomInitDone_, edm::EventSetup::get(), mps_fire::i, interpretDetId(), position, edm::ESHandle< T >::product(), rangeRorZ_, trkTool_, and z.

                                                      {
  //Retrieve tracker topology from geometry
  edm::ESHandle<TrackerTopology> tTopoHandle;
  iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
  const TrackerTopology* const tTopo = tTopoHandle.product();

  iSetup.get<TransientTrackRecord>().get("TransientTrackBuilder", trkTool_);  // Needed for vertex fits

  //
  // Note min/max radius (z) of each barrel layer (endcap disk).
  //

  geomInitDone_ = true;

  // Get Tracker geometry
  edm::ESHandle<TrackerGeometry> trackerGeometry;
  iSetup.get<TrackerDigiGeometryRecord>().get(trackerGeometry);
  const TrackingGeometry::DetContainer& dets = trackerGeometry->dets();

  // Loop over all modules in the Tracker.
  for (unsigned int i = 0; i < dets.size(); i++) {
    // Get subdet and layer of this module
    DetInfo detInfo = this->interpretDetId(dets[i]->geographicalId(), tTopo);
    uint32_t subDet = detInfo.first;

    // Note r (or z) of module if barrel (or endcap).
    double r_or_z;
    if (this->barrel(subDet)) {
      r_or_z = dets[i]->position().perp();
    } else {
      r_or_z = fabs(dets[i]->position().z());
    }

    // Recover min/max r/z value of this layer/disk found so far.
    double minRZ = 999.;
    double maxRZ = 0.;
    if (rangeRorZ_.find(detInfo) != rangeRorZ_.end()) {
      minRZ = rangeRorZ_[detInfo].first;
      maxRZ = rangeRorZ_[detInfo].second;
    }

    // Update with those of this module if it exceeds them.
    if (minRZ > r_or_z)
      minRZ = r_or_z;
    if (maxRZ < r_or_z)
      maxRZ = r_or_z;
    rangeRorZ_[detInfo] = std::pair<double, double>(minRZ, maxRZ);
  }

#ifdef DEBUG_CHECKHITPATTERN
  RZrangeMap::const_iterator d;
  for (d = rangeRorZ_.begin(); d != rangeRorZ_.end(); d++) {
    DetInfo detInfo = d->first;
    std::pair<double, double> rangeRZ = d->second;
    std::std::cout << "CHECKHITPATTERN: Tracker subdetector type=" << detInfo.first << " layer=" << detInfo.second
                   << " has min r (or z) =" << rangeRZ.first << " and max r (or z) = " << rangeRZ.second
                   << std::std::endl;
  }
#endif
}

interpretDetId()

CheckHitPattern::DetInfo CheckHitPattern::interpretDetId ( DetId  detId, const TrackerTopology *  tTopo  )

static

Definition at line 82 of file CheckHitPattern.cc.

References TrackerTopology::layer(), and DetId::subdetId().

Referenced by init().

                                                                                                {
  // Convert detId to a pair<uint32, uint32> consisting of the numbers used by HitPattern
  // to identify subdetector and layer number respectively.
  return DetInfo(detId.subdetId(), tTopo->layer(detId));
}

operator()()

CheckHitPattern::Result CheckHitPattern::operator()	(	const reco::Track &	track,
		const VertexState &	vert
	)		const

Definition at line 94 of file CheckHitPattern.cc.

References barrel(), TransientTrackBuilder::build(), Exception, geomInitDone_, reco::HitPattern::getLayer(), reco::HitPattern::getSubStructure(), trackingPlots::hp, mps_fire::i, if(), phase1PixelTopology::layer, LogDebug, M_PI, reco::HitPattern::MISSING_INNER_HITS, reco::HitPattern::missingHitFilter(), PV3DBase< T, PVType, FrameType >::perp(), VertexState::position(), rangeRorZ_, mps_fire::result, edm::second(), HLT_2022v12_cff::track, reco::HitPattern::TRACK_HITS, reco::HitPattern::trackerHitFilter(), reco::TransientTrack::trajectoryStateClosestToPoint(), trkTool_, reco::HitPattern::validHitFilter(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                                       {
  // Check if hit pattern of this track is consistent with it being produced
  // at given vertex.

  // Initialise geometry info if not yet done.
  if (!geomInitDone_)
    throw cms::Exception("CheckHitPattern::operator() called before CheckHitPattern::init");

  // Optionally set vertex position to zero for debugging.
  // VertexState vertDebug( GlobalPoint(0.,0.,0.) , GlobalError(1e-8, 0., 1e-8, 0., 0., 1e-8) );

  // Evaluate track parameters at vertex.
  reco::TransientTrack t_trk = trkTool_->build(track);
  GlobalVector p3_trk = t_trk.trajectoryStateClosestToPoint(vert.position()).momentum();
  bool trkGoesInsideOut =
      fabs(reco::deltaPhi<const GlobalVector, const GlobalPoint>(p3_trk, vert.position())) < 0.5 * M_PI;

  LogDebug("CHP") << "TRACK: in-->out ? " << trkGoesInsideOut << " dxy=" << track.dxy() << " sz=" << track.dz()
                  << " eta=" << track.eta() << " barrel hits=" << track.hitPattern().numberOfValidPixelHits() << "/"
                  << track.hitPattern().numberOfValidStripTIBHits() << "/"
                  << track.hitPattern().numberOfValidStripTOBHits();
  LogDebug("CHP") << "VERT: r=" << vert.position().perp() << " z=" << vert.position().z();
  //  if (vert.position().perp() < 3.5 && fabs(vert.position().z()) < 10. && fabs(track.eta()) < 1 && fabs(track.dxy()) < 2 && fabs(track.dz()) < 2 && track.hitPattern().numberOfValidPixelHits() == 0 && track.hitPattern().numberOfValidStripTIBHits() == 0) LogDebug("CHP")<<"LOOKATTHISTRACK";
  // Get hit patterns of this track
  const reco::HitPattern& hp = track.hitPattern();

  // Count number of valid hits on track definately in front of the vertex,
  // taking into account finite depth of each layer.
  unsigned int nHitBefore = 0;
  for (int i = 0; i < hp.numberOfAllHits(reco::HitPattern::TRACK_HITS); i++) {
    uint32_t hit = hp.getHitPattern(reco::HitPattern::TRACK_HITS, i);
    if (reco::HitPattern::trackerHitFilter(hit) && reco::HitPattern::validHitFilter(hit)) {
      uint32_t subDet = reco::HitPattern::getSubStructure(hit);
      uint32_t layer = reco::HitPattern::getLayer(hit);
      DetInfo detInfo(subDet, layer);
      auto maxRZ = (*rangeRorZ_.find(detInfo)).second.second;

      if (this->barrel(subDet)) {
        // Be careful. If the track starts by going outside-->in, it is allowed to have hits before the vertex !
        if (vert.position().perp() > maxRZ && trkGoesInsideOut)
          nHitBefore++;
      } else {
        if (fabs(vert.position().z()) > maxRZ)
          nHitBefore++;
      }
    }
  }

  // Count number of missing hits before the innermost hit on the track,
  // taking into account finite depth of each layer.
  unsigned int nMissHitAfter = 0;
  for (int i = 0; i < hp.numberOfAllHits(reco::HitPattern::MISSING_INNER_HITS); i++) {
    uint32_t hit = hp.getHitPattern(reco::HitPattern::MISSING_INNER_HITS, i);
    //    if (hp.trackerHitFilter(hit)) {
    if (reco::HitPattern::trackerHitFilter(hit) && reco::HitPattern::missingHitFilter(hit)) {
      uint32_t subDet = reco::HitPattern::getSubStructure(hit);
      uint32_t layer = reco::HitPattern::getLayer(hit);
      DetInfo detInfo(subDet, layer);
      auto minRZ = (*rangeRorZ_.find(detInfo)).second.first;

      if (this->barrel(subDet)) {
        // Be careful. If the track starts by going outside-->in, then it misses hits
        // in all layers it crosses  before its innermost valid hit.
        if (vert.position().perp() < minRZ || !trkGoesInsideOut)
          nMissHitAfter++;
      } else {
        if (fabs(vert.position().z()) < minRZ)
          nMissHitAfter++;
      }
    }
  }

  Result result;
  result.hitsInFrontOfVert = nHitBefore;
  result.missHitsAfterVert = nMissHitAfter;
  return result;
}

print() [1/2]

void CheckHitPattern::print ( const reco::Track &  track )

static

Definition at line 172 of file CheckHitPattern.cc.

References gather_cfg::cout, trackingPlots::hp, reco::HitPattern::MISSING_INNER_HITS, HLT_2022v12_cff::track, and reco::HitPattern::TRACK_HITS.

                                                  {
  // Get hit patterns of this track
  const reco::HitPattern& hp = track.hitPattern();
  std::cout << "=== Hits on Track ===" << std::endl;
  print(reco::HitPattern::TRACK_HITS, hp);
  std::cout << "=== Hits before track ===" << std::endl;
  print(reco::HitPattern::MISSING_INNER_HITS, hp);
}

print() [2/2]

void CheckHitPattern::print ( const reco::HitPattern::HitCategory  category, const reco::HitPattern &  hp  )

static

Definition at line 181 of file CheckHitPattern.cc.

References taus_updatedMVAIds_cff::category, gather_cfg::cout, reco::HitPattern::getLayer(), reco::HitPattern::getSubStructure(), trackingPlots::hp, mps_fire::i, phase1PixelTopology::layer, and reco::HitPattern::trackerHitFilter().

                                                                                              {
  for (int i = 0; i < hp.numberOfAllHits(category); i++) {
    uint32_t hit = hp.getHitPattern(category, i);
    if (reco::HitPattern::trackerHitFilter(hit)) {
      uint32_t subdet = reco::HitPattern::getSubStructure(hit);
      uint32_t layer = reco::HitPattern::getLayer(hit);
      std::cout << "hit " << i << " subdet=" << subdet << " layer=" << layer << " type " << hp.getHitType(hit)
                << std::endl;
    }
  }
}

Member Data Documentation

geomInitDone_

bool CheckHitPattern::geomInitDone_ = false

private

Definition at line 56 of file CheckHitPattern.h.

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

rangeRorZ_

RZrangeMap CheckHitPattern::rangeRorZ_

private

Definition at line 61 of file CheckHitPattern.h.

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

trkTool_

edm::ESHandle<TransientTrackBuilder> CheckHitPattern::trkTool_

private

Definition at line 64 of file CheckHitPattern.h.

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