casa6

//# UtilsTVI.h: declaration of the UVContSubResult class

//#

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

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

//#  Copyright (C) European Southern Observatory, 2021, 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

//# $Id: $

​

#ifndef UVCContSubResult_H_

#define UVCContSubResult_H_

​

// casacore types and containers

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

#include <casacore/casa/Containers/Record.h>

​

#include <unordered_map>

​

namespace casa { //# NAMESPACE CASA - BEGIN

​

namespace vi { //# NAMESPACE VI - BEGIN

​

using std::unordered_map;

​

/*

 * To accumulate results from the fits calculated by the

 * UVContSubTVI. For every row and polarization a fit result should be

 * added.

 */

class UVContSubResult

{

 public:

  // add/accumulate one result, for every row-polarization fit

  void addOneFit(int field, int scan, int spw, int pol, casacore::Complex chiSquared);

​

  // Produces a record with the information accumulated throughout iteration/fits

  casacore::Record getAccumulatedResult() const;

​

 private:

  struct FitResultAcc {

    size_t count;

    casacore::Complex chiSqAvg;

    casacore::Complex chiSqMin;

    casacore::Complex chiSqMax;

  };

​

  // chi_square per field, scan, spw, polarization, real/imag

  unordered_map<int,

    unordered_map<int, unordered_map<int, unordered_map<int, FitResultAcc>>>> accum;

};