ZDCTask.cc

Go to the documentation of this file.

#include "DQM/HcalTasks/interface/ZDCTask.h"
#include <map>

using namespace hcaldqm;
using namespace hcaldqm::constants;
ZDCTask::ZDCTask(edm::ParameterSet const& ps)
{
    
    //  tags
    _tagQIE10 = ps.getUntrackedParameter<edm::InputTag>("tagQIE10",
        edm::InputTag("hcalDigis"));
    _tokQIE10 = consumes<ZDCDigiCollection>(_tagQIE10);

    //  cuts
    _cut = ps.getUntrackedParameter<double>("cut", 50.0);
    _ped = ps.getUntrackedParameter<int>("ped", 4);
}

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

    //############################## hardcode manually the zdc mapping #############################
    //############################# this follows from https://github.com/cms-sw/cmssw/blob/CMSSW_8_0_X/EventFilter/CastorRawToDigi/src/ZdcUnpacker.cc#L118
    //##############################################################################################
    //PZDC
    std::map<HcalElectronicsId,DetId> myEMap;
    HcalElectronicsId eid = HcalElectronicsId(0, 1, 0, 3);
    eid.setHTR(18, 8, 1);
    myEMap[eid]=DetId(0x54000051);//PZDC EM1

    eid = HcalElectronicsId(1, 1, 0, 3);
    eid.setHTR(18, 8, 1);
    myEMap[eid]=DetId(0x54000052);//PZDC EM2

    eid = HcalElectronicsId(2, 1, 0, 3);
    eid.setHTR(18, 8, 1);
    myEMap[eid]=DetId(0x54000053);//PZDC EM3

    eid = HcalElectronicsId(0, 2, 0, 3);
    eid.setHTR(18, 8, 1);
    myEMap[eid]=DetId(0x54000061);//PZDC HAD1

    eid = HcalElectronicsId(1, 2, 0, 3);
    eid.setHTR(18, 8, 1);
    myEMap[eid]=DetId(0x54000054);//PZDC EM4

    eid = HcalElectronicsId(2, 2, 0, 3);
    eid.setHTR(18, 8, 1);
    myEMap[eid]=DetId(0x54000055);//PZDC EM5

    eid = HcalElectronicsId(0, 3, 0, 3);
    eid.setHTR(18, 8, 1);
    myEMap[eid]=DetId(0x54000062);//PZDC HAD2

    eid = HcalElectronicsId(1, 3, 0, 3);
    eid.setHTR(18, 8, 1);
    myEMap[eid]=DetId(0x54000063);//PZDC HAD3

    eid = HcalElectronicsId(2, 3, 0, 3);
    eid.setHTR(18, 8, 1);
    myEMap[eid]=DetId(0x54000064);//PZDC HAD4

    //NZDC
    eid = HcalElectronicsId(0,1,1,3);
    eid.setHTR(18, 8, 0);
    myEMap[eid]=DetId(0x54000011);//NZDC EM1

    eid = HcalElectronicsId(1,1,1,3);
    eid.setHTR(18, 8, 0);
    myEMap[eid]=DetId(0x54000012);//NZDC EM2

    eid = HcalElectronicsId(2,1,1,3);
    eid.setHTR(18, 8, 0);
    myEMap[eid]=DetId(0x54000013);//NZDC EM3

    eid = HcalElectronicsId(0,2,1,3);
    eid.setHTR(18, 8, 0);
    myEMap[eid]=DetId(0x54000015);//NZDC EM5

    eid = HcalElectronicsId(1,2,1,3);
    eid.setHTR(18, 8, 0);
    myEMap[eid]=DetId(0x54000021);//NZDC HAD1

    eid = HcalElectronicsId(2,2,1,3);
    eid.setHTR(18, 8, 0);
    myEMap[eid]=DetId(0x54000014);//NZDC EM4

    eid = HcalElectronicsId(0,3,1,3);
    eid.setHTR(18, 8, 0);
    myEMap[eid]=DetId(0x54000022);//NZDC HAD2

    eid = HcalElectronicsId(1,3,1,3);
    eid.setHTR(18, 8, 0);
    myEMap[eid]=DetId(0x54000023);//NZDC HAD3

    eid = HcalElectronicsId(2,3,1,3);
    eid.setHTR(18, 8, 0);
    myEMap[eid]=DetId(0x54000024);//NZDC HAD4
    //##################################### end hardcoding ###################################

    ib.cd();

    //quantities for axis
    hcaldqm::quantity::ValueQuantity xAxisShape(hcaldqm::quantity::fTiming_TS);
    hcaldqm::quantity::ValueQuantity yAxisShape(hcaldqm::quantity::ffC_10000);

    hcaldqm::quantity::ValueQuantity xAxisADC(hcaldqm::quantity::fADC_128);

    //book histos per channel
    for (std::map<HcalElectronicsId,DetId>::const_iterator itr=myEMap.begin(); itr!=myEMap.end(); ++itr)
      {
        char histoname[300];
        
        sprintf(histoname,"%d_%d_%d_%d",itr->first.fiberChanId(),itr->first.fiberIndex(),itr->first.spigot(),itr->first.dccid());

        ib.setCurrentFolder("Hcal/ZDCTask/Shape_perChannel");
        _cShape_EChannel[histoname] = ib.bookProfile(histoname,histoname,xAxisShape.nbins(),xAxisShape.min(),xAxisShape.max(),yAxisShape.nbins(),yAxisShape.min(),yAxisShape.max());
        _cShape_EChannel[histoname]->setAxisTitle("Timing",1);
        _cShape_EChannel[histoname]->setAxisTitle("fC QIE8",2);

        ib.setCurrentFolder("Hcal/ZDCTask/ADC_perChannel");
        _cADC_EChannel[histoname] = ib.book1D(histoname,histoname,xAxisADC.nbins(),xAxisADC.min(),xAxisADC.max());
        _cADC_EChannel[histoname]->getRootObject()->SetBit(BIT(hcaldqm::constants::BIT_OFFSET+hcaldqm::quantity::AxisType::fYAxis));
        _cADC_EChannel[histoname]->setAxisTitle("ADC QIE8",1);

        ib.setCurrentFolder("Hcal/ZDCTask/ADC_vs_TS_perChannel");
        _cADC_vs_TS_EChannel[histoname] = ib.book2D(histoname,histoname,xAxisShape.nbins(),xAxisShape.min(),xAxisShape.max(),xAxisADC.nbins(),xAxisADC.min(),xAxisADC.max());
        _cADC_vs_TS_EChannel[histoname]->getRootObject()->SetBit(BIT(hcaldqm::constants::BIT_OFFSET+hcaldqm::quantity::AxisType::fYAxis));
        _cADC_vs_TS_EChannel[histoname]->setAxisTitle("Timing",1);
        _cADC_vs_TS_EChannel[histoname]->setAxisTitle("ADC QIE8",2);
        

      }

    //book global histos
    ib.setCurrentFolder("Hcal/ZDCTask");

    _cShape = ib.bookProfile("Shape","Shape",xAxisShape.nbins(),xAxisShape.min(),xAxisShape.max(),yAxisShape.nbins(),yAxisShape.min(),yAxisShape.max());
    _cShape->setAxisTitle("Timing",1);
    _cShape->setAxisTitle("fC QIE8",2);

    _cADC = ib.book1D("ADC","ADC",xAxisADC.nbins(),xAxisADC.min(),xAxisADC.max());
    _cADC->getRootObject()->SetBit(BIT(hcaldqm::constants::BIT_OFFSET+hcaldqm::quantity::AxisType::fYAxis));
    _cADC->setAxisTitle("ADC QIE8",1);

    _cADC_vs_TS = ib.book2D("ADC_vs_TS","ADC_vs_TS",xAxisShape.nbins(),xAxisShape.min(),xAxisShape.max(),xAxisADC.nbins(),xAxisADC.min(),xAxisADC.max());
    _cADC_vs_TS->getRootObject()->SetBit(BIT(hcaldqm::constants::BIT_OFFSET+hcaldqm::quantity::AxisType::fYAxis));
    _cADC_vs_TS->setAxisTitle("Timing",1);
    _cADC_vs_TS->setAxisTitle("ADC QIE8",2);
        
}


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


    for (uint32_t i=0; i<cqie10->size(); i++)
      {
        ZDCDataFrame frame = static_cast<ZDCDataFrame>((*cqie10)[i]);
        HcalElectronicsId eid = frame.elecId();

        char histoname[300];
            sprintf(histoname,"%d_%d_%d_%d",eid.fiberChanId(),eid.fiberIndex(),eid.spigot(),eid.dccid());
        
        //  compute the signal, ped subracted
        //double q = hcaldqm::utilities::sumQ_v10<ZDCDataFrame>(frame, constants::adc2fC[_ped], 0, frame.size()-1);

        //  iterate thru all TS and fill
        for (int j=0; j<frame.size(); j++)
          {
        _cShape_EChannel[histoname]->Fill(j, frame[j].nominal_fC());
        _cShape->Fill(j, frame[j].nominal_fC());

        _cADC_EChannel[histoname]->Fill(frame[j].adc());
        _cADC->Fill(frame[j].adc());

        _cADC_vs_TS_EChannel[histoname]->Fill(j,frame[j].adc());
        _cADC_vs_TS->Fill(j,frame[j].adc());

        
          }
      }
}


DEFINE_FWK_MODULE(ZDCTask);