From 0997a4f0834675f976f78c5df4cf4d59ed4e438e Mon Sep 17 00:00:00 2001
From: Zhenyu Ye <yezhenyu@EIC-SoLID.local>
Date: Thu, 4 Nov 2021 21:20:06 -0500
Subject: [PATCH] add BTOF benchmark
---
benchmarks/tof/multiple_tracks.sh | 1 +
benchmarks/tof/scripts/rec_multiple_tracks.cxx | 12 ++++--------
2 files changed, 5 insertions(+), 8 deletions(-)
diff --git a/benchmarks/tof/multiple_tracks.sh b/benchmarks/tof/multiple_tracks.sh
index d5790802..97a40b1e 100644
--- a/benchmarks/tof/multiple_tracks.sh
+++ b/benchmarks/tof/multiple_tracks.sh
@@ -44,6 +44,7 @@ if [[ ! -n "${JUGGLER_N_EVENTS}" ]] ; then
fi
#export JUGGLER_N_EVENTS=$(expr ${JUGGLER_N_EVENTS} \* 1)
+export JUGGLER_DETECTOR_VERSION="131-update-tof-for-canyonlands"
export JUGGLER_N_EVENTS=1
export JUGGLER_FILE_NAME_TAG="multiple_tracks"
export JUGGLER_GEN_FILE="${JUGGLER_FILE_NAME_TAG}.hepmc"
diff --git a/benchmarks/tof/scripts/rec_multiple_tracks.cxx b/benchmarks/tof/scripts/rec_multiple_tracks.cxx
index d4d3756f..972c4e24 100644
--- a/benchmarks/tof/scripts/rec_multiple_tracks.cxx
+++ b/benchmarks/tof/scripts/rec_multiple_tracks.cxx
@@ -87,8 +87,8 @@ int rec_multiple_tracks(const char* fname = "topside/rec_multiple_tracks.root")
ROOT::EnableImplicitMT();
ROOT::RDataFrame df("events", fname);
- auto df0 = df.Define("isThrown", "mcparticles2.genStatus == 1")
- .Define("thrownParticles", "mcparticles2[isThrown]")
+ auto df0 = df.Define("isThrown", "mcparticles.genStatus == 1")
+ .Define("thrownParticles", "mcparticles[isThrown]")
.Define("thrownP", fourvec, {"thrownParticles"})
.Define("nThrown", "thrownParticles.size()")
.Define("nProto", "outputProtoTracks.size()")
@@ -102,9 +102,7 @@ int rec_multiple_tracks(const char* fname = "topside/rec_multiple_tracks.root")
.Define("N_Hits", [](std::vector<eic::TrackerHitData> hits) { return hits.size();}, {"trackingHits"})
.Define("N_BarrelHits", [](std::vector<eic::TrackerHitData> hits) { return hits.size();}, {"TrackerBarrelRecHits"})
.Define("N_EndcapHits", [](std::vector<eic::TrackerHitData> hits) { return hits.size();}, {"TrackerEndcapRecHits"})
- .Define("N_BarrelTOFHits", [](std::vector<eic::TrackerHitData> hits) { return hits.size();}, {"BarrelTOFRecHits"})
- .Define("N_ForwardTOFHits", [](std::vector<eic::TrackerHitData> hits) { return hits.size();}, {"ForwardTOFRecHits"})
- .Define("N_BackwardTOFHits", [](std::vector<eic::TrackerHitData> hits) { return hits.size();}, {"BackwardTOFRecHits"})
+ //.Define("N_TOFBarrelHits", [](std::vector<eic::TrackerHitData> hits) { return hits.size();}, {"TOFBarrelRecHits"})
;
auto h_nTracks_vs_theta = df0.Histo2D({"h_nTracks_vs_theta", "; #theta; N tracks ", 40,0,180,10, 0, 10}, "theta0","nTracks");
@@ -123,9 +121,7 @@ int rec_multiple_tracks(const char* fname = "topside/rec_multiple_tracks.root")
auto hNhits_vs_theta = df0.Histo1D({"hNhits_vs_theta", "; #theta [deg.]", 40, 0, 180 }, "theta0", "N_Hits");
auto hBarrel_N_vs_theta = df0.Histo1D({"hBarrel_N_vs_theta", "; #theta [deg.]", 40, 0, 180 }, "theta0", "N_BarrelHits");
auto hEndcap_N_vs_theta = df0.Histo1D({"hEndcap_N_vs_theta", "; #theta [deg.]", 40, 0, 180 }, "theta0", "N_EndcapHits");
- auto hBarrelTOF_N_vs_theta = df0.Histo1D({"hBarrelTOF_N_vs_theta", "; #theta [deg.]", 40, 0, 180 }, "theta0", "N_BarrelTOFHits");
- auto hForwardTOF_N_vs_theta = df0.Histo1D({"hForwardTOF_N_vs_theta", "; #theta [deg.]", 40, 0, 180 }, "theta0", "N_ForwardTOFHits");
- auto hBackwardTOF_N_vs_theta = df0.Histo1D({"hBackwardTOF_N_vs_theta", "; #theta [deg.]", 40, 0, 180 }, "theta0", "N_BackwardTOFHits");
+ //auto hTOFBarrel_N_vs_theta = df0.Histo1D({"hTOFBarrel_N_vs_theta", "; #theta [deg.]", 40, 0, 180 }, "theta0", "N_TOFBarrelHits");
auto hHits_Nhits = df0.Histo1D({"hHits_Nhits", "; #theta [deg.]", 20, 0, 20 }, "N_Hits");
auto hBarrel_Nhits = df0.Histo1D({"hBarrel_Nhits", "; #theta [deg.]", 20, 0, 20 }, "N_BarrelHits");
--
GitLab