#!/usr/bin/env python3
from LDMX.Framework import ldmxcfg
p = ldmxcfg.Process('test')

p.maxTriesPerEvent = 10000

from LDMX.Biasing import ecal, util
from LDMX.SimCore import generators as gen
mySim = ecal.photo_nuclear('ldmx-det-v14',gen.single_4gev_e_upstream_tagger())
mySim.beamSpotSmear = [20.,80.,0.]
mySim.description = 'ECal PN Test Simulation'
# step = util.StepPrinter(process_name='CoulombScat', track_id=-1)
# step.depth=3
# mySim.actions.extend([step])

p.sequence = [ mySim ]

##################################################################
# Below should be the same for all sim scenarios

import os
import sys

if 'LDMX_NUM_EVENTS' in os.environ:
    p.maxEvents = int(os.environ['LDMX_NUM_EVENTS'])
else:
    p.maxEvents = int(sys.argv[1])
if 'LDMX_RUN_NUMBER' in os.environ:
    p.run = int(os.environ['LDMX_RUN_NUMBER'])
else:
    p.run = 1




rectype = sys.argv[2]

output_base = f'swan/data/type_ecalpn_N_{p.maxEvents}_rectype_{rectype}'
output_base = f'data/type_hcalrectest_N_{p.maxEvents}_rectype_{rectype}'

p.histogramFile = f'{output_base}_kind_hist.root'
p.outputFiles = [f'{output_base}_kind_events.root']

import LDMX.Ecal.EcalGeometry
import LDMX.Ecal.ecal_hardcoded_conditions
import LDMX.Hcal.HcalGeometry
import LDMX.Hcal.hcal_hardcoded_conditions
import LDMX.Hcal.digi as digi


from LDMX.DQM import dqm

if rectype == 'simple':
    p.sequence.extend([digi.HcalSimpleDigiAndRecProducer()])
else:
    p.sequence.extend([
        digi.HcalDigiProducer(),
        digi.HcalRecProducer(),
    ])
p.sequence.extend(
                   dqm.hcal_dqm
                  )
![hcal_dqm_back_along_x](https://github.com/LDMX-Software/Hcal/assets/14233914/9f8d189f-1986-46f4-8543-2f797adf9517)