HcalZSAlgoRealistic.cc

Go to the documentation of this file.

#include "SimCalorimetry/HcalZeroSuppressionProducers/interface/HcalZSAlgoRealistic.h"
#include <iostream>

HcalZSAlgoRealistic::HcalZSAlgoRealistic(bool mp,
                                         bool use1ts,
                                         int levelHB,
                                         int levelHE,
                                         int levelHO,
                                         int levelHF,
                                         std::pair<int, int> HBsearchTS,
                                         std::pair<int, int> HEsearchTS,
                                         std::pair<int, int> HOsearchTS,
                                         std::pair<int, int> HFsearchTS)
    : HcalZeroSuppressionAlgo(mp),
      use1ts_(use1ts),
      thresholdHB_(levelHB),
      thresholdHE_(levelHE),
      thresholdHO_(levelHO),
      thresholdHF_(levelHF),
      HBsearchTS_(HBsearchTS),
      HEsearchTS_(HEsearchTS),
      HOsearchTS_(HOsearchTS),
      HFsearchTS_(HFsearchTS) {
  usingDBvalues = false;
}

HcalZSAlgoRealistic::HcalZSAlgoRealistic(bool mp,
                                         bool use1ts,
                                         std::pair<int, int> HBsearchTS,
                                         std::pair<int, int> HEsearchTS,
                                         std::pair<int, int> HOsearchTS,
                                         std::pair<int, int> HFsearchTS)
    : HcalZeroSuppressionAlgo(mp),
      use1ts_(use1ts),
      HBsearchTS_(HBsearchTS),
      HEsearchTS_(HEsearchTS),
      HOsearchTS_(HOsearchTS),
      HFsearchTS_(HFsearchTS) {
  thresholdHB_ = -1;
  thresholdHE_ = -1;
  thresholdHO_ = -1;
  thresholdHF_ = -1;
  usingDBvalues = true;
}

template <class Digi>
bool HcalZSAlgoRealistic::keepMe(const Digi &inp, int start, int finish, int threshold, uint32_t zsmask) const {
  if ((usingDBvalues) && (threshold < 0) && (m_dbService != nullptr)) {
    threshold = (m_dbService->getHcalZSThreshold(inp.id()))->getValue();
  }

  // determine the sum of 2 timeslices
  for (int i = start; i < finish; i++) {
    if ((zsmask & (1 << i)) != 0)
      continue;
    if ((inp[i].adc() + inp[i + 1].adc()) >= threshold)
      return true;
  }
  return false;
}

// zs mask not used for QIE10,11

template <>
bool HcalZSAlgoRealistic::keepMe<QIE10DataFrame>(
    const QIE10DataFrame &inp, int start, int finish, int threshold, uint32_t zsmask) const {
  if ((usingDBvalues) && (threshold < 0) && (m_dbService != nullptr)) {
    threshold = (m_dbService->getHcalZSThreshold(inp.id()))->getValue();
  }

  // determine the sum of 2 timeslices
  for (int i = start; i < finish; i++) {
    if ((inp[i].adc() + inp[i + 1].adc()) >= threshold)
      return true;
  }
  return false;
}

template <>
bool HcalZSAlgoRealistic::keepMe<QIE11DataFrame>(
    const QIE11DataFrame &inp, int start, int finish, int threshold, uint32_t zsmask) const {
  if ((usingDBvalues) && (threshold < 0) && (m_dbService != nullptr)) {
    threshold = (m_dbService->getHcalZSThreshold(inp.id()))->getValue();
  }

  // sum of 2 timeslices or just single one
  for (int i = start; i < finish; i++) {
    int sum = inp[i].adc();
    if (!use1ts_)
      sum += inp[i + 1].adc();
    if (sum >= threshold)
      return true;
  }
  return false;
}

bool HcalZSAlgoRealistic::shouldKeep(const HBHEDataFrame &digi) const {
  if (digi.id().subdet() == HcalBarrel) {
    int start = std::max(0, HBsearchTS_.first);
    int finish = std::min(digi.size() - 1, HBsearchTS_.second);
    return keepMe(digi, start, finish, thresholdHB_, digi.zsCrossingMask());
  } else {
    int start = std::max(0, HEsearchTS_.first);
    int finish = std::min(digi.size() - 1, HEsearchTS_.second);
    return keepMe(digi, start, finish, thresholdHE_, digi.zsCrossingMask());
  }
}

bool HcalZSAlgoRealistic::shouldKeep(const HODataFrame &digi) const {
  int start = std::max(0, HOsearchTS_.first);
  int finish = std::min(digi.size() - 1, HOsearchTS_.second);
  return keepMe(digi, start, finish, thresholdHO_, digi.zsCrossingMask());
}

bool HcalZSAlgoRealistic::shouldKeep(const HFDataFrame &digi) const {
  int start = std::max(0, HFsearchTS_.first);
  int finish = std::min(digi.size() - 1, HFsearchTS_.second);
  return keepMe(digi, start, finish, thresholdHF_, digi.zsCrossingMask());
}

bool HcalZSAlgoRealistic::shouldKeep(const QIE10DataFrame &digi) const {
  int start = std::max(0, HFsearchTS_.first);
  int finish = std::min((int)digi.samples() - 1, HFsearchTS_.second);
  return keepMe(digi, start, finish, thresholdHF_, 0);
}

bool HcalZSAlgoRealistic::shouldKeep(const QIE11DataFrame &digi) const {
  HcalDetId hid(digi.id());
  if (hid.subdet() == HcalBarrel) {
    int start = std::max(0, HBsearchTS_.first);
    int finish = std::min(digi.samples() - 1, HBsearchTS_.second);
    return keepMe(digi, start, finish, thresholdHB_, 0);
  } else {
    int start = std::max(0, HEsearchTS_.first);
    int finish = std::min(digi.samples() - 1, HEsearchTS_.second);
    return keepMe(digi, start, finish, thresholdHE_, 0);
  }
}