InOutConversionSeedFinder.h

Go to the documentation of this file.

#ifndef RecoEGAMMA_ConversionSeed_InOutConversionSeedFinder_h
#define RecoEGAMMA_ConversionSeed_InOutConversionSeedFinder_h

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/CaloRecHit/interface/CaloClusterFwd.h"
#include "DataFormats/Common/interface/View.h"

#include "DataFormats/TrajectorySeed/interface/TrajectorySeedCollection.h"
#include "DataFormats/TrackCandidate/interface/TrackCandidateCollection.h"

#include "RecoEgamma/EgammaPhotonAlgos/interface/ConversionSeedFinder.h"
#include "TrackingTools/DetLayers/interface/DetLayer.h"
#include "RecoTracker/TkDetLayers/interface/GeometricSearchTracker.h"
#include "RecoTracker/MeasurementDet/interface/MeasurementTracker.h"
#include "TrackingTools/DetLayers/interface/MeasurementEstimator.h"
#include "TrackingTools/PatternTools/interface/Trajectory.h"
#include "TrackingTools/MeasurementDet/interface/LayerMeasurements.h"
#include "TrackingTools/PatternTools/interface/TrajectoryMeasurement.h" 
#include "Geometry/CommonDetUnit/interface/TrackingGeometry.h"
#include <string>
#include <vector>



class MagneticField;
class FreeTrajectoryState;
class TrajectoryStateOnSurface;
class TrajectoryMeasurement;

class InOutConversionSeedFinder : public ConversionSeedFinder {
  
  
 private:
  
  typedef FreeTrajectoryState FTS;
  typedef TrajectoryStateOnSurface TSOS;
  
  
  public :
    
    
    
    InOutConversionSeedFinder( const edm::ParameterSet& config );
  
  
  
  virtual ~InOutConversionSeedFinder();
  
  
  virtual void  makeSeeds(  const edm::Handle<edm::View<reco::CaloCluster> > & allBc) const;
  
  
  
  void setTracks(std::vector<Trajectory> const & in) { theOutInTracks_ = in;}
  
  
  private :

  edm::ParameterSet conf_;
  virtual void fillClusterSeeds(  ) const ;
  void startSeed(FreeTrajectoryState * fts, const TrajectoryStateOnSurface & stateAtPreviousLayer, int charge, int layer) const ;
  virtual void findSeeds(const TrajectoryStateOnSurface & startingState,
             float signedpt, unsigned int startingLayer) const ;
  
  std::vector<const reco::CaloCluster*> getSecondCaloClusters(const GlobalPoint & conversionPosition, float charge) const;
  void completeSeed(const TrajectoryMeasurement & m1,FreeTrajectoryState & fts, const Propagator* propagator, int ilayer) const;
  void createSeed(const TrajectoryMeasurement & m1,  const TrajectoryMeasurement & m2) const ;


 private:
  float  the2ndHitdphi_;
  float  the2ndHitdzConst_;    
  float  the2ndHitdznSigma_;
  mutable int track2Charge_;
  mutable GlobalVector track2InitialMomentum_; 
  mutable int nSeedsPerInputTrack_;
  int maxNumberOfInOutSeedsPerInputTrack_;
     

  mutable TrajectoryMeasurement* myPointer;

  mutable std::vector<Trajectory> inputTracks_;
  mutable std::vector<Trajectory> theOutInTracks_;
  mutable std::vector<TrajectoryMeasurement> theFirstMeasurements_;

  mutable reco::CaloCluster theSecondBC_;
  mutable edm::Handle<edm::View<reco::CaloCluster> >  bcCollection_;

  
};

#endif