SiStripRegionCabling.h

Go to the documentation of this file.

#ifndef CalibTracker_SiStripConnectivity_SiStripRegionCabling_H
#define CalibTracker_SiStripConnectivity_SiStripRegionCabling_H
#define _USE_MATH_DEFINES

#include <boost/cstdint.hpp>
#include "CondFormats/SiStripObjects/interface/FedChannelConnection.h"
#include "DataFormats/Common/interface/Handle.h"
#include <sstream>
#include <vector>
#include <map>
#include <cmath>

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