BPHDecayVertex.h

Go to the documentation of this file.

#ifndef HeavyFlavorAnalysis_RecoDecay_BPHDecayVertex_h
#define HeavyFlavorAnalysis_RecoDecay_BPHDecayVertex_h

//----------------------
// Base Class Headers --
//----------------------
#include "HeavyFlavorAnalysis/RecoDecay/interface/BPHDecayMomentum.h"

namespace edm {
  class EventSetup;
}

namespace reco {
  class TransientTrack;
  class Vertex;
}  // namespace reco

//------------------------------------
// Collaborating Class Declarations --
//------------------------------------
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "RecoVertex/VertexPrimitives/interface/VertexFitter.h"

//---------------
// C++ Headers --
//---------------
#include <vector>
#include <map>
#include <string>

//              ---------------------
//              -- Class Interface --
//              ---------------------

class BPHDecayVertex : public virtual BPHDecayMomentum {
public:
  // deleted copy constructor and assignment operator
  BPHDecayVertex(const BPHDecayVertex& x) = delete;
  BPHDecayVertex& operator=(const BPHDecayVertex& x) = delete;

  ~BPHDecayVertex() override;

  virtual bool validTracks() const;
  virtual bool validVertex() const;

  virtual const reco::Vertex& vertex(VertexFitter<5>* fitter = nullptr,
                                     const reco::BeamSpot* bs = nullptr,
                                     const GlobalPoint* priorPos = nullptr,
                                     const GlobalError* priorError = nullptr) const;

  const std::vector<const reco::Track*>& tracks() const;

  const reco::Track* getTrack(const reco::Candidate* cand) const;

  const std::vector<reco::TransientTrack>& transientTracks() const;

  reco::TransientTrack* getTransientTrack(const reco::Candidate* cand) const;

  const edm::EventSetup* getEventSetup() const;

  const std::string& getTrackSearchList(const reco::Candidate* cand) const;

protected:
  // constructor
  BPHDecayVertex(const edm::EventSetup* es);
  // pointer used to retrieve informations from other bases
  BPHDecayVertex(const BPHDecayVertex* ptr, const edm::EventSetup* es);

  virtual void addV(const std::string& name, const reco::Candidate* daug, const std::string& searchList, double mass);
  virtual void addV(const std::string& name, const BPHRecoConstCandPtr& comp);

  // utility function used to cash reconstruction results
  void setNotUpdated() const override;

private:
  // EventSetup needed to build TransientTrack
  const edm::EventSetup* evSetup;

  // map linking particles to associated track search list
  std::map<const reco::Candidate*, std::string> searchMap;

  // reconstruction results cache
  mutable bool oldTracks;
  mutable bool oldVertex;
  mutable bool validTks;
  mutable std::vector<const reco::Track*> rTracks;
  mutable std::vector<reco::TransientTrack> trTracks;
  mutable std::map<const reco::Candidate*, const reco::Track*> tkMap;
  mutable std::map<const reco::Candidate*, reco::TransientTrack*> ttMap;
  mutable reco::Vertex fittedVertex;
  mutable VertexFitter<5>* savedFitter;
  mutable reco::BeamSpot const* savedBS;
  mutable GlobalPoint const* savedPP;
  mutable GlobalError const* savedPE;

  // create TransientTrack and fit vertex
  virtual void tTracks() const;
  virtual void fitVertex(VertexFitter<5>* fitter,
                         const reco::BeamSpot* bs,
                         const GlobalPoint* priorPos,
                         const GlobalError* priorError) const;
};

#endif