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 <boost/cstdint.hpp>
#include <cmath>
#include <map>
#include <sstream>
#include <vector>

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