OutInConversionSeedFinder.h

Go to the documentation of this file.

#ifndef RecoEGAMMA_ConversionSeed_OutInConversionSeedFinder_h
#define RecoEGAMMA_ConversionSeed_OutInConversionSeedFinder_h
 
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/CaloRecHit/interface/CaloCluster.h"
#include "DataFormats/Common/interface/View.h"
#include "DataFormats/TrajectorySeed/interface/TrajectorySeedCollection.h"
#include "RecoEgamma/EgammaPhotonAlgos/interface/ConversionSeedFinder.h"
#include "RecoTracker/TkDetLayers/interface/GeometricSearchTracker.h"
 
#include "TrackingTools/DetLayers/interface/MeasurementEstimator.h"
 
#include "RecoTracker/MeasurementDet/interface/MeasurementTracker.h"
 
#include <string>
#include <vector>
 
class MagneticField;
class FreeTrajectoryState;
class TrajectoryStateOnSurface;
class LayerMeasurements;
 
class OutInConversionSeedFinder : public ConversionSeedFinder {
private:
  typedef FreeTrajectoryState FTS;
  typedef TrajectoryStateOnSurface TSOS;
 
public:
  OutInConversionSeedFinder(const edm::ParameterSet &config, edm::ConsumesCollector &&iC);
 
  ~OutInConversionSeedFinder() override;
 
  void makeSeeds(const edm::Handle<edm::View<reco::CaloCluster> > &allBc) override;
  virtual void makeSeeds(const reco::CaloClusterPtr &aBC);
 
private:
  edm::ParameterSet conf_;
  std::pair<FreeTrajectoryState, bool> makeTrackState(int charge) const;
 
  void fillClusterSeeds(const reco::CaloClusterPtr &bc);
 
  void startSeed(const FreeTrajectoryState &);
  void completeSeed(const TrajectoryMeasurement &m1, const FreeTrajectoryState &fts, const Propagator *, int layer);
  void createSeed(const TrajectoryMeasurement &m1, const TrajectoryMeasurement &m2);
  FreeTrajectoryState createSeedFTS(const TrajectoryMeasurement &m1, const TrajectoryMeasurement &m2) const;
  GlobalPoint fixPointRadius(const TrajectoryMeasurement &) const;
 
  MeasurementEstimator *makeEstimator(const DetLayer *, float dphi) const;
 
private:
  float the2ndHitdphi_;
  float the2ndHitdzConst_;
  float the2ndHitdznSigma_;
  std::vector<TrajectoryMeasurement> theFirstMeasurements_;
  int nSeedsPerBC_;
  int maxNumberOfOutInSeedsPerBC_;
  float bcEtcut_;
  float bcEcut_;
  bool useEtCut_;
};
 
#endif