RealisticSimClusterMapper.h

Go to the documentation of this file.

#ifndef __RecoParticleFlow_PFClusterProducer_RealisticSimClusterMapper_H__
#define __RecoParticleFlow_PFClusterProducer_RealisticSimClusterMapper_H__
// Author: Felice Pantaleo
// Date:   30/06/2017
// Email: felice@cern.ch
#include "RecoParticleFlow/PFClusterProducer/interface/InitialClusteringStepBase.h"
#include "DataFormats/ParticleFlowReco/interface/PFRecHitFraction.h"
#include "RecoLocalCalo/HGCalRecAlgos/interface/RecHitTools.h"

#include "SimDataFormats/CaloAnalysis/interface/SimClusterFwd.h"

class RealisticSimClusterMapper : public InitialClusteringStepBase {
public:
  RealisticSimClusterMapper(const edm::ParameterSet& conf, edm::ConsumesCollector& sumes)
      : InitialClusteringStepBase(conf, sumes),
        invisibleFraction_(conf.getParameter<double>("invisibleFraction")),
        exclusiveFraction_(conf.getParameter<double>("exclusiveFraction")),
        maxDistanceFilter_(conf.getParameter<bool>("maxDistanceFilter")),
        maxDistance_(conf.getParameter<double>("maxDistance")),
        maxDforTimingSquared_(conf.getParameter<double>("maxDforTimingSquared")),
        timeOffset_(conf.getParameter<double>("timeOffset")),
        minNHitsforTiming_(conf.getParameter<unsigned int>("minNHitsforTiming")),
        useMCFractionsForExclEnergy_(conf.getParameter<bool>("useMCFractionsForExclEnergy")),
        calibMinEta_(conf.getParameter<double>("calibMinEta")),
        calibMaxEta_(conf.getParameter<double>("calibMaxEta")) {
    simClusterToken_ = sumes.consumes<SimClusterCollection>(conf.getParameter<edm::InputTag>("simClusterSrc"));
    hadronCalib_ = conf.getParameter<std::vector<double> >("hadronCalib");
    egammaCalib_ = conf.getParameter<std::vector<double> >("egammaCalib");
  }

  ~RealisticSimClusterMapper() override {}
  RealisticSimClusterMapper(const RealisticSimClusterMapper&) = delete;
  RealisticSimClusterMapper& operator=(const RealisticSimClusterMapper&) = delete;

  void updateEvent(const edm::Event&) final;
  void update(const edm::EventSetup&) final;

  void buildClusters(const edm::Handle<reco::PFRecHitCollection>&,
                     const std::vector<bool>&,
                     const std::vector<bool>&,
                     reco::PFClusterCollection&) override;

private:
  hgcal::RecHitTools rhtools_;
  const float invisibleFraction_ = 0.3f;
  const float exclusiveFraction_ = 0.7f;
  const bool maxDistanceFilter_ = false;
  const float maxDistance_ = 10.f;
  const float maxDforTimingSquared_ = 4.0f;
  const float timeOffset_;
  const unsigned int minNHitsforTiming_ = 3;
  const bool useMCFractionsForExclEnergy_ = false;
  const float calibMinEta_ = 1.4;
  const float calibMaxEta_ = 3.0;
  std::vector<double> hadronCalib_;
  std::vector<double> egammaCalib_;

  edm::EDGetTokenT<SimClusterCollection> simClusterToken_;
  edm::Handle<SimClusterCollection> simClusterH_;
};

DEFINE_EDM_PLUGIN(InitialClusteringStepFactory, RealisticSimClusterMapper, "RealisticSimClusterMapper");

#endif