Skip to content

Commit d7b61f7

Browse files
nzardoshnzardosh
authored
[PWGJE] Adding the option of skipping gap events in MC (#9361)
Please fix the linter in a subsequent PR.
1 parent bc3af27 commit d7b61f7

17 files changed

+135
-87
lines changed

‎PWGJE/Core/JetBkgSubUtils.cxx‎

Lines changed: 0 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -159,7 +159,6 @@ fastjet::PseudoJet JetBkgSubUtils::doRhoAreaSub(fastjet::PseudoJet& jet, double
159159

160160
std::vector<fastjet::PseudoJet> JetBkgSubUtils::doEventConstSub(std::vector<fastjet::PseudoJet>& inputParticles, double rhoParam, double rhoMParam)
161161
{
162-
163162
JetBkgSubUtils::initialise();
164163
fastjet::contrib::ConstituentSubtractor constituentSub(rhoParam, rhoMParam);
165164
constituentSub.set_distance_type(fastjet::contrib::ConstituentSubtractor::deltaR); /// deltaR=sqrt((y_i-y_j)^2+(phi_i-phi_j)^2)), longitudinal Lorentz invariant

‎PWGJE/Core/JetDerivedDataUtilities.h‎

Lines changed: 9 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -40,9 +40,17 @@ enum JCollisionSel {
4040
sel7KINT7 = 9
4141
};
4242

43+
enum JCollisionSubGeneratorId {
44+
none = -1,
45+
mbGap = 0
46+
};
47+
4348
template <typename T>
44-
bool selectCollision(T const& collision, int eventSelection = -1)
49+
bool selectCollision(T const& collision, int eventSelection = -1, bool skipMBGapEvents = true)
4550
{
51+
if (skipMBGapEvents && collision.subGeneratorId() == JCollisionSubGeneratorId::mbGap) {
52+
return false;
53+
}
4654
if (eventSelection == -1) {
4755
return true;
4856
}

‎PWGJE/DataModel/Jet.h‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -196,7 +196,7 @@ DECLARE_JET_DUPLICATE_TABLES_LEVELS(Charged, JTrackSub, HfD0Bases, HfD0PBases, "
196196
#undef STRINGIFY
197197
#undef DECLARE_JET_DUPLICATE_TABLES_LEVELS
198198

199-
using JetCollisions = JCollisions;
199+
using JetCollisions = o2::soa::Join<JCollisions, JCollisionMcInfos>;
200200
using JetCollision = JetCollisions::iterator;
201201
using JetCollisionsMCD = o2::soa::Join<JetCollisions, JMcCollisionLbs>;
202202
using JetCollisionMCD = o2::soa::Join<JetCollisions, JMcCollisionLbs>::iterator;

‎PWGJE/DataModel/JetReducedData.h‎

Lines changed: 9 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -69,6 +69,8 @@ DECLARE_SOA_COLUMN(PosY, posY, float);
6969
DECLARE_SOA_COLUMN(PosZ, posZ, float);
7070
DECLARE_SOA_COLUMN(Multiplicity, multiplicity, float);
7171
DECLARE_SOA_COLUMN(Centrality, centrality, float);
72+
DECLARE_SOA_COLUMN(Weight, weight, float);
73+
DECLARE_SOA_COLUMN(SubGeneratorId, subGeneratorId, int);
7274
DECLARE_SOA_COLUMN(EventSel, eventSel, uint16_t);
7375
DECLARE_SOA_BITMAP_COLUMN(Alias, alias, 32);
7476
DECLARE_SOA_COLUMN(TrackOccupancyInTimeRange, trackOccupancyInTimeRange, int);
@@ -107,6 +109,10 @@ DECLARE_SOA_TABLE_STAGED(JCollisions, "JCOLLISION",
107109
using JCollision = JCollisions::iterator;
108110
using StoredJCollision = StoredJCollisions::iterator;
109111

112+
DECLARE_SOA_TABLE_STAGED(JCollisionMcInfos, "JCOLLISIONMCINFO",
113+
jcollision::Weight,
114+
jcollision::SubGeneratorId);
115+
110116
DECLARE_SOA_TABLE_STAGED(JEMCCollisionLbs, "JEMCCOLLISIONLB",
111117
jcollision::IsAmbiguous,
112118
jcollision::IsEMCALReadout);
@@ -149,13 +155,15 @@ DECLARE_SOA_COLUMN(PosX, posX, float);
149155
DECLARE_SOA_COLUMN(PosY, posY, float);
150156
DECLARE_SOA_COLUMN(PosZ, posZ, float);
151157
DECLARE_SOA_COLUMN(Weight, weight, float);
158+
DECLARE_SOA_COLUMN(SubGeneratorId, subGeneratorId, int);
152159
} // namespace jmccollision
153160
DECLARE_SOA_TABLE_STAGED(JMcCollisions, "JMCCOLLISION",
154161
o2::soa::Index<>,
155162
jmccollision::PosX,
156163
jmccollision::PosY,
157164
jmccollision::PosZ,
158-
jmccollision::Weight);
165+
jmccollision::Weight,
166+
jmccollision::SubGeneratorId);
159167

160168
using JMcCollision = JMcCollisions::iterator;
161169
using StoredJMcCollision = StoredJMcCollisions::iterator;

‎PWGJE/JetFinders/jetFinder.cxx‎

Lines changed: 8 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -18,6 +18,7 @@
1818
#include <vector>
1919

2020
#include "PWGJE/Core/JetFindingUtilities.h"
21+
#include "PWGJE/Core/JetDerivedDataUtilities.h"
2122
#include "Framework/runDataProcessing.h"
2223

2324
using namespace o2;
@@ -41,6 +42,7 @@ struct JetFinderTask {
4142
Configurable<int> trackOccupancyInTimeRangeMax{"trackOccupancyInTimeRangeMax", 999999, "maximum occupancy of tracks in neighbouring collisions in a given time range"};
4243
Configurable<std::string> eventSelections{"eventSelections", "sel8", "choose event selection"};
4344
Configurable<std::string> triggerMasks{"triggerMasks", "", "possible JE Trigger masks: fJetChLowPt,fJetChHighPt,fTrackLowPt,fTrackHighPt,fJetD0ChLowPt,fJetD0ChHighPt,fJetLcChLowPt,fJetLcChHighPt,fEMCALReadout,fJetFullHighPt,fJetFullLowPt,fJetNeutralHighPt,fJetNeutralLowPt,fGammaVeryHighPtEMCAL,fGammaVeryHighPtDCAL,fGammaHighPtEMCAL,fGammaHighPtDCAL,fGammaLowPtEMCAL,fGammaLowPtDCAL,fGammaVeryLowPtEMCAL,fGammaVeryLowPtDCAL"};
45+
Configurable<bool> skipMBGapEvents{"skipMBGapEvents", true, "decide to run over MB gap events or not"};
4446

4547
// track level configurables
4648
Configurable<float> trackPtMin{"trackPtMin", 0.15, "minimum track pT"};
@@ -139,9 +141,9 @@ struct JetFinderTask {
139141
}
140142

141143
aod::EMCALClusterDefinition clusterDefinition = aod::emcalcluster::getClusterDefinitionFromString(clusterDefinitionS.value);
142-
Filter collisionFilter = (nabs(aod::jcollision::posZ) < vertexZCut && aod::jcollision::centrality >= centralityMin && aod::jcollision::centrality < centralityMax && aod::jcollision::trackOccupancyInTimeRange <= trackOccupancyInTimeRangeMax);
144+
Filter collisionFilter = (nabs(aod::jcollision::posZ) < vertexZCut && aod::jcollision::centrality >= centralityMin && aod::jcollision::centrality < centralityMax && aod::jcollision::trackOccupancyInTimeRange <= trackOccupancyInTimeRangeMax && (!skipMBGapEvents || aod::jcollision::subGeneratorId != static_cast<int>(jetderiveddatautilities::JCollisionSubGeneratorId::mbGap)));
145+
Filter mcCollisionFilter = (!skipMBGapEvents || aod::jmccollision::subGeneratorId != static_cast<int>(jetderiveddatautilities::JCollisionSubGeneratorId::mbGap)); // should we add a posZ vtx cut here or leave it to analysers?
143146
Filter trackCuts = (aod::jtrack::pt >= trackPtMin && aod::jtrack::pt < trackPtMax && aod::jtrack::eta >= trackEtaMin && aod::jtrack::eta <= trackEtaMax && aod::jtrack::phi >= trackPhiMin && aod::jtrack::phi <= trackPhiMax); // do we need eta cut both here and in globalselection?
144-
145147
Filter partCuts = (aod::jmcparticle::pt >= trackPtMin && aod::jmcparticle::pt < trackPtMax && aod::jmcparticle::eta >= trackEtaMin && aod::jmcparticle::eta <= trackEtaMax && aod::jmcparticle::phi >= trackPhiMin && aod::jmcparticle::phi <= trackPhiMax);
146148
Filter clusterFilter = (aod::jcluster::definition == static_cast<int>(clusterDefinition) && aod::jcluster::eta >= clusterEtaMin && aod::jcluster::eta <= clusterEtaMax && aod::jcluster::phi >= clusterPhiMin && aod::jcluster::phi <= clusterPhiMax && aod::jcluster::energy >= clusterEnergyMin && aod::jcluster::time > clusterTimeMin && aod::jcluster::time < clusterTimeMax && (clusterRejectExotics && aod::jcluster::isExotic != true));
147149

@@ -197,7 +199,7 @@ struct JetFinderTask {
197199
}
198200
PROCESS_SWITCH(JetFinderTask, processFullJets, "Data and reco level jet finding for full and neutral jets", false);
199201

200-
void processParticleLevelChargedJets(aod::JetMcCollision const& collision, soa::Filtered<aod::JetParticles> const& particles)
202+
void processParticleLevelChargedJets(soa::Filtered<aod::JetMcCollisions>::iterator const& collision, soa::Filtered<aod::JetParticles> const& particles)
201203
{
202204
// TODO: MC event selection?
203205
inputParticles.clear();
@@ -206,7 +208,7 @@ struct JetFinderTask {
206208
}
207209
PROCESS_SWITCH(JetFinderTask, processParticleLevelChargedJets, "Particle level charged jet finding", false);
208210

209-
void processParticleLevelChargedEvtWiseSubJets(aod::JetMcCollision const& collision, soa::Filtered<aod::JetParticlesSub> const& particles)
211+
void processParticleLevelChargedEvtWiseSubJets(soa::Filtered<aod::JetMcCollisions>::iterator const& collision, soa::Filtered<aod::JetParticlesSub> const& particles)
210212
{
211213
// TODO: MC event selection?
212214
inputParticles.clear();
@@ -215,7 +217,7 @@ struct JetFinderTask {
215217
}
216218
PROCESS_SWITCH(JetFinderTask, processParticleLevelChargedEvtWiseSubJets, "Particle level charged with event-wise constituent subtraction jet finding", false);
217219

218-
void processParticleLevelNeutralJets(aod::JetMcCollision const& collision, soa::Filtered<aod::JetParticles> const& particles)
220+
void processParticleLevelNeutralJets(soa::Filtered<aod::JetMcCollisions>::iterator const& collision, soa::Filtered<aod::JetParticles> const& particles)
219221
{
220222
// TODO: MC event selection?
221223
inputParticles.clear();
@@ -224,7 +226,7 @@ struct JetFinderTask {
224226
}
225227
PROCESS_SWITCH(JetFinderTask, processParticleLevelNeutralJets, "Particle level neutral jet finding", false);
226228

227-
void processParticleLevelFullJets(aod::JetMcCollision const& collision, soa::Filtered<aod::JetParticles> const& particles)
229+
void processParticleLevelFullJets(soa::Filtered<aod::JetMcCollisions>::iterator const& collision, soa::Filtered<aod::JetParticles> const& particles)
228230
{
229231
// TODO: MC event selection?
230232
inputParticles.clear();

‎PWGJE/JetFinders/jetFinderHF.cxx‎

Lines changed: 6 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -19,6 +19,7 @@
1919
#include "CommonConstants/PhysicsConstants.h"
2020

2121
#include "PWGJE/Core/JetFindingUtilities.h"
22+
#include "PWGJE/Core/JetDerivedDataUtilities.h"
2223
#include "Common/Core/RecoDecay.h"
2324

2425
using namespace o2;
@@ -45,6 +46,7 @@ struct JetFinderHFTask {
4546
Configurable<int> trackOccupancyInTimeRangeMax{"trackOccupancyInTimeRangeMax", 999999, "maximum occupancy of tracks in neighbouring collisions in a given time range"};
4647
Configurable<std::string> eventSelections{"eventSelections", "sel8", "choose event selection"};
4748
Configurable<std::string> triggerMasks{"triggerMasks", "", "possible JE Trigger masks: fJetChLowPt,fJetChHighPt,fTrackLowPt,fTrackHighPt,fJetD0ChLowPt,fJetD0ChHighPt,fJetLcChLowPt,fJetLcChHighPt,fEMCALReadout,fJetFullHighPt,fJetFullLowPt,fJetNeutralHighPt,fJetNeutralLowPt,fGammaVeryHighPtEMCAL,fGammaVeryHighPtDCAL,fGammaHighPtEMCAL,fGammaHighPtDCAL,fGammaLowPtEMCAL,fGammaLowPtDCAL,fGammaVeryLowPtEMCAL,fGammaVeryLowPtDCAL"};
49+
Configurable<bool> skipMBGapEvents{"skipMBGapEvents", true, "decide to run over MB gap events or not"};
4850

4951
// track level configurables
5052
Configurable<float> trackPtMin{"trackPtMin", 0.15, "minimum track pT"};
@@ -151,7 +153,8 @@ struct JetFinderHFTask {
151153
}
152154

153155
aod::EMCALClusterDefinition clusterDefinition = aod::emcalcluster::getClusterDefinitionFromString(clusterDefinitionS.value);
154-
Filter collisionFilter = (nabs(aod::jcollision::posZ) < vertexZCut && aod::jcollision::centrality >= centralityMin && aod::jcollision::centrality < centralityMax && aod::jcollision::trackOccupancyInTimeRange <= trackOccupancyInTimeRangeMax);
156+
Filter collisionFilter = (nabs(aod::jcollision::posZ) < vertexZCut && aod::jcollision::centrality >= centralityMin && aod::jcollision::centrality < centralityMax && aod::jcollision::trackOccupancyInTimeRange <= trackOccupancyInTimeRangeMax && (!skipMBGapEvents || aod::jcollision::subGeneratorId != static_cast<int>(jetderiveddatautilities::JCollisionSubGeneratorId::mbGap)));
157+
Filter mcCollisionFilter = (!skipMBGapEvents || aod::jmccollision::subGeneratorId != static_cast<int>(jetderiveddatautilities::JCollisionSubGeneratorId::mbGap)); // should we add a posZ vtx cut here or leave it to analysers?
155158
Filter trackCuts = (aod::jtrack::pt >= trackPtMin && aod::jtrack::pt < trackPtMax && aod::jtrack::eta >= trackEtaMin && aod::jtrack::eta <= trackEtaMax && aod::jtrack::phi >= trackPhiMin && aod::jtrack::phi <= trackPhiMax);
156159
Filter partCuts = (aod::jmcparticle::pt >= trackPtMin && aod::jmcparticle::pt < trackPtMax && aod::jmcparticle::eta >= trackEtaMin && aod::jmcparticle::eta <= trackEtaMax && aod::jmcparticle::phi >= trackPhiMin && aod::jmcparticle::phi <= trackPhiMax);
157160
Filter clusterFilter = (aod::jcluster::definition == static_cast<int>(clusterDefinition) && aod::jcluster::eta >= clusterEtaMin && aod::jcluster::eta <= clusterEtaMax && aod::jcluster::phi >= clusterPhiMin && aod::jcluster::phi <= clusterPhiMax && aod::jcluster::energy >= clusterEnergyMin && aod::jcluster::time > clusterTimeMin && aod::jcluster::time < clusterTimeMax && (clusterRejectExotics && aod::jcluster::isExotic != true));
@@ -251,7 +254,7 @@ struct JetFinderHFTask {
251254
}
252255
PROCESS_SWITCH(JetFinderHFTask, processChargedEvtWiseSubJetsMCD, "charged hf jet finding on MC detector level with event-wise constituent subtraction", false);
253256

254-
void processChargedJetsMCP(aod::JetMcCollision const& collision,
257+
void processChargedJetsMCP(soa::Filtered<aod::JetMcCollisions>::iterator const& collision,
255258
soa::Filtered<aod::JetParticles> const& particles,
256259
CandidateTableMCP const& candidates)
257260
{
@@ -261,7 +264,7 @@ struct JetFinderHFTask {
261264
}
262265
PROCESS_SWITCH(JetFinderHFTask, processChargedJetsMCP, "hf jet finding on MC particle level", false);
263266

264-
void processChargedEvtWiseSubJetsMCP(aod::JetMcCollision const& collision,
267+
void processChargedEvtWiseSubJetsMCP(soa::Filtered<aod::JetMcCollisions>::iterator const& collision,
265268
soa::Filtered<JetParticlesSubTable> const& particles,
266269
CandidateTableMCP const& candidates)
267270
{

‎PWGJE/JetFinders/jetFinderV0.cxx‎

Lines changed: 5 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -18,6 +18,7 @@
1818
#include "CommonConstants/PhysicsConstants.h"
1919

2020
#include "PWGJE/Core/JetFindingUtilities.h"
21+
#include "PWGJE/Core/JetDerivedDataUtilities.h"
2122
#include "Common/Core/RecoDecay.h"
2223

2324
using namespace o2;
@@ -43,6 +44,7 @@ struct JetFinderV0Task {
4344
Configurable<int> trackOccupancyInTimeRangeMax{"trackOccupancyInTimeRangeMax", 999999, "maximum occupancy of tracks in neighbouring collisions in a given time range"};
4445
Configurable<std::string> eventSelections{"eventSelections", "sel8", "choose event selection"};
4546
Configurable<std::string> triggerMasks{"triggerMasks", "", "possible JE Trigger masks: fJetChLowPt,fJetChHighPt,fTrackLowPt,fTrackHighPt,fJetD0ChLowPt,fJetD0ChHighPt,fJetLcChLowPt,fJetLcChHighPt,fEMCALReadout,fJetFullHighPt,fJetFullLowPt,fJetNeutralHighPt,fJetNeutralLowPt,fGammaVeryHighPtEMCAL,fGammaVeryHighPtDCAL,fGammaHighPtEMCAL,fGammaHighPtDCAL,fGammaLowPtEMCAL,fGammaLowPtDCAL,fGammaVeryLowPtEMCAL,fGammaVeryLowPtDCAL"};
47+
Configurable<bool> skipMBGapEvents{"skipMBGapEvents", true, "decide to run over MB gap events or not"};
4648

4749
// track level configurables
4850
Configurable<float> trackPtMin{"trackPtMin", 0.15, "minimum track pT"};
@@ -142,7 +144,8 @@ struct JetFinderV0Task {
142144
registry.add("hJetMCP", "sparse for mcp jets", {HistType::kTHnC, {{jetRadiiBins, ""}, {jetPtBinNumber, jetPtMinDouble, jetPtMaxDouble}, {40, -1.0, 1.0}, {18, 0.0, 7.0}}});
143145
}
144146

145-
Filter collisionFilter = (nabs(aod::jcollision::posZ) < vertexZCut && aod::jcollision::centrality >= centralityMin && aod::jcollision::centrality < centralityMax && aod::jcollision::trackOccupancyInTimeRange <= trackOccupancyInTimeRangeMax);
147+
Filter collisionFilter = (nabs(aod::jcollision::posZ) < vertexZCut && aod::jcollision::centrality >= centralityMin && aod::jcollision::centrality < centralityMax && aod::jcollision::trackOccupancyInTimeRange <= trackOccupancyInTimeRangeMax && (!skipMBGapEvents || aod::jcollision::subGeneratorId != static_cast<int>(jetderiveddatautilities::JCollisionSubGeneratorId::mbGap)));
148+
Filter mcCollisionFilter = (!skipMBGapEvents || aod::jmccollision::subGeneratorId != static_cast<int>(jetderiveddatautilities::JCollisionSubGeneratorId::mbGap)); // should we add a posZ vtx cut here or leave it to analysers?
146149
Filter trackCuts = (aod::jtrack::pt >= trackPtMin && aod::jtrack::pt < trackPtMax && aod::jtrack::eta >= trackEtaMin && aod::jtrack::eta <= trackEtaMax && aod::jtrack::phi >= trackPhiMin && aod::jtrack::phi <= trackPhiMax);
147150
Filter partCuts = (aod::jmcparticle::pt >= trackPtMin && aod::jmcparticle::pt < trackPtMax && aod::jmcparticle::eta >= trackEtaMin && aod::jmcparticle::eta <= trackEtaMax && aod::jmcparticle::phi >= trackPhiMin && aod::jmcparticle::phi <= trackPhiMax);
148151
// Filter candidateCuts = (aod::hfcand::pt >= candPtMin && aod::hfcand::pt < candPtMax && aod::hfcand::y >= candYMin && aod::hfcand::y < candYMax);
@@ -200,7 +203,7 @@ struct JetFinderV0Task {
200203
}
201204
PROCESS_SWITCH(JetFinderV0Task, processChargedJetsMCD, "charged hf jet finding on MC detector level", false);
202205

203-
void processChargedJetsMCP(aod::JetMcCollision const& collision,
206+
void processChargedJetsMCP(soa::Filtered<aod::JetMcCollisions>::iterator const& collision,
204207
soa::Filtered<aod::JetParticles> const& particles,
205208
CandidateTableMCP const& candidates)
206209
{

‎PWGJE/TableProducer/derivedDataProducer.cxx‎

Lines changed: 10 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -64,6 +64,7 @@ struct JetDerivedDataProducerTask {
6464
Produces<aod::JBCs> jBCsTable;
6565
Produces<aod::JBCPIs> jBCParentIndexTable;
6666
Produces<aod::JCollisions> jCollisionsTable;
67+
Produces<aod::JCollisionMcInfos> jCollisionMcInfosTable;
6768
Produces<aod::JCollisionPIs> jCollisionsParentIndexTable;
6869
Produces<aod::JCollisionBCs> jCollisionsBunchCrossingIndexTable;
6970
Produces<aod::JEMCCollisionLbs> jCollisionsEMCalLabelTable;
@@ -146,6 +147,9 @@ struct JetDerivedDataProducerTask {
146147
void processClearMaps(aod::Collisions const&)
147148
{
148149
trackCollisionMapping.clear();
150+
if (!doprocessMcCollisionLabels) {
151+
jCollisionMcInfosTable(-1.0, jetderiveddatautilities::JCollisionSubGeneratorId::none); // fill a dummy weights table if not MC
152+
}
149153
}
150154
PROCESS_SWITCH(JetDerivedDataProducerTask, processClearMaps, "clears all maps", true);
151155

@@ -212,21 +216,22 @@ struct JetDerivedDataProducerTask {
212216
}
213217
PROCESS_SWITCH(JetDerivedDataProducerTask, processEMCalCollisionLabels, "produces derived collision labels for EMCal", false);
214218

215-
void processMcCollisionLabels(soa::Join<aod::Collisions, aod::McCollisionLabels>::iterator const& collision)
219+
void processMcCollisionLabels(soa::Join<aod::Collisions, aod::McCollisionLabels>::iterator const& collision, aod::McCollisions const&)
216220
{
217-
218221
if (collision.has_mcCollision()) {
219222
jMcCollisionsLabelTable(collision.mcCollisionId());
223+
jCollisionMcInfosTable(collision.mcCollision().weight(), collision.mcCollision().getSubGeneratorId());
220224
} else {
221225
jMcCollisionsLabelTable(-1);
226+
jCollisionMcInfosTable(-1.0, jetderiveddatautilities::JCollisionSubGeneratorId::none);
222227
}
223228
}
224229
PROCESS_SWITCH(JetDerivedDataProducerTask, processMcCollisionLabels, "produces derived MC collision labels table", false);
225230

226-
void processMcCollisions(aod::McCollision const& McCollision)
231+
void processMcCollisions(aod::McCollision const& mcCollision)
227232
{
228-
jMcCollisionsTable(McCollision.posX(), McCollision.posY(), McCollision.posZ(), McCollision.weight());
229-
jMcCollisionsParentIndexTable(McCollision.globalIndex());
233+
jMcCollisionsTable(mcCollision.posX(), mcCollision.posY(), mcCollision.posZ(), mcCollision.weight(), mcCollision.getSubGeneratorId());
234+
jMcCollisionsParentIndexTable(mcCollision.globalIndex());
230235
}
231236
PROCESS_SWITCH(JetDerivedDataProducerTask, processMcCollisions, "produces derived MC collision table", false);
232237

0 commit comments

Comments
 (0)