PFCheckHitPattern Class Reference

PFCheckHitPatter. More...

#include <PFCheckHitPattern.h>

Public Types
typedef std::pair< PFTrackHitInfo, PFTrackHitInfo >	PFTrackHitFullInfo
typedef std::pair< unsigned int, unsigned int >	PFTrackHitInfo

Public Member Functions
PFTrackHitFullInfo	analyze (const TrackerTopology *tkerTopo, const TrackerGeometry *tkerGeom, const reco::TrackBaseRef track, const TransientVertex &vert)
	PFCheckHitPattern ()
void	print (const reco::TrackBaseRef track) const
	Print hit pattern on track. More...
	~PFCheckHitPattern ()

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

Private Member Functions
void	init (const TrackerTopology *, const TrackerGeometry *)
	Create map indicating r/z values of all layers/disks. More...
void	print (const reco::HitPattern::HitCategory, const reco::HitPattern &hp) const

Static Private Member Functions
static bool	barrel (uint32_t subDet)
	Return a bool indicating if a given subdetector is in the barrel. More...

Private Attributes
bool	geomInitDone_
	Note if geometry info is already initialized. More...
RZrangeMap	rangeRorZ_

Detailed Description

PFCheckHitPatter.

Author

Ian Tomalin, modified by Maxime Gouzevitch

Date

October 2009

Definition at line 34 of file PFCheckHitPattern.h.

Member Typedef Documentation

DetInfo

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

private

a pair<uint32, uint32> consisting of the numbers used by HitPattern to identify subdetector and layer number respectively.

Definition at line 62 of file PFCheckHitPattern.h.

PFTrackHitFullInfo

typedef std::pair<PFTrackHitInfo, PFTrackHitInfo> PFCheckHitPattern::PFTrackHitFullInfo

Definition at line 41 of file PFCheckHitPattern.h.

PFTrackHitInfo

typedef std::pair<unsigned int, unsigned int> PFCheckHitPattern::PFTrackHitInfo

Definition at line 40 of file PFCheckHitPattern.h.

RZrangeMap

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

private

For a given subdetector & layer number, this stores the minimum and maximum r (or z) values if it is barrel (or endcap) respectively.

Definition at line 75 of file PFCheckHitPattern.h.

Constructor & Destructor Documentation

PFCheckHitPattern()

PFCheckHitPattern::PFCheckHitPattern ( )

inline

Definition at line 36 of file PFCheckHitPattern.h.

: geomInitDone_(false) {}

~PFCheckHitPattern()

PFCheckHitPattern::~PFCheckHitPattern ( )

inline

Definition at line 38 of file PFCheckHitPattern.h.

{}

Member Function Documentation

analyze()

pair< PFCheckHitPattern::PFTrackHitInfo, PFCheckHitPattern::PFTrackHitInfo > PFCheckHitPattern::analyze	(	const TrackerTopology *	tkerTopo,
		const TrackerGeometry *	tkerGeom,
		const reco::TrackBaseRef	track,
		const TransientVertex &	vert
	)

PFCheck if hit pattern of this track is consistent with it being produced at given vertex. Pair.first gives number of hits on track in front of vertex. Pair.second gives number of missing hits between vertex and innermost hit on track.

Definition at line 77 of file PFCheckHitPattern.cc.

                                 {
  // PFCheck if hit pattern of this track is consistent with it being produced
  // at given vertex. Pair.first gives number of hits on track in front of vertex.
  // Pair.second gives number of missing hits between vertex and innermost hit
  // on track.
 
  // Initialise geometry info if not yet done.
  if (!geomInitDone_)
    this->init(tkerTopo, tkerGeom);
 
  // Get hit patterns of this track
  const reco::HitPattern& hp = track.get()->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;
  unsigned int nHitAfter = 0;
 
  for (int i = 0; i < hp.numberOfAllHits(HitPattern::TRACK_HITS); i++) {
    uint32_t hit = hp.getHitPattern(HitPattern::TRACK_HITS, i);
    if (hp.trackerHitFilter(hit) && hp.validHitFilter(hit)) {
      uint32_t subDet = hp.getSubStructure(hit);
      uint32_t layer = hp.getLayer(hit);
      DetInfo detInfo(subDet, layer);
      double maxRZ = rangeRorZ_[detInfo].second;
 
      if (this->barrel(subDet)) {
        if (vert.position().perp() > maxRZ)
          nHitBefore++;
        else
          nHitAfter++;
      } else {
        if (fabs(vert.position().z()) > maxRZ)
          nHitBefore++;
        else
          nHitAfter++;
      }
    }
  }
 
  // Count number of missing hits before the innermost hit on the track,
  // taking into account finite depth of each layer.
  unsigned int nMissHitAfter = 0;
  unsigned int nMissHitBefore = 0;
 
  for (int i = 0; i < hp.numberOfAllHits(HitPattern::MISSING_INNER_HITS); i++) {
    uint32_t hit = hp.getHitPattern(HitPattern::MISSING_INNER_HITS, i);
    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);
      double minRZ = rangeRorZ_[detInfo].first;
 
      //      cout << "subDet = " << subDet << " layer = " << layer << " minRZ = " << minRZ << endl;
 
      if (this->barrel(subDet)) {
        if (vert.position().perp() < minRZ)
          nMissHitAfter++;
        else
          nMissHitBefore++;
      } else {
        if (fabs(vert.position().z()) < minRZ)
          nMissHitAfter++;
        else
          nMissHitBefore++;
      }
    }
  }
 
  PFTrackHitInfo trackToVertex(nHitBefore, nMissHitBefore);
  PFTrackHitInfo trackFromVertex(nHitAfter, nMissHitAfter);
 
  return pair<PFTrackHitInfo, PFTrackHitInfo>(trackToVertex, trackFromVertex);
}

References Reference_intrackfit_cff::barrel, reco::HitPattern::getLayer(), reco::HitPattern::getSubStructure(), trackingPlots::hp, mps_fire::i, init, reco::HitPattern::MISSING_INNER_HITS, reco::HitPattern::missingHitFilter(), PV3DBase< T, PVType, FrameType >::perp(), TransientVertex::position(), HLT_FULL_cff::track, reco::HitPattern::TRACK_HITS, reco::HitPattern::trackerHitFilter(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by PFDisplacedVertexFinder::fitVertexFromSeed().

barrel()

bool PFCheckHitPattern::barrel ( uint32_t  subDet )

staticprivate

Return a bool indicating if a given subdetector is in the barrel.

Definition at line 71 of file PFCheckHitPattern.cc.

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

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

init()

void PFCheckHitPattern::init ( const TrackerTopology *  tkerTopo, const TrackerGeometry *  tkerGeom  )

private

Create map indicating r/z values of all layers/disks.

Definition at line 20 of file PFCheckHitPattern.cc.

                                                                                             {
  //
  // Note min/max radius (z) of each barrel layer (endcap disk).
  //
 
  geomInitDone_ = true;
 
  // Get Tracker geometry
  const TrackingGeometry::DetContainer& dets = tkerGeom->dets();
 
  // Loop over all modules in the Tracker.
  for (unsigned int i = 0; i < dets.size(); i++) {
    // Get subdet and layer of this module
    auto detId = dets[i]->geographicalId();
    auto detInfo = DetInfo(detId.subdetId(), tkerTopo->layer(detId));
    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] = pair<double, double>(minRZ, maxRZ);
  }
 
#if 0
  //def DEBUG_CHECKHITPATTERN
  RZrangeMap::const_iterator d;
  for (d = rangeRorZ_.begin(); d != rangeRorZ_.end(); d++) {
    DetInfo detInfo = d->first;
    pair<double, double> rangeRZ = d->second;
  }
#endif
}

References Reference_intrackfit_cff::barrel, ztail::d, TrackerGeometry::dets(), mps_fire::i, TrackerTopology::layer(), and position.

print() [1/2]

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

private

Definition at line 166 of file PFCheckHitPattern.cc.

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

References taus_updatedMVAIds_cff::category, gather_cfg::cout, trackingPlots::hp, and mps_fire::i.

print() [2/2]

void PFCheckHitPattern::print

(

const reco::TrackBaseRef

track

)

const

Print hit pattern on track.

Definition at line 156 of file PFCheckHitPattern.cc.

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

References gather_cfg::cout, trackingPlots::hp, reco::HitPattern::MISSING_INNER_HITS, print(), HLT_FULL_cff::track, and reco::HitPattern::TRACK_HITS.

Member Data Documentation

geomInitDone_

bool PFCheckHitPattern::geomInitDone_

private

Note if geometry info is already initialized.

Definition at line 71 of file PFCheckHitPattern.h.

rangeRorZ_

RZrangeMap PFCheckHitPattern::rangeRorZ_

private

Definition at line 76 of file PFCheckHitPattern.h.