Commit 3c82469e authored by Chao Peng's avatar Chao Peng
Browse files

adjust analysis script and fix typos

parent 6c003ef3
Pipeline #8039 passed with stages
in 8 minutes and 44 seconds
......@@ -12,7 +12,7 @@ from Configurables import Jug__Reco__CalorimeterIslandCluster as IslandCluster
from Configurables import Jug__Reco__ClusterRecoCoG as RecoCoG
geo_service = GeoSvc("GeoSvc", detectors=["../topside/test.xml"])
podioevent = EICDataSvc("EventDataSvc", inputs=["barrel_electrons.root"], OutputLevel=DEBUG)
podioevent = EICDataSvc("EventDataSvc", inputs=["../topside/test.root"], OutputLevel=DEBUG)
podioinput = PodioInput("PodioReader", collections=["mcparticles", "EcalBarrelHits"], OutputLevel=DEBUG)
emcaldigi = EcalTungstenSamplingDigi("ecal_digi",
......@@ -20,28 +20,32 @@ emcaldigi = EcalTungstenSamplingDigi("ecal_digi",
outputHitCollection="DigiEcalBarrelHits",
inputEnergyUnit=units.GeV,
inputTimeUnit=units.ns,
energyResolutions=[0., 0.02, 0.],
dynamicRangeADC=700*units.keV,
pedestalSigma=50,
OutputLevel=DEBUG)
emcalreco = EcalTungstenSamplingReco("ecal_reco",
inputHitCollection="DigiEcalBarrelHits",
outputHitCollection="RecoEcalBarrelHits",
dynamicRangeADC=700*units.keV,
OutputLevel=DEBUG)
# readout id definition for barrel ecal
# <id>system:8,barrel:3,module:4,layer:6,slice:5,x:32:-16,y:-16</id>
# xy_merger sum layers/slices, masking (8+3+4, 8+3+4+5+6-1)
# <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)
xymerger = SamplingECalHitsMerger("ecal_xy_merger",
cellIDMaskRanges=[(15, 25)],
cellIDMaskRanges=[(15, 29)],
inputHitCollection="RecoEcalBarrelHits",
outputHitCollection="RecoEcalBarrelHitsXY")
# xy_merger sum modules, masking (8+3+4+5+6, 8+3+4+5+6+32-1)
# xy_merger sum modules, masking (8+3+4+5+10, 8+3+4+5+10+32-1)
zmerger = SamplingECalHitsMerger("ecal_z_merger",
cellIDMaskRanges=[(26, 57)],
cellIDMaskRanges=[(30, 61)],
inputHitCollection="RecoEcalBarrelHits",
outputHitCollection="RecoEcalBarrelHitsZ")
emcalcluster = IslandCluster(inputHitCollection="RecoEcalBarrelHitsXY",
outputClusterCollection="EcalBarrelClusters",
minClusterCenterEdep=5.0*units.MeV,
minClusterCenterEdep=0.5*units.MeV,
splitCluster=False,
groupRange=5.0)
groupRange=10.0)
clusterreco = RecoCoG(clusterCollection="EcalBarrelClusters", logWeightBase=6.2, OutputLevel=DEBUG)
......
......@@ -33,7 +33,7 @@ namespace Jug {
Gaudi::Property<double> m_eUnit{this, "inputEnergyUnit", GeV};
Gaudi::Property<double> m_tUnit{this, "inputTimeUnit", ns};
Gaudi::Property<int> m_capADC{this, "capacityADC", 8096};
Gaudi::Property<double> m_dyRangeADC{this, "DynamicRangeADC", 100*MeV};
Gaudi::Property<double> m_dyRangeADC{this, "dynamicRangeADC", 100*MeV};
Gaudi::Property<int> m_pedMeanADC{this, "pedestalMean", 400};
Gaudi::Property<double> m_pedSigmaADC{this, "pedestalSigma", 3.2};
Rndm::Numbers m_normDist;
......
......@@ -31,7 +31,7 @@ namespace Jug::Reco {
public:
Gaudi::Property<double> m_lUnit{this, "lengthUnit", dd4hep::mm};
Gaudi::Property<int> m_capADC{this, "capacityADC", 8096};
Gaudi::Property<double> m_dyRangeADC{this, "DynamicRangeADC", 100*MeV};
Gaudi::Property<double> m_dyRangeADC{this, "dynamicRangeADC", 100*MeV};
Gaudi::Property<int> m_pedMeanADC{this, "pedestalMean", 400};
Gaudi::Property<double> m_pedSigmaADC{this, "pedestalSigma", 3.2};
Gaudi::Property<double> m_thresholdADC{this, "thresholdFactor", 3.0};
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment