ConversionSeedFinder.h

Go to the documentation of this file.

#ifndef RecoEGAMMA_ConversionSeed_ConversionSeedFinder_h
#define RecoEGAMMA_ConversionSeed_ConversionSeedFinder_h

#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/TrajectorySeed/interface/PropagationDirection.h"
#include "DataFormats/TrajectorySeed/interface/TrajectorySeedCollection.h"
#include "DataFormats/CaloRecHit/interface/CaloCluster.h"
#include "DataFormats/CaloRecHit/interface/CaloClusterFwd.h"
#include "DataFormats/Common/interface/View.h"
#include "DataFormats/Common/interface/Handle.h"
#include "TrackingTools/MaterialEffects/interface/PropagatorWithMaterial.h"
#include "TrackingTools/GeomPropagators/interface/StraightLinePropagator.h"
#include "TrackingTools/KalmanUpdators/interface/KFUpdator.h"
#include "TrackingTools/TrajectoryState/interface/FreeTrajectoryState.h"
#include "TrackingTools/PatternTools/interface/TrajectoryMeasurement.h"
#include "RecoTracker/TkDetLayers/interface/GeometricSearchTracker.h"
#include "RecoTracker/MeasurementDet/interface/MeasurementTracker.h"
#include "RecoTracker/MeasurementDet/interface/MeasurementTrackerEvent.h"
#include "DataFormats/BeamSpot/interface/BeamSpot.h"

#include "DataFormats/Math/interface/Vector3D.h"
#include "DataFormats/Math/interface/Point3D.h"

#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "DataFormats/GeometryVector/interface/GlobalVector.h"

// C/C++ headers
#include <string>
#include <vector>

//
namespace edm {
  class ConsumesCollector;
}

class DetLayer;
class FreeTrajectoryState;
class TrajectoryStateOnSurface;
class NavigationSchool;

class ConversionSeedFinder {
public:
  ConversionSeedFinder();
  ConversionSeedFinder(const edm::ParameterSet& config, edm::ConsumesCollector& iC);

  virtual ~ConversionSeedFinder() {}

  virtual void makeSeeds(const edm::Handle<edm::View<reco::CaloCluster> >& allBc) = 0;

  TrajectorySeedCollection& seeds() { return theSeeds_; }
  virtual void setCandidate(float e, GlobalPoint pos) {
    theSCenergy_ = e;
    theSCPosition_ = pos;
  }
  std::vector<const DetLayer*> const& layerList() const { return theLayerList_; }

  void setMeasurementTracker(const MeasurementTracker* tracker) const { ; }
  const MeasurementTracker* getMeasurementTracker() const { return theMeasurementTracker_; }

  void setEventSetup(const edm::EventSetup& es);
  void setNavigationSchool(const NavigationSchool* navigation) { theNavigationSchool_ = navigation; }
  void setEvent(const edm::Event& e);

  void clear() { theSeeds_.clear(); }

protected:
  //edm::ParameterSet conf_; found this to be completely unused
  void findLayers();
  void findLayers(const FreeTrajectoryState& fts);

  FreeTrajectoryState trackStateFromClusters(int aCharge,
                                             const GlobalPoint& gpOrigine,
                                             PropagationDirection dir,
                                             float scaleFactor) const;

  void printLayer(int i) const;

  TrajectorySeedCollection theSeeds_;
  GlobalPoint theSCPosition_;

  std::string theMeasurementTrackerName_;
  const MeasurementTracker* theMeasurementTracker_;
  const TrackingGeometry* theTrackerGeom_;
  const NavigationSchool* theNavigationSchool_ = nullptr;

  edm::ESHandle<MagneticField> theMF_;
  edm::ESHandle<GeometricSearchTracker> theGeomSearchTracker_;

  edm::EDGetTokenT<reco::BeamSpot> beamSpotToken_;
  edm::EDGetTokenT<MeasurementTrackerEvent> measurementTrkToken_;

  KFUpdator theUpdator_;
  PropagationDirection dir_;
  reco::CaloCluster* theSC_;
  float theSCenergy_;

  std::vector<const DetLayer*> theLayerList_;

  GlobalPoint theBCPosition_;
  float theBCEnergy_;

  const Propagator* thePropagatorAlongMomentum_;
  const Propagator* thePropagatorOppositeToMomentum_;

  reco::BeamSpot theBeamSpot_;
  edm::Handle<MeasurementTrackerEvent> theTrackerData_;
};

#endif