NjettinessAdder.h

Go to the documentation of this file.

#ifndef NjettinessAdder_h
#define NjettinessAdder_h

#include <memory>
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/stream/EDProducer.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "DataFormats/JetReco/interface/Jet.h"
#include "DataFormats/Common/interface/ValueMap.h"
#include "fastjet/contrib/Njettiness.hh"

class NjettinessAdder : public edm::stream::EDProducer<> {
public:
  enum MeasureDefinition_t {
    NormalizedMeasure = 0,      // (beta,R0)
    UnnormalizedMeasure,        // (beta)
    OriginalGeometricMeasure,   // (beta)
    NormalizedCutoffMeasure,    // (beta,R0,Rcutoff)
    UnnormalizedCutoffMeasure,  // (beta,Rcutoff)
    GeometricCutoffMeasure,     // (beta,Rcutoff)
    N_MEASURE_DEFINITIONS
  };
  enum AxesDefinition_t {
    KT_Axes = 0,
    CA_Axes,
    AntiKT_Axes,  // (axAxesR0)
    WTA_KT_Axes,
    WTA_CA_Axes,
    Manual_Axes,
    OnePass_KT_Axes,
    OnePass_CA_Axes,
    OnePass_AntiKT_Axes,  // (axAxesR0)
    OnePass_WTA_KT_Axes,
    OnePass_WTA_CA_Axes,
    OnePass_Manual_Axes,
    MultiPass_Axes,
    N_AXES_DEFINITIONS
  };

  explicit NjettinessAdder(const edm::ParameterSet& iConfig);

  ~NjettinessAdder() override {}

  void produce(edm::Event& iEvent, const edm::EventSetup& iSetup) override;
  float getTau(unsigned num, const edm::Ptr<reco::Jet>& object) const;

private:
  edm::InputTag src_;
  edm::EDGetTokenT<edm::View<reco::Jet>> src_token_;
  std::vector<unsigned> Njets_;

  // Measure definition :
  unsigned measureDefinition_;
  double beta_;
  double R0_;
  double Rcutoff_;

  // Axes definition :
  unsigned axesDefinition_;
  int nPass_;
  double akAxesR0_;

  edm::EDGetTokenT<edm::ValueMap<float>> input_weights_token_;
  edm::ValueMap<float> const* weightsHandle_;

  std::unique_ptr<fastjet::contrib::Njettiness> routine_;
};

#endif