HcalFrontEndId.cc

Go to the documentation of this file.

#include "DataFormats/HcalDetId/interface/HcalFrontEndId.h"
#include <iomanip>
#include <sstream>
#include <iostream>
#include <cstdlib>


HcalFrontEndId::HcalFrontEndId(const std::string& rbx,int rm,int pixel,int rmfiber,int fiberchannel,int qie,int adc)
{
  hcalFrontEndId_=0;

  if (rbx.size()<5) return;
  if (rm<1 || rm>5) return; //changed to 5 to incorporate CALIB channels which define RM = 5
  if (pixel<0 || pixel>19) return;
  if (rmfiber<1 || rmfiber>8) return;
  if (fiberchannel<0 || fiberchannel>2) return;
  if (qie<1 || qie>4) return;
  if (adc<0 || adc>5) return;

  int num=-1;
  if      (!rbx.compare(0,3,"HBM" ))  { num=0         + atoi(rbx.substr(3,2).c_str())-1; }
  else if (!rbx.compare(0,3,"HBP" ))  { num=18        + atoi(rbx.substr(3,2).c_str())-1; }
  else if (!rbx.compare(0,3,"HEM" ))  { num=18*2      + atoi(rbx.substr(3,2).c_str())-1; }
  else if (!rbx.compare(0,3,"HEP" ))  { num=18*3      + atoi(rbx.substr(3,2).c_str())-1; }
  else if (!rbx.compare(0,4,"HO2M"))  { num=18*4      + atoi(rbx.substr(4,2).c_str())-1; }
  else if (!rbx.compare(0,4,"HO1M"))  { num=18*4+12   + atoi(rbx.substr(4,2).c_str())-1; }
  else if (!rbx.compare(0,3,"HO0" ))  { num=18*4+12*2 + atoi(rbx.substr(3,2).c_str())-1; }
  else if (!rbx.compare(0,4,"HO1P"))  { num=18*4+12*3 + atoi(rbx.substr(4,2).c_str())-1; }
  else if (!rbx.compare(0,4,"HO2P"))  { num=18*4+12*4 + atoi(rbx.substr(4,2).c_str())-1; }
  else if (!rbx.compare(0,3,"HFM" ))  { num=18*4+12*5 + atoi(rbx.substr(3,2).c_str())-1; }
  else if (!rbx.compare(0,3,"HFP" ))  { num=18*4+12*6 + atoi(rbx.substr(3,2).c_str())-1; }
  else return;

  hcalFrontEndId_|=((adc+1)&0x7);
  hcalFrontEndId_|=((qie-1)&0x3)<<3;
  hcalFrontEndId_|=(fiberchannel&0x3)<<5;
  hcalFrontEndId_|=((rmfiber-1)&0x7)<<7;
  hcalFrontEndId_|=(pixel&0x1F)<<10;
  hcalFrontEndId_|=((rm-1)&0x7)<<15;
  hcalFrontEndId_|=(num&0xFF)<<18;
}

HcalFrontEndId::~HcalFrontEndId()
{
}

std::string HcalFrontEndId::rbx() const
{
  std::string subdets[11]={"HBM","HBP","HEM","HEP","HO2M","HO1M","HO0","HO1P","HO2P","HFM","HFP"};

  int box=hcalFrontEndId_>>18;
  int num=-1;
  int subdet_index=-1;
  if (box<18*4) {
    num=box%18;
    subdet_index=(box-num)/18;
  }
  else {
    num=(box-18*4)%12;
    subdet_index=4+(box-18*4-num)/12;
  }
  std::stringstream tempss;
  tempss << std::setw(2) << std::setfill('0') << num+1;
  return subdets[subdet_index]+tempss.str();
}

std::ostream& operator<<(std::ostream& s,const HcalFrontEndId& id) {
  return s << id.rbx() << id.rm() << '[' << id.rmFiber() << '/' << id.fiberChannel() 
       << "] pix=" << id.pixel()
       << " qiecard=" << id.qieCard() << " adc=" << id.adc();
}