feat(physmon): compare GX2F vs KF truth tracking

CI/physmon/config/info_only.yml

@@ -0,0 +1,7 @@
+checks:
+  "*":
+    Chi2Test: null
+    KolmogorovTest: null
+    RatioCheck: null
+    ResidualCheck: null
+    IntegralCheck: null

CI/physmon/phys_perf_mon.sh

@@ -20,8 +20,8 @@ shopt -s extglob
 
 
 mode=${1:-all}
-if ! [[ $mode = @(all|kf|gsf|gx2f|refit_kf|refit_gsf|fullchains|simulation) ]]; then
-    echo "Usage: $0 <all|kf|gsf|gx2f|refit_kf|refit_gsf|fullchains|simulation> (outdir)"
+if ! [[ $mode = @(all|kf|gsf|gx2f|refit_kf|refit_gsf|fullchains|simulation|gx2f_vs_kf) ]]; then
+    echo "Usage: $0 <all|kf|gsf|gx2f|refit_kf|refit_gsf|fullchains|simulation|gx2f_vs_kf> (outdir)"
     exit 1
 fi
 
@@ -163,6 +163,9 @@ if [[ "$mode" == "all" || "$mode" == "fullchains" ]]; then
     run_physmon_gen "CKF muon 50" "trackfinding_4muon_50vertices"
     run_physmon_gen "CKF ttbar 200" "trackfinding_ttbar_pu200"
 fi
+if [[ "$mode" == "all" || "$mode" == "gx2f_vs_kf" ]]; then
+    run_physmon_gen "Comparison - Truth Tracking GX2F vs KF" "trackfitting_gx2f_vs_kf"
+fi
 echo "::endgroup::"
 
 
@@ -468,6 +471,18 @@ if [[ "$mode" == "all" || "$mode" == "fullchains" ]]; then
     vertexing "trackfinding | ttbar with 200 pileup | default seeding" trackfinding_ttbar_pu200 CI/physmon/config/vertexing_ttbar_pu200.yml
 fi
 
+if [[ "$mode" == "all" || "$mode" == "gx2f_vs_kf" ]]; then
+    run_histcmp \
+        $outdir/data/trackfitting_gx2f_vs_kf/performance_trackfitting_gx2f.root \
+        $outdir/data/trackfitting_gx2f_vs_kf/performance_trackfitting_kf.root \
+        "Comparison - Truth tracking (GX2F vs KF)" \
+        trackfitting_gx2f_vs_kf/performance_trackfitting.html \
+        trackfitting_gx2f_vs_kf/performance_trackfitting_plots \
+        --config CI/physmon/config/info_only.yml \
+        --label-reference=KF \
+        --label-monitored=GX2F
+fi
+
 run CI/physmon/summary.py $histcmp_results \
   --md $outdir/summary.md \
   --html $outdir/summary.html

CI/physmon/summary.py

@@ -48,10 +48,16 @@
         )
 
         for s in summary:
-            f.write(
-                f"""
+            if s["title"].startswith("Comparison"):
+                f.write(
+                    f"""
+        <li>🔵 <a href="{s["path"]}">{s["title"]}</a></li>"""
+                )
+            else:
+                f.write(
+                    f"""
         <li>{"✅" if s["total"] else "🔴"} <a href="{s["path"]}">{s["title"]}</a></li>"""
-            )
+                )
 
         f.write(
             """
@@ -74,4 +80,8 @@
                 )
             else:
                 url = s["path"]
-            f.write(f"  - {'✅' if s['total'] else '🔴'} [{s['title']}]({url})\n")
+
+            if s["title"].startswith("Comparison"):
+                f.write(f"  - 🔵️ [{s['title']}]({url})\n")
+            else:
+                f.write(f"  - {'✅' if s['total'] else '🔴'} [{s['title']}]({url})\n")

CI/physmon/workflows/physmon_trackfitting_gx2f_vs_kf.py

@@ -0,0 +1,157 @@
+#!/usr/bin/env python3
+
+import tempfile
+
+from pathlib import Path
+
+import acts
+import acts.examples
+
+from acts.examples.simulation import (
+    addParticleGun,
+    ParticleConfig,
+    EtaConfig,
+    PhiConfig,
+    MomentumConfig,
+    ParticleSelectorConfig,
+    addFatras,
+    addDigitization,
+)
+from acts.examples.reconstruction import (
+    addSeeding,
+    SeedingAlgorithm,
+    addGx2fTracks,
+    addKalmanTracks,
+)
+
+from physmon_common import makeSetup
+
+u = acts.UnitConstants
+
+setup = makeSetup()
+
+trackingGeometry = setup.trackingGeometry
+field = setup.field
+digiConfigFile = setup.digiConfig
+
+with tempfile.TemporaryDirectory() as temp:
+    s = acts.examples.Sequencer(
+        events=100000,
+        numThreads=-1,
+        logLevel=acts.logging.INFO,
+        skip=0,
+        trackFpes=False,
+    )
+
+    rnd = acts.examples.RandomNumbers(seed=42)
+    outputDirTemp = Path(temp)
+    outputDir = setup.outdir
+
+    addParticleGun(
+        s,
+        ParticleConfig(num=1, pdg=acts.PdgParticle.eMuon, randomizeCharge=True),
+        EtaConfig(-3.0, 3.0, uniform=True),
+        MomentumConfig(1.0 * u.GeV, 100.0 * u.GeV, transverse=True),
+        PhiConfig(0.0, 360.0 * u.degree),
+        vtxGen=acts.examples.GaussianVertexGenerator(
+            mean=acts.Vector4(0, 0, 0, 0),
+            stddev=acts.Vector4(0, 0, 0, 0),
+        ),
+        multiplicity=1,
+        rnd=rnd,
+        outputDirRoot=outputDirTemp,
+    )
+
+    addFatras(
+        s,
+        trackingGeometry,
+        field,
+        rnd=rnd,
+        enableInteractions=True,
+        postSelectParticles=ParticleSelectorConfig(
+            pt=(0.9 * u.GeV, None),
+            measurements=(7, None),
+            removeNeutral=True,
+            removeSecondaries=True,
+        ),
+    )
+
+    addDigitization(
+        s,
+        trackingGeometry,
+        field,
+        digiConfigFile=digiConfigFile,
+        rnd=rnd,
+    )
+
+    addSeeding(
+        s,
+        trackingGeometry,
+        field,
+        rnd=rnd,
+        inputParticles="particles_input",
+        seedingAlgorithm=SeedingAlgorithm.TruthSmeared,
+        particleHypothesis=acts.ParticleHypothesis.muon,
+    )
+
+    addGx2fTracks(
+        s,
+        trackingGeometry,
+        field,
+        nUpdateMax=17,
+        relChi2changeCutOff=1e-7,
+        multipleScattering=True,
+    )
+
+    s.addAlgorithm(
+        acts.examples.TrackSelectorAlgorithm(
+            level=acts.logging.INFO,
+            inputTracks="tracks",
+            outputTracks="selected-tracks-gx2f",
+            selectorConfig=acts.TrackSelector.Config(
+                minMeasurements=7,
+            ),
+        )
+    )
+
+    s.addWriter(
+        acts.examples.TrackFitterPerformanceWriter(
+            level=acts.logging.INFO,
+            inputTracks="selected-tracks-gx2f",
+            inputParticles="particles_selected",
+            inputTrackParticleMatching="track_particle_matching",
+            filePath=str(outputDir / "performance_trackfitting_gx2f.root"),
+        )
+    )
+
+    addKalmanTracks(
+        s,
+        trackingGeometry,
+        field,
+        directNavigation=False,
+        reverseFilteringMomThreshold=0 * u.GeV,
+    )
+
+    s.addAlgorithm(
+        acts.examples.TrackSelectorAlgorithm(
+            level=acts.logging.INFO,
+            inputTracks="tracks",
+            outputTracks="selected-tracks-kf",
+            selectorConfig=acts.TrackSelector.Config(
+                minMeasurements=7,
+            ),
+        )
+    )
+
+    s.addWriter(
+        acts.examples.TrackFitterPerformanceWriter(
+            level=acts.logging.INFO,
+            inputTracks="selected-tracks-kf",
+            inputParticles="particles_selected",
+            inputTrackParticleMatching="track_particle_matching",
+            filePath=str(outputDir / "performance_trackfitting_kf.root"),
+        )
+    )
+
+    s.run()
+    del s