Skip to content
Snippets Groups Projects
Commit 1dec81c6 authored by Whitney Armstrong's avatar Whitney Armstrong
Browse files

modified: scripts/example_hit_recon.cxx

	modified:   scripts/example_tracking.cxx
parent 40bec3a4
Branches master
No related tags found
No related merge requests found
......@@ -47,7 +47,7 @@ void example_hit_recon(const char* fname = "gem_tracker_digi.root"){
auto nhits = [] (const std::vector<dd4hep::sim::Geant4Tracker::Hit*>& hits){ return hits.size(); };
dd4hep::Detector& detector = dd4hep::Detector::getInstance();
detector.fromCompact("gem_tracker_disc.xml");
detector.fromCompact("gem_tracker.xml");
dd4hep::rec::CellIDPositionConverter cellid_converter(detector);
dd4hep::rec::SurfaceManager& surfMan = *detector.extension<dd4hep::rec::SurfaceManager>() ;
......
R__LOAD_LIBRARY(libGenFind)
R__LOAD_LIBRARY(libgenfit2)
R__LOAD_LIBRARY(libDDG4IO.so)
#include "ConstField.h"
#include "Exception.h"
......@@ -76,7 +77,8 @@ R__LOAD_LIBRARY(libgenfit2)
void example_tracking(const char* fname = "gem_tracker_recon.root"){
ROOT::EnableImplicitMT(1);
//ROOT::EnableImplicitMT(0);
ROOT::DisableImplicitMT(); // Tell ROOT you want to go parallel
using namespace ROOT::Math;
using namespace genfind;
......@@ -98,7 +100,7 @@ void example_tracking(const char* fname = "gem_tracker_recon.root"){
auto nhits = [] (const std::vector<dd4hep::sim::Geant4Tracker::Hit*>& hits){ return hits.size(); };
dd4hep::Detector& detector = dd4hep::Detector::getInstance();
detector.fromCompact("gem_tracker_disc.xml");
detector.fromCompact("gem_tracker.xml");
dd4hep::rec::CellIDPositionConverter cellid_converter(detector);
ROOT::Math::XYZVector position = {0,0,0}; // position to calculate magnetic field at (the origin in this case)
......@@ -236,7 +238,7 @@ void example_tracking(const char* fname = "gem_tracker_recon.root"){
}
display->addEvent(all_tracks);
//std::cout << " |P| = " << (*mc_particles)[0]->psx << ", " << (*mc_particles)[0]->psy << ", " << (*mc_particles)[0]->psz << std::endl ;
//display->open();
display->open();
//std::string dummy;
//std::cout << "Enter to continue..." << std::endl;
//std::getline(std::cin, dummy);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment