SiStripRegionCabling.h

Go to the documentation of this file.

#ifndef CalibTracker_SiStripConnectivity_SiStripRegionCabling_H
#define CalibTracker_SiStripConnectivity_SiStripRegionCabling_H
#define _USE_MATH_DEFINES
 
#include "CondFormats/SiStripObjects/interface/FedChannelConnection.h"
#include "DataFormats/Common/interface/Handle.h"
#include <cmath>
#include <map>
#include <sstream>
#include <vector>
#include <cstdint>
 
class SiStripRegionCabling {
public:
  enum SubDet { TIB = 0, TOB = 1, TID = 2, TEC = 3, ALLSUBDETS = 4 };
  enum Layer { TIBLAYERS = 4, TOBLAYERS = 6, TIDLAYERS = 3, TECLAYERS = 9, ALLLAYERS = 10 };
 
  typedef std::pair<uint32_t, std::vector<FedChannelConnection>> Element;
  typedef std::vector<Element> ElementCabling;
  typedef std::vector<ElementCabling> WedgeCabling;
  typedef std::vector<WedgeCabling> RegionCabling;
  typedef std::vector<RegionCabling> Cabling;
 
  typedef std::pair<double, double> Position;
  typedef std::pair<uint32_t, uint32_t> PositionIndex;
 
  typedef uint32_t ElementIndex;
 
  SiStripRegionCabling(const uint32_t, const uint32_t, const double);
 
  ~SiStripRegionCabling() {}
 
  inline void setRegionCabling(const Cabling &);
 
  inline const Cabling &getRegionCabling() const;
 
  inline const uint32_t etadivisions() const;
 
  inline const uint32_t phidivisions() const;
 
  inline const std::pair<double, double> regionDimensions() const;
 
  inline const Position position(const uint32_t) const;
 
  inline const Position position(const PositionIndex) const;
 
  inline const PositionIndex positionIndex(const uint32_t) const;
 
  const PositionIndex positionIndex(const Position) const;
 
  const uint32_t region(const Position) const;
 
  inline const uint32_t region(const PositionIndex) const;
 
  PositionIndex increment(const PositionIndex, int, int) const;
 
  inline static const ElementIndex elementIndex(const uint32_t region, const SubDet, const uint32_t layer);
 
  inline const ElementIndex elementIndex(const PositionIndex, const SubDet, const uint32_t layer) const;
 
  inline const ElementIndex elementIndex(const Position, const SubDet, const uint32_t layer) const;
 
  inline static const uint32_t layer(const ElementIndex);
 
  inline static const SubDet subdet(const ElementIndex);
 
  inline static const uint32_t region(const ElementIndex);
 
  //  /** Methods for extracting det-id information */
  //
  static const SubDet subdetFromDetId(const uint32_t detid);
 
  void print(std::stringstream &) const;
 
private:
  SiStripRegionCabling() { ; }
 
  int etadivisions_;
  int phidivisions_;
 
  double etamax_;
 
  Cabling regioncabling_;
};
 
inline void SiStripRegionCabling::setRegionCabling(const Cabling &regioncabling) { regioncabling_ = regioncabling; }
 
inline const SiStripRegionCabling::Cabling &SiStripRegionCabling::getRegionCabling() const { return regioncabling_; }
 
inline const uint32_t SiStripRegionCabling::etadivisions() const { return static_cast<uint32_t>(etadivisions_); }
 
inline const uint32_t SiStripRegionCabling::phidivisions() const { return static_cast<uint32_t>(phidivisions_); }
 
inline const std::pair<double, double> SiStripRegionCabling::regionDimensions() const {
  return std::pair<double, double>((2. * etamax_) / etadivisions_, 2. * M_PI / phidivisions_);
}
 
inline const SiStripRegionCabling::Position SiStripRegionCabling::position(const uint32_t region) const {
  PositionIndex index = positionIndex(region);
  return position(index);
}
 
inline const SiStripRegionCabling::Position SiStripRegionCabling::position(const PositionIndex index) const {
  return Position(regionDimensions().first * (index.first + .5) - etamax_,
                  regionDimensions().second * (index.second + .5) - M_PI);
}
 
inline const SiStripRegionCabling::PositionIndex SiStripRegionCabling::positionIndex(const uint32_t region) const {
  return PositionIndex(region / phidivisions_, region % phidivisions_);
}
 
inline const uint32_t SiStripRegionCabling::region(const PositionIndex index) const {
  return index.first * phidivisions_ + index.second;
}
 
inline const uint32_t SiStripRegionCabling::elementIndex(const uint32_t region,
                                                         const SubDet subdet,
                                                         const uint32_t layer) {
  return region * ALLSUBDETS * ALLLAYERS + subdet * ALLLAYERS + layer;
}
 
inline const uint32_t SiStripRegionCabling::elementIndex(const PositionIndex index,
                                                         const SubDet subdet,
                                                         const uint32_t layer) const {
  return elementIndex(region(index), subdet, layer);
}
 
inline const uint32_t SiStripRegionCabling::elementIndex(const Position position,
                                                         const SubDet subdet,
                                                         const uint32_t layer) const {
  return elementIndex(region(position), subdet, layer);
}
 
inline const uint32_t SiStripRegionCabling::layer(const uint32_t index) { return index % ALLLAYERS; }
 
inline const SiStripRegionCabling::SubDet SiStripRegionCabling::subdet(const uint32_t index) {
  return static_cast<SiStripRegionCabling::SubDet>((index / ALLLAYERS) % ALLSUBDETS);
}
 
inline const uint32_t SiStripRegionCabling::region(const uint32_t index) { return index / (ALLSUBDETS * ALLLAYERS); }
 
#endif