ZDCQIE10Task.cc

Go to the documentation of this file.

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

ZDCQIE10Task::ZDCQIE10Task(edm::ParameterSet const& ps) {
  //    tags
  _tagQIE10 = ps.getUntrackedParameter<edm::InputTag>("tagQIE10", edm::InputTag("hcalDigis", "ZDC"));
  _tokQIE10 = consumes<QIE10DigiCollection>(_tagQIE10);
}

/* virtual */ void ZDCQIE10Task::bookHistograms(DQMStore::IBooker& ib, edm::Run const& r, edm::EventSetup const& es) {
  ib.cd();

  //book histos per channel
  std::string histoname;
  for (int channel = 1; channel < 6; channel++) {
    // EM Pos
    HcalZDCDetId didp(HcalZDCDetId::EM, true, channel);
    histoname = "EM_P_" + std::to_string(channel) + "_1";
    ib.setCurrentFolder("Hcal/ZDCQIE10Task/ADC_perChannel");
    _cADC_EChannel[didp()] = ib.book1DD(histoname.c_str(), histoname.c_str(), 256, 0, 256);
    _cADC_EChannel[didp()]->setAxisTitle("ADC", 1);
    _cADC_EChannel[didp()]->setAxisTitle("N", 2);
    ib.setCurrentFolder("Hcal/ZDCQIE10Task/ADC_vs_TS_perChannel");
    _cADC_vs_TS_EChannel[didp()] = ib.book1DD(histoname.c_str(), histoname.c_str(), 10, 0, 10);
    _cADC_vs_TS_EChannel[didp()]->setAxisTitle("TS", 1);
    _cADC_vs_TS_EChannel[didp()]->setAxisTitle("sum ADC", 2);

    // EM Minus
    HcalZDCDetId didm(HcalZDCDetId::EM, false, channel);
    histoname = "EM_M_" + std::to_string(channel) + "_1";
    ib.setCurrentFolder("Hcal/ZDCQIE10Task/ADC_perChannel");
    _cADC_EChannel[didm()] = ib.book1DD(histoname.c_str(), histoname.c_str(), 256, 0, 256);
    _cADC_EChannel[didm()]->setAxisTitle("ADC", 1);
    _cADC_EChannel[didm()]->setAxisTitle("N", 2);
    ib.setCurrentFolder("Hcal/ZDCQIE10Task/ADC_vs_TS_perChannel");
    _cADC_vs_TS_EChannel[didm()] = ib.book1DD(histoname.c_str(), histoname.c_str(), 10, 0, 10);
    _cADC_vs_TS_EChannel[didm()]->setAxisTitle("TS", 1);
    _cADC_vs_TS_EChannel[didm()]->setAxisTitle("sum ADC", 2);
  }

  for (int channel = 1; channel < 5; channel++) {
    // HAD Pos
    HcalZDCDetId didp(HcalZDCDetId::HAD, true, channel);
    histoname = "HAD_P_" + std::to_string(channel) + "_" + std::to_string(channel + 2);
    ib.setCurrentFolder("Hcal/ZDCQIE10Task/ADC_perChannel");
    _cADC_EChannel[didp()] = ib.book1DD(histoname.c_str(), histoname.c_str(), 256, 0, 256);
    _cADC_EChannel[didp()]->setAxisTitle("ADC", 1);
    _cADC_EChannel[didp()]->setAxisTitle("N", 2);
    ib.setCurrentFolder("Hcal/ZDCQIE10Task/ADC_vs_TS_perChannel");
    _cADC_vs_TS_EChannel[didp()] = ib.book1DD(histoname.c_str(), histoname.c_str(), 10, 0, 10);
    _cADC_vs_TS_EChannel[didp()]->setAxisTitle("TS", 1);
    _cADC_vs_TS_EChannel[didp()]->setAxisTitle("sum ADC", 2);

    // HAD Minus
    HcalZDCDetId didm(HcalZDCDetId::HAD, false, channel);
    histoname = "HAD_M_" + std::to_string(channel) + "_" + std::to_string(channel + 2);
    ib.setCurrentFolder("Hcal/ZDCQIE10Task/ADC_perChannel");
    _cADC_EChannel[didm()] = ib.book1DD(histoname.c_str(), histoname.c_str(), 256, 0, 256);
    _cADC_EChannel[didm()]->setAxisTitle("ADC", 1);
    _cADC_EChannel[didm()]->setAxisTitle("N", 2);
    ib.setCurrentFolder("Hcal/ZDCQIE10Task/ADC_vs_TS_perChannel");
    _cADC_vs_TS_EChannel[didm()] = ib.book1DD(histoname.c_str(), histoname.c_str(), 10, 0, 10);
    _cADC_vs_TS_EChannel[didm()]->setAxisTitle("TS", 1);
    _cADC_vs_TS_EChannel[didm()]->setAxisTitle("sum ADC", 2);
  }

  for (int channel = 1; channel < 17; channel++) {
    // RPD Pos
    HcalZDCDetId didp(HcalZDCDetId::RPD, true, channel);
    histoname = "RPD_P_" + std::to_string(channel) + "_2";
    ib.setCurrentFolder("Hcal/ZDCQIE10Task/ADC_perChannel");
    _cADC_EChannel[didp()] = ib.book1DD(histoname.c_str(), histoname.c_str(), 256, 0, 256);
    _cADC_EChannel[didp()]->setAxisTitle("ADC", 1);
    _cADC_EChannel[didp()]->setAxisTitle("N", 2);
    ib.setCurrentFolder("Hcal/ZDCQIE10Task/ADC_vs_TS_perChannel");
    _cADC_vs_TS_EChannel[didp()] = ib.book1DD(histoname.c_str(), histoname.c_str(), 10, 0, 10);
    _cADC_vs_TS_EChannel[didp()]->setAxisTitle("TS", 1);
    _cADC_vs_TS_EChannel[didp()]->setAxisTitle("sum ADC", 2);

    // RPD Minus
    HcalZDCDetId didm(HcalZDCDetId::RPD, false, channel);
    histoname = "RPD_M_" + std::to_string(channel) + "_2";
    ib.setCurrentFolder("Hcal/ZDCQIE10Task/ADC_perChannel");
    _cADC_EChannel[didm()] = ib.book1DD(histoname.c_str(), histoname.c_str(), 256, 0, 256);
    _cADC_EChannel[didm()]->setAxisTitle("ADC", 1);
    _cADC_EChannel[didm()]->setAxisTitle("N", 2);
    ib.setCurrentFolder("Hcal/ZDCQIE10Task/ADC_vs_TS_perChannel");
    _cADC_vs_TS_EChannel[didm()] = ib.book1DD(histoname.c_str(), histoname.c_str(), 10, 0, 10);
    _cADC_vs_TS_EChannel[didm()]->setAxisTitle("TS", 1);
    _cADC_vs_TS_EChannel[didm()]->setAxisTitle("sum ADC", 2);
  }
}

/* virtual */ void ZDCQIE10Task::analyze(edm::Event const& e, edm::EventSetup const&) {
  edm::Handle<QIE10DigiCollection> digis;
  if (!e.getByToken(_tokQIE10, digis))
    edm::LogError("Collection QIE10DigiCollection for ZDC isn't available" + _tagQIE10.label() + " " +
                  _tagQIE10.instance());

  for (auto it = digis->begin(); it != digis->end(); it++) {
    const QIE10DataFrame digi = static_cast<const QIE10DataFrame>(*it);
    HcalZDCDetId const& did = digi.detid();

    for (int i = 0; i < digi.samples(); i++) {
      // iter over all samples
      if (_cADC_EChannel.find(did()) != _cADC_EChannel.end()) {
        _cADC_EChannel[did()]->Fill(digi[i].adc());
      }
      if (_cADC_vs_TS_EChannel.find(did()) != _cADC_vs_TS_EChannel.end()) {
        _cADC_vs_TS_EChannel[did()]->Fill(i, digi[i].adc());
      }
    }
  }
}

DEFINE_FWK_MODULE(ZDCQIE10Task);