HFRaddamTask.cc

Go to the documentation of this file.

#include "DQM/HcalTasks/interface/HFRaddamTask.h"

using namespace hcaldqm;
using namespace hcaldqm::constants;
using namespace hcaldqm::filter;

HFRaddamTask::HFRaddamTask(edm::ParameterSet const& ps) : DQTask(ps) {
  //    List all the DetIds
  _vDetIds.push_back(HcalDetId(HcalForward, -30, 35, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, -30, 71, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, -32, 15, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, -32, 51, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, -34, 35, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, -34, 71, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, -36, 15, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, -36, 51, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, -38, 35, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, -38, 71, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, -40, 15, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, -40, 51, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, -41, 35, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, -41, 71, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, -30, 15, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, -30, 51, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, -32, 35, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, -32, 71, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, -34, 15, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, -34, 51, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, -36, 35, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, -36, 71, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, -38, 15, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, -38, 51, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, -40, 35, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, -40, 71, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, -41, 15, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, -41, 51, 2));

  _vDetIds.push_back(HcalDetId(HcalForward, 30, 21, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, 30, 57, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, 32, 1, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, 32, 37, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, 34, 21, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, 34, 57, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, 36, 1, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, 36, 37, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, 38, 21, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, 38, 57, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, 40, 35, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, 40, 71, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, 41, 19, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, 41, 55, 1));
  _vDetIds.push_back(HcalDetId(HcalForward, 30, 1, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, 30, 37, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, 32, 21, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, 32, 57, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, 34, 1, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, 34, 37, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, 36, 21, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, 36, 57, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, 38, 1, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, 38, 37, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, 40, 19, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, 40, 55, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, 41, 35, 2));
  _vDetIds.push_back(HcalDetId(HcalForward, 41, 71, 2));

  //    tags
  _tagHF = ps.getUntrackedParameter<edm::InputTag>("tagHF", edm::InputTag("hcalDigis"));
  _taguMN = ps.getUntrackedParameter<edm::InputTag>("taguMN", edm::InputTag("hcalDigis"));
  _tokHF = consumes<HFDigiCollection>(_tagHF);
  _tokuMN = consumes<HcalUMNioDigi>(_taguMN);
}

/* virtual */ void HFRaddamTask::bookHistograms(DQMStore::IBooker& ib, edm::Run const& r, edm::EventSetup const& es) {
  //    Initialize all the Single Containers
  for (std::vector<HcalDetId>::const_iterator it = _vDetIds.begin(); it != _vDetIds.end(); ++it) {
    _vcShape.push_back(ContainerSingle1D(_name,
                                         "Shape",
                                         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
                                         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000)));
  }

  DQTask::bookHistograms(ib, r, es);
  char aux[200];
  for (unsigned int i = 0; i < _vDetIds.size(); i++) {
    sprintf(aux, "ieta%diphi%dd%d", _vDetIds[i].ieta(), _vDetIds[i].iphi(), _vDetIds[i].depth());
    _vcShape[i].book(ib, _subsystem, aux);
  }
}

/* virtual */ void HFRaddamTask::_process(edm::Event const& e, edm::EventSetup const& es) {
  edm::Handle<HFDigiCollection> chf;
  if (!e.getByToken(_tokHF, chf))
    _logger.dqmthrow("Collection HFDigiCollection isn't avalaible" + _tagHF.label() + " " + _tagHF.instance());

  for (HFDigiCollection::const_iterator it = chf->begin(); it != chf->end(); ++it) {
    const HFDataFrame digi = (const HFDataFrame)(*it);
    for (unsigned int i = 0; i < _vDetIds.size(); i++)
      if (digi.id() == _vDetIds[i]) {
        for (int j = 0; j < digi.size(); j++)
          _vcShape[i].fill(j, digi.sample(j).nominal_fC() - 2.5);
      }
  }
}

/* virtual */ bool HFRaddamTask::_isApplicable(edm::Event const& e) {
  if (_ptype == fOnline) {
    edm::Handle<HcalUMNioDigi> cumn;
    if (!e.getByToken(_tokuMN, cumn))
      return false;

    //  event type check first
    uint8_t eventType = cumn->eventType();
    if (eventType != constants::EVENTTYPE_LASER)
      return false;

    //  check if this analysis task is of the right laser type
    uint32_t laserType = cumn->valueUserWord(0);
    if (laserType == constants::tHFRaddam)
      return true;
  } else {
    //  local, just return true as all the settings will be done in cfg
    return true;
  }

  return false;
}

DEFINE_FWK_MODULE(HFRaddamTask);