All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Functions
sbnpca Namespace Reference

Functions

float VecAngle (std::array< float, 2 > A, std::array< float, 2 > B)
 
std::array< float, 2 > HitVector (const recob::Hit &A, const geo::GeometryCore *geo, const detinfo::DetectorPropertiesData &dprop)
 
std::array< float, 2 > HitVector (const recob::Hit &A, const recob::Hit &B, const geo::GeometryCore *geo, const detinfo::DetectorPropertiesData &dprop)
 
float HitDistance (const recob::Hit &A, const recob::Hit &B, const geo::GeometryCore *geo, const detinfo::DetectorPropertiesData &dprop)
 
std::tuple< std::vector
< art::Ptr< recob::Hit >
>, std::vector< art::Ptr
< recob::Hit > >, bool > 		GetNearestHits (const std::vector< art::Ptr< recob::Hit >> &hits, int ihit, float distance, const geo::GeometryCore *geo, const detinfo::DetectorPropertiesData &dprop)
 
std::array< float, 2 > HitPCAVec (const std::vector< art::Ptr< recob::Hit >> &hits, const recob::Hit &center, const geo::GeometryCore *geo, const detinfo::DetectorPropertiesData &dprop)
 
std::array< float, 2 > HitPCAVec (const std::vector< art::Ptr< recob::Hit >> &hits, const recob::Vertex &center, const geo::GeometryCore *geo, const detinfo::DetectorPropertiesData &dprop)
 
std::array< float, 2 > HitPCAVec (const std::vector< art::Ptr< recob::Hit >> &hits, std::array< float, 2 > center, const geo::GeometryCore *geo, const detinfo::DetectorPropertiesData &dprop)
 
std::array< float, 2 > HitPCAEigen (const std::vector< art::Ptr< recob::Hit >> &hits, const art::Ptr< recob::Hit > &center, const geo::GeometryCore *geo, const detinfo::DetectorPropertiesData &dprop)
 
float Vert2HitDistance (const recob::Hit &hit, const recob::Vertex &vert, const geo::GeometryCore *geo, const detinfo::DetectorPropertiesData &dprop)
 
std::array< float, 2 > Vert2HitCoord (const recob::Vertex &vert, const geo::PlaneID &planeID, const geo::GeometryCore *geo, const detinfo::DetectorPropertiesData &dprop)
 

Function Documentation

std::tuple< std::vector< art::Ptr< recob::Hit > >, std::vector< art::Ptr< recob::Hit > >, bool > sbnpca::GetNearestHits ( const std::vector< art::Ptr< recob::Hit >> &  hits,
int  ihit,
float  distance,
const geo::GeometryCore *  geo,
const detinfo::DetectorPropertiesData &  dprop 
)

Definition at line 39 of file PCA.cc.

42  {
43  std::vector<art::Ptr<recob::Hit>> retlo;
44  std::vector<art::Ptr<recob::Hit>> rethi;
45 
46  bool lo_complete = false;
47  // pull in smaller ones
48  for (int j = ihit-1; j >= 0; j--) {
49  if (HitDistance(*hits[j], *hits[ihit], geo, dprop) > distance) {
50  lo_complete = true;
51  break;
52  }
53  retlo.push_back(hits[j]);
54  }
55 
56  bool hi_complete = false;
57  // pull in larger ones
58  for (unsigned j = ihit+1; j < hits.size(); j++) {
59  if (HitDistance(*hits[j], *hits[ihit], geo, dprop) > distance) {
60  hi_complete = true;
61  break;
62  }
63  rethi.push_back(hits[j]);
64  }
65 
66  return {retlo, rethi, lo_complete && hi_complete};
67 }
icarus::distance
double distance(geo::Point_t const &point, CathodeDesc_t const &cathode)
Returns the distance of a point from the cathode.
Definition: CathodeCrossingUtils.cxx:35
sbnpca::HitDistance
float HitDistance(const recob::Hit &A, const recob::Hit &B, const geo::GeometryCore *geo, const detinfo::DetectorPropertiesData &dprop)
Definition: PCA.cc:34
float sbnpca::HitDistance ( const recob::Hit &  A,
const recob::Hit &  B,
const geo::GeometryCore *  geo,
const detinfo::DetectorPropertiesData &  dprop 
)

Definition at line 34 of file PCA.cc.

34  {
35  std::array<float, 2> vec = HitVector(A, B, geo, dprop);
36  return sqrt(vec[0]*vec[0] + vec[1]*vec[1]);
37 }
HitVector
std::array< float, 2 > HitVector(const recob::Hit &hit, const geo::GeometryCore *geo, const detinfo::DetectorPropertiesData &dprop)
Definition: VertexChargeVacuum_module.cc:115
std::array< float, 2 > sbnpca::HitPCAEigen ( const std::vector< art::Ptr< recob::Hit >> &  hits,
const art::Ptr< recob::Hit > &  center,
const geo::GeometryCore *  geo,
const detinfo::DetectorPropertiesData &  dprop 
)

Definition at line 131 of file PCA.cc.

132  {
133 
134  std::array<float, 2> sum {};
135  for (const art::Ptr<recob::Hit> &h: hits) {
136  std::array<float, 2> vec = HitVector(*h, *center, geo, dprop);
137  sum[0] += vec[0];
138  sum[1] += vec[1];
139  }
140  sum[0] = sum[0] / hits.size();
141  sum[1] = sum[1] / hits.size();
142 
143  std::array<std::array<float, 2>, 2> scatter {};
144  for (const art::Ptr<recob::Hit> &h: hits) {
145  std::array<float, 2> vec = HitVector(*h, *center, geo, dprop);
146  vec[0] -= sum[0];
147  vec[1] -= sum[1];
148 
149  scatter[0][0] += vec[0] * vec[0];
150  scatter[0][1] += vec[0] * vec[1];
151  scatter[1][0] += vec[1] * vec[0];
152  scatter[1][1] += vec[1] * vec[1];
153  }
154 
155  // first get the eigenvalues of the matrix
156  float trace = scatter[0][0] + scatter[1][1];
157  float det = scatter[0][0] * scatter[1][1] - scatter[0][1] * scatter[1][0];
158 
159  float eigenP = (1. / 2.) * (trace + sqrt(trace*trace - 4 * det));
160  float eigenM = (1. / 2.) * (trace - sqrt(trace*trace - 4 * det));
161 
162  return {eigenP, eigenM};
163 }
HitVector
std::array< float, 2 > HitVector(const recob::Hit &hit, const geo::GeometryCore *geo, const detinfo::DetectorPropertiesData &dprop)
Definition: VertexChargeVacuum_module.cc:115
h
while getopts h
Definition: getStatistics.sh:13
trkf::trace
M::value_type trace(const M &m)
Definition: KalmanLinearAlgebra.h:231
std::array< float, 2 > sbnpca::HitPCAVec ( const std::vector< art::Ptr< recob::Hit >> &  hits,
const recob::Hit &  center,
const geo::GeometryCore *  geo,
const detinfo::DetectorPropertiesData &  dprop 
)

Definition at line 69 of file PCA.cc.

70  {
71 
72  return HitPCAVec(hits, HitVector(center, geo, dprop), geo, dprop);
73 }
sbnpca::HitPCAVec
std::array< float, 2 > HitPCAVec(const std::vector< art::Ptr< recob::Hit >> &hits, const recob::Hit &center, const geo::GeometryCore *geo, const detinfo::DetectorPropertiesData &dprop)
Definition: PCA.cc:69
HitVector
std::array< float, 2 > HitVector(const recob::Hit &hit, const geo::GeometryCore *geo, const detinfo::DetectorPropertiesData &dprop)
Definition: VertexChargeVacuum_module.cc:115
std::array< float, 2 > sbnpca::HitPCAVec ( const std::vector< art::Ptr< recob::Hit >> &  hits,
const recob::Vertex &  center,
const geo::GeometryCore *  geo,
const detinfo::DetectorPropertiesData &  dprop 
)

Definition at line 75 of file PCA.cc.

76  {
77 
78  if (hits.size() == 0) return {-100., -100.};
79 
80  return HitPCAVec(hits, Vert2HitCoord(center, hits[0]->WireID(), geo, dprop), geo, dprop);
81 }
sbnpca::HitPCAVec
std::array< float, 2 > HitPCAVec(const std::vector< art::Ptr< recob::Hit >> &hits, const recob::Hit &center, const geo::GeometryCore *geo, const detinfo::DetectorPropertiesData &dprop)
Definition: PCA.cc:69
geo::fhicl::WireID
IDparameter< geo::WireID > WireID
Member type of validated geo::WireID parameter.
Definition: geo_types_fhicl.h:407
sbnpca::Vert2HitCoord
std::array< float, 2 > Vert2HitCoord(const recob::Vertex &vert, const geo::PlaneID &planeID, const geo::GeometryCore *geo, const detinfo::DetectorPropertiesData &dprop)
Definition: PCA.cc:165
std::array< float, 2 > sbnpca::HitPCAVec ( const std::vector< art::Ptr< recob::Hit >> &  hits,
std::array< float, 2 >  center,
const geo::GeometryCore *  geo,
const detinfo::DetectorPropertiesData &  dprop 
)

Definition at line 83 of file PCA.cc.

84  {
85 
86  if (hits.size() < 2) return {-100., -100.};
87 
88  std::array<float, 2> sum {};
89  for (const art::Ptr<recob::Hit> &h: hits) {
90  std::array<float, 2> vec = HitVector(*h, geo, dprop);
91  sum[0] += vec[0] - center[0];
92  sum[1] += vec[1] - center[1];
93  }
94  sum[0] = sum[0] / hits.size();
95  sum[1] = sum[1] / hits.size();
96 
97  std::array<std::array<float, 2>, 2> scatter {};
98  for (const art::Ptr<recob::Hit> &h: hits) {
99  std::array<float, 2> vec = HitVector(*h, geo, dprop);
100  vec[0] -= sum[0] - center[0];
101  vec[1] -= sum[1] - center[1];
102 
103  scatter[0][0] += vec[0] * vec[0];
104  scatter[0][1] += vec[0] * vec[1];
105  scatter[1][0] += vec[1] * vec[0];
106  scatter[1][1] += vec[1] * vec[1];
107  }
108 
109  // first get the eigenvalues of the matrix
110  float trace = scatter[0][0] + scatter[1][1];
111  float det = scatter[0][0] * scatter[1][1] - scatter[0][1] * scatter[1][0];
112 
113  // this is always the max-eigenvalue
114  float eigenP = (1. / 2.) * (trace + sqrt(trace*trace - 4 * det));
115  // float eigenM = (1. / 2.) * (trace - sqrt(trace*trace - 4 * det));
116 
117  // and then the eigenvectors
118  std::array<float, 2> ret {scatter[0][1], eigenP - scatter[0][0]};
119  // std::array<float, 2> eigenVM {scatter[0][1], eigenM - scatter[0][0]};
120 
121  // make sure the sign is right
122  if (sum[0] * ret[0] + sum[1] * ret[1] < 0.) {
123  ret[0] = -ret[0];
124  ret[1] = -ret[1];
125  }
126 
127  return ret;
128 
129 }
HitVector
std::array< float, 2 > HitVector(const recob::Hit &hit, const geo::GeometryCore *geo, const detinfo::DetectorPropertiesData &dprop)
Definition: VertexChargeVacuum_module.cc:115
h
while getopts h
Definition: getStatistics.sh:13
trkf::trace
M::value_type trace(const M &m)
Definition: KalmanLinearAlgebra.h:231
std::array< float, 2 > sbnpca::HitVector ( const recob::Hit &  A,
const geo::GeometryCore *  geo,
const detinfo::DetectorPropertiesData &  dprop 
)

Definition at line 12 of file PCA.cc.

12  {
13  // get the wire distance between A and B
14  float wire_distance = A.WireID().Wire;
15  // convert to cm
16  float wire_distance_cm = wire_distance * geo->WirePitch();
17 
18  // and the time difference
19  float time_distance = A.PeakTime();
20  // convert to cm
21  float time_distance_cm = dprop.ConvertTicksToX(time_distance, A.WireID());
22 
23  return {wire_distance_cm, time_distance_cm};
24 }
recob::Hit::WireID
geo::WireID WireID() const
Definition: Hit.h:233
geo::WireID::Wire
WireID_t Wire
Index of the wire within its plane.
Definition: geo_types.h:580
geo::GeometryCore::WirePitch
geo::Length_t WirePitch(geo::PlaneID const &planeid) const
Returns the distance between two consecutive wires.
Definition: GeometryCore.cxx:873
detinfo::DetectorPropertiesData::ConvertTicksToX
double ConvertTicksToX(double ticks, int p, int t, int c) const
Definition: DetectorPropertiesData.cc:115
recob::Hit::PeakTime
float PeakTime() const
Time of the signal peak, in tick units.
Definition: Hit.h:218
std::array< float, 2 > sbnpca::HitVector ( const recob::Hit &  A,
const recob::Hit &  B,
const geo::GeometryCore *  geo,
const detinfo::DetectorPropertiesData &  dprop 
)

Definition at line 26 of file PCA.cc.

26  {
27  // get each individual vec
28  std::array<float, 2> vecA = HitVector(A, geo, dprop);
29  std::array<float, 2> vecB = HitVector(B, geo, dprop);
30 
31  return {vecA[0] - vecB[0], vecA[1] - vecB[1]};
32 }
HitVector
std::array< float, 2 > HitVector(const recob::Hit &hit, const geo::GeometryCore *geo, const detinfo::DetectorPropertiesData &dprop)
Definition: VertexChargeVacuum_module.cc:115
float sbnpca::VecAngle ( std::array< float, 2 >  A,
std::array< float, 2 >  B 
)

Definition at line 5 of file PCA.cc.

5  {
6  float costh = (A[0] * B[0] + A[1] * B[1]) \
7  / (sqrt(A[0]*A[0] + A[1] * A[1]) * sqrt(B[0]*B[0] + B[1] * B[1]));
8 
9  return acos(costh);
10 }
std::array< float, 2 > sbnpca::Vert2HitCoord ( const recob::Vertex &  vert,
const geo::PlaneID &  planeID,
const geo::GeometryCore *  geo,
const detinfo::DetectorPropertiesData &  dprop 
)

Definition at line 165 of file PCA.cc.

166  {
167 
168  float vert_wire_coord = geo->WireCoordinate(vert.position(), planeID) * geo->WirePitch();
169  float vert_time_coord = vert.position().X();
170 
171  return {vert_wire_coord, vert_time_coord};
172 }
geo::GeometryCore::WireCoordinate
geo::Length_t WireCoordinate(double YPos, double ZPos, geo::PlaneID const &planeid) const
Returns the index of the nearest wire to the specified position.
Definition: GeometryCore.cxx:1217
geo::GeometryCore::WirePitch
geo::Length_t WirePitch(geo::PlaneID const &planeid) const
Returns the distance between two consecutive wires.
Definition: GeometryCore.cxx:873
recob::Vertex::position
const Point_t & position() const
Return vertex 3D position.
Definition: Vertex.h:60
float sbnpca::Vert2HitDistance ( const recob::Hit &  hit,
const recob::Vertex &  vert,
const geo::GeometryCore *  geo,
const detinfo::DetectorPropertiesData &  dprop 
)

Definition at line 174 of file PCA.cc.

174  {
175  float vert_wire_coord = geo->WireCoordinate(vert.position(), hit.WireID()) * geo->WirePitch();
176  float hit_wire_coord = hit.WireID().Wire * geo->WirePitch();
177 
178  float vert_time_coord = vert.position().X();
179  float hit_time_coord = dprop.ConvertTicksToX(hit.PeakTime(), hit.WireID());
180 
181  return sqrt((vert_wire_coord - hit_wire_coord) * (vert_wire_coord - hit_wire_coord) +
182  (vert_time_coord - hit_time_coord) * (vert_time_coord - hit_time_coord));
183 }
geo::GeometryCore::WireCoordinate
geo::Length_t WireCoordinate(double YPos, double ZPos, geo::PlaneID const &planeid) const
Returns the index of the nearest wire to the specified position.
Definition: GeometryCore.cxx:1217
recob::Hit::WireID
geo::WireID WireID() const
Definition: Hit.h:233
geo::WireID::Wire
WireID_t Wire
Index of the wire within its plane.
Definition: geo_types.h:580
geo::GeometryCore::WirePitch
geo::Length_t WirePitch(geo::PlaneID const &planeid) const
Returns the distance between two consecutive wires.
Definition: GeometryCore.cxx:873
detinfo::DetectorPropertiesData::ConvertTicksToX
double ConvertTicksToX(double ticks, int p, int t, int c) const
Definition: DetectorPropertiesData.cc:115
recob::Hit::PeakTime
float PeakTime() const
Time of the signal peak, in tick units.
Definition: Hit.h:218
recob::Vertex::position
const Point_t & position() const
Return vertex 3D position.
Definition: Vertex.h:60
Generated by   doxygen 1.8.5