1 #ifndef TRACKRECO_PHRTREESEEDING_H
2 #define TRACKRECO_PHRTREESEEDING_H
17 #include <boost/geometry/geometries/box.hpp>
18 #include <boost/geometry/geometries/point.hpp>
19 #include <boost/geometry/index/rtree.hpp>
33 #define LogDebug(exp) std::cout << "DEBUG: " << __FILE__ << ": " << __LINE__ << ": " << exp
34 #define LogError(exp) std::cout << "ERROR: " << __FILE__ << ": " << __LINE__ << ": " << exp
35 #define LogWarning(exp) std::cout << "WARNING: " << __FILE__ << ": " << __LINE__ << ": " << exp
40 #define _USE_ALAN_TRACK_REFITTING_
51 namespace bg = boost::geometry;
52 namespace bgi = boost::geometry::index;
53 typedef bg::model::point<float, 3, bg::cs::cartesian>
point;
54 typedef bg::model::box<point>
box;
55 typedef std::pair<point, TrkrDefs::cluskey>
pointKey;
63 const std::string &
name =
"PHRTreeSeeding",
64 unsigned int nlayers_maps = 3,
65 unsigned int nlayers_intt = 4,
66 unsigned int nlayers_tpc = 48,
67 unsigned int start_layer = 53);
105 double phiadd(
double phi1,
double phi2);
106 double phidiff(
double phi1,
double phi2);
111 void QueryTree(
const bgi::rtree<
pointKey, bgi::quadratic<16>> &rtree,
double phimin,
double etamin,
double lmin,
double phimax,
double etamax,
double lmax, std::vector<pointKey> &returned_values);
129 bgi::rtree<pointKey, bgi::quadratic<16>>
_rtree;