Debug TOF

benchmarks/tof/scripts/gen_tof_hits.cxx

@@ -34,12 +34,12 @@ void gen_tof_hits(int n_events = 100,
     // pdgid 11 - electron
     // pdgid 111 - pi0
     // pdgid 2212 - proton
+    Double_t p     = r1->Uniform(0.2, 2);
     for (int ip = 0; ip < n_parts; ip++) {
       GenParticlePtr p1 = std::make_shared<GenParticle>(FourVector(0.0, 0.0, 10.0, 10.0), 11, 4);
       GenParticlePtr p2 = std::make_shared<GenParticle>(FourVector(0.0, 0.0, 0.0, 0.938), 2212, 4);
 
       // Define momentum
-      Double_t p     = r1->Uniform(0.2, 2);
       Double_t phi   = r1->Uniform(0.0, 2.0 * M_PI);
       Double_t costh = 0; //r1->Uniform(cos_theta_min, cos_theta_max);
       Double_t th    = std::acos(costh);

benchmarks/tof/scripts/sim_tof_hits.cxx

@@ -96,7 +96,6 @@ int sim_tof_hits(const char* fname = "sim_tof_hits.root")
                  .Define("N_VertexBarrelHits", [](std::vector<dd4pod::TrackerHitData> hits) { return hits.size();}, {"VertexBarrelHits"})
                  .Define("N_VertexEndcapHits", [](std::vector<dd4pod::TrackerHitData> hits) { return hits.size();}, {"VertexEndcapHits"})
                  .Define("N_BarrelHits", [](std::vector<dd4pod::TrackerHitData> hits) { return hits.size();}, {"TrackerBarrelHits"})
-                 .Define("p_barrelhits", momentum, {"BarrelTOFHits.momentum"})
                  .Define("N_EndcapHits", [](std::vector<dd4pod::TrackerHitData> hits) { return hits.size();}, {"TrackerEndcapHits"})
                  .Define("N_BarrelTOFHits", [](std::vector<dd4pod::TrackerHitData> hits) { return hits.size();}, {"BarrelTOFHits"})
                  .Define("N_ForwardTOFHits", [](std::vector<dd4pod::TrackerHitData> hits) { return hits.size();}, {"ForwardTOFHits"})
@@ -148,9 +147,9 @@ int sim_tof_hits(const char* fname = "sim_tof_hits.root")
   auto hVertexEndcap_N_vs_theta = df0.Histo1D({"hVertexEndcap_N_vs_theta", "; #theta [deg.]",   20, 0, 180 }, "theta0", "N_VertexEndcapHits");
     
 
-    auto hBarrelTof_time_vs_momentum = df0.Histo2D(
-         {"hBarrelTof_time_vs_momentum", "; time ; momentum ", 100, 0, 3, 100, 0, 5},
-         "BarrelTOFHits.truth.time", "p_barrelhits");
+    auto hBarrelTof_time_vs_length = df0.Histo2D(
+         {"hBarrelTof_time_vs_length", "; time ; length ", 100, 0, 10, 100, 0, 1},
+         "BarrelTOFHits.truth.time", "BarrelTOFHits.length");
 
     auto hBarrelTof_x_vs_y = df0.Histo2D(
          {"hBarrelTof_x_vs_y", "; x ; y ", 100, -600, 600, 100, -600, 600},
@@ -291,9 +290,9 @@ int sim_tof_hits(const char* fname = "sim_tof_hits.root")
     
     
     c = new TCanvas();
-    hBarrelTof_time_vs_momentum->DrawCopy("colz");
-    c->SaveAs("results/tof/sim_tof_hits_tofBarrel_tp.png");
-    c->SaveAs("results/tof/sim_tof_hits_tofBarrel_tp.pdf");
+    hBarrelTof_time_vs_length->DrawCopy("colz");
+    c->SaveAs("results/tof/sim_tof_hits_tofBarrel_tl.png");
+    c->SaveAs("results/tof/sim_tof_hits_tofBarrel_tl.pdf");
     c = new TCanvas();
     hBarrelTof_x_vs_y->DrawCopy("colz");
     c->SaveAs("results/tof/sim_tof_hits_tofBarrel_xy.png");

benchmarks/tof/tof_hits.sh

@@ -57,7 +57,7 @@ if [[ ! -n  "${JUGGLER_N_EVENTS}" ]] ; then
   export JUGGLER_N_EVENTS=100
 fi
 
-export JUGGLER_N_EVENTS=5
+export JUGGLER_N_EVENTS=100
 export JUGGLER_FILE_NAME_TAG="tof_hits"
 export JUGGLER_GEN_FILE="${LOCAL_DATA_PATH}/${JUGGLER_FILE_NAME_TAG}.hepmc"
 
@@ -111,7 +111,7 @@ fi
 root_filesize=$(stat --format=%s "${JUGGLER_SIM_FILE}")
 if [[ "${JUGGLER_N_EVENTS}" -lt "500" ]] ; then 
   # file must be less than 10 MB to upload
-  if [[ "${root_filesize}" -lt "10000000" ]] ; then 
+  if [[ "${root_filesize}" -lt "50000000" ]] ; then 
     cp ${JUGGLER_SIM_FILE} results/.
   fi
 fi