QjetsAdder.h

Go to the documentation of this file.

#ifndef QjetsAdder_h
#define QjetsAdder_h

#include <memory>
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/EDProducer.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "DataFormats/PatCandidates/interface/Jet.h"
#include "RecoJets/JetAlgorithms/interface/QjetsPlugin.h"

class QjetsAdder : public edm::EDProducer {
public:
  explicit QjetsAdder(const edm::ParameterSet& iConfig)
      : src_(iConfig.getParameter<edm::InputTag>("src")),
        src_token_(consumes<edm::View<reco::Jet>>(src_)),
        qjetsAlgo_(iConfig.getParameter<double>("zcut"),
                   iConfig.getParameter<double>("dcutfctr"),
                   iConfig.getParameter<double>("expmin"),
                   iConfig.getParameter<double>("expmax"),
                   iConfig.getParameter<double>("rigidity")),
        ntrial_(iConfig.getParameter<int>("ntrial")),
        cutoff_(iConfig.getParameter<double>("cutoff")),
        jetRad_(iConfig.getParameter<double>("jetRad")),
        mJetAlgo_(iConfig.getParameter<std::string>("jetAlgo")),
        QjetsPreclustering_(iConfig.getParameter<int>("preclustering")) {
    produces<edm::ValueMap<float>>("QjetsVolatility");
  }

  ~QjetsAdder() override {}

  void produce(edm::Event& iEvent, const edm::EventSetup& iSetup) override;

private:
  edm::InputTag src_;
  edm::EDGetTokenT<edm::View<reco::Jet>> src_token_;
  QjetsPlugin qjetsAlgo_;
  int ntrial_;
  double cutoff_;
  double jetRad_;
  std::string mJetAlgo_;
  int QjetsPreclustering_;
  edm::Service<edm::RandomNumberGenerator> rng_;
};

#endif