Draft: Faster tracker reconstruction turn-around by using truth init params

benchmarks/tracking/options/tracker_reconstruction.py

@@ -47,7 +47,12 @@ from Configurables import Jug__Reco__SimpleClustering as SimpleClustering
 
 
 podioinput = PodioInput("PodioReader", 
-                        collections=["mcparticles","TrackerEndcapHits","TrackerBarrelHits","VertexBarrelHits","VertexEndcapHits","EcalBarrelHits"])#, OutputLevel=DEBUG)
+                        collections=["mcparticles","TrackerEndcapHits","TrackerBarrelHits","VertexBarrelHits","VertexEndcapHits"])#, OutputLevel=DEBUG)
+
+## Track param init
+truth_trk_init = TrackParamTruthInit("truth_trk_init",
+        inputMCParticles="mcparticles",
+        outputInitialTrackParameters="InitTrackParamsFromTruth")
 
 ## copiers to get around input --> output copy bug. Note the "2" appended to the output collection.
 copier = MCCopier("MCCopier", 
@@ -57,10 +62,6 @@ trkcopier = TrkCopier("TrkCopier",
         inputCollection="TrackerBarrelHits", 
         outputCollection="TrackerBarrelHits2") 
 
-ecal_digi = EMCalorimeterDigi("ecal_digi", 
-        inputHitCollection="EcalBarrelHits", 
-        outputHitCollection="RawEcalBarrelHits")
-
 trk_b_digi = TrackerDigi("trk_b_digi", 
         inputHitCollection="TrackerBarrelHits",
         outputHitCollection="TrackerBarrelRawHits",
@@ -80,19 +81,6 @@ vtx_ec_digi = TrackerDigi("vtx_ec_digi",
         outputHitCollection="VertexEndcapRawHits",
         timeResolution=8)
 
-ecal_reco = EMCalReconstruction("ecal_reco", 
-        inputHitCollection="RawEcalBarrelHits", 
-        outputHitCollection="RecEcalBarrelHits",
-        minModuleEdep=0.0*units.MeV,
-        OutputLevel=WARNING)
-
-simple_cluster = SimpleClustering("simple_cluster", 
-        inputHitCollection="RecEcalBarrelHits", 
-        outputClusters="SimpleClusters",
-        minModuleEdep=1.0*units.MeV,
-        maxDistance=50.0*units.cm,
-        OutputLevel=WARNING)
-
 # Tracker and vertex reconstruction
 trk_b_reco = TrackerHitReconstruction("trk_b_reco",
         inputHitCollection = trk_b_digi.outputHitCollection,
@@ -147,7 +135,7 @@ clust_trk_init = TrackParamClusterInit("clust_trk_init",
 trk_find_alg = TrackFindingAlgorithm("trk_find_alg",
         inputSourceLinks = sourcelinker.outputSourceLinks,
         inputMeasurements = sourcelinker.outputMeasurements,
-        inputInitialTrackParameters= "InitTrackParams",#"InitTrackParamsFromClusters", 
+        inputInitialTrackParameters= "InitTrackParamsFromTruth",
         outputTrajectories="trajectories",
         OutputLevel=DEBUG)
 
@@ -214,10 +202,8 @@ out.outputCommands = ["keep *",
 ApplicationMgr(
     TopAlg = [podioinput, 
               copier, trkcopier,
-              ecal_digi, 
               trk_b_digi, trk_ec_digi, vtx_b_digi, vtx_ec_digi, 
               trk_b_reco, trk_ec_reco, vtx_b_reco, vtx_ec_reco, 
-              ecal_reco, simple_cluster,
               sourcelinker, #trk_hits_srclnkr,
               clust_trk_init, truth_trk_init, #vtxcluster_trk_init, 
               trk_find_alg, parts_from_fit,