PMTChannelMapDumper.cxx File Reference

Utility dumping the content of PMT channel mapping on screen. More...

#include "icaruscode/Decode/ChannelMapping/ICARUSChannelMapProvider.h"
#include "larcorealg/CoreUtils/enumerate.h"
#include "larcorealg/TestUtils/unit_test_base.h"
#include "messagefacility/MessageLogger/MessageLogger.h"
#include <iomanip>
#include <iostream>
#include <algorithm>
#include <numeric>
#include <array>

Go to the source code of this file.

Classes
struct	SortByElement< KeyNo >

Functions
int	main (int argc, char **argv)

Detailed Description

Utility dumping the content of PMT channel mapping on screen.

Author

Gianluca Petrillo (petri.nosp@m.llo@.nosp@m.slac..nosp@m.stan.nosp@m.ford..nosp@m.edu)

This utility can be run with any configuration file including a configuration for IICARUSChannel service.

It is using art facilities for tool loading, but it does not run in art environment. So it may break without warning and without solution.

Definition in file PMTChannelMapDumper.cxx.

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 51 of file PMTChannelMapDumper.cxx.

                                {
  
  using Environment
    = testing::TesterEnvironment<testing::BasicEnvironmentConfiguration>;
  
  testing::BasicEnvironmentConfiguration config("PMTchannelMappingDumper");

  //
  // parameter parsing
  //
  int iParam = 0;

  // first argument: configuration file (mandatory)
  if (++iParam < argc)
    config.SetConfigurationPath(argv[iParam]);
  else {
    std::cerr << "FHiCL configuration file path required as first argument!"
      << std::endl;
    return 1;
  }

  Environment const Env { config };
  
  icarusDB::ICARUSChannelMapProvider const channelMapping
    { Env.ServiceParameters("IICARUSChannelMap") };
  
  // hard-coded list of fragment ID; don't like it?
  // ask for an extension of the channel mapping features.
  std::array<unsigned int, 24U> FragmentIDs;
  std::iota(FragmentIDs.begin(), FragmentIDs.end(), 0x2000);
  
  mf::LogVerbatim("PMTchannelMappingDumper") << "Fragment IDs:";
  for (auto const [ iFragment, fragmentID]: util::enumerate(FragmentIDs)) {
    unsigned int const effFragmentID = fragmentID & 0xFF;
    
    if (!channelMapping.hasPMTDigitizerID(effFragmentID)) {
      mf::LogVerbatim("PMTchannelMappingDumper")
        << "[" << iFragment << "] " << std::hex << fragmentID << std::dec
        << " not found in the database (as "
        << std::hex << effFragmentID << std::dec << ")";
      continue;
    }
    
    icarusDB::DigitizerChannelChannelIDPairVec digitizerChannels
      = channelMapping.getChannelIDPairVec(effFragmentID);
    
    std::sort
      (digitizerChannels.begin(), digitizerChannels.end(), SortByElement<1U>{});
    
    
    mf::LogVerbatim log("PMTchannelMappingDumper");
    log
      << "[" << iFragment << "] " << std::hex << fragmentID << std::dec
      << " includes " << digitizerChannels.size()
      << " LArSoft channels between " << digitizerChannels.front().second
      << " and " << digitizerChannels.back().second
      << " [board channel index in brackets]:";
    constexpr unsigned int Cols = 8U;
    unsigned int n = 0;
    for(auto const [ digitizerChannel, channelID ]: digitizerChannels) {
      if (n-- == 0) { log << "\n     "; n = Cols - 1U; }
      log << " "  << std::setw(3) << channelID
        << " [" << std::setw(3) << digitizerChannel << "]";
    } // for channel
    
  } // for fragment
  
  return 0;
} // main()