casa6

//#  CASA - Common Astronomy Software Applications (http://casa.nrao.edu/)

//#  Copyright (C) Associated Universities, Inc. Washington DC, USA 2011, All rights reserved.

//#  Copyright (C) European Southern Observatory, 2011, All rights reserved.

//#

//#  This library is free software; you can redistribute it and/or

//#  modify it under the terms of the GNU Lesser General Public

//#  License as published by the Free software Foundation; either

//#  version 2.1 of the License, or (at your option) any later version.

//#

//#  This library is distributed in the hope that it will be useful,

//#  but WITHOUT ANY WARRANTY, without even the implied warranty of

//#  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU

//#  Lesser General Public License for more details.

//#

//#  You should have received a copy of the GNU Lesser General Public

//#  License along with this library; if not, write to the Free Software

//#  Foundation, Inc., 59 Temple Place, Suite 330, Boston,

//#  MA 02111-1307  USA

​

#ifndef STATWTTVI_H_

#define STATWTTVI_H_

​

#include <msvis/MSVis/TransformingVi2.h>

​

#include <casacore/casa/aipstype.h>

#include <casacore/casa/Arrays/ArrayIter.h>

#include <casacore/casa/BasicSL/Complex.h>

#include <casacore/ms/MSSel/MSSelection.h>

#include <casacore/scimath/StatsFramework/StatisticsAlgorithmFactory.h>

​

#include <msvis/MSVis/VisBuffer2.h>

#include <msvis/MSVis/VisibilityIterator2.h>

#include <mstransform/TVI/StatWtClassicalDataAggregator.h>

#include <mstransform/TVI/StatWtTypes.h>

#include <mstransform/TVI/UtilsTVI.h>

#include <stdcasa/variant.h>

#include <stdcasa/StdCasa/CasacSupport.h>

​

#include <map>

#include <vector>

#include <utility>

#include <set>

#include <memory>

​

namespace casa {

​

class StatWtVarianceAndWeightCalculator;

​

namespace vi {

​

class StatWtTVI : public TransformingVi2 {

​

public:

​

    static const casacore::String CHANBIN;

​

    // The following fields are supported in the input configuration record

    // combine           String, if contains "corr", data will be aggregated

    //                   across correlations.

    // value in CHANBIN: Int or Quantity String, describes channel bin widths

    //                   in which to aggregate data within spectral windows

    //                   (spw boundaries are not crossed). If not supplied,

    //                   data for all channels in each spectral window are

    //                   aggregated.

    // minsamp:          Int, minimum number of samples required in an

    //                   aggregated set, if less than that, stats are not

    //                   computed and the data in the sample are flagged. If

    //                   not supplied, 2 is used.

    // statalg           String representing what statistics algorithm to use.

    //                   "cl", "ch", "f", "h".

    // maxiter           Int max number of iterations for Chauvenet algorithm

    // zscore            Double zscore for Chauvenet algorithm

    // center            String center for FitToHalf algorithm, "mean",

    //                   "median", or "zero"

    // lside             Bool side to use for FitToHalf algorithm, True means

    //                   <= center side.

    // fence             Double fence value for HingesFences algorithm

    // wtrange           Zero or two element Array<Double>. Specifies the range

    //                   of "good"

    //                   weight values. Data with weights computed to be outside

    //                   this range will be flagged. Both elements must be

    //                   non-negative. If zero length, all weights are

    //                   acceptable.

    // fitspw            String. MSSelection string representing channels to

    //                   exclude from weight computation.

    // datacolumn        String. Data column to use for computing weights.

    //                   Supports 'data' or 'corrected'. Minimum match, case

    //                   insensitive. If not provided. 'corrected' is used.

    // slidetimebin      Bool. If true, use a sliding window for binning in

    //                   time.

    // timebin           Double. Width of sliding time window. Not used if