Skip to content
Snippets Groups Projects

Draft: Resolve "Add BTOF Benchmark"

Open Zhenyu Ye requested to merge 80-add-btof-benchmark into master
2 files
+ 5
8
Compare changes
  • Side-by-side
  • Inline
Files
2
@@ -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");
Loading