Line data Source code
1 : //# RFChunkStats.cc: this defines RFChunkStats 2 : //# Copyright (C) 2000,2001,2002,2003 3 : //# Associated Universities, Inc. Washington DC, USA. 4 : //# 5 : //# This library is free software; you can redistribute it and/or modify it 6 : //# under the terms of the GNU Library General Public License as published by 7 : //# the Free Software Foundation; either version 2 of the License, or (at your 8 : //# option) any later version. 9 : //# 10 : //# This library is distributed in the hope that it will be useful, but WITHOUT 11 : //# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or 12 : //# FITNESS FOR A PARTICULAR PURPOSE. See the GNU Library General Public 13 : //# License for more details. 14 : //# 15 : //# You should have received a copy of the GNU Library General Public License 16 : //# along with this library; if not, write to the Free Software Foundation, 17 : //# Inc., 675 Massachusetts Ave, Cambridge, MA 02139, USA. 18 : //# 19 : //# Correspondence concerning AIPS++ should be addressed as follows: 20 : //# Internet email: casa-feedback@nrao.edu. 21 : //# Postal address: AIPS++ Project Office 22 : //# National Radio Astronomy Observatory 23 : //# 520 Edgemont Road 24 : //# Charlottesville, VA 22903-2475 USA 25 : //# 26 : //# $Id$ 27 : #include <casacore/casa/Arrays/ArrayMath.h> 28 : #include <casacore/casa/Exceptions/Error.h> 29 : #include <flagging/Flagging/Flagger.h> 30 : #include <flagging/Flagging/RFChunkStats.h> 31 : #include <msvis/MSVis/VisibilityIterator.h> 32 : #include <msvis/MSVis/VisBuffer.h> 33 : #include <casacore/casa/System/PGPlotter.h> 34 : 35 : 36 : namespace casa { //# NAMESPACE CASA - BEGIN 37 : 38 0 : template<class T> RFlagWord RFChunkStats::getCorrMask ( const casacore::Vector<T> &corrspec ) 39 : { 40 0 : RFlagWord mask=0; 41 : // loop over polzn spec 42 0 : for( casacore::uInt i=0; i<corrspec.nelements(); i++ ) 43 : { 44 : // convert element of polspec to casacore::Stokes type 45 0 : casacore::Stokes::StokesTypes type = casacore::Stokes::type( corrspec(i) ); 46 0 : if( type == casacore::Stokes::Undefined ){ 47 0 : std::ostringstream oss; 48 0 : oss << corrspec(i); 49 0 : throw(casacore::AipsError( casacore::String("Unknown correlation type: ")+ casacore::String(oss))); 50 0 : } 51 : // find this type in current corrarizations 52 0 : casacore::Int icorr = findCorrType(type,corrtypes); 53 0 : if( icorr>=0 ) 54 0 : mask |= (1<<icorr); 55 : } 56 0 : return mask; 57 : } 58 : 59 : } //# NAMESPACE CASA - END 60 :