From f232feafddb7d2227949056f460631c0fcb43bf4 Mon Sep 17 00:00:00 2001
From: Wouter Deconinck <wouter.deconinck@umanitoba.ca>
Date: Tue, 5 Oct 2021 21:53:29 +0000
Subject: [PATCH] Resolve "Add neutrals to options file"

---
 options/reconstruction.py | 23 +++++++++++++++++++++--
 1 file changed, 21 insertions(+), 2 deletions(-)

diff --git a/options/reconstruction.py b/options/reconstruction.py
index 90b32b89..e716d514 100644
--- a/options/reconstruction.py
+++ b/options/reconstruction.py
@@ -59,6 +59,7 @@ from Configurables import Jug__Base__InputCopier_dd4pod__PhotoMultiplierHitColle
 from Configurables import Jug__Fast__MC2SmearedParticle as MC2DummyParticle
 from Configurables import Jug__Fast__ParticlesWithTruthPID as ParticlesWithTruthPID
 from Configurables import Jug__Fast__SmearedFarForwardParticles as SmearedFarForwardParticles
+from Configurables import Jug__Fast__MatchClusters as MatchClusters
 
 from Configurables import Jug__Digi__PhotoMultiplierDigi as PhotoMultiplierDigi
 from Configurables import Jug__Digi__CalorimeterHitDigi as CalHitDigi
@@ -543,16 +544,34 @@ algorithms.append(parts_from_fit)
 trajs_from_fit = TrajectoryFromTrackFit("trajs_from_fit",
         inputTrajectories = trk_find_alg.outputTrajectories,
         outputTrajectoryParameters = "outputTrajectoryParameters")
-        #OutputLevel=DEBUG)
 algorithms.append(trajs_from_fit)
 
 # Event building
 parts_with_truth_pid = ParticlesWithTruthPID("parts_with_truth_pid",
         inputMCParticles = "mcparticles",
         inputTrackParameters = parts_from_fit.outputTrackParameters,
+        outputParticles = "ReconstructedChargedParticles",
+        outputRelations = "ReconstructedChargedParticleRelations")
+algorithms.append(parts_with_truth_pid)
+
+match_clusters = MatchClusters("match_clusters",
+        inputMCParticles = "mcparticles",
+        inputParticles = parts_with_truth_pid.outputParticles,
+        inputRelations = parts_with_truth_pid.outputRelations,
+        inputEcalClusters = [
+                str(ce_ecal_clreco.outputClusterCollection),
+                str(img_barrel_clreco.outputClusterCollection),
+                str(scfi_barrel_clreco.outputClusterCollection),
+                str(ci_ecal_clreco.outputClusterCollection)
+        ],
+        inputHcalClusters = [
+                str(ce_hcal_clreco.outputClusterCollection),
+                str(cb_hcal_clreco.outputClusterCollection),
+                str(ci_hcal_clreco.outputClusterCollection)
+        ],
         outputParticles = "ReconstructedParticles",
         outputRelations = "ReconstructedParticleRelations")
-algorithms.append(parts_with_truth_pid)
+algorithms.append(match_clusters)
 
 # DRICH
 drich_digi = PhotoMultiplierDigi("drich_digi",
-- 
GitLab