/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_5/src/FastSimulation/EgammaElectronAlgos/src/FastPixelHitMatcher.h

Go to the documentation of this file.

#ifndef FastPixelHitMatcher_H
#define FastPixelHitMatcher_H

// -*- C++ -*-
//
// Package:    EgammaElectronAlgos
// Class:      FastPixelHitMatcher
//
//
// Original Author:  Patrick Janot.
//
//

#include "TrackingTools/TransientTrackingRecHit/interface/GenericTransientTrackingRecHit.h"
#include <vector>

class TrackerGeometry;
class GeometricSearchTracker;
class TrackerInteractionGeometry;
class TrackerLayer;
class TrackerRecHit;
class TrackingRecHit;
class ParticlePropagator;
class MagneticFieldMap;

class FastPixelHitMatcher{

 public:
  //RC
  typedef TransientTrackingRecHit::ConstRecHitPointer   ConstRecHitPointer;
  typedef TransientTrackingRecHit::RecHitPointer        RecHitPointer;
  typedef TransientTrackingRecHit::RecHitContainer      RecHitContainer;


  FastPixelHitMatcher(float,float,float,float,float,float,float,float,float,float,float,float,bool);

  virtual ~FastPixelHitMatcher();

  void setES(const MagneticFieldMap* aFieldMap,
             const TrackerGeometry* aTrackerGeometry,
             const GeometricSearchTracker* geomSearchTracker,
             const TrackerInteractionGeometry* interactionGeometry);

  std::vector< std::pair<ConstRecHitPointer,ConstRecHitPointer> >
    compatibleHits(const GlobalPoint& xmeas,
                   const GlobalPoint& vprim,
                   float energy,
                   std::vector<TrackerRecHit>& theHits);

  inline double getVertex() { return vertex; }

  void set1stLayer (float ephimin, float ephimax,
                    float pphimin, float pphimax) {

    ephi1min = ephimin;
    ephi1max = ephimax;
    pphi1min = pphimin;
    pphi1max = pphimax;

  }

  void set2ndLayer (float phimin, float phimax) {
    phi2min = phimin;
    phi2max = phimax;
  }

  bool isASeed(const ParticlePropagator& myElec,
               const ParticlePropagator& myPosi,
               const GlobalPoint& theVertex,
               double rCluster,
               double zCluster,
               const TrackerRecHit& hit1,
               const TrackerRecHit& hit2);

  bool propagateToLayer(ParticlePropagator& myPart,
                        const GlobalPoint& theVertex,
                        GlobalPoint& theHit,
                        double phimin,
                        double phimax,
                        unsigned layer);

  double zVertex(double zCluster,
                 double rCluster,
                 GlobalPoint& theHit);

  bool checkRZCompatibility(double zCluster,double rCluster,
                            double zVertex,
                            float rzMin, float rzMax,
                            GlobalPoint& theHit,
                            bool forward);


  void set1stLayerZRange(double zmin1, double zmax1) {
    z1max = zmax1;
    z1min = zmin1;
  }

 private:

  RecHitContainer hitsInTrack;

  float ephi1min, ephi1max;
  float pphi1min, pphi1max;
  float phi2min, phi2max;
  float z1min, z1max;
  float z2minB, z2maxB;
  float r2minF, r2maxF;
  float rMinI, rMaxI;
  bool searchInTIDTEC;

  const TrackerGeometry* theTrackerGeometry;
  const MagneticField* theMagneticField;
  const GeometricSearchTracker* theGeomSearchTracker;
  const TrackerInteractionGeometry* _theGeometry;
  const MagneticFieldMap* theFieldMap;
  std::vector<const TrackerLayer*> thePixelLayers;
  double vertex;


};

#endif