HcalFiberPattern.cc

Go to the documentation of this file.

#include "IORawData/CaloPatterns/interface/HcalFiberPattern.h"
#include "FWCore/Utilities/interface/Exception.h"

static inline int setIf(const std::string& name,const std::map<std::string, std::string>& params) {
  std::map<std::string, std::string>::const_iterator j=params.find(name);
  if (j==params.end()) throw cms::Exception("InvalidFormat") << "Missing parameter '" << name << "'";
  else return strtol(j->second.c_str(),0,0);
} 

HcalFiberPattern::HcalFiberPattern(const std::map<std::string, std::string>& params, const std::vector<uint32_t>& data) : pattern_(data) {
  crate_=setIf("CRATE",params);
  slot_=setIf("SLOT",params);
  fiber_=setIf("FIBER",params);
  dcc_=setIf("DCC",params);
  spigot_=setIf("SPIGOT",params);
  tb_=setIf("TOPBOTTOM",params);
}

HcalQIESample HcalFiberPattern::unpack(int bc, int fc) {
  uint32_t w1=pattern_[bc*2]; // lsw
  uint32_t w2=pattern_[bc*2+1]; // msw

  int adc=0, capid=0;
  bool dv=(w1&0x10000)!=0;
  bool er=(w1&0x20000)!=0;
  
  switch (fc) {
  case (0):
    adc=(w2&0xFE00)>>9;
    capid=(w1&0x0180)>>7;
    break;
  case (1):
    adc=(w2&0xFE)>>1;
    capid=(w1&0x0060)>>5;
    break;
  case (2):
    adc=(w1&0xFE00)>>9;
    capid=(w1&0x0018)>>3;
    break;
  default:
    break;
  }
  return HcalQIESample(adc,capid,fiber_,fc,dv,er);
}

std::vector<HcalQIESample> HcalFiberPattern::getSamples(int bunch, int npresamples, int nsamples, int fiberChan) {
  if (bunch<npresamples) throw cms::Exception("InvalidArgument") << "Asked for " << npresamples << " presamples with event at bunch " << bunch;
  if (nsamples-npresamples+bunch>=(int)(pattern_.size()/2)) throw cms::Exception("InvalidArgument") << "Asked for " << nsamples << " with event at " << bunch << " and " << npresamples << " presamples, but only " << pattern_.size()/2 << " bunches are available";

  std::vector<HcalQIESample> retval;
  retval.reserve(nsamples);

  for (int i=0; i<nsamples; i++) {
    int bc=bunch+i-npresamples;
    retval.push_back(unpack(bc,fiberChan));
  }
  return retval;
}

HcalElectronicsId HcalFiberPattern::getId(int fiberChan) {
  HcalElectronicsId retval(fiberChan,fiber_,spigot_,dcc_);
  retval.setHTR(crate_,slot_,tb_);
  return retval;
}