/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_14/src/RecoParticleFlow/PFTracking/interface/PFCheckHitPattern.h

Go to the documentation of this file.

#ifndef PFCheckHitPattern_H
#define PFCheckHitPattern_H

// standard EDAnalyser include files
#include <memory>
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.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(edm::ESHandle<TrackerGeometry>, const reco::TrackBaseRef track, 
            const TransientVertex& vert);

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



private:
  void init (edm::ESHandle<TrackerGeometry>);

  typedef std::pair<uint32_t, uint32_t> DetInfo;
  static DetInfo interpretDetId(DetId detId);

  static bool barrel(uint32_t subDet);

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

private:
  bool geomInitDone_;

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


};

#endif