ViImplementation2.h
Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028 #if ! defined (MSVIS_ViImplementation2_H_121115_0950)
00029 #define MSVIS_ViImplementation2_H_121115_0950
00030
00031 #include <casa/aips.h>
00032 #include <casa/BasicSL.h>
00033 #include <msvis/MSVis/VisBufferComponents2.h>
00034 #include <measures/Measures/MFrequency.h>
00035 #include <measures/Measures/Stokes.h>
00036
00037 #include <map>
00038 #include <vector>
00039
00040 namespace casa {
00041
00042 template <typename T> class Array;
00043 template <typename T> class Block;
00044
00045 template <typename T> class Cube;
00046 template <typename T> class Matrix;
00047 class MDirection;
00048 class MeasurementSet;
00049 class MEpoch;
00050 class MPosition;
00051 class MSDerivedValues;
00052 class RecordInterface;
00053 template <typename T, Int n> class RigidVector;
00054 class Slice;
00055 class String;
00056 template <typename T, Int n> class SquareMatrix;
00057 template <typename T> class Vector;
00058 class VisImagingWeight;
00059
00060 namespace vi {
00061
00062
00063
00064
00065 class ChannelSelector;
00066 class ChannelSelectorCache;
00067 class FrequencySelections;
00068 class SortColumns;
00069 class SpectralWindowChannels;
00070 class SpectralWindowChannelsCache;
00071 class Subchunk;
00072 class SubtableColumns;
00073 class VisBuffer2;
00074 class VisibilityIterator2;
00075 class WeightScaling;
00076 enum VisBufferType : int;
00077 enum VisBufferOptions : int;
00078
00079
00080
00081
00082
00083
00084
00085
00086
00087
00088
00089
00090
00091
00092
00093
00094
00095
00096
00097
00098
00099
00100
00101
00102
00103
00104
00105
00106
00107
00108
00109
00110
00111
00112
00113
00114
00115
00116
00117
00118
00119
00120
00121
00122
00123
00124
00125
00126
00127
00128
00129 class ViImplementation2 {
00130
00131 friend class VisibilityIterator2;
00132
00133 public:
00134
00135
00136 ViImplementation2( const ViImplementation2& ) = delete;
00137 ViImplementation2& operator=( const ViImplementation2& ) = delete;
00138 ViImplementation2 () {}
00139
00140
00141
00142 virtual ~ViImplementation2 () {}
00143
00144
00145 virtual String ViiType() const = 0;
00146
00147
00148
00149
00150
00151
00152
00153
00154
00155 virtual void origin () = 0;
00156 virtual Bool more () const = 0;
00157 virtual void next () = 0;
00158 virtual Subchunk getSubchunkId () const = 0;
00159
00160
00161
00162 virtual void originChunks (Bool forceRewind = False) = 0;
00163 virtual Bool moreChunks () const = 0;
00164 virtual void nextChunk () = 0;
00165
00166 virtual Bool isWritable () const = 0;
00167
00168
00169
00170
00171
00172
00173 virtual Double getInterval() const = 0;
00174 virtual void setInterval (Double timeInterval) = 0;
00175
00176
00177
00178
00179 virtual void setFrequencySelections (const FrequencySelections & selection) = 0;
00180
00181
00182
00183
00184
00185
00186
00187 virtual void setRowBlocking (Int nRows = 0) = 0;
00188
00189 virtual Bool existsColumn (VisBufferComponent2 id) const = 0;
00190
00191 virtual const SortColumns & getSortColumns() const = 0;
00192
00193 virtual Bool isNewArrayId () const = 0;
00194 virtual Bool isNewFieldId () const = 0;
00195 virtual Bool isNewMs () const = 0;
00196 virtual Bool isNewSpectralWindow () const = 0;
00197
00198
00199
00200 virtual Int nRows () const = 0;
00201
00202
00203
00204
00205
00206 virtual void getRowIds (Vector<uInt> & rowids) const = 0;
00207
00208 virtual VisBuffer2 * getVisBuffer (const VisibilityIterator2 *) = 0;
00209 virtual VisBuffer2 * getVisBuffer () = 0;
00210
00211
00212
00213
00214
00215
00216
00217
00218
00219
00220 virtual void antenna1 (Vector<Int> & ant1) const = 0;
00221
00222
00223
00224 virtual void antenna2 (Vector<Int> & ant2) const = 0;
00225
00226
00227
00228 virtual void corrType (Vector<Int> & corrTypes) const = 0;
00229
00230
00231
00232 virtual Int dataDescriptionId () const = 0;
00233
00234 virtual void dataDescriptionIds (Vector<Int> & ddis) const = 0;
00235
00236
00237
00238 virtual void exposure (Vector<Double> & expo) const = 0;
00239
00240
00241
00242 virtual void feed1 (Vector<Int> & fd1) const = 0;
00243
00244
00245
00246 virtual void feed2 (Vector<Int> & fd2) const = 0;
00247
00248
00249
00250 virtual void fieldIds (Vector<Int>&) const = 0;
00251
00252
00253
00254 virtual void arrayIds (Vector<Int>&) const = 0;
00255
00256
00257
00258 virtual String fieldName () const = 0;
00259
00260
00261
00262 virtual void flag (Cube<Bool> & flags) const = 0;
00263
00264
00265
00266 virtual void flag (Matrix<Bool> & flags) const = 0;
00267
00268
00269
00270 virtual Bool flagCategoryExists () const = 0;
00271
00272
00273
00274 virtual void flagCategory (Array<Bool> & flagCategories) const = 0;
00275
00276
00277
00278 virtual void flagRow (Vector<Bool> & rowflags) const = 0;
00279
00280
00281
00282 virtual void observationId (Vector<Int> & obsids) const = 0;
00283
00284
00285
00286 virtual Int polarizationId () const = 0;
00287
00288
00289
00290 virtual void processorId (Vector<Int> & procids) const = 0;
00291
00292
00293
00294 virtual void scan (Vector<Int> & scans) const = 0;
00295
00296
00297
00298 virtual String sourceName () const = 0;
00299
00300
00301
00302 virtual void stateId (Vector<Int> & stateids) const = 0;
00303
00304
00305
00306
00307 virtual void jonesC (Vector<SquareMatrix<Complex, 2> > & cjones) const = 0;
00308
00309
00310
00311 virtual Int polFrame () const = 0;
00312
00313
00314
00315 virtual void sigma (Matrix<Float> & sigmat) const = 0;
00316
00317
00318
00319 virtual Int spectralWindow () const = 0;
00320
00321 virtual void spectralWindows (Vector<Int> & spws) const = 0;
00322
00323
00324
00325 virtual void time (Vector<Double> & t) const = 0;
00326
00327
00328
00329 virtual void timeCentroid (Vector<Double> & t) const = 0;
00330
00331
00332
00333 virtual void timeInterval (Vector<Double> & ti) const = 0;
00334
00335
00336
00337 virtual void uvw (Matrix<Double> & uvwmat) const = 0;
00338
00339
00340
00341 virtual void visibilityCorrected (Cube<Complex> & vis) const = 0;
00342 virtual void visibilityModel (Cube<Complex> & vis) const = 0;
00343 virtual void visibilityObserved (Cube<Complex> & vis) const = 0;
00344
00345
00346
00347 virtual void floatData (Cube<Float> & fcube) const = 0;
00348
00349
00350
00351
00352
00353
00354
00355
00356
00357
00358
00359 virtual IPosition visibilityShape () const = 0;
00360
00361
00362
00363 virtual void weight (Matrix<Float> & wtmat) const = 0;
00364
00365
00366
00367 virtual Bool weightSpectrumExists () const = 0;
00368
00369
00370
00371 virtual Bool sigmaSpectrumExists () const = 0;
00372
00373
00374
00375 virtual void weightSpectrum (Cube<Float> & wtsp) const = 0;
00376
00377
00378
00379 virtual void sigmaSpectrum (Cube<Float> & wtsp) const = 0;
00380
00381
00382 virtual void setWeightScaling (CountedPtr<WeightScaling> weightscaling) = 0;
00383 virtual Bool hasWeightScaling () const = 0;
00384 virtual CountedPtr<WeightScaling> getWeightScaling () const = 0;
00385
00386
00387
00388
00389
00390
00391
00392
00393
00394
00395
00396 virtual Bool allBeamOffsetsZero () const = 0;
00397
00398
00399
00400 virtual MDirection azel0 (Double time) const = 0;
00401 static void azel0Calculate (Double time, MSDerivedValues & msd,
00402 MDirection & azel0, const MEpoch & mEpoch0);
00403
00404 virtual const Vector<MDirection> & azel (Double time) const = 0;
00405 static void azelCalculate (Double time, MSDerivedValues & msd, Vector<MDirection> & azel,
00406 Int nAnt, const MEpoch & mEpoch0);
00407
00408
00409
00410 virtual const Vector<Float> & feed_pa (Double time) const = 0;
00411 static Vector<Float> feed_paCalculate (Double time, MSDerivedValues & msd,
00412 Int nAntennas, const MEpoch & mEpoch0,
00413 const Vector<Float> & receptor0Angle);
00414
00415
00416
00417
00418
00419
00420 virtual const Cube<RigidVector<Double, 2> > & getBeamOffsets () const = 0;
00421
00422
00423
00424 virtual Double hourang (Double time) const = 0;
00425 static Double hourangCalculate (Double time, MSDerivedValues & msd, const MEpoch & mEpoch0);
00426
00427
00428
00429
00430
00431 virtual const Float & parang0 (Double time) const = 0;
00432 static Float parang0Calculate (Double time, MSDerivedValues & msd, const MEpoch & epoch0);
00433
00434
00435
00436 virtual const Vector<Float> & parang (Double time) const = 0;
00437 static Vector<Float> parangCalculate (Double time, MSDerivedValues & msd,
00438 int nAntennas, const MEpoch mEpoch0);
00439
00440
00441
00442 virtual const MDirection & phaseCenter () const = 0;
00443
00444
00445
00446
00447
00448
00449
00450
00451 virtual const Cube<Double> & receptorAngles () const = 0;
00452
00453
00454
00455
00456
00457
00458
00459
00460
00461 virtual const Vector<String> & antennaMounts () const = 0;
00462
00463 virtual MEpoch getEpoch () const = 0;
00464
00465
00466
00467
00468 virtual const VisImagingWeight & getImagingWeightGenerator () const = 0;
00469
00470 virtual MFrequency::Types getObservatoryFrequencyType () const = 0;
00471 virtual MPosition getObservatoryPosition () const = 0;
00472 virtual Vector<Float> getReceptor0Angle () = 0;
00473
00474 virtual Int getReportingFrameOfReference () const = 0;
00475 virtual void setReportingFrameOfReference (Int frame) = 0;
00476
00477 virtual Vector<Int> getCorrelations () const = 0;
00478 virtual Vector<Stokes::StokesTypes> getCorrelationTypesDefined () const = 0;
00479 virtual Vector<Stokes::StokesTypes> getCorrelationTypesSelected () const = 0;
00480
00481 virtual Vector<Int> getChannels (Double time, Int frameOfReference, Int spectralWndow = -1,
00482 Int msId = -1) const = 0;
00483 virtual Vector<Double> getFrequencies (Double time, Int frameOfReference, Int spectralWndow = -1,
00484 Int msId = -1) const = 0;
00485
00486
00487
00488
00489 virtual Int msId () const = 0;
00490 virtual const MeasurementSet & ms () const = 0;
00491 virtual Int getNMs () const = 0;
00492
00493
00494 virtual String msName() const = 0;
00495
00496
00497
00498
00499
00500
00501
00502
00503
00504
00505
00506 virtual void slurp () const = 0;
00507
00508
00509
00510 virtual const vi::SubtableColumns & subtableColumns () const = 0;
00511
00512
00513
00514
00515 virtual const SpectralWindowChannels & getSpectralWindowChannels (Int msId, Int spectralWindowId) const = 0;
00516
00517
00518
00519 virtual void useImagingWeight (const VisImagingWeight & imWgt) = 0;
00520
00521
00522
00523 virtual Int nAntennas () const = 0;
00524 virtual Int nDataDescriptionIds () const = 0;
00525 virtual Int nPolarizationIds () const = 0;
00526 virtual Int nRowsInChunk () const = 0;
00527 virtual Int nRowsViWillSweep () const = 0;
00528 virtual Int nSpectralWindows () const = 0;
00529
00530
00531
00532
00533
00534
00535
00536
00537
00538
00539 virtual void writeBackChanges (VisBuffer2 * vb) = 0;
00540
00541
00542
00543
00544
00545
00546
00547
00548
00549
00550
00551 virtual void writeFlag (const Cube<Bool> & flag) = 0;
00552
00553
00554 virtual void writeFlagRow (const Vector<Bool> & rowflags) = 0;
00555
00556 virtual void writeFlagCategory(const Array<Bool>& fc) = 0;
00557
00558
00559
00560
00561
00562
00563
00564
00565
00566
00567
00568
00569
00570 virtual void writeVisCorrected (const Cube<Complex> & vis) = 0;
00571 virtual void writeVisModel (const Cube<Complex> & vis) = 0;
00572 virtual void writeVisObserved (const Cube<Complex> & vis) = 0;
00573
00574
00575 virtual void writeWeight (const Matrix<Float> & wt) = 0;
00576
00577
00578
00579
00580
00581 virtual void writeWeightSpectrum (const Cube<Float> & wtsp) = 0;
00582
00583
00584 virtual void initWeightSpectrum (const Cube<Float> & ) {};
00585
00586
00587 virtual void writeSigmaSpectrum (const Cube<Float> & sigsp) = 0;
00588
00589
00590 virtual void initSigmaSpectrum (const Cube<Float> & ) {};
00591
00592
00593 virtual void writeSigma (const Matrix<Float> & sig) = 0;
00594
00595
00596
00597
00598
00599
00600 virtual void writeModel(const RecordInterface& rec, Bool iscomponentlist=True,
00601 Bool incremental=False) = 0;
00602
00603
00604 protected:
00605
00606 void associateVbWithVi2 (VisBuffer2 * vb, const VisibilityIterator2 * vi);
00607 VisBuffer2 * createAttachedVisBuffer (VisBufferType t, VisBufferOptions options);
00608
00609
00610 static void doWeightScaling (Bool hasWeightScaling,
00611 WeightScaling * scaling,
00612 const Array<Float>& unscaled,
00613 Array<Float>& scaled);
00614
00615 };
00616
00617 }
00618
00619 }
00620
00621 #endif // ! defined (MSVIS_ViImplementation2_H_121115_0950)
00622
00623