HGCalValidator.h

Go to the documentation of this file.

#ifndef HGCalValidator_h
#define HGCalValidator_h

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/EDGetToken.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/ParameterSet/interface/FileInPath.h"

#include "RecoLocalCalo/HGCalRecAlgos/interface/RecHitTools.h"
#include "DataFormats/ParticleFlowReco/interface/PFCluster.h"
#include "DataFormats/ParticleFlowCandidate/interface/PFCandidate.h"
#include "SimDataFormats/CaloAnalysis/interface/CaloParticle.h"
#include "SimDataFormats/CaloAnalysis/interface/SimCluster.h"
#include "SimDataFormats/Vertex/interface/SimVertex.h"

#include "DQMServices/Core/interface/DQMGlobalEDAnalyzer.h"

#include "Validation/HGCalValidation/interface/HGVHistoProducerAlgo.h"
#include "Validation/HGCalValidation/interface/CaloParticleSelector.h"
#include "RecoLocalCalo/HGCalRecProducers/interface/HGCalClusteringAlgoBase.h"

class PileupSummaryInfo;

struct HGCalValidatorHistograms {
  HGVHistoProducerAlgoHistograms histoProducerAlgo;
  std::vector<ConcurrentMonitorElement> h_layerclusters_coll;
};

class HGCalValidator : public DQMGlobalEDAnalyzer<HGCalValidatorHistograms> {
 public:
  using Histograms = HGCalValidatorHistograms;

  HGCalValidator(const edm::ParameterSet& pset);

  ~HGCalValidator() override;


  void dqmAnalyze(const edm::Event&, const edm::EventSetup&, const Histograms& ) const override;
  void bookHistograms(DQMStore::ConcurrentBooker&, edm::Run const&, edm::EventSetup const&, Histograms&) const override;

  void cpParametersAndSelection(const Histograms& histograms, std::vector<CaloParticle> const & cPeff, std::vector<SimVertex> const & simVertices, std::vector<size_t>& selected_cPeff) const;


 protected:

  std::vector<edm::InputTag> label;
  const bool SaveGeneralInfo_;
  const bool doCaloParticlePlots_;
  const bool dolayerclustersPlots_;
  const edm::FileInPath cummatbudinxo_;

  std::vector<edm::EDGetTokenT<reco::CaloClusterCollection> > labelToken;
  edm::EDGetTokenT<std::vector<CaloParticle> > label_cp_effic;
  edm::EDGetTokenT<std::vector<CaloParticle> > label_cp_fake;
  edm::EDGetTokenT<std::vector<SimVertex> > simVertices_;
  edm::EDGetTokenT<HGCRecHitCollection> recHitsEE_;
  edm::EDGetTokenT<HGCRecHitCollection> recHitsFH_;
  edm::EDGetTokenT<HGCRecHitCollection> recHitsBH_;
  edm::EDGetTokenT<Density> density_;
  std::unique_ptr<HGVHistoProducerAlgo> histoProducerAlgo_;



 private:

  void fillHitMap(std::map<DetId, const HGCRecHit * > &,
      const HGCRecHitCollection &,
      const HGCRecHitCollection &,
      const HGCRecHitCollection &) const;
  CaloParticleSelector cpSelector;
  std::shared_ptr<hgcal::RecHitTools> tools_;
  std::map<double, double> cummatbudg;
  std::vector<int> particles_to_monitor_;
  unsigned totallayers_to_monitor_;
  std::vector<int> thicknesses_to_monitor_;
  std::string dirName_;
};


#endif