Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
R__LOAD_LIBRARY(libActsDD4hepPlugin.so)
R__LOAD_LIBRARY(libDDG4IO.so)
R__LOAD_LIBRARY(libGenDetectors.so)
#include "DD4hep/Detector.h"
#include "DDG4/Geant4Data.h"
#include "DDRec/CellIDPositionConverter.h"
#include "DDRec/SurfaceManager.h"
#include "DDRec/Surface.h"
#include "ROOT/RDataFrame.hxx"
#include "TCanvas.h"
#include "TChain.h"
#include "Acts/Geometry/TrackingGeometry.hpp"
#include "Acts/Geometry/TrackingVolume.hpp"
#include "Acts/Plugins/DD4hep/ConvertDD4hepDetector.hpp"
/** Example loading ACTs.
*
*
*/
void example_Acts(const char* fname = "test_tracker_disc.root"){
using namespace ROOT::Math;
ROOT::EnableImplicitMT(4);
TChain* t = new TChain("EVENT");
t->Add(fname);
ROOT::RDataFrame d0(*t, {"GEMTrackerHits","MCParticles"});
//How to get the type of the initial branch: (vector<dd4hep::sim::Geant4Tracker::Hit*>)
//std::cout << t->GetBranch("GEMTrackerHits")->GetClassName() << std::endl;
// -------------------------
// Get the DD4hep instance
// Load the compact XML file
// Initialize the position converter tool
dd4hep::Detector& detector = dd4hep::Detector::getInstance();
detector.fromCompact("gem_tracker_disc.xml");
dd4hep::rec::CellIDPositionConverter cellid_converter(detector);
// -------------------------
// Get the surfaces map
dd4hep::rec::SurfaceManager& surfMan = *detector.extension<dd4hep::rec::SurfaceManager>() ;
auto surfMap = surfMan.map( "world" ) ;
//std::unique_ptr<const Acts::TrackingGeometry>
auto acts_tracking_geometry = Acts::convertDD4hepDetector (detector.world(),Acts::Logging::Level::VERBOSE);
if(acts_tracking_geometry) {
if(acts_tracking_geometry->highestTrackingVolume()) {
std::cout << " volume name \n ";
std::cout << acts_tracking_geometry->highestTrackingVolume()->volumeName() << std::endl;
} else {
std::cout << "derp\n";
}
}
}