L1TdeRCT.h

Go to the documentation of this file.

#ifndef L1TdeRCT_H
#define L1TdeRCT_H

// system include files
#include <memory>
#include <unistd.h>

#include <iostream>
#include <fstream>
#include <vector>
#include <bitset>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

// DQM
#include "DQMServices/Core/interface/DQMStore.h"

// GCT and RCT data formats
#include "DataFormats/L1CaloTrigger/interface/L1CaloCollections.h"
#include "CondFormats/L1TObjects/interface/L1GtTriggerMenu.h"
#include "CondFormats/DataRecord/interface/L1GtTriggerMenuRcd.h"
#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutSetupFwd.h"
#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutSetup.h"
#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutRecord.h"
// TPGs
#include "DataFormats/EcalDigi/interface/EcalDigiCollections.h"
#include "DataFormats/HcalDigi/interface/HcalDigiCollections.h"

#include "DQMServices/Core/interface/DQMOneEDAnalyzer.h"

class RunInfoRcd;
class RunInfo;
// Trigger Headers
//
// class declaration
//
namespace l1tderct {
  struct Empty {};
}  // namespace l1tderct

class L1TdeRCT : public DQMOneEDAnalyzer<edm::LuminosityBlockCache<l1tderct::Empty>> {
public:
  // Constructor
  L1TdeRCT(const edm::ParameterSet& ps);

  // Destructor
  ~L1TdeRCT() override;

protected:
  // Analyze
  void analyze(const edm::Event& e, const edm::EventSetup& c) override;

  //For FED vector monitoring
  void bookHistograms(DQMStore::IBooker& ibooker, const edm::Run&, const edm::EventSetup&) override;
  std::shared_ptr<l1tderct::Empty> globalBeginLuminosityBlock(const edm::LuminosityBlock&,
                                                              const edm::EventSetup&) const override;
  void globalEndLuminosityBlock(const edm::LuminosityBlock&, const edm::EventSetup&) final {}
  void readFEDVector(MonitorElement*, const edm::EventSetup&, const bool isLumitransition = true) const;

private:
  // ----------member data ---------------------------

  // begin GT decision information
  MonitorElement* triggerAlgoNumbers_;

  // trigger type information
  MonitorElement* triggerType_;

  // begin region information
  MonitorElement* rctRegDataOcc1D_;
  MonitorElement* rctRegEmulOcc1D_;
  MonitorElement* rctRegMatchedOcc1D_;
  MonitorElement* rctRegUnmatchedDataOcc1D_;
  MonitorElement* rctRegUnmatchedEmulOcc1D_;
  MonitorElement* rctRegSpEffOcc1D_;
  MonitorElement* rctRegSpIneffOcc1D_;

  MonitorElement* rctRegEff1D_;
  MonitorElement* rctRegIneff1D_;
  MonitorElement* rctRegOvereff1D_;
  MonitorElement* rctRegSpEff1D_;
  MonitorElement* rctRegSpIneff1D_;

  MonitorElement* rctRegDataOcc2D_;
  MonitorElement* rctRegEmulOcc2D_;
  MonitorElement* rctRegMatchedOcc2D_;
  MonitorElement* rctRegUnmatchedDataOcc2D_;
  MonitorElement* rctRegUnmatchedEmulOcc2D_;
  //  MonitorElement *rctRegDeltaEt2D_;
  MonitorElement* rctRegSpEffOcc2D_;
  MonitorElement* rctRegSpIneffOcc2D_;

  MonitorElement* rctRegEff2D_;
  MonitorElement* rctRegIneff2D_;
  MonitorElement* rctRegOvereff2D_;
  MonitorElement* rctRegSpEff2D_;
  MonitorElement* rctRegSpIneff2D_;

  MonitorElement* rctRegBitOn_;
  MonitorElement* rctRegBitOff_;
  MonitorElement* rctRegBitDiff_;

  // end region information

  // begin bit information
  MonitorElement* rctBitEmulOverFlow2D_;
  MonitorElement* rctBitDataOverFlow2D_;
  MonitorElement* rctBitMatchedOverFlow2D_;
  MonitorElement* rctBitUnmatchedEmulOverFlow2D_;
  MonitorElement* rctBitUnmatchedDataOverFlow2D_;
  MonitorElement* rctBitOverFlowEff2D_;
  MonitorElement* rctBitOverFlowIneff2D_;
  MonitorElement* rctBitOverFlowOvereff2D_;
  MonitorElement* rctBitEmulTauVeto2D_;
  MonitorElement* rctBitDataTauVeto2D_;
  MonitorElement* rctBitMatchedTauVeto2D_;
  MonitorElement* rctBitUnmatchedEmulTauVeto2D_;
  MonitorElement* rctBitUnmatchedDataTauVeto2D_;
  MonitorElement* rctBitTauVetoEff2D_;
  MonitorElement* rctBitTauVetoIneff2D_;
  MonitorElement* rctBitTauVetoOvereff2D_;
  MonitorElement* rctBitEmulMip2D_;
  MonitorElement* rctBitDataMip2D_;
  MonitorElement* rctBitMatchedMip2D_;
  MonitorElement* rctBitUnmatchedEmulMip2D_;
  MonitorElement* rctBitUnmatchedDataMip2D_;
  MonitorElement* rctBitMipEff2D_;
  MonitorElement* rctBitMipIneff2D_;
  MonitorElement* rctBitMipOvereff2D_;
  MonitorElement* rctBitEmulQuiet2D_;
  MonitorElement* rctBitDataQuiet2D_;
  MonitorElement* rctBitMatchedQuiet2D_;
  MonitorElement* rctBitUnmatchedEmulQuiet2D_;
  MonitorElement* rctBitUnmatchedDataQuiet2D_;
  // QUIETBIT: To add quiet bit information, uncomment following 3 lines:
  // MonitorElement *rctBitQuietEff2D_;
  // MonitorElement *rctBitQuietIneff2D_;
  // MonitorElement *rctBitQuietOvereff2D_;
  MonitorElement* rctBitEmulHfPlusTau2D_;
  MonitorElement* rctBitDataHfPlusTau2D_;
  MonitorElement* rctBitMatchedHfPlusTau2D_;
  MonitorElement* rctBitUnmatchedEmulHfPlusTau2D_;
  MonitorElement* rctBitUnmatchedDataHfPlusTau2D_;
  MonitorElement* rctBitHfPlusTauEff2D_;
  MonitorElement* rctBitHfPlusTauIneff2D_;
  MonitorElement* rctBitHfPlusTauOvereff2D_;

  // end bit information

  MonitorElement* rctInputTPGEcalOcc_;
  MonitorElement* rctInputTPGEcalOccNoCut_;
  MonitorElement* rctInputTPGEcalRank_;
  MonitorElement* rctInputTPGHcalOcc_;
  MonitorElement* rctInputTPGHcalRank_;
  MonitorElement* rctInputTPGHcalSample_;

  MonitorElement* rctIsoEmDataOcc_;
  MonitorElement* rctIsoEmEmulOcc_;
  MonitorElement* rctIsoEmEff1Occ_;
  MonitorElement* rctIsoEmEff2Occ_;
  MonitorElement* rctIsoEmIneff2Occ_;
  MonitorElement* rctIsoEmIneffOcc_;
  MonitorElement* rctIsoEmOvereffOcc_;
  MonitorElement* rctIsoEmEff1_;
  MonitorElement* rctIsoEmEff2_;
  MonitorElement* rctIsoEmIneff2_;
  MonitorElement* rctIsoEmIneff_;
  MonitorElement* rctIsoEmOvereff_;

  MonitorElement* rctIsoEmDataOcc1D_;
  MonitorElement* rctIsoEmEmulOcc1D_;
  MonitorElement* rctIsoEmEff1Occ1D_;
  MonitorElement* rctIsoEmEff2Occ1D_;
  MonitorElement* rctIsoEmIneff2Occ1D_;
  MonitorElement* rctIsoEmIneffOcc1D_;
  MonitorElement* rctIsoEmOvereffOcc1D_;
  MonitorElement* rctIsoEmEff1oneD_;
  MonitorElement* rctIsoEmEff2oneD_;
  MonitorElement* rctIsoEmIneff2oneD_;
  MonitorElement* rctIsoEmIneff1D_;
  MonitorElement* rctIsoEmOvereff1D_;

  MonitorElement* rctIsoEmBitOn_;
  MonitorElement* rctIsoEmBitOff_;
  MonitorElement* rctIsoEmBitDiff_;

  MonitorElement* rctNisoEmDataOcc_;
  MonitorElement* rctNisoEmEmulOcc_;
  MonitorElement* rctNisoEmEff1Occ_;
  MonitorElement* rctNisoEmEff2Occ_;
  MonitorElement* rctNisoEmIneff2Occ_;
  MonitorElement* rctNisoEmIneffOcc_;
  MonitorElement* rctNisoEmOvereffOcc_;
  MonitorElement* rctNisoEmEff1_;
  MonitorElement* rctNisoEmEff2_;
  MonitorElement* rctNisoEmIneff2_;
  MonitorElement* rctNisoEmIneff_;
  MonitorElement* rctNisoEmOvereff_;

  MonitorElement* rctNisoEmDataOcc1D_;
  MonitorElement* rctNisoEmEmulOcc1D_;
  MonitorElement* rctNisoEmEff1Occ1D_;
  MonitorElement* rctNisoEmEff2Occ1D_;
  MonitorElement* rctNisoEmIneff2Occ1D_;
  MonitorElement* rctNisoEmIneffOcc1D_;
  MonitorElement* rctNisoEmOvereffOcc1D_;
  MonitorElement* rctNisoEmEff1oneD_;
  MonitorElement* rctNisoEmEff2oneD_;
  MonitorElement* rctNisoEmIneff2oneD_;
  MonitorElement* rctNisoEmIneff1D_;
  MonitorElement* rctNisoEmOvereff1D_;

  MonitorElement* rctNIsoEmBitOn_;
  MonitorElement* rctNIsoEmBitOff_;
  MonitorElement* rctNIsoEmBitDiff_;

  MonitorElement* rctIsoEffChannel_[396];
  MonitorElement* rctIsoIneffChannel_[396];
  MonitorElement* rctIsoOvereffChannel_[396];

  MonitorElement* rctNisoEffChannel_[396];
  MonitorElement* rctNisoIneffChannel_[396];
  MonitorElement* rctNisoOvereffChannel_[396];

  // begin region channel information
  MonitorElement* rctRegEffChannel_[396];
  MonitorElement* rctRegIneffChannel_[396];
  MonitorElement* rctRegOvereffChannel_[396];

  //efficiency
  MonitorElement* trigEffThresh_;
  MonitorElement* trigEffThreshOcc_;
  MonitorElement* trigEffTriggThreshOcc_;
  MonitorElement* trigEff_[396];
  MonitorElement* trigEffOcc_[396];
  MonitorElement* trigEffTriggOcc_[396];

  // end region channel information

  //begin fed vector information
  static const int crateFED[108];
  MonitorElement* fedVectorMonitorRUN_;
  MonitorElement* fedVectorMonitorLS_;

  int nev_;                 // Number of events processed
  std::string histFolder_;  // base dqm folder
  bool verbose_;
  bool singlechannelhistos_;
  bool perLSsaving_;  //to avoid nanoDQMIO crashing, driven by  DQMServices/Core/python/DQMStore_cfi.py

  edm::EDGetTokenT<L1CaloRegionCollection> rctSourceEmul_rgnEmul_;
  edm::EDGetTokenT<L1CaloEmCollection> rctSourceEmul_emEmul_;
  edm::EDGetTokenT<L1CaloRegionCollection> rctSourceData_rgnData_;
  edm::EDGetTokenT<L1CaloEmCollection> rctSourceData_emData_;
  edm::EDGetTokenT<L1CaloRegionCollection> gctSourceData_rgnData_;
  edm::EDGetTokenT<L1CaloEmCollection> gctSourceData_emData_;
  edm::EDGetTokenT<EcalTrigPrimDigiCollection> ecalTPGData_;
  edm::EDGetTokenT<HcalTrigPrimDigiCollection> hcalTPGData_;
  edm::EDGetTokenT<L1GlobalTriggerReadoutRecord> gtDigisLabel_;
  edm::ESGetToken<RunInfo, RunInfoRcd> runInfoToken_;
  edm::ESGetToken<RunInfo, RunInfoRcd> runInfolumiToken_;
  std::string gtEGAlgoName_;  // name of algo to determine EG trigger threshold
  int doubleThreshold_;       // value of ET at which to make 2-D eff plot

  int filterTriggerType_;
  int selectBX_;

  std::string dataInputTagName_;

  int trigCount, notrigCount;

protected:
  void DivideME1D(MonitorElement* numerator, MonitorElement* denominator, MonitorElement* result);
  void DivideME2D(MonitorElement* numerator, MonitorElement* denominator, MonitorElement* result);
};

#endif