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.book1D( 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.book1D( 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.book1D( 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.book1D( 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.book1D( 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.book1D( 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.book1D( 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.book1D( 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.book1D( 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.book1D( 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.book1D( 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.book1D( 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);