RecoTauPiZeroTrivialPlugin.cc

Go to the documentation of this file.

/*
 * RecoTauPiZeroTrivialPlugin
 *
 * Author: Evan K. Friis (UC Davis)
 *
 * Given an input PFJet, produces collection of trivial 'un-merged' PiZero
 * RecoTauPiZeros.  Each PiZero is composed of only one photon from
 * the jet.
 *
 *
 */

#include "RecoTauTag/RecoTau/interface/RecoTauPiZeroPlugins.h"
#include "DataFormats/Candidate/interface/CandidateFwd.h"
#include "DataFormats/Candidate/interface/Candidate.h"
#include "DataFormats/JetReco/interface/PFJet.h"
#include "DataFormats/TauReco/interface/RecoTauPiZero.h"

#include "RecoTauTag/RecoTau/interface/RecoTauCommonUtilities.h"
#include "RecoTauTag/RecoTau/interface/RecoTauQualityCuts.h"

namespace reco::tau {

  class RecoTauPiZeroTrivialPlugin : public RecoTauPiZeroBuilderPlugin {
  public:
    explicit RecoTauPiZeroTrivialPlugin(const edm::ParameterSet& pset, edm::ConsumesCollector&& iC);
    ~RecoTauPiZeroTrivialPlugin() override {}
    return_type operator()(const reco::Jet& jet) const override;

  private:
    RecoTauQualityCuts qcuts_;
  };

  RecoTauPiZeroTrivialPlugin::RecoTauPiZeroTrivialPlugin(const edm::ParameterSet& pset, edm::ConsumesCollector&& iC)
      : RecoTauPiZeroBuilderPlugin(pset, std::move(iC)),
        qcuts_(pset.getParameterSet("qualityCuts").getParameterSet("signalQualityCuts")) {}

  RecoTauPiZeroBuilderPlugin::return_type RecoTauPiZeroTrivialPlugin::operator()(const reco::Jet& jet) const {
    std::vector<CandidatePtr> pfGammaCands = qcuts_.filterCandRefs(pfGammas(jet));
    PiZeroVector output;
    output.reserve(pfGammaCands.size());

    for (auto const& gamma : pfGammaCands) {
      std::unique_ptr<RecoTauPiZero> piZero(
          new RecoTauPiZero(0, (*gamma).p4(), (*gamma).vertex(), 22, 1000, true, RecoTauPiZero::kTrivial));
      piZero->addDaughter(gamma);
      output.push_back(std::move(piZero));
    }
    return output.release();
  }

}  // namespace reco::tau

#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_EDM_PLUGIN(RecoTauPiZeroBuilderPluginFactory,
                  reco::tau::RecoTauPiZeroTrivialPlugin,
                  "RecoTauPiZeroTrivialPlugin");