modified: multiple_tracks.sh

	modified:   scripts/rec_multiple_tracks.cxx

benchmarks/track_finding/multiple_tracks.sh

@@ -54,38 +54,40 @@ export JUGGLER_REC_FILE="rec_${JUGGLER_FILE_NAME_TAG}.root"
 echo "JUGGLER_N_EVENTS = ${JUGGLER_N_EVENTS}"
 echo "JUGGLER_DETECTOR = ${JUGGLER_DETECTOR}"
 
-## generate the input events
-root -b -q "benchmarks/track_finding/scripts/gen_multiple_tracks.cxx(${JUGGLER_N_EVENTS}, \"${JUGGLER_FILE_NAME_TAG}.hepmc\")"
-if [[ "$?" -ne "0" ]] ; then
-  echo "ERROR running script"
-  exit 1
-fi
-
-echo "Running geant4 simulation"
-## run geant4 simulations
-npsim --runType batch \
-  --part.minimalKineticEnergy 1000*GeV  \
-  -v WARNING \
-  --numberOfEvents ${JUGGLER_N_EVENTS} \
-  --compactFile ${DETECTOR_PATH}/${JUGGLER_DETECTOR}.xml \
-  --inputFiles  ${JUGGLER_FILE_NAME_TAG}.hepmc \
-  --outputFile  ${JUGGLER_SIM_FILE}
-if [[ "$?" -ne "0" ]] ; then
-  echo "ERROR running script"
-  exit 1
-fi
-
-rootls -t ${JUGGLER_SIM_FILE}
-
-if [[ -z "${ANALYSIS_ONLY}" ]] ;
-then
-  # Need to figure out how to pass file name to juggler from the commandline
-  xenv -x ${JUGGLER_INSTALL_PREFIX}/Juggler.xenv gaudirun.py benchmarks/track_finding/options/track_reconstruction.py
-  if [[ "$?" -ne "0" ]] ; then
-    echo "ERROR running juggler"
-    exit 1
-  fi
-fi
+### generate the input events
+#root -b -q "benchmarks/track_finding/scripts/gen_multiple_tracks.cxx(${JUGGLER_N_EVENTS}, \"${JUGGLER_FILE_NAME_TAG}.hepmc\")"
+#if [[ "$?" -ne "0" ]] ; then
+#  echo "ERROR running script"
+#  exit 1
+#fi
+#
+#echo "Running geant4 simulation"
+### run geant4 simulations
+#npsim --runType batch \
+#  --part.minimalKineticEnergy 1000*GeV  \
+#  -v WARNING \
+#  --numberOfEvents ${JUGGLER_N_EVENTS} \
+#  --compactFile ${DETECTOR_PATH}/${JUGGLER_DETECTOR}.xml \
+#  --inputFiles  ${JUGGLER_FILE_NAME_TAG}.hepmc \
+#  --outputFile  ${JUGGLER_SIM_FILE}
+#if [[ "$?" -ne "0" ]] ; then
+#  echo "ERROR running script"
+#  exit 1
+#fi
+
+#rootls -t ${JUGGLER_SIM_FILE}
+#
+#if [[ -z "${ANALYSIS_ONLY}" ]] ;
+#then
+#  # Need to figure out how to pass file name to juggler from the commandline
+#  xenv -x ${JUGGLER_INSTALL_PREFIX}/Juggler.xenv gaudirun.py benchmarks/track_finding/options/track_reconstruction.py
+#  if [[ "$?" -ne "0" ]] ; then
+#    echo "ERROR running juggler"
+#    exit 1
+#  fi
+#fi
+
+rootls -t ${JUGGLER_REC_FILE}
 
 mkdir -p results/track_finding
 

benchmarks/track_finding/scripts/rec_multiple_tracks.cxx

@@ -5,7 +5,7 @@
 #include "TProfile.h"
 
 #include <iostream>
-
+R__LOAD_LIBRARY(libJugBase.so)
 R__LOAD_LIBRARY(libeicd.so)
 R__LOAD_LIBRARY(libDD4pod.so)
 #include "dd4pod/Geant4ParticleCollection.h"
@@ -77,6 +77,9 @@ auto delta_p_over_p = [](const std::vector<double>& tracks, const std::vector<do
   }
   return res;
 };
+//gInterpreter->GenerateDictionary("vector<vector<float> >", "vector")
+//gInterpreter->GenerateDictionary("vector<unsigned long>", "vector")
+//gInterpreter->GenerateDictionary("vector<vector<unsigned long> >", "vector")
 
 int rec_multiple_tracks(const char* fname = "topside/rec_multiple_tracks.root")
 {
@@ -109,7 +112,7 @@ int rec_multiple_tracks(const char* fname = "topside/rec_multiple_tracks.root")
   auto h_delta_p0  = df0.Histo1D({"h_delta_p0", "Truth Track Init; GeV/c ",  100, -10, 10}, "delta_p0");
 
   auto h_nProtoTracks = df0.Histo1D({"h_nProtoTracks", "; n ", 10, 0, 10}, "nProto");
-  auto h_nProtoTracks2 = df0.Histo1D({"h_nProtoTracks2", "; n ", 10, 0, 10}, "nProtoTrakcs");
+  auto h_nProtoTracks2 = df0.Histo1D({"h_nProtoTracks2", "; n ", 10, 0, 10}, "nProtoTracks");
   auto h_nThrown      = df0.Histo1D({"h_nThrown",      "; n ", 10, 0, 10}, "nThrown");
 
   auto h_delta_p0_over_p = df0.Histo1D({"h_delta_p0_over_p",  "Truth Track Init; delta p/p ",  100, -0.1, 0.1}, "delta_p_over_p0");