casa6

//# SimpleSubMS.cc 

//# Copyright (C) 2010

//# Associated Universities, Inc. Washington DC, USA.

//#

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

//# it under the terms of the GNU General Public License as published by

//# the Free Software Foundation; either version 2 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 General Public License for more details.

//# 

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

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

//# Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.

//#

//# Correspondence concerning AIPS++ should be addressed as follows:

//#        Internet email: aips2-request@nrao.edu.

//#        Postal address: AIPS++ Project Office

//#                        National Radio Astronomy Observatory

//#                        520 Edgemont Road

//#                        Charlottesville, VA 22903-2475 USA

//#

//# $Id: $

#include <casacore/casa/System/AppInfo.h>

#include <casacore/casa/Utilities/GenSort.h>

#include <msvis/MSVis/VisBuffer.h>

#include <msvis/MSVis/VisibilityIterator.h>

#include <msvis/MSVis/SimpleSubMS.h>

​

using namespace casacore;

namespace casa {

​

  SimpleSubMS::SimpleSubMS(MeasurementSet& ms) : SubMS(ms){

​

  }

​

  SimpleSubMS::~SimpleSubMS(){

​

  }

​

  MeasurementSet* SimpleSubMS::makeMemSubMS(const MS::PredefinedColumns& whichcol, const String& name){

    LogIO os(LogOrigin("SimpleSubMS", "makeMemSubMS()"));

    if(max(fieldid_p) >= Int(ms_p.field().nrow())){

      os << LogIO::SEVERE 

     << "Field selection contains elements that do not exist in "

     << "this MS"

     << LogIO::POST;

      ms_p=MeasurementSet();

      return 0;  

    }

    if(max(spw_p) >= Int(ms_p.spectralWindow().nrow())){

      os << LogIO::SEVERE 

         << "SpectralWindow selection contains elements that do not exist in "

         << "this MS"

         << LogIO::POST;

      ms_p=MeasurementSet();

      return 0;

    }