Warning, /detector_benchmarks/benchmarks/femc_electron/Snakefile is written in an unsupported language. File is not indexed.
0001 def get_n_events(wildcards):
0002 energy = float(wildcards.P)
0003 n_events = 1000
0004 n_events = int(n_events // ((energy / 20) ** 0.5))
0005 return n_events
0006
0007
0008 rule femc_electron_generate:
0009 input:
0010 script="benchmarks/femc_electron/analysis/gen_particles.cxx",
0011 output:
0012 GEN_FILE="sim_output/femc_electron/e-_{P}GeV.hepmc",
0013 params:
0014 N_EVENTS=get_n_events,
0015 th_max=28,
0016 th_min=2.0,
0017 P=lambda wildcards: wildcards.P,
0018 shell:
0019 """
0020 mkdir -p sim_output/femc_electron
0021 root -l -b -q '{input.script}({params.N_EVENTS},"{output.GEN_FILE}", "e-", {params.th_min}, {params.th_max}, 0., 360., {params.P})'
0022 """
0023
0024 rule femc_electron_simulate:
0025 input:
0026 GEN_FILE="sim_output/femc_electron/e-_{P}GeV.hepmc",
0027 warmup="warmup/{DETECTOR_CONFIG}.edm4hep.root",
0028 geometry_lib=find_epic_libraries(),
0029 output:
0030 SIM_FILE="sim_output/femc_electron/{DETECTOR_CONFIG}_sim_e-_{P}GeV.edm4hep.root",
0031 params:
0032 N_EVENTS=get_n_events,
0033 PHYSICS_LIST="FTFP_BERT",
0034 DETECTOR_PATH=os.environ["DETECTOR_PATH"],
0035 DETECTOR_CONFIG=lambda wildcards: wildcards.DETECTOR_CONFIG,
0036 DD4HEP_HASH=get_spack_package_hash("dd4hep"),
0037 NPSIM_HASH=get_spack_package_hash("npsim"),
0038 cache: True
0039 shell:
0040 """
0041 # Running simulation
0042 npsim \
0043 --compactFile {params.DETECTOR_PATH}/{params.DETECTOR_CONFIG}.xml \
0044 --numberOfEvents {params.N_EVENTS} \
0045 --physicsList {params.PHYSICS_LIST} \
0046 --inputFiles {input.GEN_FILE} \
0047 --outputFile {output.SIM_FILE}
0048 """
0049
0050 rule femc_electron_recon:
0051 input:
0052 SIM_FILE="sim_output/femc_electron/{DETECTOR_CONFIG}_sim_e-_{P}GeV.edm4hep.root",
0053 output:
0054 REC_FILE="sim_output/femc_electron/{DETECTOR_CONFIG}_rec_e-_{P}GeV.edm4eic.root",
0055 params:
0056 N_EVENTS=get_n_events,
0057 DETECTOR_CONFIG=lambda wildcards: wildcards.DETECTOR_CONFIG,
0058 EICRECON_HASH=get_spack_package_hash("eicrecon"),
0059 cache: True
0060 shell:
0061 """
0062 env DETECTOR_CONFIG={params.DETECTOR_CONFIG} \
0063 eicrecon {input.SIM_FILE} -Ppodio:output_file={output.REC_FILE} -Ppodio:output_collections=MCParticles,EcalEndcapPClusters -Pjana:nevents={params.N_EVENTS}
0064 """
0065
0066 rule femc_electron_analysis:
0067 input:
0068 expand("sim_output/femc_electron/{DETECTOR_CONFIG}_rec_e-_{P}GeV.edm4eic.root",
0069 P=[10, 20, 30, 40, 50, 60, 70, 80],
0070 DETECTOR_CONFIG=["{DETECTOR_CONFIG}"]),
0071 script="benchmarks/femc_electron/analysis/femc_electron_plots.py",
0072 output:
0073 results_dir=directory("results/{DETECTOR_CONFIG}/femc_electron"),
0074 shell:
0075 """
0076 mkdir -p {output.results_dir}
0077 python {input.script} {output.results_dir}
0078 """
