Warning, /detector_benchmarks/benchmarks/backgrounds/Snakefile is written in an unsupported language. File is not indexed.
0001 import os
0002 import shutil
0003
0004
0005 rule backgrounds_sim:
0006 input:
0007 warmup="warmup.edm4hep.root",
0008 geometry_lib=find_epic_libraries(),
0009 output:
0010 "sim_output/{DETECTOR_CONFIG}/backgrounds/{PATH}.edm4hep.root",
0011 log:
0012 "sim_output/{DETECTOR_CONFIG}/backgrounds/{PATH}.edm4hep.root.log",
0013 params:
0014 N_EVENTS=100,
0015 hepmc=lambda wildcards: get_remote_path(f"{wildcards.PATH}.hepmc3.tree.root"),
0016 DD4HEP_HASH=get_spack_package_hash("dd4hep"),
0017 NPSIM_HASH=get_spack_package_hash("npsim"),
0018 cache: True
0019 singularity: EIC_SINGULARITY_CONTAINER,
0020 shell:
0021 """
0022 set -m # monitor mode to prevent lingering processes
0023 exec npsim \
0024 --runType batch \
0025 --part.minimalKineticEnergy 100*GeV \
0026 --filter.tracker edep0 \
0027 -v WARNING \
0028 --numberOfEvents {params.N_EVENTS} \
0029 --compactFile $DETECTOR_PATH/{wildcards.DETECTOR_CONFIG}.xml \
0030 --inputFiles {params.hepmc} \
0031 --outputFile {output}
0032 """
0033
0034
0035 DETECTOR_CONFIG=os.environ["DETECTOR_CONFIG"]
0036
0037 rule backgrounds_ecal_backwards:
0038 input:
0039 matplotlibrc=".matplotlibrc",
0040 script="benchmarks/backgrounds/ecal_backwards.org2py.py",
0041 electron_beam_gas_sim="sim_output/" + DETECTOR_CONFIG + "/backgrounds/EPIC/EVGEN/BACKGROUNDS/BEAMGAS/electron/GETaLM1.0.0-1.0/10GeV/GETaLM1.0.0-1.0_ElectronBeamGas_10GeV_foam_emin10keV_run001.edm4hep.root",
0042 physics_signal_sim="sim_output/" + DETECTOR_CONFIG + "/backgrounds/EPIC/EVGEN/DIS/NC/10x100/minQ2=1/pythia8NCDIS_10x100_minQ2=1_beamEffects_xAngle=-0.025_hiDiv_1.edm4hep.root",
0043 proton_beam_gas_sim="sim_output/" + DETECTOR_CONFIG + "/backgrounds/EPIC/EVGEN/BACKGROUNDS/BEAMGAS/proton/pythia8.306-1.0/100GeV/pythia8.306-1.0_ProtonBeamGas_100GeV_run001.edm4hep.root",
0044 output:
0045 directory("results/backgrounds/backwards_ecal")
0046 log:
0047 scheduler=".logs/results/backgrounds/backwards_ecal/scheduler.log",
0048 worker=".logs/results/backgrounds/backwards_ecal/worker.log",
0049 threads: workflow.cores
0050 singularity: EIC_SINGULARITY_CONTAINER,
0051 shell:
0052 """
0053 set -m # monitor mode to prevent lingering processes
0054 cleanup() {{
0055 echo Cleaning up
0056 kill $WORKER_PID $SCHEDULER_PID
0057 }}
0058 trap cleanup EXIT
0059
0060 PORT=$RANDOM
0061 dask scheduler --port $PORT 2>{log.scheduler} &
0062 export DASK_SCHEDULER=localhost:$PORT
0063 SCHEDULER_PID=$!
0064 dask worker tcp://$DASK_SCHEDULER --nworkers {threads} --nthreads 1 2>{log.worker} &
0065 WORKER_PID=$!
0066 env \
0067 MATPLOTLIBRC={input.matplotlibrc} \
0068 ELECTRON_BEAM_GAS_GEN=root://dtn-eic.jlab.org//volatile/eic/EPIC/EVGEN/BACKGROUNDS/BEAMGAS/electron/GETaLM1.0.0-1.0/10GeV/GETaLM1.0.0-1.0_ElectronBeamGas_10GeV_foam_emin10keV_run001.hepmc3.tree.root \
0069 ELECTRON_BEAM_GAS_SIM=$(realpath {input.electron_beam_gas_sim}) \
0070 PHYSICS_PROCESS_SIM=$(realpath {input.physics_signal_sim}) \
0071 PROTON_BEAM_GAS_GEN=root://dtn-eic.jlab.org//volatile/eic/EPIC/EVGEN/BACKGROUNDS/BEAMGAS/proton/pythia8.306-1.0/100GeV/pythia8.306-1.0_ProtonBeamGas_100GeV_run001.hepmc3.tree.root \
0072 PROTON_BEAM_GAS_SIM=$(realpath {input.proton_beam_gas_sim}) \
0073 OUTPUT_DIR={output} \
0074 python {input.script}
0075 """
