Calol2 online dqm 91x

DQM/Integration/python/clients/l1tstage2emulator_dqm_sourceclient-live_cfg.py

@@ -29,8 +29,8 @@
 
 process.load("DQM.Integration.config.environment_cfi")
 
-process.dqmEnv.subSystemFolder = "L1T2016EMU"
-process.dqmSaver.tag = "L1T2016EMU"
+process.dqmEnv.subSystemFolder = "L1TEMU"
+process.dqmSaver.tag = "L1TEMU"
 process.DQMStore.referenceFileName = "/dqmdata/dqm/reference/l1temu_reference.root"
 
 process.dqmEndPath = cms.EndPath(

DQM/L1TMonitor/interface/L1TStage2CaloLayer2.h

@@ -102,6 +102,18 @@ class L1TStage2CaloLayer2 : public DQMEDAnalyzer {
   MonitorElement* stage2CaloLayer2MHTRank_;
   MonitorElement* stage2CaloLayer2MHTPhi_;
   MonitorElement* stage2CaloLayer2HTTRank_;
+  MonitorElement* stage2CaloLayer2METHFRank_;
+  MonitorElement* stage2CaloLayer2METHFPhi_;
+  // MonitorElement* stage2CaloLayer2ETTHFRank_;
+  MonitorElement* stage2CaloLayer2MHTHFRank_;
+  MonitorElement* stage2CaloLayer2MHTHFPhi_;
+  // MonitorElement* stage2CaloLayer2HTTHFRank_;
+  MonitorElement* stage2CaloLayer2MinBiasHFP0_;
+  MonitorElement* stage2CaloLayer2MinBiasHFM0_;
+  MonitorElement* stage2CaloLayer2MinBiasHFP1_;
+  MonitorElement* stage2CaloLayer2MinBiasHFM1_;
+  MonitorElement* stage2CaloLayer2ETTEMRank_;
+  MonitorElement* stage2CaloLayer2TowCount_;
 
   MonitorElement* timingStage2CaloLayer2CenJetBxOcc_;
   MonitorElement* timingStage2CaloLayer2ForJetBxOcc_;

DQM/L1TMonitor/interface/L1TdeStage2CaloLayer2.h

@@ -0,0 +1,369 @@
+#ifndef DQM_L1TMonitor_L1TdeStage2CaloLayer2
+#define DQM_L1TMonitor_L1TdeStage2CaloLayer2
+
+#include "FWCore/Framework/interface/MakerMacros.h"
+
+#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
+#include "DQMServices/Core/interface/MonitorElement.h"
+
+#include "FWCore/Framework/interface/Event.h"
+#include "FWCore/MessageLogger/interface/MessageLogger.h"
+#include "FWCore/ParameterSet/interface/ParameterSet.h"
+
+#include "DataFormats/L1Trigger/interface/EGamma.h"
+#include "DataFormats/L1Trigger/interface/Jet.h"
+#include "DataFormats/L1Trigger/interface/EtSum.h"
+#include "DataFormats/L1Trigger/interface/Tau.h"
+
+/**
+ * Short class description.
+ *
+ * Longer class description...
+ * ... desc continued.
+ */
+class L1TdeStage2CaloLayer2 : public DQMEDAnalyzer {
+
+ public:
+  /**
+   * Class constructor
+   *
+   * Receives the set of parameters, specified by the python configuration file
+   * used to initialise the module as a part of a sequence. The contents of the
+   * set is used to configure the internal state of the objects of this class.
+   * Values from the parameter set are extracted and used to initialise
+   * bxcollections for jet, e/g, tau and sum objects reconstructed by firmware
+   * (data) and emulator. These collections are the basis of the comparisons
+   * performed by this module.
+   *
+   * @param edm::ParamterSet & ps A pointer to the parameter set used
+   */
+  L1TdeStage2CaloLayer2 (const edm::ParameterSet & ps);
+
+
+ protected:
+
+  /**
+   * Method to declare or "book" all histograms that will be part of module
+   *
+   * Histograms that are to be visualised as part of the DQM module should be
+   * registered with the IBooker object any additional configuration such as
+   * title or axis labels and ranges. A good rule of thumb for the amount of
+   * configuration is that it should be possible to understnand the contents of
+   * the histogram using the configuration received from this method since the
+   * plots generated by this module would later be stored into ROOT files for
+   * transfer to the DQM system and it should be possible to ...
+   *
+   * @param DQMStore::IBooker& ibooker Object that handles the creation of plots
+   * @param edm::Run const &           Reference to run object
+   * @param edm::EventSetup const &    Reference to event configuration object
+   *
+   * @return void
+   */
+  virtual void bookHistograms (DQMStore::IBooker&,
+			       const edm::Run&,
+			       const edm::EventSetup&) override;
+
+  /**
+   * Main method where the analysis code resides, executed once for each run
+   *
+   * The main body of the module code is contained in this method. The different
+   * object collections are extracted from the run and are passed to the
+   * respective comparison methods for processing of each object type.
+   *
+   * @param edm::Event const &         Reference to event object
+   * @param edm::EventSetup const &    Reference to event configuration object
+   *
+   * @return void
+   */
+  virtual void analyze (const edm::Event&, const edm::EventSetup&) override;
+
+ private:
+
+  /**
+   * Encapsulates the code required for performing a comparison of
+   * the jets contained in a given event.
+   *
+   * Method is called once per each event with the jet collections associated
+   * with the event being extracted for bx = 0 as the RAW information required
+   * to run the emulator is only available for bx 0. The implementation checks
+   * if the size of collections is the same and when so, compares the jets in
+   * the same positions within the data/emul collections. When a discrepancy is
+   * found, the properties (Et, eta, phi) of the problematic data/emul objects
+   * are stored in dedicated histograms. In cases of differences, a distinction
+   * is made between objects whose energy or position are in disagremeent.When
+   * the length of data/emul collections is different, all of the objects are
+   * "stored" in the histograms dedicated to problematic objects.
+   *
+   * @param edm::Handle<l1t::JetBXCollection>& dataCol Reference to jet
+   *    collection from data
+   * @param edm::Handle<l1t::JetBXCollection>& emulCol Reference to jet
+   *    collection from emulation
+   *
+   * @return bool Flag of whether the agreement was perfect
+   */
+  bool compareJets(const edm::Handle<l1t::JetBxCollection> & dataCol,
+                   const edm::Handle<l1t::JetBxCollection> & emulCol);
+
+  /**
+   * Encapsulates the code required for performing a comparison of
+   * the e/gs contained in a given event.
+   *
+   * Method is called once per each event with the e/g collections associated
+   * with the event being extracted for bx = 0 as the RAW information required
+   * to run the emulator is only available for bx 0. The implementation checks
+   * if the size of collections is the same and when so, compares the e/gs in
+   * the same positions within the data/emul collections. When a discrepancy is
+   * found, the properties (Et, eta, phi) of the problematic data/emul objects
+   * are stored in dedicated histograms. In cases of differences, a distinction
+   * is made between objects whose energy or position are in disagremeent.
+   * Another distinction is made between isolated and non-isolated e/g
+   * candidates and problematic objects are handled accordingly. When the length
+   * of data/emul collections is different, all of the objects are "stored" in
+   * the histograms dedicated to problematic objects.
+   *
+   * @param edm::Handle<l1t::EGammaBXCollection>& dataCol Reference to e/gamma
+   *    collection from data
+   * @param edm::Handle<l1t::EGammaBXCollection>& emulCol Reference to e/gamma
+   *    collection from emulation
+   *
+   * @return bool Flag of whether the agreement was perfect
+   */
+  bool compareEGs(const edm::Handle<l1t::EGammaBxCollection> & dataCol,
+                  const edm::Handle<l1t::EGammaBxCollection> & emulCol);
+
+  /**
+   * Encapsulates the code required for performing a comparison of
+   * the taus contained in a given event.
+   *
+   * Method is called once per each event with the tau collections associated
+   * with the event being extracted for bx = 0 as the RAW information required
+   * to run the emulator is only available for bx 0. The implementation checks
+   * if the size of collections is the same and when so, compares the taus in
+   * the same positions within the data/emul collections. When a discrepancy is
+   * found, the properties (Et, eta, phi) of the problematic data/emul objects
+   * are stored in dedicated histograms. In cases of differences, a distinction
+   * is made between objects whose energy or position are in disagremeent.
+   * Another distinction is made between isolated and non-isolated tau
+   * candidates and problematic objects are handled accordingly. When the length
+   * of data/emul collections is different, all of the objects are "stored" in
+   * the histograms dedicated to problematic objects.
+   *
+   * @param edm::Handle<l1t::TauBXCollection>& dataCol Reference to tau
+   *    collection from data
+   * @param edm::Handle<l1t::TauBXCollection>& emulCol Reference to tau
+   *    collection from emulation
+   *
+   * @return bool Flag of whether the agreement was perfect
+   */
+  bool compareTaus(const edm::Handle<l1t::TauBxCollection> & dataCol,
+                   const edm::Handle<l1t::TauBxCollection> & emulCol);
+
+  /**
+   * Encapsulates the code required for performing a comparison of
+   * the taus contained in a given event.
+   *
+   * Method is called once per each event with the sum collections associated
+   * with the event being extracted for bx = 0 as the RAW information required
+   * to run the emulator is only available for bx 0. The implementation loops
+   * over the collection and depending of the sum type, each variant is compared
+   * independently. If any disagreement is found, the event is marked a bad and
+   * the properties of the sum are stored in the data/emulator problematic
+   * histograms.
+   *
+   * @param edm::Handle<l1t::TauBXCollection>& dataCol Reference to tau
+   *    collection from data
+   * @param edm::Handle<l1t::TauBXCollection>& emulCol Reference to tau
+   *    collection from emulation
+   *
+   * @return bool Flag of whether the agreement was perfect
+   */
+  bool compareSums(const edm::Handle<l1t::EtSumBxCollection> & dataCol,
+                   const edm::Handle<l1t::EtSumBxCollection> & emulCol);
+
+  // Holds the name of directory in DQM where module hostograms will be shown.
+  // Value is taken from python configuration file (passed in class constructor)
+  std::string monitorDir;
+
+  // collections to hold entities reconstructed from data and emulation
+  edm::EDGetTokenT<l1t::JetBxCollection> calol2JetCollectionData;
+  edm::EDGetTokenT<l1t::JetBxCollection> calol2JetCollectionEmul;
+  edm::EDGetTokenT<l1t::EGammaBxCollection> calol2EGammaCollectionData;
+  edm::EDGetTokenT<l1t::EGammaBxCollection> calol2EGammaCollectionEmul;
+  edm::EDGetTokenT<l1t::TauBxCollection> calol2TauCollectionData;
+  edm::EDGetTokenT<l1t::TauBxCollection> calol2TauCollectionEmul;
+  edm::EDGetTokenT<l1t::EtSumBxCollection> calol2EtSumCollectionData;
+  edm::EDGetTokenT<l1t::EtSumBxCollection> calol2EtSumCollectionEmul;
+
+  enum summaryBins {
+    NEVENTS = 1,      // total number of events
+    EVENTGOOD,        // number of good events (100% agreement)
+    NJETS_S,          // total number of jets objects found
+    JETGOOD_S,        // number of jets in agreement (energy and pos)
+    NEGS_S,           // total number of e/g objects found
+    EGGOOD_S,         // number of e/g in agremeent (energy and pos)
+    NTAUS_S,          // total number of tau objects found
+    TAUGOOD_S,        // number of taus in agremenet (energy and pos)
+    NSUMS_S,          // total number of sums
+    SUMGOOD_S         // number of good sums across all events
+  };
+
+  enum problemCauses {
+    NEVENTS_P = 1,    // total number of events
+    JETCOLLSIZE,      // no. events with different data/emul obj. in jet coll.
+    EGCOLLSIZE,       // no. events with different data/emul obj. in eg coll.
+    TAUCOLLSIZE,      // no. events with different data/emul obj. in tau coll.
+    JETMISMATCH,      // no. events failed due to a jet mismatch
+    EGMISMATCH,       // no. events failed due to an e/g mismatch
+    TAUMISMATCH,      // no. events failed due to a tau mismatch
+    SUMMISMATCH       // no. events failed due to a sum mismatch
+  };
+
+  enum jetVars {
+    NJETS = 1,
+    JETGOOD,
+    JETPOSOFF,
+    JETETOFF
+  };
+
+  enum egVars {
+    NEGS = 1,
+    EGGOOD,
+    EGPOSOFF,
+    EGETOFF,
+    NISOEGS,
+    ISOEGGOOD,
+    ISOEGPOSOFF,
+    ISOEGETOFF
+  };
+
+  enum tauVars {
+    NTAUS = 1,
+    TAUGOOD,
+    TAUPOSOFF,
+    TAUETOFF,
+    NISOTAUS,
+    ISOTAUGOOD,
+    ISOTAUPOSOFF,
+    ISOTAUETOFF
+  };
+
+  enum sumVars {
+    NSUMS = 1,
+    SUMGOOD,
+    NETTSUMS,
+    ETTSUMGOOD,
+    NHTTSUMS,
+    HTTSUMGOOD,
+    NMETSUMS,
+    METSUMGOOD,
+    NMHTSUMS,
+    MHTSUMGOOD,
+    NMBHFSUMS,
+    MBHFSUMGOOD,
+    NTOWCOUNTS,
+    TOWCOUNTGOOD
+  };
+
+  // objects to represent individual plots shown in DQM
+  MonitorElement * agreementSummary;
+  MonitorElement * jetSummary;
+  MonitorElement * tauSummary;
+  MonitorElement * egSummary;
+  MonitorElement * sumSummary;
+  MonitorElement * problemSummary;
+
+  // histograms to store the properties of mismatched jets
+  MonitorElement * jetEtData;
+  MonitorElement * jetEtaData;
+  MonitorElement * jetPhiData;
+  MonitorElement * jetEtEmul;
+  MonitorElement * jetEtaEmul;
+  MonitorElement * jetPhiEmul;
+
+  // histograms to store the properties of mismatched non-isolated e/g
+  MonitorElement * egEtData;
+  MonitorElement * egEtaData;
+  MonitorElement * egPhiData;
+  MonitorElement * egEtEmul;
+  MonitorElement * egEtaEmul;
+  MonitorElement * egPhiEmul;
+
+  // histograms to store the properties of mismatched isolated e/g
+  MonitorElement * isoEgEtData;
+  MonitorElement * isoEgEtaData;
+  MonitorElement * isoEgPhiData;
+  MonitorElement * isoEgEtEmul;
+  MonitorElement * isoEgEtaEmul;
+  MonitorElement * isoEgPhiEmul;
+
+  // histograms to store the properties of mismatched non-isolated taus
+  MonitorElement * tauEtData;
+  MonitorElement * tauEtaData;
+  MonitorElement * tauPhiData;
+  MonitorElement * tauEtEmul;
+  MonitorElement * tauEtaEmul;
+  MonitorElement * tauPhiEmul;
+
+  // histograms to store the properties of mismatched isolated taus
+  MonitorElement * isoTauEtData;
+  MonitorElement * isoTauEtaData;
+  MonitorElement * isoTauPhiData;
+  MonitorElement * isoTauEtEmul;
+  MonitorElement * isoTauEtaEmul;
+  MonitorElement * isoTauPhiEmul;
+
+  // histograms for mismatched ett sums
+  MonitorElement * ettData;
+  MonitorElement * ettEmul;
+  MonitorElement * ettHFData;
+  MonitorElement * ettHFEmul;
+  MonitorElement * ettEmData;
+  MonitorElement * ettEmEmul;
+
+  // mismatched htt sums
+  MonitorElement * httData;
+  MonitorElement * httEmul;
+  MonitorElement * httHFData;
+  MonitorElement * httHFEmul;
+
+  // mismatched met sums
+  MonitorElement * metEtData;
+  MonitorElement * metEtEmul;
+  MonitorElement * metPhiData;
+  MonitorElement * metPhiEmul;
+  MonitorElement * metHFEtData;
+  MonitorElement * metHFEtEmul;
+  MonitorElement * metHFPhiData;
+  MonitorElement * metHFPhiEmul;
+
+  // mismatched mht sums
+  MonitorElement * mhtEtData;
+  MonitorElement * mhtEtEmul;
+  MonitorElement * mhtPhiData;
+  MonitorElement * mhtPhiEmul;
+  MonitorElement * mhtHFEtData;
+  MonitorElement * mhtHFEtEmul;
+  MonitorElement * mhtHFPhiData;
+  MonitorElement * mhtHFPhiEmul;
+
+  // mismatched min bias sums
+  MonitorElement * mbhfp0Data;
+  MonitorElement * mbhfp0Emul;
+  MonitorElement * mbhfm0Data;
+  MonitorElement * mbhfm0Emul;
+  MonitorElement * mbhfp1Data;
+  MonitorElement * mbhfp1Emul;
+  MonitorElement * mbhfm1Data;
+  MonitorElement * mbhfm1Emul;
+
+  // mismatched towercount sum
+  MonitorElement * towCountData;
+  MonitorElement * towCountEmul;
+
+  bool verbose;
+
+  // use only bx = 0 since it only contains RAW data (needed for emulator)
+  const unsigned int currBx = 0;
+};
+
+#endif