StartingLayerFinder.h

Go to the documentation of this file.

#ifndef TkNavigation_StartingLayerFinder_H_
#define TkNavigation_StartingLayerFinder_H_


#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/Event.h"


#include "TrackingTools/DetLayers/interface/DetLayer.h"
#include "TrackingTools/DetLayers/interface/BarrelDetLayer.h"
#include "TrackingTools/DetLayers/interface/ForwardDetLayer.h"
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
#include "TrackingTools/TrajectoryState/interface/FreeTrajectoryState.h"
#include "TrackingTools/GeomPropagators/interface/Propagator.h"
#include "DataFormats/TrajectorySeed/interface/TrajectorySeedCollection.h"
#include "RecoTracker/MeasurementDet/interface/MeasurementTracker.h"
#include "RecoTracker/TkDetLayers/interface/GeometricSearchTracker.h"

#include <vector>

class PTrajectoryStateOnDet;

class StartingLayerFinder {

private:

  typedef FreeTrajectoryState FTS;
  typedef TrajectoryStateOnSurface TSOS;
  typedef std::pair<float, float> Range;

public: 

  StartingLayerFinder(const Propagator* aPropagator, const MeasurementTracker*  tracker ) :

    thePropagator(aPropagator),
    theMeasurementTracker(tracker),
    thePixelLayersValid(false),
    theFirstPixelBarrelLayer(nullptr),
    theFirstNegPixelFwdLayer(0),
    theFirstPosPixelFwdLayer(0) { }

  ~StartingLayerFinder() {}

  std::vector<const DetLayer*> startingLayers(const FTS& aFts, float dr, float dz) const;
  

  std::vector<const DetLayer*> startingLayers(const TrajectorySeed& aSeed) const;

  const BarrelDetLayer* firstPixelBarrelLayer() const;
  const std::vector<const ForwardDetLayer*> firstNegPixelFwdLayer() const;
  const std::vector<const ForwardDetLayer*> firstPosPixelFwdLayer() const;

  const Propagator* propagator() const {return thePropagator;}



private:

  const Propagator* thePropagator;
  const MeasurementTracker*     theMeasurementTracker;
  mutable bool thePixelLayersValid;
  mutable const BarrelDetLayer* theFirstPixelBarrelLayer;
  mutable std::vector<const ForwardDetLayer*> theFirstNegPixelFwdLayer;
  mutable std::vector<const ForwardDetLayer*> theFirstPosPixelFwdLayer;


  void checkPixelLayers() const;
  
  
  
  inline bool rangesIntersect( const Range& a, const Range& b) const {
    if ( a.first > b.second || b.first > a.second) return false;
    else return true;
  }




};
#endif //TR_StartingLayerFinder_H_