HcalDetectorMap.cxx

#include "Hcal/HcalDetectorMap.h"
 
#include <sstream>
 
#include "Framework/ConditionsObjectProvider.h"
#include "Framework/EventHeader.h"
 
namespace hcal {
 
class HcalDetectorMapLoader : public framework::ConditionsObjectProvider {
 public:
  HcalDetectorMapLoader(const std::string& name, const std::string& tagname,
                        const framework::config::Parameters& parameters,
                        framework::Process& process)
      : ConditionsObjectProvider(HcalDetectorMap::CONDITIONS_OBJECT_NAME,
                                 tagname, parameters, process),
        the_map_{nullptr} {
    want_d2e_ = parameters.getParameter<bool>("want_d2e");
    connections_table_ =
        parameters.getParameter<std::string>("connections_table");
  }
 
  virtual std::pair<const framework::ConditionsObject*,
                    framework::ConditionsIOV>
  getCondition(const ldmx::EventHeader& context) {
    if (!the_map_) {
      the_map_ = new HcalDetectorMap(connections_table_, want_d2e_);
    }
 
    return std::make_pair(
        the_map_, framework::ConditionsIOV(context.getRun(), context.getRun(),
                                           true, true));
  }
 
  virtual void releaseConditionsObject(const framework::ConditionsObject* co) {}
 
 private:
  HcalDetectorMap* the_map_;
  std::string connections_table_;
  bool want_d2e_;
};
 
HcalDetectorMap::HcalDetectorMap(const std::string& connections_table,
                                 bool want_d2e)
    : framework::ConditionsObject(CONDITIONS_OBJECT_NAME),
      ldmx::ElectronicsMap<ldmx::HcalElectronicsID, ldmx::HcalDigiID>(
          want_d2e) {
  this->clear();
  conditions::StreamCSVLoader csv(connections_table);
  while (csv.nextRow()) {
    int polarfire = csv.getInteger("Polarfire");
    int hgcroc = csv.getInteger("HGCROC");
    int chan = csv.getInteger("Channel");
 
    int link = 2 * hgcroc;  // link count starts from 0
    if (chan >= 36) {
      link += 1;
      chan -= 36;
    }
 
    // one quadbar groups 4 strips and each quadbar is connected to 2 CMBs - one
    // in each end
    int end = csv.getInteger("CMB") % 2;
    if (end == 0)
      end = 1;  // negative end (?)
    else
      end = 0;  // positive end (?)
 
    // The Bar field in the map represents a numbering *within a quadbar* not
    // the strip number.
    //
    // This gives correct range for the strip numbers [0,7] and [0,11] but the
    // actual order of the quadbars has to be checked
 
    auto QuadBar{csv.getInteger("QuadBar")};
    auto BarNumber{csv.getInteger("Bar")};
    // Strip = (4 - Bar) + 4 * (QuadBar - 1)
    // QuadBar starts at 1 so we need to subtract to get first value
    // to be zero
    auto strip{(BarNumber - 1) +
               4 * (QuadBar - 1)};  // This is flipped from the original
    if (csv.getInteger("Plane") % 2 == 1)
      strip = (4 - BarNumber) +
              4 * (QuadBar - 1);  // Plane 2 is backward by mistake
    ldmx::HcalDigiID detid(0 /*section - only one section during test beam*/,
                           csv.getInteger("Plane") /*layer*/, strip /*strip*/,
                           end /*end*/);
    ldmx::HcalElectronicsID eleid(polarfire /*polarfire fpga*/, link /*elink*/,
                                  chan /*channel*/);
 
    if (this->exists(eleid)) {
      std::stringstream ss;
      ss << "Two different mappings for electronics channel " << eleid;
      EXCEPTION_RAISE("DuplicateMapping", ss.str());
    }
    this->addEntry(eleid, detid);
  }
}
 
}  // namespace hcal
DECLARE_CONDITIONS_PROVIDER_NS(hcal, HcalDetectorMapLoader);