TPCChannelMapSP.cpp

Go to the documentation of this file.

// Class:       TPCChannelMapSP
// Module type: standalone algorithm
// File:        TPCChannelMapSP.cxx
// Author:      Tom Junk, Jan 2025
//
// Implementation of hardware-offline channel mapping reading from a file.
 
#include "TPCChannelMapSP.h"
 
#include <fstream>
#include <iostream>
#include <sstream>
#include <stdexcept>
 
dune::TPCChannelMapSP::TPCChannelMapSP() {
  fSubstituteCrate = 1;
  fNChans = 0;
}
 
void dune::TPCChannelMapSP::ReadMapFromFile(std::string& fullname)
{
  fNChans = 0;
  
  std::ifstream inFile(fullname, std::ios::in);
  std::string line;
 
  while (std::getline(inFile, line)) {
    std::stringstream linestream(line);
 
    TPCChanInfo_t chanInfo;
    linestream >>
      chanInfo.offlchan >>
      chanInfo.detid >>
      chanInfo.detelement >>
      chanInfo.crate >>
      chanInfo.slot >>
      chanInfo.stream >>
      chanInfo.streamchan >>
      chanInfo.plane >>
      chanInfo.chan_in_plane >>
      chanInfo.femb >>
      chanInfo.asic >> 
      chanInfo.asicchan;
 
    chanInfo.valid = true;
    ++fNChans;
 
    // fill maps.
 
    check_offline_channel(chanInfo.offlchan);
 
    auto hashval = make_hash(chanInfo.detid,chanInfo.crate,chanInfo.slot,chanInfo.stream,chanInfo.streamchan);
    if ( DetToChanInfo.find(hashval) != DetToChanInfo.end() ){
      std::cout << "TPCChannelMapSP: duplicate electronics channel found.  detid, crate, slot, stream, streamchan: " <<
    chanInfo.detid << " " <<
    chanInfo.crate << " " <<
    chanInfo.slot << " " <<
    chanInfo.stream << " " <<
    chanInfo.streamchan << std::endl;
      throw std::range_error("Duplicate Electronics ID");
    }
    DetToChanInfo[hashval] = chanInfo;
 
    if ( OfflToChanInfo.find(chanInfo.offlchan) != OfflToChanInfo.end() ) {
      std::cout << "TPCChannelMapSP: duplicate offline channel found: " <<
    chanInfo.offlchan << std::endl;
      throw std::range_error("Duplicate Offline TPC Channel ID");      
    }
    OfflToChanInfo[chanInfo.offlchan] = chanInfo;
  }
  inFile.close();
}
 
dune::TPCChannelMapSP::TPCChanInfo_t dune::TPCChannelMapSP::GetChanInfoFromElectronicsIDs(
                                                                                          unsigned int detid,
                                                                                          unsigned int crate,
                                                                                          unsigned int slot,
                                                                                          unsigned int stream,
                                                                                          unsigned int streamchan) const
{
  size_t h = make_hash(detid,crate,slot,stream,streamchan);
  auto ret = DetToChanInfo.find(h); 
  if (ret == DetToChanInfo.end())
    {
      h = make_hash(detid,fSubstituteCrate,slot,stream,streamchan);
      ret = DetToChanInfo.find(h);
      if (ret == DetToChanInfo.end())
    {
          TPCChanInfo_t badInfo = {};
          badInfo.valid = false;
          return badInfo;
    }
    }
  return ret->second;
}
 
dune::TPCChannelMapSP::TPCChanInfo_t dune::TPCChannelMapSP::GetChanInfoFromOfflChan(
                                                                                    unsigned int offlineChannel) const
{
  auto ci = OfflToChanInfo.find(offlineChannel);
  if (ci == OfflToChanInfo.end()) {
    TPCChanInfo_t badInfo = {};
    badInfo.valid = false;
    return badInfo;
  }
  return ci->second;
}
 
size_t dune::TPCChannelMapSP::make_hash( unsigned int detid,   //6 bits
                  unsigned int crate,   //10 bits
                  unsigned int slot,    //4 bits
                  unsigned int stream,  //8 bits
                  unsigned int streamch //12 bits
                     ) const
{
  size_t hashval = (((size_t)detid & 0x3f)<<34);
  hashval ^=((crate & 0x3ff)<<24);
  hashval ^=((slot & 0xf)<<20);
  hashval ^=((stream & 0xff)<<12);
  hashval ^=(streamch & 0xfff);
  return hashval;
}