Skip to content
Snippets Groups Projects

Resolve "Move sampling detector benchmarks to reconstruction_benchmarks"

Files
6
@@ -37,13 +37,6 @@ from Configurables import PodioInput
from Configurables import Jug__Base__InputCopier_dd4pod__Geant4ParticleCollection_dd4pod__Geant4ParticleCollection_ as MCCopier
from Configurables import Jug__Base__InputCopier_dd4pod__CalorimeterHitCollection_dd4pod__CalorimeterHitCollection_ as CalCopier
from Configurables import Jug__Digi__EcalTungstenSamplingDigi as EcalTungstenSamplingDigi
from Configurables import Jug__Reco__EcalTungstenSamplingReco as EcalTungstenSamplingReco
from Configurables import Jug__Reco__SamplingECalHitsMerger as SamplingECalHitsMerger
from Configurables import Jug__Reco__CalorimeterIslandCluster as IslandCluster
from Configurables import Jug__Reco__ClusterRecoCoG as RecoCoG
podioinput = PodioInput("PodioReader", collections=["mcparticles","EcalBarrelAstroPixHits"], OutputLevel=DEBUG)
# Thrown Information
@@ -56,55 +49,13 @@ embarrelcopier = CalCopier("CalBarrelCopier",
inputCollection="EcalBarrelAstroPixHits",
outputCollection="EcalBarrelAstroPixHits2",
OutputLevel=DEBUG)
# Digitization
embarreldigi = EcalTungstenSamplingDigi("ecal_barrel_digi",
inputHitCollection="EcalBarrelAstroPixHits",
outputHitCollection="RawEcalBarrelAstroPixHits",
inputEnergyUnit=units.GeV,
inputTimeUnit=units.ns,
dynamicRangeADC=700*units.keV,
energyResolutions=[0., 0.02, 0.],
pedestalSigma=40,
OutputLevel=DEBUG)
# Reconstruction
embarrelreco = EcalTungstenSamplingReco("ecal_barrel_reco",
inputHitCollection="RawEcalBarrelAstroPixHits",
outputHitCollection="RecoEcalBarrelAstroPixHits",
dynamicRangeADC=700*units.keV,
pedestalSigma=40,
OutputLevel=DEBUG)
# 2D+1 Clusterings
# readout id definition for barrel ecal
# <id>system:8,barrel:3,module:4,layer:10,slice:5,x:32:-16,y:-16</id>
# xy_merger sum layers/slices, masking (8+3+4, 8+3+4+5+10-1)
embarrelxymerger = SamplingECalHitsMerger("ecal_barrel_xy_merger",
cellIDMaskRanges=[(15, 29)],
inputHitCollection="RecoEcalBarrelAstroPixHits",
outputHitCollection="RecoEcalBarrelAstroPixHitsXY")
# xy_merger sum modules, masking (8+3+4+5+10, 8+3+4+5+10+32-1)
embarrelzmerger = SamplingECalHitsMerger("ecal_barrel_z_merger",
cellIDMaskRanges=[(30, 61)],
inputHitCollection="RecoEcalBarrelAstroPixHits",
outputHitCollection="RecoEcalBarrelAstroPixHitsZ")
# Clustering
embarrelcluster = IslandCluster("ecal_barrel_cluster",
inputHitCollection="RecoEcalBarrelAstroPixHitsXY",
outputClusterCollection="EcalBarrelAstroPixClusters",
minClusterCenterEdep=0.5*units.MeV,
splitCluster=False,
groupRanges=[2.0*units.cm, 2.0*units.cm, 2.0*units.cm])
# Reconstruct the cluster with Center of Gravity method
embarrelclusterreco = RecoCoG("ecal_barrel_clusterreco",
clusterCollection="EcalBarrelAstroPixClusters",
logWeightBase=6.2)
out = PodioOutput("out", filename=output_rec_file)
out.outputCommands = ["keep *"]
ApplicationMgr(
TopAlg = [podioinput, copier, embarrelcopier, embarreldigi,
embarrelreco, embarrelxymerger, embarrelzmerger, embarrelcluster, embarrelclusterreco, out],
TopAlg = [podioinput, copier, embarrelcopier, out],
EvtSel = 'NONE',
EvtMax = n_events,
ExtSvc = [podioevent],
Loading