Missing CaloParticle From PU in CMSSW_15_0_0_pre2

[waredjeb]

Hello,

I suspect there's a problem with the Caloparticle starting from CMSSW_15_0_0_pre2. I was trying to generate some CloseByParticleGun events in PU 200 [1] and I noticed that the number of CaloParticle is always 1 (== number of particles generated), suggesting that the CaloParticle coming from the PU are not present.

For a quick check I generated 5 events with PU 50 and I used the TiclDebugger that runs some sort of debugging on the CaloParticle, adding a simple cout on the number of CP.
The same test have been done in CMSSW_15_0_0_pre1.
All the files / logs can be found here

Unfortunately the TiclDebugger is not up to date so it requires some changes in case one wants to use it [2]
But even inspecting the step2.root file from pre2 one can see how the number of CP represents only the CP from the generated particle. The SimCluster seem to be fine instead

I've also tried to look at recent PR tests but I couldn't find any Standard Mixing PU samples for Phase-2. But also in the Premixing the CP are not present, and only the SC are there.
Am I missing something?

Pre2

CaloParticle Plot

SimCluster Plot

Pre1

CaloParticle Plot

SimCluster Plot

PR47327 step2.root from 29834.999

CaloParticle Plot

SimCluster Plot

[1] WF = 29896.203
[2]

diff --git a/RecoHGCal/TICL/test/TiclDebugger.cc b/RecoHGCal/TICL/test/TiclDebugger.cc
index beaf4bbc872..2ed355ae345 100644
--- a/RecoHGCal/TICL/test/TiclDebugger.cc
+++ b/RecoHGCal/TICL/test/TiclDebugger.cc
@@ -105,6 +105,7 @@ void TiclDebugger::analyze(const edm::Event& iEvent, const edm::EventSetup& iSet
   edm::Handle<std::vector<CaloParticle>> caloParticlesH;
   iEvent.getByToken(caloParticlesToken_, caloParticlesH);
   auto const& caloParticles = *caloParticlesH.product();
+  LogVerbatim("TICLDebugger") << "Caloparticle size " << caloParticles.size() << std::endl;
   std::vector<std::pair<int, float>> bestCPMatches;
 
   auto bestCaloParticleMatches = [&](const ticl::Trackster& t) -> void {
diff --git a/RecoHGCal/TICL/test/ticlDebugger_cfg.py b/RecoHGCal/TICL/test/ticlDebugger_cfg.py
index 76cef2d34e4..1d92873b39b 100644
--- a/RecoHGCal/TICL/test/ticlDebugger_cfg.py
+++ b/RecoHGCal/TICL/test/ticlDebugger_cfg.py
@@ -3,21 +3,26 @@ import FWCore.ParameterSet.Config as cms
 process = cms.Process("TICLDEBUG")
 
 process.load("FWCore.MessageService.MessageLogger_cfi")
-process.load('Configuration.Geometry.GeometryExtendedRun4D49Reco_cff')
+process.load('Configuration.Geometry.GeometryExtendedRun4D110Reco_cff')
 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
 from Configuration.AlCa.GlobalTag import GlobalTag
-process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:phase2_realistic_T15', '')
+process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:phase2_realistic_T33', '')
 
-process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(1) )
+process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(10) )
 
 process.source = cms.Source("PoolSource",
     # replace 'myfile.root' with the source file you want to use
     fileNames = cms.untracked.vstring(
-        'file:step3.root'
+        'file:/eos/user/w/wredjeb/HGCAL/TestCPFromPU/step3_pre1.root'
     )
 )
 
 process.load("RecoHGCal.TICL.ticlDebugger_cfi")
+from RecoHGCal.TICL.ticlDebugger_cfi import ticlDebugger as ticlDebugger_
+process.ticlDebugger = ticlDebugger_.clone(
+    trackstersMerge = cms.InputTag("ticlCandidate")
+)
+
 process.load("SimGeneral.Debugging.caloParticleDebugger_cfi")
 # MessageLogger customizations

@AuroraPerego @felicepantaleo @rovere FYI

[cmsbuild]

cms-bot internal usage

[cmsbuild]

A new Issue was created by @waredjeb.

@Dr15Jones, @antoniovilela, @makortel, @mandrenguyen, @rappoccio, @sextonkennedy, @smuzaffar can you please review it and eventually sign/assign? Thanks.

cms-bot commands are listed here

[makortel]

assign simulation, upgrade

[cmsbuild]

New categories assigned: simulation,upgrade

@civanch,@kpedro88,@mdhildreth,@Moanwar,@srimanob,@subirsarkar you have been requested to review this Pull request/Issue and eventually sign? Thanks

[kpedro88]

@waredjeb can you please provide the exact cmsDriver commands that you used?

[waredjeb]

@waredjeb can you please provide the exact cmsDriver commands that you used?

I just used runTheMatrix.py -w upgrade -l 29896.203 -j 4

The cmsDriver commands are here:

Step1

cmsDriver CE_E_Front_120um_cfi -s GEN,SIM -n 10 --conditions auto:phase2_realistic_T33 --beamspot HGCALCloseBy --datatier GEN-SIM --eventcontent FEVTDEBUG --geometry ExtendedRun4D110 --era Phase2C17I13M9 --relval 9000,100  --fileout file:step1.root

step2

cmsDriver.py step2 -s DIGI:pdigi_valid,L1TrackTrigger,L1,L1P2GT,DIGI2RAW,HLT:@relvalRun4 --conditions auto:phase2_realistic_T33 --datatier GEN-SIM-DIGI-RAW -n 10 --eventcontent FEVTDEBUGHLT --geometry ExtendedRun4D110 --era Phase2C17I13M9 --procModifiers ticl_v5 --pileup AVE_200_BX_25ns --pileup_input das:/RelValMinBias_14TeV/CMSSW_14_1_0-141X_mcRun4_realistic_v1_STD_RegeneratedGS_2026D110_noPU-v1/GEN-SIM --filein file:step1.root --fileout file:step2.root

step3

cmsDriver.py step3 -s RAW2DIGI,RECO,RECOSIM,PAT,VALIDATION:@phase2Validation+@miniAODValidation,DQM:@phase2+@miniAODDQM --conditions auto:phase2_realistic_T33 --datatier GEN-SIM-RECO,MINIAODSIM,DQMIO -n 10 --eventcontent FEVTDEBUGHLT,MINIAODSIM,DQM --geometry ExtendedRun4D110 --era Phase2C17I13M9 --procModifiers ticl_v5 --pileup AVE_200_BX_25ns --pileup_input das:/RelValMinBias_14TeV/CMSSW_14_1_0-141X_mcRun4_realistic_v1_STD_RegeneratedGS_2026D110_noPU-v1/GEN-SIM  --filein file:step2.root --fileout file:step3.root

For the tests I modify the averageNumber of PU from 200 to 50 in both step2.py and step3.py

 process.mix.input.nbPileupEvents.averageNumber = cms.double(50.000000)

[felicepantaleo]

@civanch @kpedro88 how is your investigation going? We won't be able to go on with old samples for long for our studies

[kpedro88]

Resolved by #47682 and #47765 (backport: #47771)

[kpedro88]

+simulation