Resolve "Move sampling detector benchmarks to reconstruction_benchmarks"

Closes #12 (closed)

calorimeters/options/emcal_barrel_reco.py

@@ -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],