File indexing completed on 2024-04-06 12:28:36
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 }
0040 }
0041
0042 int main() { return test::PixelTriplets_InvPrbl_prec::test(); }