Skip to content
Snippets Groups Projects
Commit 7abd8ea3 authored by Dmitry Kalinkin's avatar Dmitry Kalinkin
Browse files

insert_neutron: use spaces

parent 3f0f317e
No related branches found
No related tags found
No related merge requests found
rule insert_neutron_generate: rule insert_neutron_generate:
input: input:
script="benchmarks/insert_neutron/analysis/gen_particles.cxx", script="benchmarks/insert_neutron/analysis/gen_particles.cxx",
params: params:
NEVENTS_GEN=1000, NEVENTS_GEN=1000,
th_max=5.7, th_max=5.7,
th_min=2.0 th_min=2.0
output: output:
GEN_FILE="sim_output/insert_neutron/neutron_{P}GeV.hepmc" GEN_FILE="sim_output/insert_neutron/neutron_{P}GeV.hepmc"
shell: shell:
""" """
mkdir -p sim_output/insert_neutron mkdir -p sim_output/insert_neutron
root -l -b -q '{input.script}({params.NEVENTS_GEN},"{output.GEN_FILE}", "neutron", {params.th_min}, {params.th_max}, 0., 360., {wildcards.P})' root -l -b -q '{input.script}({params.NEVENTS_GEN},"{output.GEN_FILE}", "neutron", {params.th_min}, {params.th_max}, 0., 360., {wildcards.P})'
""" """
rule insert_neutron_simulate: rule insert_neutron_simulate:
input: input:
GEN_FILE="sim_output/insert_neutron/neutron_{P}GeV.hepmc", GEN_FILE="sim_output/insert_neutron/neutron_{P}GeV.hepmc",
warmup="warmup/{DETECTOR_CONFIG}.edm4hep.root", warmup="warmup/{DETECTOR_CONFIG}.edm4hep.root",
params: params:
PHYSICS_LIST="FTFP_BERT" PHYSICS_LIST="FTFP_BERT"
output: output:
SIM_FILE="sim_output/insert_neutron/{DETECTOR_CONFIG}_sim_neutron_{P}GeV_{INDEX}.edm4hep.root" SIM_FILE="sim_output/insert_neutron/{DETECTOR_CONFIG}_sim_neutron_{P}GeV_{INDEX}.edm4hep.root"
shell: shell:
""" """
NEVENTS_SIM=200 NEVENTS_SIM=200
# Running simulation # Running simulation
npsim \ npsim \
...@@ -35,28 +35,28 @@ npsim \ ...@@ -35,28 +35,28 @@ npsim \
""" """
rule insert_neutron_recon: rule insert_neutron_recon:
input: input:
SIM_FILE="sim_output/insert_neutron/{DETECTOR_CONFIG}_sim_neutron_{P}GeV_{INDEX}.edm4hep.root" SIM_FILE="sim_output/insert_neutron/{DETECTOR_CONFIG}_sim_neutron_{P}GeV_{INDEX}.edm4hep.root"
output: output:
REC_FILE="sim_output/insert_neutron/{DETECTOR_CONFIG}_rec_neutron_{P}GeV_{INDEX}.edm4eic.root" REC_FILE="sim_output/insert_neutron/{DETECTOR_CONFIG}_rec_neutron_{P}GeV_{INDEX}.edm4eic.root"
shell: shell:
""" """
NEVENTS_REC=200 NEVENTS_REC=200
eicrecon {input.SIM_FILE} -Ppodio:output_file={output.REC_FILE} -Pdd4hep:xml_files=$DETECTOR_PATH/{wildcards.DETECTOR_CONFIG}.xml -Ppodio:output_collections=MCParticles,HcalEndcapPInsertRecHits,HcalEndcapPInsertClusters,HcalEndcapPInsertSubcellHits,EcalEndcapPInsertRecHits,EcalEndcapPInsertClusters -Pjana:nevents=$NEVENTS_REC eicrecon {input.SIM_FILE} -Ppodio:output_file={output.REC_FILE} -Pdd4hep:xml_files=$DETECTOR_PATH/{wildcards.DETECTOR_CONFIG}.xml -Ppodio:output_collections=MCParticles,HcalEndcapPInsertRecHits,HcalEndcapPInsertClusters,HcalEndcapPInsertSubcellHits,EcalEndcapPInsertRecHits,EcalEndcapPInsertClusters -Pjana:nevents=$NEVENTS_REC
""" """
rule insert_neutron_analysis: rule insert_neutron_analysis:
input: input:
expand("sim_output/insert_neutron/{DETECTOR_CONFIG}_rec_neutron_{P}GeV_{INDEX}.edm4eic.root", expand("sim_output/insert_neutron/{DETECTOR_CONFIG}_rec_neutron_{P}GeV_{INDEX}.edm4eic.root",
P=[20, 30, 40, 50, 60, 70, 80], P=[20, 30, 40, 50, 60, 70, 80],
DETECTOR_CONFIG=["{DETECTOR_CONFIG}"], DETECTOR_CONFIG=["{DETECTOR_CONFIG}"],
INDEX=range(5), INDEX=range(5),
), ),
script="benchmarks/insert_neutron/analysis/neutron_plots.py", script="benchmarks/insert_neutron/analysis/neutron_plots.py",
output: output:
results_dir=directory("results/{DETECTOR_CONFIG}/insert_neutron"), results_dir=directory("results/{DETECTOR_CONFIG}/insert_neutron"),
shell: shell:
""" """
mkdir -p {output.results_dir} mkdir -p {output.results_dir}
python {input.script} {output.results_dir} python {input.script} {output.results_dir}
""" """
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment