Commits

Renaud Miel authored 2a37507e746
Fixed functional bug when plotting RA(t)

* Refactored PlotMSIndexer constructor * Added PlotMSIndexerFactory class
No tags

code/plotms/Data/MSCache.cc

Modified
1918 1918 //auto = interpolator.pointingDir(antId,vbTime);
1919 1919 dir1Vec[k] = (piTvi_->getPointingAngle(antId,vbTime)).second;
1920 1920 }
1921 1921
1922 1922 Matrix<Double> raDecMat;
1923 1923 raDecMat.resize(2, dir1Vec.nelements());
1924 1924 auto nAnts = dir1Vec.nelements();
1925 1925 for (decltype(nAnts) iant = 0; iant < nAnts; ++iant) {
1926 1926 raDecMat.column(iant) = dir1Vec(iant).getAngle("deg").getValue();
1927 1927 }
1928 - *ra_[vbnum] = raDecMat.row(0);
1929 - *dec_[vbnum] = raDecMat.row(1);
1928 + //*ra_[vbnum] = raDecMat.row(0);
1929 + //*dec_[vbnum] = raDecMat.row(1);
1930 1930 *((*loadRa_)[vbnum]) = raDecMat.row(0);
1931 1931 *((*loadDec_)[vbnum]) = raDecMat.row(1);
1932 1932 break;
1933 1933 }
1934 1934 case PMS::RADIAL_VELOCITY: {
1935 1935 Int fieldId = vb->fieldId()(0);
1936 1936 const ROMSFieldColumns& fieldColumns = vi_p->subtableColumns().field();
1937 1937 MRadialVelocity radVelocity = fieldColumns.radVelMeas(fieldId, vb->time()(0));
1938 1938 radialVelocity_(vbnum) = radVelocity.get("AU/d").getValue( "km/s");
1939 1939 break;

Everything looks good. We'll let you know here if there's anything you should know about.

Add shortcut