SiStripRegionConnectivity.cc

Go to the documentation of this file.

#include "FWCore/Framework/interface/ESProducer.h"
#include "FWCore/Framework/interface/ModuleFactory.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
#include "Geometry/TrackerGeometryBuilder/interface/StripGeomDetUnit.h"
#include "Geometry/CommonDetUnit/interface/GeomDet.h"
#include "CalibTracker/Records/interface/SiStripDetCablingRcd.h"
#include "CalibTracker/Records/interface/SiStripRegionCablingRcd.h"
#include "CalibFormats/SiStripObjects/interface/SiStripDetCabling.h"
#include "CalibFormats/SiStripObjects/interface/SiStripRegionCabling.h"
#include "DataFormats/SiStripCommon/interface/SiStripConstants.h"
#include "Geometry/Records/interface/TrackerTopologyRcd.h"
 
using namespace sistrip;
 
class SiStripRegionConnectivity : public edm::ESProducer {
public:
  SiStripRegionConnectivity(const edm::ParameterSet&);
  ~SiStripRegionConnectivity() override;
 
  std::unique_ptr<SiStripRegionCabling> produceRegionCabling(const SiStripRegionCablingRcd&);
 
private:
  edm::ESGetToken<SiStripDetCabling, SiStripDetCablingRcd> detcablingToken_;
  edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> tkgeomToken_;
  edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> tTopoToken_;
 
  uint32_t etadivisions_;
  uint32_t phidivisions_;
 
  double etamax_;
};
 
SiStripRegionConnectivity::SiStripRegionConnectivity(const edm::ParameterSet& pset)
    : etadivisions_(pset.getUntrackedParameter<unsigned int>("EtaDivisions", 10)),
      phidivisions_(pset.getUntrackedParameter<unsigned int>("PhiDivisions", 10)),
      etamax_(pset.getUntrackedParameter<double>("EtaMax", 2.4))
 
{
  auto cc = setWhatProduced(this, &SiStripRegionConnectivity::produceRegionCabling);
  detcablingToken_ = cc.consumes();
  tkgeomToken_ = cc.consumes();
  tTopoToken_ = cc.consumes();
}
 
SiStripRegionConnectivity::~SiStripRegionConnectivity() {}
 
std::unique_ptr<SiStripRegionCabling> SiStripRegionConnectivity::produceRegionCabling(
    const SiStripRegionCablingRcd& iRecord) {
  const auto& detcabling = iRecord.get(detcablingToken_);
  const auto& tkgeom = iRecord.get(tkgeomToken_);
  const auto& tTopo = iRecord.get(tTopoToken_);
 
  //here build an object of type SiStripRegionCabling using the information from class SiStripDetCabling **PLUS** the geometry.
 
  //Construct region cabling object
  auto RegionConnections = std::make_unique<SiStripRegionCabling>(etadivisions_, phidivisions_, etamax_);
 
  //Construct region cabling map
  SiStripRegionCabling::Cabling regioncabling(
      etadivisions_ * phidivisions_,
      SiStripRegionCabling::RegionCabling(
          SiStripRegionCabling::ALLSUBDETS,
          SiStripRegionCabling::WedgeCabling(SiStripRegionCabling::ALLLAYERS, SiStripRegionCabling::ElementCabling())));
 
  //Loop det cabling
  for (const auto& idet : detcabling.getDetCabling()) {
    if (!idet.first || (idet.first == sistrip::invalid32_))
      continue;
 
    // Check if geom det unit exists
    auto geom_det = tkgeom.idToDetUnit(DetId(idet.first));
    auto strip_det = dynamic_cast<StripGeomDetUnit const*>(geom_det);
    if (!strip_det) {
      continue;
    }
 
    //Calculate region from geometry
    double eta = tkgeom.idToDet(DetId(idet.first))->position().eta();
    double phi = tkgeom.idToDet(DetId(idet.first))->position().phi().value();
    uint32_t reg = RegionConnections->region(SiStripRegionCabling::Position(eta, phi));
 
    //Find subdet from det-id
    uint32_t subdet = static_cast<uint32_t>(SiStripRegionCabling::subdetFromDetId(idet.first));
 
    //Find layer from det-id
    uint32_t layer = tTopo.layer(idet.first);
 
    //@@ BELOW IS TEMP FIX TO HANDLE BUG IN DET CABLING
    const std::vector<const FedChannelConnection*>& conns = idet.second;
 
    //Update region cabling map
    regioncabling[reg][subdet][layer].push_back(SiStripRegionCabling::Element());
    auto& elem = regioncabling[reg][subdet][layer].back();
    elem.first = idet.first;
    elem.second.resize(conns.size());
    for (const auto& iconn : conns) {
      if ((iconn != nullptr) && (iconn->apvPairNumber() < conns.size())) {
        elem.second[iconn->apvPairNumber()] = *iconn;
      }
    }
  }
 
  //Add map to region cabling object
  RegionConnections->setRegionCabling(regioncabling);
 
  return RegionConnections;
}
 
DEFINE_FWK_EVENTSETUP_MODULE(SiStripRegionConnectivity);