PixelForwardLayerPhase1.h

Go to the documentation of this file.

#ifndef TkDetLayers_PixelForwardLayerPhase1_h
#define TkDetLayers_PixelForwardLayerPhase1_h

#include "TrackingTools/DetLayers/interface/ForwardDetLayer.h"
#include "Phase1PixelBlade.h"
#include "Utilities/BinningTools/interface/PeriodicBinFinderInPhi.h"

#pragma GCC visibility push(hidden)

class PixelForwardLayerPhase1 final : public ForwardDetLayer {
public:
  PixelForwardLayerPhase1(std::vector<const Phase1PixelBlade*>& blades);
  ~PixelForwardLayerPhase1() override;

  // GeometricSearchDet interface

  const std::vector<const GeomDet*>& basicComponents() const override { return theBasicComps; }

  const std::vector<const GeometricSearchDet*>& components() const override __attribute__((cold)) { return theComps; }

  void groupedCompatibleDetsV(const TrajectoryStateOnSurface& tsos,
                              const Propagator& prop,
                              const MeasurementEstimator& est,
                              std::vector<DetGroup>& result) const override __attribute__((hot));

  // DetLayer interface
  SubDetector subDetector() const override { return GeomDetEnumerators::subDetGeom[GeomDetEnumerators::P1PXEC]; }

private:
  // methods for groupedCompatibleDets implementation
  static int computeHelicity(const GeometricSearchDet* firstBlade, const GeometricSearchDet* secondBlade);

  struct SubTurbineCrossings {
    SubTurbineCrossings() : isValid(false), closestIndex(0), nextIndex(0), nextDistance(0.0f) {}
    SubTurbineCrossings(int ci, int ni, float nd) : isValid(true), closestIndex(ci), nextIndex(ni), nextDistance(nd) {}

    bool isValid;
    int closestIndex;
    int nextIndex;
    float nextDistance;
  };

  void searchNeighbors(const TrajectoryStateOnSurface& tsos,
                       const Propagator& prop,
                       const MeasurementEstimator& est,
                       const SubTurbineCrossings& crossings,
                       float window,
                       std::vector<DetGroup>& result,
                       bool innerDisk) const __attribute__((hot));

  SubTurbineCrossings computeCrossings(const TrajectoryStateOnSurface& startingState,
                                       PropagationDirection propDir,
                                       bool innerDisk) const __attribute__((hot));

  static float computeWindowSize(const GeomDet* det,
                                 const TrajectoryStateOnSurface& tsos,
                                 const MeasurementEstimator& est);

private:
  typedef PeriodicBinFinderInPhi<float> BinFinderType;
  // need separate objects for inner and outer disk
  // or a smarter bin finder class
  BinFinderType theBinFinder_inner;
  BinFinderType theBinFinder_outer;
  unsigned int _num_innerpanels;
  unsigned int _num_outerpanels;

  std::vector<float> theBinFinder_byR;
  std::vector<unsigned int> theBinFinder_byR_index;
  std::vector<unsigned int> theBinFinder_byR_nextindex;
  // bool useR;
  std::vector<const GeometricSearchDet*> theComps;
  std::vector<const GeomDet*> theBasicComps;
};

#pragma GCC visibility pop
#endif