BPHKinFitChi2Select.h

Go to the documentation of this file.

#ifndef HeavyFlavorAnalysis_SpecificDecay_BPHKinFitChi2Select_h
#define HeavyFlavorAnalysis_SpecificDecay_BPHKinFitChi2Select_h

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

//------------------------------------
// Collaborating Class Declarations --
//------------------------------------
#include "HeavyFlavorAnalysis/RecoDecay/interface/BPHKinematicFit.h"
#include "TMath.h"

//---------------
// C++ Headers --
//---------------

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

class BPHKinFitChi2Select : public BPHFitSelect {
public:
  BPHKinFitChi2Select(double prob) : probMin(prob) {}

  // deleted copy constructor and assignment operator
  BPHKinFitChi2Select(const BPHKinFitChi2Select& x) = delete;
  BPHKinFitChi2Select& operator=(const BPHKinFitChi2Select& x) = delete;

  ~BPHKinFitChi2Select() override {}

  bool accept(const BPHKinematicFit& cand) const override {
    if (probMin < 0.0)
      return true;
    const RefCountedKinematicVertex tdv = cand.topDecayVertex();
    if (tdv.get() == nullptr)
      return false;
    if (!tdv->vertexIsValid())
      return false;
    reco::Vertex v(*tdv);
    if (v.isFake())
      return false;
    if (!v.isValid())
      return false;
    return (TMath::Prob(v.chi2(), lround(v.ndof())) >= probMin);
  }

  void setProbMin(double p) {
    probMin = p;
    return;
  }

  double getProbMin() const { return probMin; }

private:
  double probMin;
};

#endif