HGCalLayerClusterHeterogeneousSoADumper.cc

Go to the documentation of this file.

#include "DataFormats/HGCalReco/interface/HGCalSoAClustersHostCollection.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/global/EDAnalyzer.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
#include "FWCore/Utilities/interface/EDGetToken.h"
#include "FWCore/Utilities/interface/Exception.h"
#include "FWCore/Utilities/interface/InputTag.h"

#include <fmt/format.h>

class HGCalLayerClusterHeterogeneousSoADumper : public edm::global::EDAnalyzer<> {
public:
  HGCalLayerClusterHeterogeneousSoADumper(edm::ParameterSet const& iConfig)
      : token_{consumes(iConfig.getParameter<edm::InputTag>("src"))} {}

  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
    edm::ParameterSetDescription desc;
    desc.add<edm::InputTag>("src", edm::InputTag("hltHgcalSoALayerClustersProducer"));
    descriptions.addWithDefaultLabel(desc);
  }

  void analyze(edm::StreamID iStream, edm::Event const& iEvent, edm::EventSetup const& iSetup) const override {
    auto const& data = iEvent.get(token_);

    auto const view = data.view();
    std::cout << fmt::format("hgcalSoALayerClustersProducer size = {}", view.metadata().size()) << std::endl;
    for (int i = 0; i < data->metadata().size(); ++i) {
      std::cout << fmt::format("CLUSTERS_SOA {}, energy = {:.{}f}, x = {:.{}f}, y = {:.{}f}, z= {:.{}f}",
                               i,
                               view.energy(i),
                               std::numeric_limits<float>::max_digits10,
                               view.x(i),
                               std::numeric_limits<float>::max_digits10,
                               view.y(i),
                               std::numeric_limits<float>::max_digits10,
                               view.z(i),
                               std::numeric_limits<float>::max_digits10)
                << std::endl;
    }
  }

private:
  edm::EDGetTokenT<HGCalSoAClustersHostCollection> const token_;
};

#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_FWK_MODULE(HGCalLayerClusterHeterogeneousSoADumper);