PixelLowPtUtilities/test/LowPtGlobalDir_t.cpp

0001 #include "RecoTracker/PixelLowPtUtilities/src/LowPtClusterShapeSeedComparitor.cc"
0002 #include <iostream>
0003
0004 #ifdef OLDCODE
0005 void go(std::vector<GlobalPoint> g) {
0006   std::vector<GlobalVector> v = LowPtClusterShapeSeedComparitor::getGlobalDirs(g);
0007   for (int i = 0; i != 3; i++)
0008     std::cout << g[i] << std::endl;
0009   if (v.empty())
0010     std::cout << "failed" << std::endl;
0011   else
0012     for (int i = 0; i != 3; i++) {
0013       std::cout << v[i];
0014       std::cout << " norm " << v[i].mag2() << std::endl;
0015     }
0016 }
0017
0018 int main() {
0019   std::vector<GlobalPoint> g1({GlobalPoint(1., 3., 0.), GlobalPoint(-1., 1., 1.), GlobalPoint(1., -1., 2.)});
0020
0021   std::vector<GlobalPoint> g2({GlobalPoint(1., -1., 0.), GlobalPoint(-1., 1., 1.), GlobalPoint(1., -3., 2.)});
0022
0023   std::vector<GlobalPoint> g3({GlobalPoint(1., 3., 0.), GlobalPoint(4., 1., -1.), GlobalPoint(1., -1., -2.)});
0024
0025   go(g1);
0026   go(g2);
0027   go(g3);
0028
0029   return 0;
0030 }
0031
0032 #else
0033 void go(GlobalPoint const* g) {
0034   GlobalVector v[3];
0035   bool ok = getGlobalDirs(g, v);
0036   for (int i = 0; i != 3; i++)
0037     std::cout << g[i] << std::endl;
0038   if (!ok)
0039     std::cout << "failed" << std::endl;
0040   else
0041     for (int i = 0; i != 3; i++) {
0042       std::cout << v[i];
0043       std::cout << " norm " << v[i].mag2() << std::endl;
0044     }
0045 }
0046
0047 int main() {
0048   GlobalPoint g1[3] = {GlobalPoint(1., 3., 0.), GlobalPoint(-1., 1., 1.), GlobalPoint(1., -1., 2.)};
0049
0050   GlobalPoint g2[3] = {GlobalPoint(1., -1., 0.), GlobalPoint(-1., 1., 1.), GlobalPoint(1., -3., 2.)};
0051
0052   GlobalPoint g3[3] = {GlobalPoint(1., 3., 0.), GlobalPoint(4., 1., -1.), GlobalPoint(1., -1., -2.)};
0053
0054   go(g1);
0055   go(g2);
0056   go(g3);
0057
0058   return 0;
0059 }
0060 #endif