BPHPlusMinusVertex.h

Go to the documentation of this file.

#ifndef HeavyFlavorAnalysis_RecoDecay_BPHPlusMinusVertex_h
#define HeavyFlavorAnalysis_RecoDecay_BPHPlusMinusVertex_h

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

//------------------------------------
// Collaborating Class Declarations --
//------------------------------------
#include "HeavyFlavorAnalysis/RecoDecay/interface/BPHPlusMinusCandidatePtr.h"
#include "TrackingTools/PatternTools/interface/ClosestApproachInRPhi.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

//---------------
// C++ Headers --
//---------------
#include <iostream>

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

class BPHPlusMinusVertex : public virtual BPHDecayVertex {
public:
  ~BPHPlusMinusVertex() override;

  virtual const ClosestApproachInRPhi& cAppInRPhi() const;

protected:
  BPHPlusMinusVertex(const edm::EventSetup* es);

  // utility functions to check/enforce the number of decay particles
  // at 2
  template <class T>
  static bool chkName(const T& cont, const std::string& name, const std::string& msg);
  template <class T>
  static bool chkSize(const T& cont, const std::string& msg);
  bool chkSize(const std::string& msg) const;

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

private:
  // reconstruction results cache
  mutable bool oldA;
  mutable ClosestApproachInRPhi* inRPhi;

  // compute closest approach distance and cache it
  virtual void computeApp() const;
};

template <class T>
bool BPHPlusMinusVertex::chkName(const T& cont, const std::string& name, const std::string& msg) {
  if (cont.find(name) != cont.end())
    return true;
  edm::LogPrint("ParticleNotFound") << msg << ", " << name << " not found";
  return false;
}

template <class T>
bool BPHPlusMinusVertex::chkSize(const T& cont, const std::string& msg) {
  int n = cont.size();
  if (n == 2)
    return true;
  edm::LogPrint("WrongDataSize") << msg << ", size = " << n;
  return false;
}

#endif