From 91a4322a7fbcc6f021d1f52ebdebafdde552fcf8 Mon Sep 17 00:00:00 2001 From: Zhenyu Ye <yezhenyu@EIC-SoLID.local> Date: Sun, 26 Sep 2021 09:19:00 -0500 Subject: [PATCH] debug TOF benchmark --- benchmarks/track_finding/scripts/gen_multiple_tracks.cxx | 2 +- benchmarks/track_finding/scripts/rec_multiple_tracks.cxx | 6 ++++++ 2 files changed, 7 insertions(+), 1 deletion(-) diff --git a/benchmarks/track_finding/scripts/gen_multiple_tracks.cxx b/benchmarks/track_finding/scripts/gen_multiple_tracks.cxx index 7960c8ca..8596e929 100644 --- a/benchmarks/track_finding/scripts/gen_multiple_tracks.cxx +++ b/benchmarks/track_finding/scripts/gen_multiple_tracks.cxx @@ -44,7 +44,7 @@ void gen_multiple_tracks(int n_events = 100, // Define momentum Double_t p = r1->Uniform(1.0, 10.0); Double_t phi = r1->Uniform(0.0, 2.0 * M_PI); - Double_t costh = r1->Uniform(-0.05, 0.05); //r1->Uniform(cos_theta_min, cos_theta_max); + Double_t costh = r1->Uniform(cos_theta_min, cos_theta_max); Double_t th = std::acos(costh); Double_t px = p * std::cos(phi) * std::sin(th); Double_t py = p * std::sin(phi) * std::sin(th); diff --git a/benchmarks/track_finding/scripts/rec_multiple_tracks.cxx b/benchmarks/track_finding/scripts/rec_multiple_tracks.cxx index 788b9e2a..4a842984 100644 --- a/benchmarks/track_finding/scripts/rec_multiple_tracks.cxx +++ b/benchmarks/track_finding/scripts/rec_multiple_tracks.cxx @@ -102,6 +102,9 @@ 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();}, {"TrackerBarrelTOFRecHits"}) + .Define("N_ForwardTOFHits", [](std::vector<eic::TrackerHitData> hits) { return hits.size();}, {"TrackerForwardTOFRecHits"}) + .Define("N_BackwardTOFHits", [](std::vector<eic::TrackerHitData> hits) { return hits.size();}, {"TrackerBackwardTOFRecHits"}) ; auto h_nTracks_vs_theta = df0.Histo2D({"h_nTracks_vs_theta", "; #theta; N tracks ", 40,0,180,10, 0, 10}, "theta0","nTracks"); @@ -120,6 +123,9 @@ 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 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