diff --git a/options/reconstruction.py b/options/reconstruction.py index 90b32b89d6e08decc6fe7cccb9c2e415aabf681e..e716d514632a4abca22da46cc39eb35343d92990 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",