Merge pull request #48341 from bsunanda/Run3-alca253

Run3-alca253 Provide run range in the IsoTrackAnalyzer code and provide cfg's for ACaReco and AlCaRaw o/p's

Author: cmsbuild

Calibration/HcalCalibAlgos/plugins/HcalIsoTrackAnalyzer.cc

@@ -45,8 +45,8 @@ class HcalIsoTrackAnalyzer : public edm::one::EDAnalyzer<edm::one::WatchRuns, ed
   double gainFactor(const HcalDbService* dbserv, const HcalDetId& id);
 
   const double pTrackLow_, pTrackHigh_;
-  const int useRaw_, dataType_, unCorrect_;
-  const bool fillInRange_;
+  const int useRaw_, dataType_, unCorrect_, runLow_, runHigh_;
+  const bool fillInRange_, fillRunRange_;
   const edm::InputTag labelIsoTkVar_, labelIsoTkEvt_;
   const std::vector<int> debEvents_;
   const edm::ESGetToken<HcalTopology, HcalRecNumberingRecord> tok_htopo_;
@@ -89,7 +89,10 @@ HcalIsoTrackAnalyzer::HcalIsoTrackAnalyzer(const edm::ParameterSet& iConfig)
       useRaw_(iConfig.getUntrackedParameter<int>("useRaw", 0)),
       dataType_(iConfig.getUntrackedParameter<int>("dataType", 0)),
       unCorrect_(iConfig.getUntrackedParameter<int>("unCorrect", 0)),
+      runLow_(iConfig.getUntrackedParameter<int>("runLow", 0)),
+      runHigh_(iConfig.getUntrackedParameter<int>("runHigh", -1)),
       fillInRange_(iConfig.getUntrackedParameter<bool>("fillInRange", false)),
+      fillRunRange_((runLow_ < runHigh_) && (runLow_ > 0)),
       labelIsoTkVar_(iConfig.getParameter<edm::InputTag>("isoTrackVarLabel")),
       labelIsoTkEvt_(iConfig.getParameter<edm::InputTag>("isoTrackEvtLabel")),
       debEvents_(iConfig.getParameter<std::vector<int>>("debugEvents")),
@@ -112,7 +115,8 @@ HcalIsoTrackAnalyzer::HcalIsoTrackAnalyzer(const edm::ParameterSet& iConfig)
   edm::LogVerbatim("HcalIsoTrack") << "Parameters read from config file \n\t momentumLow_ " << pTrackLow_
                                    << "\t momentumHigh_ " << pTrackHigh_ << "\t useRaw_ " << useRaw_
                                    << "\t dataType_      " << dataType_ << "\t unCorrect " << unCorrect_
-                                   << "\t fillInRange " << fillInRange_ << "\t and " << debEvents_.size()
+                                   << "\t fillInRange " << fillInRange_ << "\t fillRunRange " << fillRunRange_
+                                   << " for " << runLow_ << ":" << runHigh_ << "\t and " << debEvents_.size()
                                    << " events to be debugged";
 }
 
@@ -264,6 +268,10 @@ void HcalIsoTrackAnalyzer::analyze(edm::Event const& iEvent, edm::EventSetup con
         if ((t_p < pTrackLow_) || (t_p > pTrackHigh_))
           select = false;
       }
+      if (select && fillRunRange_) {
+        if ((t_Run < runLow_) || (t_Run > runHigh_))
+          select = false;
+      }
       if (select) {
         tree->Fill();
         edm::LogVerbatim("HcalIsoTrackX") << "Run " << t_Run << " Event " << t_Event << " p " << t_p;
@@ -306,7 +314,13 @@ void HcalIsoTrackAnalyzer::analyze(edm::Event const& iEvent, edm::EventSetup con
       t_ietaGood = itr->ietaGood_;
       t_trackType = itr->trackType_;
       t_hltbits = itr->hltbits_;
-      tree2->Fill();
+      bool select(true);
+      if (fillRunRange_) {
+        if ((t_RunNo < static_cast<unsigned int>(runLow_)) || (t_RunNo > static_cast<unsigned int>(runHigh_)))
+          select = false;
+      }
+      if (select)
+        tree2->Fill();
     }
   } else {
     edm::LogVerbatim("HcalIsoTrack") << "Cannot find HcalIsoTrkEventVariablesCollections";
@@ -412,6 +426,8 @@ void HcalIsoTrackAnalyzer::fillDescriptions(edm::ConfigurationDescriptions& desc
   desc.addUntracked<int>("dataType", 0);
   desc.addUntracked<int>("unCorrect", 0);
   desc.addUntracked<bool>("fillInRange", false);
+  desc.addUntracked<int>("runLow", 0);
+  desc.addUntracked<int>("runHigh", -1);
   desc.add<edm::InputTag>("isoTrackVarLabel", edm::InputTag("alcaHcalIsotrkProducer", "HcalIsoTrack"));
   desc.add<edm::InputTag>("isoTrackEvtLabel", edm::InputTag("alcaHcalIsotrkProducer", "HcalIsoTrackEvent"));
   std::vector<int> events;

Calibration/HcalCalibAlgos/test/python/isoTrackAnalysisAlCaRaw_cfg.py

@@ -0,0 +1,41 @@
+import FWCore.ParameterSet.Config as cms
+
+from Configuration.Eras.Era_Run3_2025_cff import Run3_2025
+process = cms.Process("ANALYSIS",Run3_2025)
+
+process.load("FWCore.MessageService.MessageLogger_cfi")
+process.load('Configuration.StandardSequences.Services_cff')
+process.load("Configuration.StandardSequences.GeometryRecoDB_cff")
+process.load('Configuration.StandardSequences.MagneticField_cff')
+process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
+
+from Configuration.AlCa.autoCond import autoCond
+process.GlobalTag.globaltag = autoCond['phase1_2025_realistic']
+
+if 'MessageLogger' in process.__dict__:
+    process.MessageLogger.HcalIsoTrack=dict()
+process.MessageLogger.cerr.FwkReport.reportEvery = 1
+process.options = cms.untracked.PSet( wantSummary = cms.untracked.bool(True) )
+
+process.load('Calibration.HcalCalibAlgos.hcalIsoTrackAnalyzer_cfi')
+process.source = cms.Source("PoolSource", 
+                            fileNames = cms.untracked.vstring(
+       'file:/eos/cms/store/group/dpg_hcal/comm_hcal/suman/13p6TeV/2025EraC/AlcaRaw_Commissioning/23144728-40d9-4349-af8c-f6224bacb99d.root',
+    )
+)
+
+process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )
+
+process.TFileService = cms.Service("TFileService",
+   fileName = cms.string('outputNew.root')
+)
+
+process.hcalIsoTrackAnalyzer.fillInRange = True  # fils only 40-60 GeV
+process.hcalIsoTrackAnalyzer.useRaw = 0   # 1 for Raw
+process.hcalIsoTrackAnalyzer.unCorrect = 1   # 1 for RespCorr; 2 for Gain
+process.hcalIsoTrackAnalyzer.dataType = 1   # 1 for from Jet; 0 for others
+process.hcalIsoTrackAnalyzer.isoTrackVarLabel = ["alcaHcalIsotrkFromAlCaRawProducer", "HcalIsoTrack"]
+process.hcalIsoTrackAnalyzer.isoTrackEvtLabel = ["alcaHcalIsotrkFromAlCaRawProducer", "HcalIsoTrackEvent"]
+
+process.p = cms.Path(process.hcalIsoTrackAnalyzer)
+

Calibration/HcalCalibAlgos/test/python/isoTrackAnalysis_cfg.py

@@ -1,16 +1,16 @@
 import FWCore.ParameterSet.Config as cms
 
-from Configuration.Eras.Era_Run2_2018_cff import Run2_2018
-process = cms.Process("ANALYSIS",Run2_2018)
+from Configuration.Eras.Era_Run3_2025_cff import Run3_2025
+process = cms.Process("ANALYSIS",Run3_2025)
 
 process.load("FWCore.MessageService.MessageLogger_cfi")
 process.load('Configuration.StandardSequences.Services_cff')
 process.load("Configuration.StandardSequences.GeometryRecoDB_cff")
 process.load("Configuration.StandardSequences.MagneticField_cff")
 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
 
-from Configuration.AlCa.GlobalTag import GlobalTag
-process.GlobalTag = GlobalTag(process.GlobalTag,'auto:run2_data','')
+from Configuration.AlCa.autoCond import autoCond
+process.GlobalTag.globaltag = autoCond['phase1_2025_realistic']
 
 if 'MessageLogger' in process.__dict__:
     process.MessageLogger.HcalIsoTrack=dict()
@@ -30,7 +30,10 @@
    fileName = cms.string('outputNew.root')
 )
 
+process.hcalIsoTrackAnalyzer.fillInRange = True  # fils only 40-60 GeV
 process.hcalIsoTrackAnalyzer.useRaw = 0   # 1 for Raw
+process.hcalIsoTrackAnalyzer.unCorrect = 1   # 1 for RespCorr; 2 for Gain
+process.hcalIsoTrackAnalyzer.dataType = 0   # 1 for from Jet; 0 for others
 
 process.p = cms.Path(process.hcalIsoTrackAnalyzer)