EcalFenixTcpFgvbEE.cc

Go to the documentation of this file.

#include "CondFormats/EcalObjects/interface/EcalTPGFineGrainTowerEE.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <SimCalorimetry/EcalTrigPrimAlgos/interface/EcalFenixTcpFgvbEE.h>
#include <iostream>

//---------------------------------------------------------------
EcalFenixTcpFgvbEE::EcalFenixTcpFgvbEE(int maxNrSamples) {
  indexLut_.resize(maxNrSamples);
}  //---------------------------------------------------------------
EcalFenixTcpFgvbEE::~EcalFenixTcpFgvbEE() {}
//---------------------------------------------------------------
void EcalFenixTcpFgvbEE::process(std::vector<std::vector<int>> &bypasslin_out,
                                 int nStr,
                                 int bitMask,
                                 std::vector<int> &output) {
  //  std::vector<int> indexLut(output.size());

  for (unsigned int i = 0; i < output.size(); i++) {
    output[i] = 0;
    indexLut_[i] = 0;
  }

  for (unsigned int i = 0; i < output.size(); i++) {
    for (int istrip = 0; istrip < nStr; istrip++) {
      int res = (bypasslin_out[istrip])[i];
      res = (res >> bitMask) & 1;  // res is FGVB at this stage
      indexLut_[i] = indexLut_[i] | (res << istrip);
    }
    indexLut_[i] = indexLut_[i] | (nStr << 5);

    int mask = 1 << indexLut_[i];
    output[i] = fgee_lut_ & mask;
    if (output[i] > 0)
      output[i] = 1;
  }
  return;
}

//-------------------------------------------------------------------

void EcalFenixTcpFgvbEE::setParameters(uint32_t towid, const EcalTPGFineGrainTowerEE *ecaltpgFineGrainTowerEE) {
  const EcalTPGFineGrainTowerEEMap &fgee_map = ecaltpgFineGrainTowerEE->getMap();

  EcalTPGFineGrainTowerEEMapIterator it = fgee_map.find(towid);
  if (it != fgee_map.end())
    fgee_lut_ = (*it).second;
  else
    edm::LogWarning("EcalTPG") << " could not find EcalTPGFineGrainTowerEEMap for " << towid;
}