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/ParameterSet.h"
#include "fhiclcpp/types/Atom.h"
#include "messagefacility/MessageLogger/MessageLogger.h"

// C/C++ standard libraries
#include <string>
#include <vector>
#include <iomanip> // std::setw()
#include <algorithm> // std::accumulate()
#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)

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