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

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.

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

Definition at line 41 of file PFCheckHitPattern.h.

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

Definition at line 40 of file PFCheckHitPattern.h.

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 ( )

inline

Definition at line 36 of file PFCheckHitPattern.h.

: geomInitDone_(false) {}

PFCheckHitPattern::~PFCheckHitPattern ( )

inline

Definition at line 38 of file PFCheckHitPattern.h.

{}

Member Function Documentation

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.

References Reference_intrackfit_cff::barrel, edm::RefToBase< T >::get(), reco::HitPattern::getHitPattern(), reco::HitPattern::getLayer(), reco::HitPattern::getSubStructure(), reco::TrackBase::hitPattern(), AnalysisDataFormats_SUSYBSMObjects::hp, mps_fire::i, init, reco::HitPattern::MISSING_INNER_HITS, reco::HitPattern::missingHitFilter(), reco::HitPattern::numberOfAllHits(), PV3DBase< T, PVType, FrameType >::perp(), TransientVertex::position(), reco::HitPattern::TRACK_HITS, reco::HitPattern::trackerHitFilter(), reco::HitPattern::validHitFilter(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by PFDisplacedVertexFinder::fitVertexFromSeed().

                                 {
  // 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);
}

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.

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

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

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.

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

                                                                                             {
  //
  // 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
}

void PFCheckHitPattern::print

(

const reco::TrackBaseRef

track

)

const

Print hit pattern on track.

Definition at line 156 of file PFCheckHitPattern.cc.

References gather_cfg::cout, edm::RefToBase< T >::get(), reco::TrackBase::hitPattern(), AnalysisDataFormats_SUSYBSMObjects::hp, reco::HitPattern::MISSING_INNER_HITS, edm::print(), and reco::HitPattern::TRACK_HITS.

                                                            {
  // 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);
}

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

private

Definition at line 166 of file PFCheckHitPattern.cc.

References gather_cfg::cout, reco::HitPattern::getHitPattern(), reco::HitPattern::getHitType(), reco::HitPattern::getLayer(), reco::HitPattern::getSubStructure(), mps_fire::i, reco::HitPattern::numberOfAllHits(), and reco::HitPattern::trackerHitFilter().

                                                                                                      {
  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;
    }
  }
}

Member Data Documentation

bool PFCheckHitPattern::geomInitDone_

private

Note if geometry info is already initialized.

Definition at line 71 of file PFCheckHitPattern.h.

RZrangeMap PFCheckHitPattern::rangeRorZ_

private

Definition at line 76 of file PFCheckHitPattern.h.