ValidHitPairFilter.h

Go to the documentation of this file.

#ifndef _ValidHitPairFilter_h_
#define _ValidHitPairFilter_h_

#include <vector>

#include "RecoPixelVertexing/PixelTrackFitting/interface/PixelTrackFilterBase.h"

namespace edm { class EventSetup; }
class TrackingRecHit;
class Track;
class FreeTrajectoryState;
class TrackerGeometry;
class GeometricSearchTracker;
class MagneticField;
class Propagator;
class DetLayer;
class GeomDet;
class TrackerTopology;

class ValidHitPairFilter : public PixelTrackFilterBase
{
public:
  ValidHitPairFilter(const edm::EventSetup& es);
  ~ValidHitPairFilter() override;
  bool operator()(const reco::Track * track,
                          const std::vector<const TrackingRecHit *>& recHits) const override;

private:
  int getLayer(const TrackingRecHit & recHit) const;
  std::vector<int> getMissingLayers(int a, int b) const;
  FreeTrajectoryState getTrajectory(const reco::Track & track) const;
  std::vector<const GeomDet *> getCloseDets
    (int il, float rz, const std::vector<float>& rzB,
     float ph, const std::vector<float>& phB) const;

  const TrackerGeometry * theTracker;
  const GeometricSearchTracker * theGSTracker;
  const MagneticField * theMagneticField;
  const Propagator*    thePropagator;
  const TrackerTopology *tTopo;

  const std::vector<DetLayer *> detLayers;
  std::vector<float> rzBounds[7];
  std::vector<float> phBounds[7];
};

#endif