MTDTopologyEP.cc

Go to the documentation of this file.

//#define EDM_ML_DEBUG

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
#include "FWCore/Framework/interface/ModuleFactory.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/ESProducer.h"
#include "Geometry/MTDNumberingBuilder/interface/MTDTopology.h"
#include "Geometry/MTDCommonData/interface/MTDTopologyMode.h"
#include "Geometry/Records/interface/MTDTopologyRcd.h"
#include "CondFormats/GeometryObjects/interface/PMTDParameters.h"
#include "Geometry/Records/interface/PMTDParametersRcd.h"

#include <memory>

class MTDTopologyEP : public edm::ESProducer {
public:
  MTDTopologyEP(const edm::ParameterSet&);

  using ReturnType = std::unique_ptr<MTDTopology>;

  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);

  ReturnType produce(const MTDTopologyRcd&);

private:
  void fillParameters(const PMTDParameters&, int& mtdTopologyMode, MTDTopology::ETLValues&);

  const edm::ESGetToken<PMTDParameters, PMTDParametersRcd> token_;
};

MTDTopologyEP::MTDTopologyEP(const edm::ParameterSet& conf)
    : token_{setWhatProduced(this).consumesFrom<PMTDParameters, PMTDParametersRcd>(edm::ESInputTag())} {}

void MTDTopologyEP::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
  edm::ParameterSetDescription ttc;
  descriptions.add("mtdTopology", ttc);
}

MTDTopologyEP::ReturnType MTDTopologyEP::produce(const MTDTopologyRcd& iRecord) {
  int mtdTopologyMode;
  MTDTopology::ETLValues etlVals;

  fillParameters(iRecord.get(token_), mtdTopologyMode, etlVals);

  return std::make_unique<MTDTopology>(mtdTopologyMode, etlVals);
}

void MTDTopologyEP::fillParameters(const PMTDParameters& ptp, int& mtdTopologyMode, MTDTopology::ETLValues& etlVals) {
  mtdTopologyMode = ptp.topologyMode_;

  // for legacy geometry scenarios no topology informastion is stored, only for newer ETL 2-discs layout

  if (mtdTopologyMode <= static_cast<int>(MTDTopologyMode::Mode::barphiflat)) {
    return;
  }

  // Check on the internal consistency of thr ETL layout information provided by parameters

  for (size_t it = 3; it <= 9; it++) {
    if (ptp.vitems_[it].vpars_.size() != ptp.vitems_[2].vpars_.size()) {
      throw cms::Exception("MTDTopologyEP") << "Inconsistent size of ETL structure arrays";
    }
  }

  MTDTopology::ETLfaceLayout tmpFace;

  // Front Face (0), starting with type Right (2)

  tmpFace.idDiscSide_ = 0;  // ETL front side
  tmpFace.idDetType1_ = 2;  // ETL module type right

  tmpFace.start_copy_[0] = ptp.vitems_[3].vpars_;  // start_copy_FR
  tmpFace.start_copy_[1] = ptp.vitems_[2].vpars_;  // start_copy_FL
  tmpFace.offset_[0] = ptp.vitems_[7].vpars_;      // offset_FR
  tmpFace.offset_[1] = ptp.vitems_[6].vpars_;      // offset_FL

  etlVals.emplace_back(tmpFace);

  // Back Face (1), starting with type Left (1)

  tmpFace.idDiscSide_ = 1;  // ETL back side
  tmpFace.idDetType1_ = 1;  // ETL module type left

  tmpFace.start_copy_[0] = ptp.vitems_[4].vpars_;  // start_copy_BL
  tmpFace.start_copy_[1] = ptp.vitems_[5].vpars_;  // start_copy_BR
  tmpFace.offset_[0] = ptp.vitems_[8].vpars_;      // offset_BL
  tmpFace.offset_[1] = ptp.vitems_[9].vpars_;      // offset_BR

  etlVals.emplace_back(tmpFace);

#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("MTDTopologyEP") << " Topology mode = " << mtdTopologyMode << "\n";
  auto print_array = [&](std::vector<int> vector) {
    std::stringstream ss;
    for (auto const& elem : vector) {
      ss << " " << elem;
    }
    ss << "\n";
    return ss.str();
  };

  for (auto const& ilay : etlVals) {
    edm::LogVerbatim("MTDTopologyEP") << " disc face = " << ilay.idDiscSide_ << " start det type = " << ilay.idDetType1_
                                      << "\n start_copy[0]= " << print_array(ilay.start_copy_[0])
                                      << "\n start_copy[1]= " << print_array(ilay.start_copy_[1])
                                      << "\n offset[0]= " << print_array(ilay.offset_[0])
                                      << "\n offset[1]= " << print_array(ilay.offset_[1]);
  }

#endif
}

DEFINE_FWK_EVENTSETUP_MODULE(MTDTopologyEP);