Draft: Resolve "Update pi0 reconstruction"

Closes #68

benchmarks/imaging_ecal/options/imaging_topocluster.py

@@ -41,11 +41,9 @@ out = PodioOutput("out", filename=kwargs['output'])
 podioinput = PodioInput("PodioReader", collections=["mcparticles", "EcalBarrelHits"], OutputLevel=DEBUG)
 
 copier = MCCopier("MCCopier",
-        # OutputLevel=DEBUG,
         inputCollection="mcparticles",
         outputCollection="mcparticles2")
 calcopier = CalCopier("CalCopier",
-        # OutputLevel=DEBUG,
         inputCollection="EcalBarrelHits",
         outputCollection="EcalBarrelHits2")
 
@@ -56,31 +54,30 @@ daq_setting = dict(
         pedestalMean=400,
         pedestalSigma=50)   # 50/32767*3 MeV ~ 5 keV
 imcaldigi = CalorimeterHitDigi("imcal_digi",
-        # OutputLevel=DEBUG,
         inputHitCollection="EcalBarrelHits",
         outputHitCollection="DigiEcalBarrelImagingHits",
         energyResolutions=[0., 0.02, 0.],
         **daq_setting)
 imcalreco = ImagingPixelReco("imcal_reco",
-        # OutputLevel=DEBUG,
         inputHitCollection=imcaldigi.outputHitCollection,
         outputHitCollection="RecoEcalBarrelImagingHits",
         readoutClass="EcalBarrelHits",
         layerField="layer",
         sectorField="module",
+        thresholdFactor=5.0,
         **daq_setting)
 
 imcalcluster = ImagingTopoCluster("imcal_cluster",
-        # OutputLevel=DEBUG,
         inputHitCollection=imcalreco.outputHitCollection,
         outputProtoClusterCollection="EcalBarrelImagingProtoClusters",
         localDistXY=[2.*units.mm, 2*units.mm],
         layerDistEtaPhi=[10*units.mrad, 10*units.mrad],
         neighbourLayersRange=2,
         sectorDist=3.*units.cm,
-        minClusterNhits=5)
+        minClusterNhits=15,
+        minClusterEdep=0.5*MeV)
+
 clusterreco = ImagingClusterReco("imcal_clreco",
-        # OutputLevel=DEBUG,
         inputHitCollection=imcalcluster.inputHitCollection,
         inputProtoClusterCollection=imcalcluster.outputProtoClusterCollection,
         outputLayerCollection="EcalBarrelImagingClustersLayers",