DQM/CSCMonitorModule/plugins/CSCDQM_EventProcessor_processEvent.cc - master - githubmirror/cmssw

DQM/CSCMonitorModule/plugins/CSCDQM_EventProcessor_processEvent.cc
350 строк13 KB
Cms Build
Clang-Format
10 май 2019, 18:01
10 май 2019, 18:01ee68c24
/*
 * =====================================================================================
 *
 *       Filename:  EventProcessor_processEvent.cc
 *
 *    Description:  EventProcessor Object Event entry methods
 *
 *        Version:  1.0
 *        Created:  10/03/2008 10:47:11 AM
 *       Revision:  none
 *       Compiler:  gcc
 *
 *         Author:  Valdas Rapsevicius (VR), Valdas.Rapsevicius@cern.ch
 *        Company:  CERN, CH
 *
 * =====================================================================================
 */

#include "CSCDQM_EventProcessor.h"

namespace cscdqm {

  /**
 * @brief  Common Local and Global DQM function to be called before processing Event
 */
  void EventProcessor::preProcessEvent() {
    config->incNEvents();
    EmuEventDisplayWasReset = false;
  }

#ifdef DQMLOCAL

  /**
 * @brief  Process event (Local DQM)
 * @param  data Event Data buffer
 * @param  dataSize Event Data buffer size
 * @param  errorStat Error status received by reading DAQ buffer
 * @param  nodeNumber DAQ node number
 * @return
 */
  void EventProcessor::processEvent(const char* data,
                                    const int32_t dataSize,
                                    const uint32_t errorStat,
                                    const int32_t nodeNumber) {
    preProcessEvent();

    config->incNEventsCSC();

    MonitorObject* me = 0;
    if (getEMUHisto(h::EMU_ALL_READOUT_ERRORS, me)) {
      if (errorStat != 0) {
        me->Fill(nodeNumber, 1);
        for (unsigned int i = 0; i < 16; i++) {
          if ((errorStat >> i) & 0x1) {
            me->Fill(nodeNumber, i + 2);
          }
        }
      } else {
        me->Fill(nodeNumber, 0);
      }
    }

    bool eventDenied = false;
    if (((uint32_t)errorStat & config->getDDU_CHECK_MASK()) > 0) {
      eventDenied = true;
    }

    const uint16_t* tmp = reinterpret_cast<const uint16_t*>(data);
    const uint16_t tmpSize = dataSize / sizeof(short);

    /** CSC DCC Examiner object */
    CSCDCCExaminer binChecker;
    binChecker.crcALCT(config->getBINCHECKER_CRC_ALCT());
    binChecker.crcTMB(config->getBINCHECKER_CRC_CLCT());
    binChecker.crcCFEB(config->getBINCHECKER_CRC_CFEB());
    binChecker.modeDDU(config->getBINCHECKER_MODE_DDU());

    if (config->getBINCHECKER_OUTPUT()) {
      // binChecker.output1().show();
      // binChecker.output2().show();
    } else {
      // binChecker.output1().hide();
      // binChecker.output2().hide();
    }

    binChecker.setMask(config->getBINCHECK_MASK());

    if (binChecker.check(tmp, tmpSize) < 0) {
      /** No ddu trailer found - force checker to summarize errors by adding artificial trailer */
      const uint16_t dduTrailer[4] = {0x8000, 0x8000, 0xFFFF, 0x8000};
      const uint16_t* tmp = dduTrailer;
      binChecker.check(tmp, uint32_t(4));
    }

    CSCDCCFormatStatusDigi formatStatusDigi(nodeNumber,
                                            config->getBINCHECK_MASK(),
                                            binChecker.getMask(),
                                            binChecker.errors(),
                                            binChecker.errorsDetailedDDU(),
                                            binChecker.errorsDetailed(),
                                            binChecker.payloadDetailed(),
                                            binChecker.statusDetailed());

    if (processExaminer(binChecker, formatStatusDigi)) {
      config->incNEventsGood();

      if (config->getPROCESS_DDU()) {
        CSCDDUEventData dduData((short unsigned int*)tmp, &binChecker);
        processDDU(dduData, binChecker);
      }
    }
  }

#endif

#ifdef DQMGLOBAL

  /**
 * @brief  Process event (Global DQM)
 * @param  e Event object
 * @param  inputTag Tag to search Event Data in
 * @param  standby HW standby status indicator
 * @return
 */
  void EventProcessor::processEvent(const edm::Event& e, const edm::InputTag& inputTag) {
    preProcessEvent();

    edm::Handle<FEDRawDataCollection> rawdata;
    if (!e.getByToken(frdtoken, rawdata)) {
      LOG_WARN << "No product: " << inputTag << " in stream";
      return;
    }

    bCSCEventCounted = false;
    size_t eventSize = 0;  /// Total Event Size in bytes
    cntDMBs = 0;
    cntCFEBs = 0;
    cntALCTs = 0;
    cntTMBs = 0;
    int nDDUs_out_of_sync = 0;
    int nDDUs_with_CSC_data_out_of_sync = 0;
    bool fGlobal_DCC_DDU_L1A_mismatch = false;
    bool fGlobal_DCC_DDU_L1A_mismatch_with_CSC_data = false;
    MonitorObject* mo = nullptr;

    // For new CSC readout layout, which wont include DCCs need to loop over DDU FED IDs. DCC IDs are included for backward compatibility with old data
    std::vector<unsigned int> cscFEDids;

    for (unsigned int id = FEDNumbering::MINCSCFEDID; id <= FEDNumbering::MAXCSCFEDID; ++id)  // loop over DCCs
    {
      cscFEDids.push_back(id);
    }

    for (unsigned int id = FEDNumbering::MINCSCDDUFEDID; id <= FEDNumbering::MAXCSCDDUFEDID; ++id)  // loop over DDUs
    {
      cscFEDids.push_back(id);
    }

    // run through the FED's
    for (unsigned int i = 0; i < cscFEDids.size(); i++)  // loop over all CSC FEDs (DCCs and DDUs)
    {
      unsigned int id = cscFEDids[i];
      bool isDDU_FED = ((id >= FEDNumbering::MINCSCDDUFEDID) && (id <= FEDNumbering::MAXCSCDDUFEDID)) ? true : false;
      // Take a reference to this FED's data and
      // construct the DCC data object
      const FEDRawData& fedData = rawdata->FEDData(id);

      //if fed has data then unpack it
      if (fedData.size() >= 32) {
        eventSize += fedData.size();

        // Count in CSC Event!
        if (!bCSCEventCounted) {
          config->incNEventsCSC();
          bCSCEventCounted = true;
        }

        if (!isDDU_FED) {
          // Filling in FED Entries histogram
          if (getEMUHisto(h::EMU_FED_ENTRIES, mo))
            mo->Fill(id);

          if (getEMUHisto(h::EMU_FED_BUFFER_SIZE, mo))
            mo->Fill(id, log10((double)fedData.size()));

          if (getFEDHisto(h::FED_BUFFER_SIZE, id, mo))
            mo->Fill(fedData.size());
        }

        const uint16_t* data = (uint16_t*)fedData.data();
        const size_t dataSize = fedData.size() / 2;
        const short unsigned int* udata = (short unsigned int*)fedData.data();

        /** CSC DCC Examiner object */
        CSCDCCExaminer binChecker;
        binChecker.crcALCT(config->getBINCHECKER_CRC_ALCT());
        binChecker.crcTMB(config->getBINCHECKER_CRC_CLCT());
        binChecker.crcCFEB(config->getBINCHECKER_CRC_CFEB());
        binChecker.modeDDU(config->getBINCHECKER_MODE_DDU());

        /// Check if we have DCC or DDU FED, by checking FED ID. Then set examiner to use DCC or DDU mode
        if (isDDU_FED) {
          binChecker.modeDDU(true);
        }

        if (config->getBINCHECKER_OUTPUT()) {
          // binChecker.output1().show();
          // binChecker.output2().show();
        } else {
          // binChecker.output1().hide();
          // binChecker.output2().hide();
        }

        binChecker.setMask(config->getBINCHECK_MASK());

        if (isDDU_FED) {
          if (binChecker.check(data, dataSize) >= 0) {
            CSCDCCFormatStatusDigi formatStatusDigi(id,
                                                    config->getBINCHECK_MASK(),
                                                    binChecker.getMask(),
                                                    binChecker.errors(),
                                                    binChecker.errorsDetailedDDU(),
                                                    binChecker.errorsDetailed(),
                                                    binChecker.payloadDetailed(),
                                                    binChecker.statusDetailed());

            if (processExaminer(binChecker, formatStatusDigi)) {
              config->incNEventsGood();

              if (config->getPROCESS_DDU()) {
                CSCDDUEventData dduData(const_cast<short unsigned int*>(udata), &binChecker);
                processDDU(dduData, binChecker);
              }
            }
          }
        }

        else {
          if (binChecker.check(data, dataSize) < 0) {
            if (getEMUHisto(h::EMU_FED_FATAL, mo))
              mo->Fill(id);
          } else {
            CSCDCCFormatStatusDigi formatStatusDigi(id,
                                                    config->getBINCHECK_MASK(),
                                                    binChecker.getMask(),
                                                    binChecker.errors(),
                                                    binChecker.errorsDetailedDDU(),
                                                    binChecker.errorsDetailed(),
                                                    binChecker.payloadDetailed(),
                                                    binChecker.statusDetailed());

            /// Fill summary FED/DCC Format Errors histos
            if (getEMUHisto(h::EMU_FED_FORMAT_ERRORS, mo)) {
              uint32_t errs = binChecker.errors();

              if (errs != 0) {
                for (int i = 0; i < 29; i++) {
                  if ((errs >> i) & 0x1) {
                    mo->Fill(id, i + 1);
                  }
                }
              } else {
                mo->Fill(id, 0);
              }
            }

            if (processExaminer(binChecker, formatStatusDigi)) {
              config->incNEventsGood();

              if (binChecker.warnings() != 0) {
                if (getEMUHisto(h::EMU_FED_NONFATAL, mo))
                  mo->Fill(id);
              }

              if (config->getPROCESS_DDU()) {
                CSCDCCEventData dccData(const_cast<short unsigned int*>(udata), &binChecker);
                const std::vector<CSCDDUEventData>& dduData = dccData.dduData();

                /// Check for DCC and DDUs L1A counter mismatch and fill corresponding histogram
                bool fDCC_DDU_L1A_mismatch = false;
                bool fDCC_DDU_L1A_mismatch_with_CSC_data = false;
                int DCC_L1A = dccData.dccHeader().getCDFEventNumber();
                for (int ddu = 0; ddu < (int)dduData.size(); ddu++) {
                  if (DCC_L1A != dduData[ddu].header().lvl1num()) {
                    fDCC_DDU_L1A_mismatch = true;
                    fGlobal_DCC_DDU_L1A_mismatch = true;
                    nDDUs_out_of_sync++;

                    /// Check if DDU potentially has CSC data
                    if (dduData[ddu].sizeInWords() > 24) {
                      fDCC_DDU_L1A_mismatch_with_CSC_data = true;
                      fGlobal_DCC_DDU_L1A_mismatch_with_CSC_data = true;
                      nDDUs_with_CSC_data_out_of_sync++;
                      /*
                                                        std::cout <<  "FED" << id << " L1A:" << DCC_L1A <<  ", ";
                                                        std::cout << "DDU" << (dduData[ddu].header().source_id() & 0xFF )
                                                        	<< " L1A:" << dduData[ddu].header().lvl1num() << " size:" << (dduData[ddu].sizeInWords()*2);
                                                        std::cout << " - L1A mismatch";
                                                        std::cout << std::endl;
                                      */
                    }
                  }
                  processDDU(dduData[ddu], binChecker);
                }

                if (fDCC_DDU_L1A_mismatch && getEMUHisto(h::EMU_FED_DDU_L1A_MISMATCH, mo))
                  mo->Fill(id);
                if (fDCC_DDU_L1A_mismatch_with_CSC_data && getEMUHisto(h::EMU_FED_DDU_L1A_MISMATCH_WITH_CSC_DATA, mo))
                  mo->Fill(id);
              }

            } else {
              if (getEMUHisto(h::EMU_FED_FORMAT_FATAL, mo))
                mo->Fill(id);
            }
          }
        }
      }
    }

    if (fGlobal_DCC_DDU_L1A_mismatch && getEMUHisto(h::EMU_FED_DDU_L1A_MISMATCH_CNT, mo))
      mo->Fill(nDDUs_out_of_sync);
    if (fGlobal_DCC_DDU_L1A_mismatch && getEMUHisto(h::EMU_FED_DDU_L1A_MISMATCH_WITH_CSC_DATA_CNT, mo))
      mo->Fill(nDDUs_with_CSC_data_out_of_sync);

    if (getEMUHisto(h::EMU_FED_STATS, mo)) {
      mo->Fill(0);
      if (fGlobal_DCC_DDU_L1A_mismatch)
        mo->Fill(1);
      if (fGlobal_DCC_DDU_L1A_mismatch_with_CSC_data)
        mo->Fill(2);
    }

    if (getEMUHisto(h::EMU_FED_EVENT_SIZE, mo))
      mo->Fill(eventSize / 1024.);  /// CSC Event Size in KBytes
    if (getEMUHisto(h::EMU_FED_TOTAL_CSC_NUMBER, mo))
      mo->Fill(cntDMBs);  /// Total Number of CSC/DMBs in event (DDU Header DAV)
    if (getEMUHisto(h::EMU_FED_TOTAL_CFEB_NUMBER, mo))
      mo->Fill(cntCFEBs);  /// Total Number of CFEBs in event (DMB Header DAV)
    if (getEMUHisto(h::EMU_FED_TOTAL_ALCT_NUMBER, mo))
      mo->Fill(cntALCTs);  /// Total Number of ALCTs in event (DMB Header DAV)
    if (getEMUHisto(h::EMU_FED_TOTAL_TMB_NUMBER, mo))
      mo->Fill(cntTMBs);  /// Total Number of TMBs in event (DMB Header DAV)
  }

#endif

#undef ECHO_FUNCTION

}  // namespace cscdqm