BPHPlusMinusCandidate.cc

Go to the documentation of this file.

/*
 *  See header file for a description of this class.
 *
 *  \author Paolo Ronchese INFN Padova
 *
 */

//-----------------------
// This Class' Header --
//-----------------------
#include "HeavyFlavorAnalysis/RecoDecay/interface/BPHPlusMinusCandidate.h"

//-------------------------------
// Collaborating Class Headers --
//-------------------------------
#include "HeavyFlavorAnalysis/RecoDecay/interface/BPHRecoBuilder.h"
#include "HeavyFlavorAnalysis/RecoDecay/interface/BPHRecoSelect.h"
#include "HeavyFlavorAnalysis/RecoDecay/interface/BPHMomentumSelect.h"
#include "DataFormats/PatCandidates/interface/CompositeCandidate.h"
#include "DataFormats/Math/interface/deltaPhi.h"
#include "CommonTools/CandUtils/interface/AddFourMomenta.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

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

using namespace std;

//-------------------
// Initializations --
//-------------------

//----------------
// Constructors --
//----------------
BPHPlusMinusCandidate::BPHPlusMinusCandidate(const edm::EventSetup* es)
    : BPHDecayVertex(es), BPHPlusMinusVertex(es), BPHRecoCandidate(es) {}

BPHPlusMinusCandidate::BPHPlusMinusCandidate(const edm::EventSetup* es, const BPHRecoBuilder::ComponentSet& compList)
    : BPHDecayMomentum(compList.daugMap, compList.compMap),
      BPHDecayVertex(this, es),
      BPHKinematicFit(this),
      BPHPlusMinusVertex(es),
      BPHRecoCandidate(es, compList) {}

//--------------
// Destructor --
//--------------
BPHPlusMinusCandidate::~BPHPlusMinusCandidate() {}

//--------------
// Operations --
//--------------
void BPHPlusMinusCandidate::add(const string& name, const reco::Candidate* daug, double mass, double sigma) {
  add(name, daug, "cfhpmig", mass, sigma);
  return;
}

void BPHPlusMinusCandidate::add(
    const string& name, const reco::Candidate* daug, const string& searchList, double mass, double sigma) {
  const vector<const reco::Candidate*>& dL = daughters();
  switch (dL.size()) {
    case 2:
      edm::LogPrint("TooManyParticles") << "BPHPlusMinusCandidate::add: "
                                        << "complete, add rejected";
      return;
    case 1:
      if ((daug->charge() * dL.front()->charge()) > 0) {
        edm::LogPrint("TooManyParticles") << "BPHPlusMinusCandidate::add: "
                                          << "already containing same sign particle, add rejected";
        return;
      }
      [[fallthrough]];
    case 0:
      addK(name, daug, searchList, mass, sigma);
  }
  return;
}

vector<BPHPlusMinusConstCandPtr> BPHPlusMinusCandidate::build(
    const BPHRecoBuilder& builder, const string& nPos, const string& nNeg, double mass, double msig) {
  vector<BPHPlusMinusConstCandPtr> cList;
  class ChargeSelect : public BPHRecoSelect {
  public:
    ChargeSelect(int c) : charge(c) {}
    ~ChargeSelect() override {}
    bool accept(const reco::Candidate& cand) const override { return ((charge * cand.charge()) > 0); }

  private:
    int charge;
  };
  ChargeSelect tkPos(+1);
  ChargeSelect tkNeg(-1);
  builder.filter(nPos, tkPos);
  builder.filter(nNeg, tkNeg);
  fill<BPHPlusMinusCandidate>(cList, builder, mass, msig);
  return cList;
}

const pat::CompositeCandidate& BPHPlusMinusCandidate::composite() const {
  static const pat::CompositeCandidate compCand;
  static const string msg = "BPHPlusMinusCandidate incomplete, no composite available";
  if (!chkSize(msg))
    return compCand;
  return BPHDecayMomentum::composite();
}

bool BPHPlusMinusCandidate::isCowboy() const {
  static const string msg = "BPHPlusMinusCandidate incomplete, no cowboy/sailor classification";
  return (chkSize(msg) && phiDiff());
}

bool BPHPlusMinusCandidate::isSailor() const {
  static const string msg = "BPHPlusMinusCandidate incomplete, no cowboy/sailor classification";
  return (chkSize(msg) && !phiDiff());
}

bool BPHPlusMinusCandidate::phiDiff() const {
  const vector<const reco::Candidate*>& dL = daughters();
  int idPos = (dL.front()->charge() > 0 ? 0 : 1);
  return reco::deltaPhi(dL[idPos]->phi(), dL[1 - idPos]->phi()) > 0;
}