@@ -152,6 +152,7 @@ struct StudyPnch {
152152 histos.add (" hMultiplicityMCrecMod" , " hMultiplicityMCrecMod" , kTH1F , {axisMult}, true );
153153 histos.add (" hMultiplicityMCgenMod" , " hMultiplicityMCgenMod" , kTH1F , {axisMult}, true );
154154 histos.add (" hResponseMatrixMod" , " hResponseMatrixMod" , kTH2F , {axisMult, axisMult}, true );
155+ histos.add (" hCountNTracksMod" , " hCountNTracksMod" , kTH1F , {axisCountNumberTracks}, true );
155156 }
156157 if (doprocessEvtLossSigLossMC) {
157158 histos.add (" MCEventHist" , " MCEventHist" , kTH1F , {axisEvent}, false );
@@ -255,7 +256,6 @@ struct StudyPnch {
255256 if (!isGenTrackSelected (track)) {
256257 continue ;
257258 }
258- // Verify that the track belongs to the given MC collision
259259 if (track.mcCollisionId () != McCol.globalIndex ()) {
260260 continue ;
261261 }
@@ -275,7 +275,6 @@ struct StudyPnch {
275275 if (!isTrackSelected (track)) {
276276 continue ;
277277 }
278- // Verify that the track belongs to the given MC collision
279278 if (track.has_mcParticle ()) {
280279 auto particle = track.mcParticle ();
281280 if (particle.mcCollisionId () != McCol.mcCollisionId ()) {
@@ -286,7 +285,7 @@ struct StudyPnch {
286285 }
287286 histos.fill (HIST (" PhiVsEtaHist" ), track.phi (), track.eta ());
288287 }
289- // Once all the frequencies have been counted, a loop can be made to fill the histogram
288+ // Loop to fill the histogram without cloned tracks
290289 for (const auto & [globalIndex, frequency] : recoFrequencies) {
291290 histos.fill (HIST (" hCountNTracks" ), frequency);
292291 // Fill histogram with not cloned tracks
@@ -354,6 +353,7 @@ struct StudyPnch {
354353 if (!isEventSelected (RecCol)) {
355354 continue ;
356355 }
356+ // Evaluation of reconstructed collisions with more than 1 contributor
357357 if (RecCol.globalIndex () != mcCollision.bestCollisionIndex ()) {
358358 continue ;
359359 }
0 commit comments