Draft: Resolve "Simple clustering: double free or corruption"

Closes #38 (closed)

benchmarks/full/options/full_reconstruction.py

@@ -460,9 +460,9 @@ truth_trk_init = TrackParamTruthInit("truth_trk_init",
         inputMCParticles="mcparticles",
         outputInitialTrackParameters="InitTrackParamsFromTruth")
 
-#clust_trk_init = TrackParamClusterInit("clust_trk_init",
-#        inputClusters="SimpleClusters",
-#        outputInitialTrackParameters="InitTrackParamsFromClusters")
+clust_trk_init = TrackParamClusterInit("clust_trk_init",
+        inputClusters="EcalBarrelClusters",
+        outputInitialTrackParameters="InitTrackParamsFromClusters")
 
 #vtxcluster_trk_init = TrackParamVertexClusterInit("vtxcluster_trk_init",
 #        inputVertexHits="VertexBarrelRecHits",
@@ -482,28 +482,28 @@ parts_from_fit = ParticlesFromTrackFit("parts_from_fit",
         outputParticles = "ReconstructedParticlesInitFromTruth",
         outputTrackParameters = "outputTrackParameters")
 
-#trk_find_alg1 = TrackFindingAlgorithm("trk_find_alg1",
-#        inputSourceLinks = sourcelinker.outputSourceLinks,
-#        inputMeasurements = sourcelinker.outputMeasurements,
-#        inputInitialTrackParameters= "InitTrackParamsFromClusters",
-#        outputTrajectories="trajectories1")
+trk_find_alg1 = TrackFindingAlgorithm("trk_find_alg1",
+        inputSourceLinks = sourcelinker.outputSourceLinks,
+        inputMeasurements = sourcelinker.outputMeasurements,
+        inputInitialTrackParameters= "InitTrackParamsFromClusters",
+        outputTrajectories="trajectories1")
 
-#parts_from_fit1 = ParticlesFromTrackFit("parts_from_fit1",
-#        inputTrajectories="trajectories1",
-#        outputParticles="ReconstructedParticlesInitFromClusters",
-#        outputTrackParameters="outputTrackParameters1")
+parts_from_fit1 = ParticlesFromTrackFit("parts_from_fit1",
+        inputTrajectories="trajectories1",
+        outputParticles="ReconstructedParticlesInitFromClusters",
+        outputTrackParameters="outputTrackParameters1")
 
 #trk_find_alg2 = TrackFindingAlgorithm("trk_find_alg2",
 #        inputSourceLinks = trk_hits_srclnkr.outputSourceLinks,
 #        inputMeasurements = trk_hits_srclnkr.outputMeasurements,
 #        inputInitialTrackParameters= "InitTrackParamsFromVtxClusters",
 #        outputTrajectories="trajectories2")
+
 #parts_from_fit2 = ParticlesFromTrackFit("parts_from_fit2",
 #        inputTrajectories="trajectories2",
-#        outputParticles="ReconstructedParticles2",
+#        outputParticles="ReconstructedParticlesInitFromVertex",
 #        outputTrackParameters="outputTrackParameters2")
 
-
 pmtdigi = PhotoMultiplierDigi("pmtdigi",
         inputHitCollection="ForwardRICHHits",
         outputHitCollection="ForwardRICHHitsDigi",
@@ -532,7 +532,7 @@ podout.outputCommands = ["keep *",
 ApplicationMgr(
     TopAlg = [podin, copier, trkcopier,
               dummy,
-              ecal_digi, ecal_reco, #simple_cluster,
+              #ecal_digi, #ecal_reco, #simple_cluster,
               ce_ecal_digi, ce_ecal_reco, ce_ecal_cl, ce_ecal_clreco,
               ci_ecal_digi, ci_ecal_reco, ci_ecal_merger, ci_ecal_cl, ci_ecal_clreco,
               cb_ecal_digi, cb_ecal_reco, cb_ecal_cl, cb_ecal_clreco,
@@ -543,11 +543,11 @@ ApplicationMgr(
               trk_b_digi, trk_ec_digi, vtx_b_digi, vtx_ec_digi,
               trk_b_reco, trk_ec_reco, vtx_b_reco, vtx_ec_reco,
               sourcelinker, #trk_hits_srclnkr,
-              #clust_trk_init,
+              clust_trk_init,
               truth_trk_init,
               #vtxcluster_trk_init,
               trk_find_alg, parts_from_fit,
-              #trk_find_alg1, parts_from_fit1,
+              trk_find_alg1, parts_from_fit1,
               #trk_find_alg2, parts_from_fit2,
               pmtcopier, pmtdigi, pmtreco, #richcluster,
               podout