PixelSeeding/test/PixelTriplets_InvPrbl_prec.cpp

0001 #include "RecoTracker/PixelSeeding/plugins/ThirdHitPredictionFromInvParabola.cc"
0002
0003 #include <iostream>
0004 #include <string>
0005
0006 namespace test {
0007   namespace PixelTriplets_InvPrbl_prec {
0008     int test() {
0009       std::string c("++Constr");
0010       std::string r("++R");
0011       std::string a;
0012       double par[7];
0013       double q[2];
0014       while (std::cin) {
0015         ThirdHitPredictionFromInvParabola pred;
0016         std::cin >> a;
0017         if (a == c) {
0018           for (auto& p : par)
0019             std::cin >> p;
0020           pred = ThirdHitPredictionFromInvParabola(par[0], par[1], par[2], par[3], par[4], par[5], par[6]);
0021           std::cout << "ip min, max " << pred.theIpRangePlus.min() << " " << pred.theIpRangePlus.max() << "  "
0022                     << pred.theIpRangeMinus.min() << " " << pred.theIpRangeMinus.max() << std::endl;
0023         } else if (a == r) {
0024           std::cin >> q[0] >> q[1];
0025           {
0026             auto rp = pred.rangeRPhi(q[0], 1);
0027             auto rn = pred.rangeRPhi(q[0], -1);
0028             std::cout << "range " << rp.min() << " " << rp.max() << " " << rn.min() << " " << rn.max() << std::endl;
0029           }
0030           {
0031             auto rp = pred.rangeRPhi(q[1], 1);
0032             auto rn = pred.rangeRPhi(q[1], -1);
0033             std::cout << "range " << rp.min() << " " << rp.max() << " " << rn.min() << " " << rn.max() << std::endl;
0034           }
0035         }
0036       }
0037       return 0;
0038     }
0039   }  // namespace PixelTriplets_InvPrbl_prec
0040 }  // namespace test
0041
0042 int main() { return test::PixelTriplets_InvPrbl_prec::test(); }