BTLDetId.cc

Go to the documentation of this file.

#include "DataFormats/ForwardDetId/interface/BTLDetId.h"

int BTLDetId::iphi( CrysLayout lay ) const {
  int kCrystalsInPhi = 1;
  switch ( lay ) {
    case CrysLayout::tile : {
      kCrystalsInPhi = kCrystalsInPhiTile ;
      break ;
    }
    case CrysLayout::bar  : {
      kCrystalsInPhi = kCrystalsInPhiBar ;
      break ;
    }
    case CrysLayout::barzflat  : {
      kCrystalsInPhi = kCrystalsInPhiBarZ ;
      break ;
    }
    case CrysLayout::barphiflat  : {
      kCrystalsInPhi = kCrystalsInPhiBarPhi ;
      break ;
    }
    default: {
      break ;
    }
  }
  return kCrystalsInPhi * ( ( mtdRR()%HALF_ROD > 0 ? mtdRR()%HALF_ROD : HALF_ROD ) -1 ) 
    + ( crystal()%kCrystalsInPhi > 0 ? crystal()%kCrystalsInPhi : kCrystalsInPhi ) ;
}

int BTLDetId::ietaAbs( CrysLayout lay ) const {
  int kCrystalsInEta = 1, kCrystalsInPhi = 1;
  std::array<int,4> kTypeBoundaries = {};
  switch ( lay ) {
    case CrysLayout::tile : {
      kCrystalsInEta = kCrystalsInEtaTile;
      kCrystalsInPhi = kCrystalsInPhiTile;
      kTypeBoundaries = kTypeBoundariesReference;
      break ;
    }
    case CrysLayout::bar  : {
      kCrystalsInEta = kCrystalsInEtaBar ;
      kCrystalsInPhi = kCrystalsInPhiBar ;
      kTypeBoundaries = kTypeBoundariesReference;
      break ;
    }
    case CrysLayout::barzflat  : {
      kCrystalsInEta = kCrystalsInEtaBarZ ;
      kCrystalsInPhi = kCrystalsInPhiBarZ ;
      kTypeBoundaries = kTypeBoundariesBarZflat ;
      break ;
    }
    case CrysLayout::barphiflat  : {
      kCrystalsInEta = kCrystalsInEtaBarPhi ;
      kCrystalsInPhi = kCrystalsInPhiBarPhi ;
      kTypeBoundaries = kTypeBoundariesBarPhiFlat ;
      break ;
    }
    default: {
      break ;
    }
  }
  int etaRowInModule = zside() > 0 ? (crystal()-1)/kCrystalsInPhi + 1 : kCrystalsInEta - (crystal()-1)/kCrystalsInPhi;
  return kCrystalsInEta * ( module() -1 ) 
    + kCrystalsInEta * kTypeBoundaries[(modType() -1)] 
    + etaRowInModule; 
}

int BTLDetId::hashedIndex( CrysLayout lay ) const { 
  int max_iphi = 1, max_ieta = 1;
  switch ( lay ) {
    case CrysLayout::tile : {
      max_iphi = MAX_IPHI_TILE;
      max_ieta = MAX_IETA_TILE;
      break ;
    }
    case CrysLayout::bar : {
      max_iphi = MAX_IPHI_BAR;
      max_ieta = MAX_IETA_BAR;
      break ;
    }
    case CrysLayout::barzflat : {
      max_iphi = MAX_IPHI_BARZFLAT;
      max_ieta = MAX_IETA_BARZFLAT;
      break ;
    }
    case CrysLayout::barphiflat : {
      max_iphi = MAX_IPHI_BARPHIFLAT;
      max_ieta = MAX_IETA_BARPHIFLAT;
      break ;
    }
    default: {
      break ;
    }
  }
  return (max_ieta + ( zside() > 0 ? ietaAbs( lay ) - 1 : -ietaAbs( lay ) ) )*max_iphi+ iphi( lay ) - 1;
}

BTLDetId BTLDetId::getUnhashedIndex( int hi, CrysLayout lay ) const {
  int max_iphi =1 ,max_ieta = 1, nphi = 0, keta = 0, tmphi = hi + 1;
  std::array<int,4> kTypeBoundaries = {};
  switch ( lay ) {
    case CrysLayout::tile : {
      max_iphi = MAX_IPHI_TILE;
      max_ieta = MAX_IETA_TILE;
      nphi = kCrystalsInPhiTile;
      keta = kCrystalsInEtaTile;
      kTypeBoundaries = kTypeBoundariesReference;
      break ;
    }
    case CrysLayout::bar : {
      max_iphi = MAX_IPHI_BAR;
      max_ieta = MAX_IETA_BAR;
      nphi = kCrystalsInPhiBar;
      keta = kCrystalsInEtaBar;
      kTypeBoundaries = kTypeBoundariesReference;
      break ;
    }
    case CrysLayout::barzflat : {
      max_iphi = MAX_IPHI_BARZFLAT;
      max_ieta = MAX_IETA_BARZFLAT;
      nphi = kCrystalsInPhiBarZ;
      keta = kCrystalsInEtaBarZ;
      kTypeBoundaries = kTypeBoundariesBarZflat;
      break ;
    }
    case CrysLayout::barphiflat : {
      max_iphi = MAX_IPHI_BARPHIFLAT;
      max_ieta = MAX_IETA_BARPHIFLAT;
      nphi = kCrystalsInPhiBarPhi;
      keta = kCrystalsInEtaBarPhi;
      kTypeBoundaries = kTypeBoundariesBarPhiFlat;
      break ;
    }
    default: {
      break ;
    }
  }
  int zside = 0, rod = 0, module = 0, modtype = 1, crystal = 0;
  if ( tmphi > max_ieta*max_iphi ) { zside = 1; }
  int ip = (tmphi-1)%max_iphi+1;
  int ie = (tmphi-1)/max_iphi - max_ieta;
  ie = ( zside == 1 ? ie + 1 : -ie ) ; 
  rod = (ip-1)/nphi+1;
  module = (ie-1)/keta+1 ; 
  if ( module > kTypeBoundaries[1] ) { modtype = (module > kTypeBoundaries[2] ? 3 : 2 ) ;  }
  if ( modtype > 1 ) { module = module - kTypeBoundaries[modtype-1]; }
  crystal = zside == 1 ? ((ip-1)%nphi+1)+((ie-1)%keta)*nphi : ((ip-1)%nphi+1)+(keta - 1 - (ie-1)%keta)*nphi ;
  return  BTLDetId( zside, rod, module, modtype, crystal);
}

int BTLDetId::modulesPerType( CrysLayout lay ) const {
  int mod = kTypeBoundariesReference[1];
  if ( lay == CrysLayout::barzflat ) { mod = kTypeBoundariesBarZflat[1]; }
  else if ( lay == CrysLayout::barphiflat ) { mod = kTypeBoundariesBarPhiFlat[1]; }
  return mod;
}

BTLDetId BTLDetId::geographicalId( CrysLayout lay ) const {
  // reorganize the modules to count from 0 to 54
  //    (0 to 42 in the case of BarZflat geometry)
  // remove module type
  // remove crystal index

  int boundRef = modulesPerType(lay);

  return BTLDetId(mtdSide(),mtdRR(),module()+boundRef*(modType()-1),0,1);
}

#include <iomanip>

std::ostream& operator<< ( std::ostream& os, const BTLDetId& id ) {
  os << ( MTDDetId& ) id;
  os << " BTL " << std::endl
     << " Side        : " << id.mtdSide() << std::endl
     << " Rod         : " << id.mtdRR() << std::endl
     << " Module      : " << id.module() << std::endl
     << " Crystal type: " << id.modType() << std::endl
     << " Crystal     : " << id.crystal() << std::endl;
  return os;
}