Skip to content
Snippets Groups Projects
Commit b11372ea authored by Whitney Armstrong's avatar Whitney Armstrong
Browse files

modified: options/tracker_reconstruction.py

parent 56e69f0e
Branches
Tags
1 merge request!100Improve tracking setup
...@@ -26,7 +26,7 @@ from Configurables import Jug__Base__InputCopier_dd4pod__CalorimeterHitCollectio ...@@ -26,7 +26,7 @@ from Configurables import Jug__Base__InputCopier_dd4pod__CalorimeterHitCollectio
from Configurables import Jug__Base__InputCopier_dd4pod__TrackerHitCollection_dd4pod__TrackerHitCollection_ as TrkCopier from Configurables import Jug__Base__InputCopier_dd4pod__TrackerHitCollection_dd4pod__TrackerHitCollection_ as TrkCopier
from Configurables import Jug__Digi__ExampleCaloDigi as ExampleCaloDigi from Configurables import Jug__Digi__ExampleCaloDigi as ExampleCaloDigi
from Configurables import Jug__Digi__UFSDTrackerDigi as UFSDTrackerDigi from Configurables import Jug__Digi__UFSDTrackerDigi as TrackerDigi
from Configurables import Jug__Digi__EMCalorimeterDigi as EMCalorimeterDigi from Configurables import Jug__Digi__EMCalorimeterDigi as EMCalorimeterDigi
from Configurables import Jug__Reco__TrackerHitReconstruction as TrackerHitReconstruction from Configurables import Jug__Reco__TrackerHitReconstruction as TrackerHitReconstruction
...@@ -47,7 +47,7 @@ from Configurables import Jug__Reco__SimpleClustering as SimpleClustering ...@@ -47,7 +47,7 @@ from Configurables import Jug__Reco__SimpleClustering as SimpleClustering
podioinput = PodioInput("PodioReader", podioinput = PodioInput("PodioReader",
collections=["mcparticles","TrackerEndcapHits","TrackerBarrelHits","EcalBarrelHits"])#, OutputLevel=DEBUG)"SiVertexBarrelHits", collections=["mcparticles","TrackerEndcapHits","TrackerBarrelHits","VertexBarrelHits","VertexEndcapHits","EcalBarrelHits"])#, OutputLevel=DEBUG)
## copiers to get around input --> output copy bug. Note the "2" appended to the output collection. ## copiers to get around input --> output copy bug. Note the "2" appended to the output collection.
copier = MCCopier("MCCopier", copier = MCCopier("MCCopier",
...@@ -61,20 +61,23 @@ ecal_digi = EMCalorimeterDigi("ecal_digi", ...@@ -61,20 +61,23 @@ ecal_digi = EMCalorimeterDigi("ecal_digi",
inputHitCollection="EcalBarrelHits", inputHitCollection="EcalBarrelHits",
outputHitCollection="RawEcalBarrelHits") outputHitCollection="RawEcalBarrelHits")
ufsd_digi = UFSDTrackerDigi("ufsd_digi", trk_b_digi = TrackerDigi("trk_b_digi",
inputHitCollection="TrackerBarrelHits", inputHitCollection="TrackerBarrelHits",
outputHitCollection="TrackerBarrelRawHits", outputHitCollection="TrackerBarrelRawHits",
timeResolution=8) timeResolution=8)
ufsd_digi2 = UFSDTrackerDigi("ufsd_digi2", trk_ec_digi = TrackerDigi("trk_ec_digi",
inputHitCollection="TrackerEndcapHits", inputHitCollection="TrackerEndcapHits",
outputHitCollection="TrackerEndcapRawHits", outputHitCollection="TrackerEndcapRawHits",
timeResolution=8) timeResolution=8)
#vtx_digi = UFSDTrackerDigi("vtx_digi", vtx_b_digi = UFSDTrackerDigi("vtx_b_digi",
# inputHitCollection="SiVertexBarrelHits", inputHitCollection="VertexBarrelHits",
# outputHitCollection="SiVertexBarrelRawHits", outputHitCollection="VertexBarrelRawHits",
# timeResolution=8) timeResolution=8)
vtx_ec_digi = UFSDTrackerDigi("vtx_ec_digi",
inputHitCollection="VertexEndcapHits",
outputHitCollection="VertexEndcapRawHits",
timeResolution=8)
ecal_reco = EMCalReconstruction("ecal_reco", ecal_reco = EMCalReconstruction("ecal_reco",
inputHitCollection="RawEcalBarrelHits", inputHitCollection="RawEcalBarrelHits",
...@@ -90,32 +93,36 @@ simple_cluster = SimpleClustering("simple_cluster", ...@@ -90,32 +93,36 @@ simple_cluster = SimpleClustering("simple_cluster",
OutputLevel=WARNING) OutputLevel=WARNING)
# Tracker and vertex reconstruction # Tracker and vertex reconstruction
trk_barrel_reco = TrackerHitReconstruction("trk_barrel_reco", trk_b_reco = TrackerHitReconstruction("trk_b_reco",
inputHitCollection = ufsd_digi.outputHitCollection, inputHitCollection = trk_b_digi.outputHitCollection,
outputHitCollection="TrackerBarrelRecHits") outputHitCollection="TrackerBarrelRecHits")
trk_endcap_reco = TrackerHitReconstruction("trk_endcap_reco", trk_ec_reco = TrackerHitReconstruction("trk_ec_reco",
inputHitCollection = ufsd_digi2.outputHitCollection, inputHitCollection = trk_ec_digi.outputHitCollection,
outputHitCollection="TrackerEndcapRecHits") outputHitCollection="TrackerEndcapRecHits")
#vtx_barrel_reco = TrackerHitReconstruction("vtx_barrel_reco", vtx_b_reco = TrackerHitReconstruction("vtx_b_digi",
# inputHitCollection = vtx_digi.outputHitCollection, inputHitCollection = vtx_b_digi.outputHitCollection,
# outputHitCollection="VertexBarrelRecHits") outputHitCollection="VertexBarrelRecHits")
vtx_ec_reco = TrackerHitReconstruction("vtx_ec_digi",
inputHitCollection = vtx_ec_digi.outputHitCollection,
outputHitCollection="VertexEndcapRecHits")
# Hit Source linker # Hit Source linker
sourcelinker = TrackerSourceLinker("sourcelinker", sourcelinker = TrackerSourceLinker("sourcelinker",
inputHitCollection="TrackerBarrelRecHits", inputHitCollection=trk_b_reco.outputHitCollection,
outputSourceLinks="BarrelTrackSourceLinks", outputSourceLinks="BarrelTrackSourceLinks",
outputMeasurements="BarrelTrackMeasurements", outputMeasurements="BarrelTrackMeasurements",
OutputLevel=DEBUG) OutputLevel=DEBUG)
#trk_hits_srclnkr = TrackerSourcesLinker("trk_srcslnkr", trk_hits_srclnkr = TrackerSourcesLinker("trk_srcslnkr",
# ITrackerBarrelHits = "TrackerBarrelRecHits", ITrackerBarrelHits = vtx_b_reco.outputHitCollection,
# ITrackerEndcapHits = "TrackerEndcapRecHits", ITrackerEndcapHits = vtx_ec_reco.outputHitCollection,
# OTrackerBarrelHits = "TrackerBarrelRecHits", OTrackerBarrelHits = trk_b_reco.outputHitCollection,
# OTrackerEndcapHits = "TrackerEndcapRecHits", OTrackerEndcapHits = trk_ec_reco.outputHitCollection,
# outputSourceLinks="TrackerMeasurements", outputSourceLinks="TrackerMeasurements",
# OutputLevel=DEBUG) OutputLevel=DEBUG)
## Track param init ## Track param init
truth_trk_init = TrackParamTruthInit("truth_trk_init", truth_trk_init = TrackParamTruthInit("truth_trk_init",
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment