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