HitQuadrupletGeneratorFromLayerPairForPhotonConversion.h

Go to the documentation of this file.

#ifndef HitQuadrupletGeneratorFromLayerPairForPhotonConversion_h
#define HitQuadrupletGeneratorFromLayerPairForPhotonConversion_h

#include "RecoTracker/TkHitPairs/interface/OrderedHitPairs.h"
#include "RecoTracker/TkHitPairs/interface/LayerHitMapCache.h"
#include "FWCore/Framework/interface/EventSetup.h"

#include "ConversionRegion.h"

class DetLayer;
class TrackingRegion;

class HitQuadrupletGeneratorFromLayerPairForPhotonConversion {

public:

  typedef LayerHitMapCache LayerCacheType;
  typedef SeedingLayerSetsHits::SeedingLayerSet Layers;
  typedef SeedingLayerSetsHits::SeedingLayer Layer;
 
  HitQuadrupletGeneratorFromLayerPairForPhotonConversion(unsigned int inner,
                                unsigned int outer,
                LayerCacheType* layerCache,
                unsigned int max=0);

  ~HitQuadrupletGeneratorFromLayerPairForPhotonConversion() { }

  void hitPairs( const TrackingRegion& reg, OrderedHitPairs & prs,
                 const Layers& layers,
                 const edm::Event & ev,  const edm::EventSetup& es);

  bool failCheckRZCompatibility(const RecHitsSortedInPhi::Hit & hit, const DetLayer& layer, const HitRZCompatibility *checkRZ, const TrackingRegion & region);
  //void checkPhiRange(double phi1, double phi2);

  bool failCheckSlopeTest(const RecHitsSortedInPhi::Hit & ohit, const RecHitsSortedInPhi::Hit & nohit, const RecHitsSortedInPhi::Hit & ihit, const RecHitsSortedInPhi::Hit & nihit, const TrackingRegion & region);
  void bubbleSortVsR(int n, double* ax, double* ay, double* aey);
  bool failCheckSegmentZCompatibility(double &rInn, double &zInnMin, double &zInnMax,
                      double &rInt, double &zIntMin, double &zIntMax,
                      double &rOut, double &zOutMin, double &zOutMax);
  double getZAtR(double &rInn, double &zInn, double &r, double &rOut, double &zOut);
  double verySimpleFit(int size, double* ax, double* ay, double* e2y, double& p0, double& e2p0, double& p1);
  double getSqrEffectiveErrorOnZ(const RecHitsSortedInPhi::Hit & hit, const TrackingRegion & region);
  double getEffectiveErrorOnZ(const RecHitsSortedInPhi::Hit & hit, const TrackingRegion & region);

private:
  
  LayerCacheType & theLayerCache;
  const unsigned int theOuterLayer;
  const unsigned int theInnerLayer;
  const unsigned int theMaxElement;

  std::stringstream *ss;

};

#endif