SiStripRegionCabling.cc

Go to the documentation of this file.

#include "CalibFormats/SiStripObjects/interface/SiStripRegionCabling.h"
#include "FWCore/Utilities/interface/typelookup.h"
#include "DataFormats/SiStripDetId/interface/TIBDetId.h"
#include "DataFormats/SiStripDetId/interface/TOBDetId.h"
#include "DataFormats/SiStripDetId/interface/TECDetId.h"
#include "DataFormats/SiStripDetId/interface/TIDDetId.h"

using namespace sistrip;

SiStripRegionCabling::SiStripRegionCabling(const uint32_t etadivisions, const uint32_t phidivisions, const double etamax) :

  etadivisions_(static_cast<int>(etadivisions)),
  phidivisions_(static_cast<int>(phidivisions)),
  etamax_(etamax),
  regioncabling_()

{}

const SiStripRegionCabling::PositionIndex SiStripRegionCabling::positionIndex(const Position position) const {
  int eta = static_cast<int>((position.first+etamax_)*etadivisions_/(2.*etamax_));
  int phi = static_cast<int>((position.second+M_PI)*phidivisions_/(2.*M_PI));
  return PositionIndex(static_cast<uint32_t>(eta),static_cast<uint32_t>(phi));
}

const uint32_t SiStripRegionCabling::region(const Position position) const {
  PositionIndex index = positionIndex(position); 
  return region(index);
}

SiStripRegionCabling::PositionIndex SiStripRegionCabling::increment(const PositionIndex index, int deta, int dphi) const {
  
  int eta = static_cast<int>(index.first)+deta;
  if (eta > etadivisions_-1) eta = etadivisions_-1;
  else if (eta < 0) eta = 0;
 
  int phi = static_cast<int>(index.second)+dphi;
  while (phi<0) phi+=phidivisions_;
  while (phi>phidivisions_-1) phi-=phidivisions_;

  return PositionIndex(static_cast<uint32_t>(eta),static_cast<uint32_t>(phi));
}

const SiStripRegionCabling::SubDet SiStripRegionCabling::subdetFromDetId(const uint32_t detid) {

  SiStripDetId::SubDetector subdet = SiStripDetId(detid).subDetector();
  if (subdet == 3) return TIB;
  else if (subdet == 4) return TID;
  else if (subdet == 5) return TOB;
  else if (subdet == 6) return TEC;
  else return ALLSUBDETS;
}

const uint32_t SiStripRegionCabling::layerFromDetId(const uint32_t detid) {
 
  SiStripRegionCabling::SubDet subdet = subdetFromDetId(detid);
  if (subdet == TIB) return TIBDetId(detid).layer();
  else if (subdet == TID) return TIDDetId(detid).wheel();
  else if (subdet == TOB) return TOBDetId(detid).layer(); 
  else if (subdet == TEC) return TECDetId(detid).wheel();
  else return ALLLAYERS;
}

const uint32_t SiStripRegionCabling::physicalLayerFromDetId(const uint32_t detid) {
  return physicalLayer(subdetFromDetId(detid),layerFromDetId(detid));
}

const uint32_t SiStripRegionCabling::physicalLayer(const SubDet subdet, const uint32_t layer) {

  if (subdet == TIB) return layer;
  else if (subdet == TOB) return TIBLAYERS + layer;
  else if (subdet == TID) return layer;
  else if (subdet == TEC) return TIDLAYERS + layer;
  else return ALLLAYERS;
}

// -----------------------------------------------------------------------------
//
void SiStripRegionCabling::print( std::stringstream& ss ) const {
  uint32_t valid = 0;
  uint32_t total = 0;
  ss << "[SiStripRegionCabling::" << __func__ << "] Printing REGION cabling:" << std::endl;
  ss << "Printing cabling for " << regioncabling_.size() << " regions" << std::endl;
  Cabling::const_iterator id = regioncabling_.begin();
  Cabling::const_iterator jd = regioncabling_.end();
  for ( ; id != jd; ++id ) {
    ss << "Printing cabling for " << id->size()
       << " regions for partition " << static_cast<int32_t>( id - regioncabling_.begin() )
       << std::endl;
    RegionCabling::const_iterator ir = id->begin();
    RegionCabling::const_iterator jr = id->end();
    for ( ; ir != jr; ++ir ) {
      ss << "Printing cabling for " << ir->size()
     << " wedges for region " << static_cast<int32_t>( ir - id->begin() )
     << std::endl;
      WedgeCabling::const_iterator iw = ir->begin();
      WedgeCabling::const_iterator jw = ir->end();
      for ( ; iw != jw; ++iw ) {
    ss << "Printing cabling for " << iw->size()
       << " elements for wedge " << static_cast<int32_t>( iw - ir->begin() )
       << std::endl;
    ElementCabling::const_iterator ie = iw->begin();
    ElementCabling::const_iterator je = iw->end();
    for ( ; ie != je; ++ie ) {
      ss << "Printing cabling for " << ie->second.size()
         << " connections for element (DetId) " << ie->first
         << std::endl;
      std::vector<FedChannelConnection>::const_iterator ic = ie->second.begin();
      std::vector<FedChannelConnection>::const_iterator jc = ie->second.end();
      for ( ; ic != jc; ++ic ) {
        if ( ic->isConnected() ) { valid++; }
        total++;
        ic->print(ss);
        ss << std::endl;
      }
    }
      }
    }
  }
  ss << "Number of connected:   " << valid << std::endl
     << "Number of connections: " << total << std::endl;
}