hgcalUtils::DumpClusters Class Reference

#include <DumpClustersDetails.h>

Public Member Functions
	DumpClusters ()=default
template<typename T >
void	dumpInfos (const T &clusters, const std::string &moduleLabel, edm::RunNumber_t run, edm::LuminosityBlockNumber_t lumi, edm::EventNumber_t event, bool dumpCellsDetId=false) const

Detailed Description

Definition at line 22 of file DumpClustersDetails.h.

Constructor & Destructor Documentation

DumpClusters()

hgcalUtils::DumpClusters::DumpClusters ( )

default

Member Function Documentation

dumpInfos()

template<typename T >

void hgcalUtils::DumpClusters::dumpInfos ( const T &  clusters, const std::string &  moduleLabel, edm::RunNumber_t  run, edm::LuminosityBlockNumber_t  lumi, edm::EventNumber_t  event, bool  dumpCellsDetId = false  ) const

inline

Definition at line 27 of file DumpClustersDetails.h.

References DummyCfis::c, bsc_activity_cfg::clusters, submitPVResolutionJobs::count, corrVsCorr::filename, nano_mu_digi_cff::float, dqm-mbProfile::format, mps_fire::i, HerwigMaxPtPartonFilter_cfi::moduleLabel, das-up-to-nevents::outfile, writedatasetfile::run, and hgcalUtils::sortByDetId().

Referenced by HGCalLayerClustersFromSoAProducer::produce(), and HGCalLayerClusterProducer::produce().

                                                      {
      // Get the process ID
      pid_t pid = getpid();

      // Create the filename using the PID
      std::ostringstream filename;
      filename << "CLUSTERS_" << pid << "_" << moduleLabel << "_" << run << "_" << lumi << "_" << event << ".txt";
      // Open the file
      std::ofstream outfile(filename.str());
      int count = 0;
      for (auto const& i : clusters) {
        outfile << fmt::format(
            "Seed: {}, Idx: {}, energy: {:.{}f}, x: {:.{}f}, y: {:.{}f}, z: {:.{}f}, eta: {:.{}f}, phi: {:.{}f}",
            i.seed().rawId(),
            count++,
            (float)i.energy(),
            std::numeric_limits<float>::max_digits10,
            i.x(),
            std::numeric_limits<float>::max_digits10,
            i.y(),
            std::numeric_limits<float>::max_digits10,
            i.z(),
            std::numeric_limits<float>::max_digits10,
            i.eta(),
            std::numeric_limits<float>::max_digits10,
            i.phi(),
            std::numeric_limits<float>::max_digits10);
        if (dumpCellsDetId) {
          auto sorted = i.hitsAndFractions();  // copy...
          std::stable_sort(std::begin(sorted), std::end(sorted), sortByDetId);
          for (auto const& c : sorted) {
            outfile << fmt::format(" ({}, {:.{}f})", c.first, c.second, std::numeric_limits<float>::max_digits10);
          }  // loop on hits and fractions
        } else {
          outfile << fmt::format(" ({} cells)", i.hitsAndFractions().size());
        }
        outfile << std::endl;
      }  // loop on clusters
      outfile.close();
    }