Skip to content
Snippets Groups Projects

Draft: Resolve "Update pi0 reconstruction"

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