Skip to content

Updating the TICL offline NanoAOD #51725

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
cmsbuild merged 3 commits into from
Aug 31, 2026
Merged

Updating the TICL offline NanoAOD #51725

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
b23905a
update offline TICL nanoAOD
Moanwar Aug 18, 2026
c87ee4c
code format check
Moanwar Aug 18, 2026
9b64388
small update
Moanwar Aug 21, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
3 changes: 3 additions & 0 deletions DPGAnalysis/HGCalNanoAOD/plugins/BuildFile.xml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -5,6 +5,9 @@
<use name="DataFormats/NanoAOD"/>
<use name="DataFormats/ParticleFlowCandidate"/>
<use name="DataFormats/StdDictionaries"/>
<use name="DataFormats/TrackReco"/>
<use name="DataFormats/MuonReco"/>
<use name="RecoParticleFlow/PFProducer"/>
<use name="SimDataFormats/CaloAnalysis"/>
<use name="FWCore/Common"/>
<use name="FWCore/Framework"/>
Expand Down
28 changes: 22 additions & 6 deletions DPGAnalysis/HGCalNanoAOD/plugins/SimTracksterTableProducer.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -48,6 +48,7 @@ class SimTracksterTableProducer : public edm::global::EDProducer<> {
const size_t nSimTracksters = simTrackstersHandle.isValid() ? simTrackstersHandle->size() : 0;

static constexpr float default_value = std::numeric_limits<float>::quiet_NaN();
static constexpr int default_int_value = -1;

std::vector<float> boundaryX(nSimTracksters, default_value);
std::vector<float> boundaryY(nSimTracksters, default_value);
Expand All @@ -61,6 +62,8 @@ class SimTracksterTableProducer : public edm::global::EDProducer<> {
std::vector<float> simTime(nSimTracksters, default_value);
std::vector<float> genPt(nSimTracksters, default_value);
std::vector<float> mass(nSimTracksters, default_value);
std::vector<int> caloParticleIdx(nSimTracksters, default_int_value);
std::vector<int8_t> isPU(nSimTracksters, default_int_value);

if ((simTrackstersHandle.isValid() && caloParticlesHandle.isValid() && simClustersHandle.isValid() &&
cpToSCMapHandle.isValid()) ||
Expand All @@ -71,7 +74,7 @@ class SimTracksterTableProducer : public edm::global::EDProducer<> {
const auto& cpToSCMap = *cpToSCMapHandle;

//utility lambda for filling vectors
auto fillVectors = [&](const auto& obj, size_t iSim, float time) {
auto fillVectors = [&](const auto& obj, size_t iSim, float time, int cpIdx) {
const auto& simTrack = obj.g4Tracks()[0];
const auto caloPt = obj.pt();
const auto simHitSumEnergy = obj.simEnergy();
Expand All @@ -90,26 +93,35 @@ class SimTracksterTableProducer : public edm::global::EDProducer<> {
simEnergy[iSim] = simHitSumEnergy;
genPt[iSim] = caloPt;
mass[iSim] = caloMass;

// PU flag: a non-zero event/bunch-crossing on the
// seed g4Track means the CaloParticle/SimCluster came from a pileup interaction.
isPU[iSim] = (simTrack.eventId().event() != 0 || simTrack.eventId().bunchCrossing() != 0) ? 1 : 0;

caloParticleIdx[iSim] = cpIdx;
};

for (size_t iSim = 0; iSim < simTracksters.size(); ++iSim) {
const auto& simT = simTracksters[iSim];
float time = default_value;

if (simT.seedID() == caloParticlesHandle.id()) {
const int cpIdx = static_cast<int>(simT.seedIndex());
const auto& cp = caloParticles[simT.seedIndex()];
time = cp.simTime();
fillVectors(cp, iSim, time);
fillVectors(cp, iSim, time, cpIdx);
} else {
const auto& sc = simClusters[simT.seedIndex()];
int cpIdx = default_int_value;
//SCtoCP map not availalbe, use CPtoSC map instead
for (const auto& [cpIdx, scVec] : cpToSCMap) {
for (const auto& [cpIdxCandidate, scVec] : cpToSCMap) {
if (std::ranges::find(scVec, simT.seedIndex()) != scVec.end()) {
time = caloParticles[cpIdx].simTime();
cpIdx = static_cast<int>(cpIdxCandidate);
time = caloParticles[cpIdxCandidate].simTime();
break; //dont need to check further
}
}
fillVectors(sc, iSim, time);
fillVectors(sc, iSim, time, cpIdx);
}
}
}
Expand All @@ -125,7 +137,7 @@ class SimTracksterTableProducer : public edm::global::EDProducer<> {
simTrackstersTable->addColumn<float>(
"boundaryEta", boundaryEta, "CaloVolume boundary pseudorapidity of associated Simobject", precision_);
simTrackstersTable->addColumn<float>(
"boundaryPhi", boundaryEta, "CaloVolume boundary phi of associated Simobject", precision_);
"boundaryPhi", boundaryPhi, "CaloVolume boundary phi of associated Simobject", precision_);
simTrackstersTable->addColumn<float>(
"boundaryPx", boundaryPx, "X component of momentum at CaloVolume boundary of associated Simobject", precision_);
simTrackstersTable->addColumn<float>(
Expand All @@ -135,6 +147,10 @@ class SimTracksterTableProducer : public edm::global::EDProducer<> {
simTrackstersTable->addColumn<float>("simTime", simTime, "Sim-Time of simulated object [ns]", precision_);
simTrackstersTable->addColumn<float>("genPt", genPt, "Gen-pT associated with SimObject", precision_);
simTrackstersTable->addColumn<float>("mass", mass, "mass associated with SimObject", precision_);
simTrackstersTable->addColumn<int>(
"caloParticleIdx", caloParticleIdx, "Index of the parent CaloParticle in the CaloParticle collection");
simTrackstersTable->addColumn<int>(
"isPU", isPU, "PU flag of the associated Simobject: 1 = pileup (non-zero event/bx), 0 = otherwise");

event.put(std::move(simTrackstersTable), tableName_);
}
Expand Down
Loading