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