RectangularEtaPhiTrackingRegion.h

Go to the documentation of this file.

#ifndef RectangularEtaPhiTrackingRegion_H
#define RectangularEtaPhiTrackingRegion_H

#include "RecoTracker/TkTrackingRegions/interface/TrackingRegion.h"
#include "RecoTracker/TkTrackingRegions/interface/TkTrackingRegionsMargin.h"
//#include "CommonDet/TrajectoryParametrization/interface/GlobalTrajectoryParameters.h"
#include "RecoTracker/TkTrackingRegions/interface/HitRZConstraint.h"
#include "RecoTracker/TkTrackingRegions/interface/OuterHitPhiPrediction.h"
#include "FWCore/Framework/interface/EventSetup.h"
class OuterEstimator;
class BarrelDetLayer;
class ForwardDetLayer;

class RectangularEtaPhiTrackingRegion GCC11_FINAL : public TrackingRegion {
public:

  typedef TkTrackingRegionsMargin<float> Margin;

  RectangularEtaPhiTrackingRegion( const GlobalVector & dir, 
                   const GlobalPoint & vertexPos,
                                   float ptMin, float rVertex, float zVertex,
                                   float deltaEta, float deltaPhi,
                                   float whereToUseMeasurementTracker = 0.,
                                   bool precise = true,
                                   const std::string & measurementTrackerName = "",
                   bool etaPhiRegion=false) 
    : TrackingRegionBase( dir, vertexPos, Range( -1/ptMin, 1/ptMin), 
              rVertex, zVertex),
    thePhiMargin( Margin( fabs(deltaPhi),fabs(deltaPhi))),
    theMeasurementTrackerUsage(whereToUseMeasurementTracker), thePrecise(precise), theMeasurementTrackerName(measurementTrackerName),
    theUseEtaPhi(etaPhiRegion)
   { initEtaRange(dir, Margin( fabs(deltaEta),fabs(deltaEta))); }
 
  RectangularEtaPhiTrackingRegion( const GlobalVector & dir, 
                               const GlobalPoint & vertexPos,
                                   float ptMin, float rVertex, float zVertex,
                                   Margin etaMargin,
                                   Margin phiMargin,
                                   float whereToUseMeasurementTracker = 0.,
                   bool precise = true, 
                                   const std::string & measurementTrackerName = "",
                   bool etaPhiRegion=false) 
    : TrackingRegionBase( dir, vertexPos, Range( -1/ptMin, 1/ptMin), 
      rVertex, zVertex), thePhiMargin( phiMargin), theMeasurementTrackerUsage(whereToUseMeasurementTracker), thePrecise(precise),
      theMeasurementTrackerName(measurementTrackerName),
      theUseEtaPhi(etaPhiRegion)
    { initEtaRange(dir, etaMargin); }

  RectangularEtaPhiTrackingRegion( const GlobalVector & dir, 
                               const GlobalPoint & vertexPos,
                                   Range invPtRange, 
                                   float rVertex, float zVertex,
                                   Margin etaMargin,
                                   Margin phiMargin,
                                   float whereToUseMeasurementTracker = 0.,
                                   bool precise = true,
                                   const std::string & measurementTrackerName = "",
                   bool etaPhiRegion=false) 
    : TrackingRegionBase( dir, vertexPos, invPtRange, rVertex, zVertex),
      thePhiMargin( phiMargin), theMeasurementTrackerUsage(whereToUseMeasurementTracker), thePrecise(precise),
      theMeasurementTrackerName(measurementTrackerName),
      theUseEtaPhi(etaPhiRegion)
    { initEtaRange(dir, etaMargin); }


  const Range & etaRange() const { return theEtaRange; }

  const Margin & phiMargin() const { return thePhiMargin; }

  bool  isPrecise() const { return thePrecise; }

  virtual TrackingRegion::Hits hits(
      const edm::Event& ev,  
      const edm::EventSetup& es, 
      const ctfseeding::SeedingLayer* layer) const;

  virtual HitRZCompatibility * checkRZ(const DetLayer* layer,  
                       const Hit &  outerHit,
                       const edm::EventSetup& iSetup,
                       const DetLayer* outerlayer=0,
                       float lr=0, float gz=0, float dr=0, float dz=0) const
  { return checkRZOld(layer,outerHit->hit(),iSetup); }

  virtual RectangularEtaPhiTrackingRegion* clone() const { 
    return new RectangularEtaPhiTrackingRegion(*this);
  }

  virtual std::string name() const { return "RectangularEtaPhiTrackingRegion"; }
  virtual std::string print() const;

private:
  HitRZCompatibility* checkRZOld(
      const DetLayer* layer, 
      const TrackingRecHit*  outerHit,
      const edm::EventSetup& iSetup) const;

  OuterEstimator * estimator(const BarrelDetLayer* layer,const edm::EventSetup& iSetup) const;
  OuterEstimator * estimator(const ForwardDetLayer* layer,const edm::EventSetup& iSetup) const;

  OuterHitPhiPrediction phiWindow(const edm::EventSetup& iSetup) const;
  HitRZConstraint rzConstraint() const;

  void  initEtaRange( const GlobalVector & dir, const Margin& margin);

private:

  Range theEtaRange;
  Margin thePhiMargin;
  double theMeasurementTrackerUsage;
  bool thePrecise;
  std::string theMeasurementTrackerName;
  bool theUseEtaPhi;
};

#endif