File indexing completed on 2024-04-06 12:31:41
0001 #include "TrackingTools/TrajectoryState/interface/FreeTrajectoryState.h"
0002 #include <ostream>
0003
0004 using namespace std;
0005
0006 ostream& operator<<(ostream& os, const FreeTrajectoryState& fts) {
0007 os << "parameters" << endl;
0008 {
0009 const AlgebraicVector6& v = fts.parameters().vector();
0010 os << "x = ";
0011 {
0012 for (int i = 0; i < 3; i++) {
0013 os.precision(6);
0014 os.width(13);
0015 os << v[i];
0016 }
0017 }
0018 os << endl;
0019 os << "p = ";
0020 {
0021 for (int i = 3; i < 6; i++) {
0022 os.precision(6);
0023 os.width(13);
0024 os << v[i];
0025 }
0026 }
0027 os << endl;
0028 }
0029 if (fts.hasError()) {
0030 os << "error" << endl;
0031 {
0032 const AlgebraicSymMatrix55& m = fts.curvilinearError().matrix();
0033 for (int i = 0; i < 5; i++) {
0034 for (int j = 0; j < 5; j++) {
0035 os.precision(6);
0036 os.width(13);
0037 os << m(i, j);
0038 }
0039 os << endl;
0040 }
0041 }
0042 } else {
0043 os << "no error defined." << endl;
0044 }
0045 return os;
0046 }