PFDisplacedVertexFinder.h

Go to the documentation of this file.

#ifndef RecoParticleFlow_PFTracking_PFDisplacedVertexFinder_h
#define RecoParticleFlow_PFTracking_PFDisplacedVertexFinder_h 

#include "RecoParticleFlow/PFTracking/interface/PFDisplacedVertexHelper.h"
#include "FWCore/Framework/interface/ESHandle.h"

#include "DataFormats/ParticleFlowReco/interface/PFDisplacedVertexCandidate.h"
#include "DataFormats/ParticleFlowReco/interface/PFDisplacedVertexCandidateFwd.h"
#include "DataFormats/ParticleFlowReco/interface/PFDisplacedVertexSeed.h"
#include "DataFormats/ParticleFlowReco/interface/PFDisplacedVertexSeedFwd.h"
#include "DataFormats/ParticleFlowReco/interface/PFDisplacedVertex.h"
#include "DataFormats/ParticleFlowReco/interface/PFDisplacedVertexFwd.h"

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

#include "DataFormats/TrackReco/interface/TrackFwd.h"

#include "RecoParticleFlow/PFTracking/interface/PFCheckHitPattern.h"



class TrackingGeometry;
class TrackerGeometry;
class MagneticField;

class PFDisplacedVertexFinder {

 public:

  PFDisplacedVertexFinder();

  ~PFDisplacedVertexFinder();


  typedef std::set< reco::TrackBaseRef >::iterator IEset;
  typedef reco::PFDisplacedVertexCandidateCollection::iterator IDVC;
  typedef reco::PFDisplacedVertexSeedCollection::iterator IDVS;
  typedef reco::PFDisplacedVertexCollection::iterator IDV;

  typedef std::pair <unsigned int, unsigned int> PFTrackHitInfo;
  typedef std::pair <PFTrackHitInfo, PFTrackHitInfo> PFTrackHitFullInfo;

  enum FitterType {
    F_NOTDEFINED,
    F_DONOTREFIT,
    F_KALMAN,
    F_ADAPTIVE
  };



  void setParameters(double transvSize, double longSize, 
             double primaryVertexCut, double tobCut, 
             double tecCut, double minAdaptWeight, bool switchOff2TrackVertex) {
    transvSize_ = transvSize;
    longSize_   = longSize;
    primaryVertexCut_ = primaryVertexCut;
    tobCut_ = tobCut;
    tecCut_ = tecCut;
    minAdaptWeight_ = minAdaptWeight;
    switchOff2TrackVertex_ = switchOff2TrackVertex;
  }

  void setDebug( bool debug ) {debug_ = debug;}

  void setEdmParameters( const MagneticField* magField,
             edm::ESHandle<GlobalTrackingGeometry> globTkGeomHandle,
             edm::ESHandle<TrackerGeometry> tkerGeomHandle){ 
    magField_ = magField; 
    globTkGeomHandle_ = globTkGeomHandle;
    tkerGeomHandle_ = tkerGeomHandle; 
  }

  void setTracksSelector(edm::ParameterSet ps){
    helper_.setTracksSelector(ps);
  }

  void setVertexIdentifier(edm::ParameterSet ps){
    helper_.setVertexIdentifier(ps);
  }

  void setPrimaryVertex(edm::Handle< reco::VertexCollection > mainVertexHandle, 
            edm::Handle< reco::BeamSpot > beamSpotHandle){
    helper_.setPrimaryVertex(mainVertexHandle, beamSpotHandle);
  }

  void setAVFParameters(edm::ParameterSet ps){
      sigmacut_ = ps.getParameter<double>("sigmacut");
      t_ini_    = ps.getParameter<double>("Tini");
      ratio_    = ps.getParameter<double>("ratio");
  }

  void  setInput(const edm::Handle< reco::PFDisplacedVertexCandidateCollection >&); 
  
  
  std::auto_ptr< reco::PFDisplacedVertexCollection > transferDisplacedVertices() {return displacedVertices_;}

  const std::auto_ptr< reco::PFDisplacedVertexCollection >& displacedVertices() const {return displacedVertices_;}




  void findDisplacedVertices();


 private:
  

  void findSeedsFromCandidate(reco::PFDisplacedVertexCandidate&, reco::PFDisplacedVertexSeedCollection&);

  void mergeSeeds(reco::PFDisplacedVertexSeedCollection&, std::vector<bool>& bLocked);

  bool fitVertexFromSeed(reco::PFDisplacedVertexSeed&, reco::PFDisplacedVertex&);

  void selectAndLabelVertices(reco::PFDisplacedVertexCollection&,  std::vector <bool>&);

  bool rejectAndLabelVertex(reco::PFDisplacedVertex& dv);


  bool isCloseTo(const reco::PFDisplacedVertexSeed&, const reco::PFDisplacedVertexSeed&) const;

  double getTransvDiff(const GlobalPoint&, const GlobalPoint&) const;
  double getLongDiff(const GlobalPoint&, const GlobalPoint&) const;
  double getLongProj(const GlobalPoint&, const GlobalVector&) const;

  reco::PFDisplacedVertex::VertexTrackType getVertexTrackType(PFTrackHitFullInfo&) const;

  unsigned commonTracks(const reco::PFDisplacedVertex&, const reco::PFDisplacedVertex&) const;

  friend std::ostream& operator<<(std::ostream&, const PFDisplacedVertexFinder&);



  std::auto_ptr< reco::PFDisplacedVertexCandidateCollection >  displacedVertexCandidates_;
  std::auto_ptr< reco::PFDisplacedVertexCollection >    displacedVertices_;



  double transvSize_;
  double longSize_;
  double primaryVertexCut_;
  double tobCut_;
  double tecCut_;
  double minAdaptWeight_;

  bool switchOff2TrackVertex_;

  
  double sigmacut_; //= 6;
  double t_ini_; //= 256.;
  double ratio_; //= 0.25;


  bool   debug_;
  
  edm::ESHandle<GlobalTrackingGeometry> globTkGeomHandle_;

  edm::ESHandle<TrackerGeometry> tkerGeomHandle_;

  const MagneticField* magField_;

  
  PFCheckHitPattern hitPattern_;

  PFDisplacedVertexHelper helper_;

};

#endif