feat: add stand-alone IRT auxfile creator (#32)

Author: c-dilks

Makefile

@@ -5,49 +5,36 @@ FLAGS = -g -Wno-deprecated -fPIC -m64 -fno-inline -Wno-write-strings
 DEPS += $(shell root-config --cflags)
 LIBS += $(shell root-config --glibs)
 
-# PODIO
-LIBS += -L/usr/local/lib -lDD4pod -lpodio -lpodioRootIO -ledm4hep
+# image libs
+LIBS += -L/usr/local/lib
+LIBS += -lDDCore
+LIBS += -lDD4pod -lpodio -lpodioRootIO -ledm4hep
+LIBS += -lfmt
+
+# local libs
+LIBS += -L${EIC_SHELL_PREFIX}/lib
+LIBS += -lIRT
+DEPS += -I${EIC_SHELL_PREFIX}/include/IRT
 
 #--------------------------------------------
 
 INSTALL_PREFIX = bin
-SRC_MAIN := $(basename $(notdir $(wildcard src/*.cpp)))
-SRC_EXAMPLES := $(basename $(notdir $(wildcard src/examples/*.cpp)))
+EXECUTABLES := $(addprefix $(INSTALL_PREFIX)/, $(basename $(notdir $(wildcard src/*.cpp))))
 
 #--------------------------------------------
 
-all: 
-	@echo ""
-	@echo "BUILDING SOURCES ==========================================="
-	@echo "$(SRC_MAIN)"
-	@echo "============================================================"
-	make main
-	@echo ""
-	@echo "BUILDING EXAMPLES =========================================="
-	@echo "$(SRC_EXAMPLES)"
-	@echo "============================================================"
-	make examples
-
-main: $(SRC_MAIN)
-examples: $(SRC_EXAMPLES)
+all: $(EXECUTABLES)
 
 clean:
 	@echo "CLEAN ======================================================"
-	$(RM) $(addprefix $(INSTALL_PREFIX)/, $(SRC_MAIN))
-	$(RM) $(addprefix $(INSTALL_PREFIX)/, $(SRC_EXAMPLES))
+	$(RM) $(EXECUTABLES)
 
 #--------------------------------------------
 
-%: %.o
-	@echo "--- make executable $(INSTALL_PREFIX)/$@"
-	$(CXX) -o $(INSTALL_PREFIX)/$@ $< $(LIBS)
-
-%.o: src/%.cpp
-	mkdir -p $(INSTALL_PREFIX)
-	@echo "----- build $@ -----"
-	$(CXX) -c $^ -o $@ $(FLAGS) $(DEPS)
-
-%.o: src/examples/%.cpp
+$(INSTALL_PREFIX)/%: src/%.cpp
 	mkdir -p $(INSTALL_PREFIX)
-	@echo "----- build $@ -----"
-	$(CXX) -c $^ -o $@ $(FLAGS) $(DEPS)
+	@echo "----- build $@.o -----"
+	$(CXX) -c $^ -o $@.o $(FLAGS) $(DEPS)
+	@echo "--- make executable $@"
+	$(CXX) -o $@ $@.o $(LIBS)
+	$(RM) $@.o

README.md

@@ -305,10 +305,6 @@ corresponding executables and install them to `bin/`
 - `src/draw_segmentation.cpp` (run with `bin/draw_segmentation`)
   - reads simulation output and draws the hits within sensor pixels, which is
     useful for checking mapping of sensor segmentation (pixels)
-  - relies on `text/sensorLUT.dat`, which must be up-to-date
-    - you can produce a new version of this file by uncommenting relevant lines
-      in `epic/src/DRICH_geo.cpp` (search for `generate LUT`), and running
-      something like `./rebuild_all.sh && ./run_dd_web_display.sh`
   - build with `make`, execute with `bin/draw_segmentation [simulation_output_file]`
   - specific for dRICH; for pfRICH version, see `deprecated/pfrich/`
 
@@ -340,11 +336,10 @@ until it is time to migrate to the new reconstruction framework.
 
 Procedure:
 
-- Create the auxiliary IRT configuration file; this uses a temporary "backdoor"
-  dependency on `irt` in `epic` to produce a ROOT file containing `libIRT`
+- Create the auxiliary IRT configuration file; this is a ROOT file containing `libIRT`
   objects, such as optical boundaries, based on the dRICH geometry description.
 ```bash
-scripts/create_irt_auxfile.sh
+bin/create_irt_auxfile
 ```
 
 - Run the simulation, for example:

build.sh

@@ -63,7 +63,6 @@ case $module in
     genOpt EVALUATION=OFF
     ;;
   epic)
-    genOpt IRT_AUXFILE=ON
     ;;
   juggler)
     prefix=$JUGGLER_INSTALL_PREFIX

scripts/src/create_irt_auxfile.py renamed to deprecated/create_irt_auxfile_athena.py

@@ -1,5 +1,10 @@
 #!/usr/bin/env python
-# produce auxiliary ROOT file for the IRT algorithm
+# produce auxiliary ROOT file for the IRT algorithm; assumes the
+# production of `libIRT` objects is in the geometry `cpp` file
+# 
+# DEPRECATED for usage in EPIC; this script is preserved for 2 reasons:
+# 1. usage for ATHENA support
+# 2. example how to edit compact files on-the-fly with Python
 
 import shutil, os, sys, argparse
 import xml.etree.ElementTree as et

doc/athena.md

@@ -19,6 +19,11 @@ build.sh athena clean
 build.sh juggler clean
 ```
 
+For reconstruction, you may need to create an IRT auxfile, using
+```bash
+deprecated/create_irt_auxfile_athena.py
+```
+
 Revert to EPIC
 ==============
 To switch your environment back to EPIC:

ruby/variator/optics2.rb

@@ -51,7 +51,7 @@ def initialize
     # sensor_sphere_radius  = 100
     # sensor_sphere_centerz = -50
     # 
-    # corresponding mirror parameters (dumped from `scripts/create_irt_auxfile.sh`):
+    # corresponding mirror parameters (dumped from `create_irt_auxfile`):
     # SECTOR 0 MIRROR:
     #   mirror x = 113.888311 cm
     #   mirror y = 0.000000 cm

scripts/create_irt_auxfile.sh

@@ -7,8 +7,8 @@
 require 'pycall/import'
 
 ## settings
-NumEvents      = 100                 # number of events per fixed momentum
-NumPoints      = 20                  # number of momenta to sample
+NumEvents      = 50                  # number of events per fixed momentum
+NumPoints      = 10                  # number of momenta to sample
 PoolSize       = 6                   # number of parallel threads to run
 OutputDir      = 'out/momentum_scan' # output directory ( ! will be overwritten ! )
 RunSimRec      = true                # if false, do not run simulation+reconstruction, only draw the result