DumpSimPhotonsLite_module.cc

Go to the documentation of this file.

// LArSoft libraries
#include "lardataobj/Simulation/SimPhotons.h"

// framework libraries
#include "art/Framework/Core/EDAnalyzer.h"
#include "art/Framework/Core/ModuleMacros.h"
#include "art/Framework/Principal/Event.h"
#include "art/Framework/Principal/Handle.h"
#include "canvas/Utilities/InputTag.h"
#include "fhiclcpp/types/Atom.h"
#include "messagefacility/MessageLogger/MessageLogger.h"

// C/C++ standard libraries
#include <iomanip> // std::setw()
#include <numeric> // std::accumulate()
#include <string>
#include <utility> // std::forward()

namespace sim {
  class DumpSimPhotonsLite;
} // namespace sim

namespace {
  using namespace fhicl;

  struct Config {
    using Name = fhicl::Name;
    using Comment = fhicl::Comment;

    fhicl::Atom<art::InputTag> InputPhotons{
      Name("InputPhotons"),
      Comment("data product with the SimPhotonsLite to be dumped")};

    fhicl::Atom<std::string> OutputCategory{
      Name("OutputCategory"),
      Comment("name of the output stream (managed by the message facility)"),
      "DumpSimPhotonsLite" /* default value */
    };

  }; // struct Config

} // local namespace

class sim::DumpSimPhotonsLite : public art::EDAnalyzer {
public:
  // type to enable module parameters description by art
  using Parameters = art::EDAnalyzer::Table<Config>;

  explicit DumpSimPhotonsLite(Parameters const& config);

  // Plugins should not be copied or assigned.
  DumpSimPhotonsLite(DumpSimPhotonsLite const&) = delete;
  DumpSimPhotonsLite(DumpSimPhotonsLite&&) = delete;
  DumpSimPhotonsLite& operator=(DumpSimPhotonsLite const&) = delete;
  DumpSimPhotonsLite& operator=(DumpSimPhotonsLite&&) = delete;

  // Operates on the event
  void analyze(art::Event const& event) override;

  template <typename Stream>
  void DumpPhoton(Stream&& out,
                  sim::SimPhotonsLite const& photons,
                  std::string indent,
                  std::string firstIndent) const;

  template <typename Stream>
  void DumpPhoton(Stream&& out, sim::SimPhotonsLite const& photons, std::string indent = "") const
  {
    DumpPhoton(std::forward<Stream>(out), photons, indent, indent);
  }

private:
  art::InputTag fInputPhotons; 
  std::string fOutputCategory; 

}; // class sim::DumpSimPhotonsLite

//------------------------------------------------------------------------------
//---  module implementation
//---
//------------------------------------------------------------------------------
sim::DumpSimPhotonsLite::DumpSimPhotonsLite(Parameters const& config)
  : EDAnalyzer(config)
  , fInputPhotons(config().InputPhotons())
  , fOutputCategory(config().OutputCategory())
{}

//------------------------------------------------------------------------------
template <typename Stream>
void sim::DumpSimPhotonsLite::DumpPhoton(Stream&& out,
                                         sim::SimPhotonsLite const& photons,
                                         std::string indent,
                                         std::string firstIndent) const
{

  unsigned int const nPhotons =
    std::accumulate(photons.DetectedPhotons.begin(),
                    photons.DetectedPhotons.end(),
                    0U,
                    [](auto sum, auto const& entry) { return sum + entry.second; });

  out << firstIndent << "channel=" << photons.OpChannel << " has ";
  if (nPhotons) {
    out << nPhotons << " photons (format: [tick] photons):";
    constexpr unsigned int PageSize = 5;
    unsigned int pager = 0;
    for (auto const& pair : photons.DetectedPhotons) {
      if (pager-- == 0) {
        pager = PageSize - 1;
        out << "\n" << indent << " ";
      }
      out << " [" << pair.first << "] " << std::setw(6) << pair.second;
    } // for
  }
  else {
    out << "no photons";
  }

} // sim::DumpSimPhotonsLite::DumpPhoton()

//------------------------------------------------------------------------------
void sim::DumpSimPhotonsLite::analyze(art::Event const& event)
{

  // get the particles from the event
  auto const& Photons = *(event.getValidHandle<std::vector<sim::SimPhotonsLite>>(fInputPhotons));

  mf::LogVerbatim(fOutputCategory)
    << "Event " << event.id() << " : data product '" << fInputPhotons.encode() << "' contains "
    << Photons.size() << " SimPhotonsLite";

  unsigned int iChannel = 0;
  for (sim::SimPhotonsLite const& photons : Photons) {

    mf::LogVerbatim log(fOutputCategory);
    // a bit of a header
    log << "[#" << (iChannel++) << "] ";
    DumpPhoton(log, photons, "  ");

  }                                         // for
  mf::LogVerbatim(fOutputCategory) << "\n"; // just an empty line

} // sim::DumpSimPhotonsLite::analyze()

//------------------------------------------------------------------------------
DEFINE_ART_MODULE(sim::DumpSimPhotonsLite)

//------------------------------------------------------------------------------