PFCheckHitPattern.h

Go to the documentation of this file.

#ifndef RecoParticleFlow_PFTracking_PFCheckHitPattern_H
#define RecoParticleFlow_PFTracking_PFCheckHitPattern_H

// standard EDAnalyser include files
#include <memory>
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "FWCore/Framework/interface/EventSetup.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "RecoVertex/VertexPrimitives/interface/TransientVertex.h"
#include "Geometry/CommonDetUnit/interface/TrackingGeometry.h"
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"

#define DEBUG_CHECKHITPATTERN

class DetId;


/*
 * Determine if a track has hits in front of its assumed production point.
 * Also determine if it misses hits between its assumed production point and its innermost hit.
 */

class PFCheckHitPattern {
public:
  PFCheckHitPattern() : geomInitDone_(false) {}

  ~PFCheckHitPattern() {}

  typedef std::pair<unsigned int, unsigned int> PFTrackHitInfo;
  typedef std::pair<PFTrackHitInfo, PFTrackHitInfo> PFTrackHitFullInfo;


  PFTrackHitFullInfo analyze(const TrackerTopology* tkerTopo,
                             const TrackerGeometry* tkerGeom,
                             const reco::TrackBaseRef track,
                             const TransientVertex& vert);

  void print(const reco::TrackBaseRef track) const;

private:
  void init(const TrackerTopology*, const TrackerGeometry*);

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

  static bool barrel(uint32_t subDet);

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

private:
  bool geomInitDone_;

  typedef std::map<DetInfo, std::pair<double, double> > RZrangeMap;
  RZrangeMap rangeRorZ_;
};

#endif