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>

//---------------------------------------------------------------
void EcalFenixTcpFgvbEE::process(std::vector<std::vector<int>> &bypasslin_out,
                                 int nStr,
                                 int bitMask,
                                 std::vector<int> &output) {
  for (unsigned int i = 0; i < output.size(); i++) {
    output[i] = 0;
  }

  //Return if fgee_lut_ is 0
  if (fgee_lut_ == 0) {
    return;
  }

  int indexLut = 0;
  for (unsigned int i = 0; i < output.size(); i++) {
    indexLut = 0;
    for (int istrip = 0; istrip < nStr; istrip++) {
      int res = (bypasslin_out[istrip])[i];
      res = (res >> bitMask) & 1;  // res is FGVB at this stage
      indexLut = indexLut | (res << istrip);
    }
    int mask = 1 << indexLut;
    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;
}