Scouting Offline DQM for 2025 data-taking (EGM+MUO+autoDQM)

DQMOffline/Configuration/python/DQMOffline_SecondStep_cff.py

@@ -84,6 +84,7 @@
 from DQM.TrackingMonitorClient.TrackingClientConfig_Tier0_cff import *
 from DQMOffline.Muon.muonQualityTests_cff import *
 from DQMOffline.EGamma.egammaPostProcessing_cff import *
+from DQMOffline.HLTScouting.HLTScoutingPostProcessing_cff import *
 from DQMOffline.Trigger.DQMOffline_Trigger_Client_cff import *
 from DQMOffline.Trigger.DQMOffline_HLT_Client_cff import *
 from DQMOffline.RecoB.dqmCollector_cff import *
@@ -267,6 +268,9 @@
 
 DQMHarvestBTag = cms.Sequence( bTagCollectorSequenceDATA )
 
+# HLT Scouting trigger sequence
+DQMHarvestHLTScouting = cms.Sequence( hltScoutingPostProcessing )
+
 from PhysicsTools.NanoAOD.nanoDQM_cff import *
 from Validation.RecoParticleFlow.DQMForPF_MiniAOD_cff import *
 from DQMOffline.RecoB.bTagMiniDQM_cff import *

DQMOffline/Configuration/python/DQMOffline_cff.py

@@ -4,6 +4,7 @@
 from DQMServices.Components.DQMProvInfo_cfi import *
 from DQMServices.Components.DQMFastTimerService_cff import *
 
+from DQMOffline.HLTScouting.HLTScoutingDqmOffline_cff import *
 from DQMOffline.L1Trigger.L1TriggerDqmOffline_cff import *
 from DQMOffline.Ecal.ecal_dqm_source_offline_cff import *
 from DQM.EcalPreshowerMonitorModule.es_dqm_source_offline_cff import *
@@ -28,6 +29,9 @@
 dqmProvInfo.dcsRecord = cms.untracked.InputTag("onlineMetaDataDigis")
 DQMOfflineDCS = cms.Sequence( dqmProvInfo )
 
+# HLT Scouting trigger sequence
+DQMOfflineScouting = cms.Sequence( hltScoutingDqmOffline ) 
+
 # L1 trigger sequences
 DQMOfflineL1T = cms.Sequence( l1TriggerDqmOffline ) # L1 emulator is run within this sequence for real data
 

DQMOffline/Configuration/python/autoDQM.py

@@ -6,9 +6,13 @@
                         'PostDQMOffline',
                         '@dcs+@DQMMessageLogger+@hlt+@beam+@fed+dqmFastTimerServiceClient'],
 
-	   'common': ['@dcs+@DQMMessageLogger+@stripCommon+@pixel+@tracking+@hlt+@beam+@castor+@physics',
-                        'PostDQMOffline',
-                        '@dcs+@DQMMessageLogger+@stripCommon+@pixel+@tracking+@hlt+@beam+@fed+dqmFastTimerServiceClient'],
+	    'common': ['@dcs+@DQMMessageLogger+@stripCommon+@pixel+@tracking+@hlt+@beam+@castor+@physics',
+                      'PostDQMOffline',
+                      '@dcs+@DQMMessageLogger+@stripCommon+@pixel+@tracking+@hlt+@beam+@fed+dqmFastTimerServiceClient'],
+
+            'commonWithScouting': ['@dcs+@DQMMessageLogger+@stripCommon+@pixel+@tracking+@hlt+@beam+@castor+@physics+@hltScouting',
+                                   'PostDQMOffline',
+                                   '@dcs+@DQMMessageLogger+@stripCommon+@pixel+@tracking+@hlt+@beam+@fed+dqmFastTimerServiceClient'],
 
             'commonFakeHLT': ['@dcs+@DQMMessageLogger+@stripCommon+@pixel+@tracking+@beam+@castor+@physics',
                         'PostDQMOffline',
@@ -141,7 +145,7 @@
 	    'physics': ['DQMOfflinePhysics',
 			'PostDQMOffline',
 			'DQMNone'],
-
+            
             'heavyFlavor': ['DQMOfflineHeavyFlavor',
                             'PostDQMOffline',
                             'DQMNone'],
@@ -219,6 +223,10 @@
                             'PostDQMOffline',
                             'dqmHarvesting'],
 
+	    'hltScouting': ['DQMOfflineScouting',
+			    'PostDQMOffline',
+			    'DQMHarvestHLTScouting'],
+
             'standardDQMExpress': ['DQMOfflineExpress',
                                    'PostDQMOffline',
                                    'dqmHarvestingExpress'],

DQMOffline/HLTScouting/BuildFile.xml

@@ -0,0 +1,47 @@
+<use name="CalibFormats/SiStripObjects"/>
+<use name="CalibTracker/Records"/>
+<use name="CommonTools/TriggerUtils"/>
+<use name="CondFormats/DTObjects"/>
+<use name="CondFormats/DataRecord"/>
+<use name="CondFormats/SiStripObjects"/>
+<use name="DQMServices/Core"/>
+<use name="DataFormats/CaloTowers"/>
+<use name="DataFormats/CSCDigi"/>
+<use name="DataFormats/CSCRecHit"/>
+<use name="DataFormats/DTRecHit"/>
+<use name="DataFormats/DetId"/>
+<use name="DataFormats/EcalDetId"/>
+<use name="DataFormats/EgammaCandidates"/>
+<use name="DataFormats/EgammaReco"/>
+<use name="DataFormats/FEDRawData"/>
+<use name="DataFormats/HcalDetId"/>
+<use name="DataFormats/HcalRecHit"/>
+<use name="DataFormats/JetReco"/>
+<use name="DataFormats/METReco"/>
+<use name="DataFormats/Math"/>
+<use name="DataFormats/MuonReco"/>
+<use name="DataFormats/PatCandidates"/>
+<use name="DataFormats/RecoCandidate"/>
+<use name="DataFormats/SiStripCluster"/>
+<use name="DataFormats/SiStripDetId"/>
+<use name="DataFormats/TrackReco"/>
+<use name="DataFormats/TrackerRecHit2D"/>
+<use name="DataFormats/TrackingRecHit"/>
+<use name="EventFilter/CSCRawToDigi"/>
+<use name="FWCore/Framework"/>
+<use name="FWCore/PluginManager"/>
+<use name="FWCore/Utilities"/>
+<use name="Geometry/CaloGeometry"/>
+<use name="Geometry/CSCGeometry"/>
+<use name="Geometry/HcalTowerAlgo"/>
+<use name="HLTrigger/HLTcore"/>
+<use name="HLTriggerOffline/Scouting"/>
+<use name="JetMETCorrections/JetCorrector"/>
+<use name="MagneticField/Engine"/>
+<use name="MagneticField/Records"/>
+<use name="RecoJets/JetAssociationAlgorithms"/>
+<use name="RecoJets/JetProducers"/>
+<use name="RecoMuon/TrackingTools"/>
+<use name="TrackingTools/Records"/>
+<use name="TrackingTools/TransientTrack"/>
+<flags EDM_PLUGIN="1"/>

DQMOffline/HLTScouting/python/HLTScoutingDqmOffline_cff.py

@@ -0,0 +1,21 @@
+'''
+Scouting DQM sequences for offline DQM developed for 2025 pp data-taking
+and used by DQM GUI (DQMOffline/Configuration):
+currently running EGM and MUO monitoring modules.
+'''
+
+import FWCore.ParameterSet.Config as cms
+from DQMServices.Core.DQMEDHarvester import DQMEDHarvester                                
+
+from HLTriggerOffline.Scouting.ScoutingMuonTriggerAnalyzer_cfi import *
+from HLTriggerOffline.Scouting.ScoutingMuonTagProbeAnalyzer_cfi import *
+from HLTriggerOffline.Scouting.ScoutingMuonMonitoring_Client_cff import *
+
+from HLTriggerOffline.Scouting.HLTScoutingEGammaDqmOffline_cff import *
+
+hltScoutingMuonDqmOffline = cms.Sequence(scoutingMonitoringTagProbeMuonNoVtx
+                                         * scoutingMonitoringTagProbeMuonVtx                                                         
+                                         * scoutingMonitoringTriggerMuon                                                              
+) 
+
+hltScoutingDqmOffline = cms.Sequence(hltScoutingMuonDqmOffline + hltScoutingEGammaDqmOffline)

DQMOffline/HLTScouting/python/HLTScoutingPostProcessing_cff.py

@@ -0,0 +1,19 @@
+'''
+Scouting DQM sequences for offline DQM developed for 2025 pp data-taking
+and used by DQM GUI (DQMOffline/Configuration):
+currently running EGM and MUO monitoring modules.
+'''
+
+import FWCore.ParameterSet.Config as cms
+from DQMServices.Core.DQMEDHarvester import DQMEDHarvester
+from DQMOffline.HLTScouting.HLTScoutingDqmOffline_cff import *
+
+from HLTriggerOffline.Scouting.ScoutingMuonMonitoring_Client_cff import *
+from HLTriggerOffline.Scouting.HLTScoutingEGammaPostProcessing_cff import *
+
+hltScoutingMuonPostProcessing = cms.Sequence(muonEfficiencyNoVtx
+                                             * muonEfficiencyVtx  
+                                             * muonTriggerEfficiency
+                                             )
+
+hltScoutingPostProcessing = cms.Sequence(hltScoutingMuonPostProcessing + hltScoutingEGammaPostProcessing)

HLTriggerOffline/Scouting/plugins/BuildFile.xml

@@ -0,0 +1,14 @@
+<use name="FWCore/Framework"/>
+<use name="FWCore/PluginManager"/>
+<use name="FWCore/ParameterSet"/>
+<use name="DQMServices/Core"/>
+<use name="CLHEP"/>
+<use name="DataFormats/PatCandidates"/>
+<use name="DataFormats/Scouting"/>
+<use name="L1Trigger/L1TGlobal"/>
+<use name="DataFormats/L1TGlobal"/>
+<use name="FWCore/Utilities"/>
+<use name="TrackingTools/Records"/>
+<use name="HLTrigger/HLTcore"/>
+<use name="DataFormats/HLTReco"/>
+<flags EDM_PLUGIN="1"/>

HLTriggerOffline/Scouting/plugins/ElectronEfficiencyPlotter.cc

@@ -0,0 +1,101 @@
+#include "ElectronEfficiencyPlotter.h"
+
+// Framework
+#include <cmath>
+#include <string>
+
+#include "DQMServices/Core/interface/DQMStore.h"
+#include "DataFormats/Common/interface/Handle.h"
+#include "FWCore/Framework/interface/ESHandle.h"
+#include "FWCore/Framework/interface/Event.h"
+#include "FWCore/Framework/interface/EventSetup.h"
+#include "FWCore/Framework/interface/MakerMacros.h"
+#include "FWCore/MessageLogger/interface/MessageLogger.h"
+#include "FWCore/ParameterSet/interface/ParameterSet.h"
+#include "TF1.h"
+#include "TH1F.h"
+
+using namespace edm;
+using namespace std;
+
+// Constructor
+ElectronEfficiencyPlotter::ElectronEfficiencyPlotter(const edm::ParameterSet &ps) {
+  ptBin_ = ps.getParameter<int>("ptBin");
+  ptMin_ = ps.getParameter<double>("ptMin");
+  ptMax_ = ps.getParameter<double>("ptMax");
+
+  ID_ = ps.getParameter<string>("sctElectronID");
+  theFolder_ = ps.getParameter<string>("folder");
+  sourceFolder_ = ps.getParameter<string>("srcFolder");
+}
+
+void ElectronEfficiencyPlotter::dqmEndJob(DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter) {
+  ibooker.setCurrentFolder(theFolder_);
+
+  h_eff_pt_EB_doubleEG_HLT =
+      ibooker.book1D("Eff_pt_barrel_DSTdoubleEG", "DSTdoubleEG Eff. vs Pt (barrel)", ptBin_, ptMin_, ptMax_);
+  h_eff_pt_EE_doubleEG_HLT =
+      ibooker.book1D("Eff_pt_endcap_DSTdoubleEG", "DSTdoubleEG Eff. vs Pt (endcap)", ptBin_, ptMin_, ptMax_);
+  h_eff_pt_EB_singlePhoton_HLT =
+      ibooker.book1D("Eff_pt_barrel_DSTsinglePhoton", "DSTsinglePhoton Eff. vs Pt (barrel)", ptBin_, ptMin_, ptMax_);
+  h_eff_pt_EE_singlePhoton_HLT =
+      ibooker.book1D("Eff_pt_endcap_DSTsinglePhoton", "DSTsinglePhoton Eff. vs Pt (endcap)", ptBin_, ptMin_, ptMax_);
+
+  // Axis title
+  h_eff_pt_EB_singlePhoton_HLT->setAxisTitle("p_{T} (GeV)", 1);
+  h_eff_pt_EE_singlePhoton_HLT->setAxisTitle("p_{T} (GeV)", 1);
+  h_eff_pt_EB_doubleEG_HLT->setAxisTitle("p_{T} (GeV)", 1);
+  h_eff_pt_EE_doubleEG_HLT->setAxisTitle("p_{T} (GeV)", 1);
+
+  MonitorElement *Numerator_pt_barrel_doubleEG_hlt =
+      igetter.get(sourceFolder_ +
+                  "/resonanceZ_Tag_pat_Probe_sctElectron_passDoubleEG_DST_"
+                  "fireTrigObj_Pt_Barrel");
+  MonitorElement *Numerator_pt_endcap_doubleEG_hlt =
+      igetter.get(sourceFolder_ +
+                  "/resonanceZ_Tag_pat_Probe_sctElectron_passDoubleEG_DST_"
+                  "fireTrigObj_Pt_Endcap");
+  MonitorElement *Numerator_pt_barrel_singlePhoton_hlt =
+      igetter.get(sourceFolder_ +
+                  "/resonanceZ_Tag_pat_Probe_sctElectron_passSinglePhoton_DST_"
+                  "fireTrigObj_Pt_Barrel");
+  MonitorElement *Numerator_pt_endcap_singlePhoton_hlt =
+      igetter.get(sourceFolder_ +
+                  "/resonanceZ_Tag_pat_Probe_sctElectron_passSinglePhoton_DST_"
+                  "fireTrigObj_Pt_Endcap");
+  MonitorElement *Denominator_pt_barrel =
+      igetter.get(sourceFolder_ + "/resonanceZ_Tag_pat_Probe_sctElectron_Pt_Barrel");
+  MonitorElement *Denominator_pt_endcap =
+      igetter.get(sourceFolder_ + "/resonanceZ_Tag_pat_Probe_sctElectron_Pt_Endcap");
+
+  if (Numerator_pt_barrel_doubleEG_hlt && Denominator_pt_barrel)
+    calculateEfficiency(Numerator_pt_barrel_doubleEG_hlt, Denominator_pt_barrel, h_eff_pt_EB_doubleEG_HLT);
+  if (Numerator_pt_endcap_doubleEG_hlt && Denominator_pt_endcap)
+    calculateEfficiency(Numerator_pt_endcap_doubleEG_hlt, Denominator_pt_endcap, h_eff_pt_EE_doubleEG_HLT);
+  if (Numerator_pt_barrel_singlePhoton_hlt && Denominator_pt_barrel)
+    calculateEfficiency(Numerator_pt_barrel_singlePhoton_hlt, Denominator_pt_barrel, h_eff_pt_EB_singlePhoton_HLT);
+  if (Numerator_pt_endcap_singlePhoton_hlt && Denominator_pt_endcap)
+    calculateEfficiency(Numerator_pt_endcap_singlePhoton_hlt, Denominator_pt_endcap, h_eff_pt_EE_singlePhoton_HLT);
+}
+
+void ElectronEfficiencyPlotter::calculateEfficiency(MonitorElement *Numerator,
+                                                    MonitorElement *Denominator,
+                                                    MonitorElement *Efficiency) {
+  TH1F *h_numerator_pt = Numerator->getTH1F();
+  TH1F *h_denominator_pt = Denominator->getTH1F();
+  TH1F *h_eff_pt = Efficiency->getTH1F();
+  if (h_eff_pt->GetSumw2N() == 0)
+    h_eff_pt->Sumw2();
+
+  // ReBin
+  int nBins = h_eff_pt->GetNbinsX();
+  double *binEdges = new double[nBins + 1];
+  for (int i = 0; i <= nBins; i++)
+    binEdges[i] = h_eff_pt->GetBinLowEdge(i + 1);
+
+  TH1F *h_numerator_pt_rebin = (TH1F *)h_numerator_pt->Rebin(nBins, "num_pt_rebinned", binEdges);
+  TH1F *h_denominator_pt_rebin = (TH1F *)h_denominator_pt->Rebin(nBins, "num_pt_rebinned", binEdges);
+  h_eff_pt->Divide(h_numerator_pt_rebin, h_denominator_pt_rebin, 1., 1., "B");
+}
+
+DEFINE_FWK_MODULE(ElectronEfficiencyPlotter);

HLTriggerOffline/Scouting/plugins/ElectronEfficiencyPlotter.h

@@ -0,0 +1,45 @@
+#ifndef HLTriggerOffline_Scouting_ElectronEfficiencyPlotter_H
+#define HLTriggerOffline_Scouting_ElectronEfficiencyPlotter_H
+
+#include "DQMServices/Core/interface/DQMEDHarvester.h"
+#include "DQMServices/Core/interface/DQMStore.h"
+#include "DataFormats/Common/interface/Handle.h"
+#include "FWCore/Framework/interface/ESHandle.h"
+#include "FWCore/Framework/interface/Event.h"
+#include "FWCore/Framework/interface/Frameworkfwd.h"
+#include "FWCore/Framework/interface/LuminosityBlock.h"
+#include "FWCore/Framework/interface/MakerMacros.h"
+#include "FWCore/ParameterSet/interface/ParameterSet.h"
+#include "FWCore/ServiceRegistry/interface/Service.h"
+
+class ElectronEfficiencyPlotter : public DQMEDHarvester {
+public:
+  // Constructor
+  ElectronEfficiencyPlotter(const edm::ParameterSet& ps);
+  // Destructor
+  ~ElectronEfficiencyPlotter() override = default;
+
+protected:
+  // DQM Client Diagnostic
+  void dqmEndJob(DQMStore::IBooker&, DQMStore::IGetter&) override;
+
+private:
+  // counters
+
+  int ptBin_;
+  double ptMin_;
+  double ptMax_;
+
+  std::string ID_;
+
+  MonitorElement* h_eff_pt_EB_doubleEG_HLT;
+  MonitorElement* h_eff_pt_EE_doubleEG_HLT;
+  MonitorElement* h_eff_pt_EB_singlePhoton_HLT;
+  MonitorElement* h_eff_pt_EE_singlePhoton_HLT;
+
+  std::string theFolder_;
+  std::string sourceFolder_;
+  void calculateEfficiency(MonitorElement* Numerator, MonitorElement* Denominator, MonitorElement* Efficiency);
+};
+
+#endif