Skip to content
Snippets Groups Projects

Resolve "CI pipeline with AMEGO-X type EMCal Barrel"

Compare and
10 files
+ 584
13
Compare changes
  • Side-by-side
  • Inline

Files

@@ -10,9 +10,11 @@ from GaudiKernel.DataObjectHandleBase import DataObjectHandleBase
@@ -10,9 +10,11 @@ from GaudiKernel.DataObjectHandleBase import DataObjectHandleBase
from Configurables import ApplicationMgr, EICDataSvc, PodioOutput, GeoSvc
from Configurables import ApplicationMgr, EICDataSvc, PodioOutput, GeoSvc
from GaudiKernel import SystemOfUnits as units
from GaudiKernel import SystemOfUnits as units
 
directory_name = "calorimeters"
detector_name = "topside"
detector_name = "topside"
if "JUGGLER_DETECTOR" in os.environ :
if "JUGGLER_DETECTOR" in os.environ :
detector_name = str(os.environ["JUGGLER_DETECTOR"])+"/"+detector_name
#detector_name = str(os.environ["JUGGLER_DETECTOR"])+"/"+detector_name
 
detector_name = directory_name+"/"+str(os.environ["JUGGLER_DETECTOR"])
input_sim_file = "jug_input.root"
input_sim_file = "jug_input.root"
if "JUGGLER_SIM_FILE" in os.environ :
if "JUGGLER_SIM_FILE" in os.environ :
@@ -62,23 +64,28 @@ embarreldigi = EcalTungstenSamplingDigi("ecal_barrel_digi",
@@ -62,23 +64,28 @@ embarreldigi = EcalTungstenSamplingDigi("ecal_barrel_digi",
outputHitCollection="RawEcalBarrelHits",
outputHitCollection="RawEcalBarrelHits",
inputEnergyUnit=units.GeV,
inputEnergyUnit=units.GeV,
inputTimeUnit=units.ns,
inputTimeUnit=units.ns,
 
DynamicRangeADC=700*units.keV,
 
energyResolutions=[0., 0.02, 0.],
 
pedestalSigma=50,
OutputLevel=DEBUG)
OutputLevel=DEBUG)
# Reconstruction
# Reconstruction
embarrelreco = EcalTungstenSamplingReco("ecal_barrel_reco",
embarrelreco = EcalTungstenSamplingReco("ecal_barrel_reco",
inputHitCollection="RawEcalBarrelHits",
inputHitCollection="RawEcalBarrelHits",
outputHitCollection="RecoEcalBarrelHits",
outputHitCollection="RecoEcalBarrelHits",
 
DynamicRangeADC=700*units.keV,
 
pedestalSigma=50,
OutputLevel=DEBUG)
OutputLevel=DEBUG)
# 2D+1 Clusterings
# 2D+1 Clusterings
# readout id definition for barrel ecal
# readout id definition for barrel ecal
# <id>system:8,barrel:3,module:4,layer:6,slice:5,x:32:-16,y:-16</id>
# <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+6-1)
# xy_merger sum layers/slices, masking (8+3+4, 8+3+4+5+10-1)
embarrelxymerger = SamplingECalHitsMerger("ecal_barrel_xy_merger",
embarrelxymerger = SamplingECalHitsMerger("ecal_barrel_xy_merger",
cellIDMaskRanges=[(15, 25)],
cellIDMaskRanges=[(15, 29)],
inputHitCollection="RecoEcalBarrelHits",
inputHitCollection="RecoEcalBarrelHits",
outputHitCollection="RecoEcalBarrelHitsXY")
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)
embarrelzmerger = SamplingECalHitsMerger("ecal_barrel_z_merger",
embarrelzmerger = SamplingECalHitsMerger("ecal_barrel_z_merger",
cellIDMaskRanges=[(26, 57)],
cellIDMaskRanges=[(30, 61)],
inputHitCollection="RecoEcalBarrelHits",
inputHitCollection="RecoEcalBarrelHits",
outputHitCollection="RecoEcalBarrelHitsZ")
outputHitCollection="RecoEcalBarrelHitsZ")
# Clustering
# Clustering
Loading