My Project
 All Classes Files Functions Variables Enumerations Enumerator Friends Macros Pages
Public Member Functions | Static Public Member Functions | List of all members
gcp::util::NvssReader Class Reference
Inheritance diagram for gcp::util::NvssReader:
gcp::util::PtSrcFitsReader gcp::util::PtSrcReader

Public Member Functions

 NvssReader (std::string catalogFile)
 
 NvssReader ()
 
virtual ~NvssReader ()
 
void setNvssPsFn (PtSrcReader::Source &src)
 
void setMajAxisError (PtSrcReader::Source &src)
 
void setMinAxisError (PtSrcReader::Source &src)
 
void setPositionAngleError (PtSrcReader::Source &src)
 
void correctForPositionBias (PtSrcReader::Source &src)
 
void setPositionErrors (PtSrcReader::Source &src)
 
void correctForConfusionBias (PtSrcReader::Source &src)
 
void setPeakError (PtSrcReader::Source &src)
 
void setIntegratedFlux (PtSrcReader::Source &src)
 
void setUnresolvedIntegratedFlux (PtSrcReader::Source &src)
 
void setUnresolvedIntegratedFluxError (PtSrcReader::Source &src)
 
void setResolvedIntegratedFlux (PtSrcReader::Source &src)
 
void setResolvedIntegratedFluxError (PtSrcReader::Source &src)
 
void setPartiallyResolvedIntegratedFlux (PtSrcReader::Source &src)
 
void setPartiallyResolvedIntegratedFluxError (PtSrcReader::Source &src, double pseudoPeakJy)
 
double getBiasErrInJy (PtSrcReader::Source &src)
 
void applyCorrections (PtSrcReader::Source &src)
 
void deconvolve (PtSrcReader::Source &src)
 
PtSrcReader::Source parseData ()
 
void readFitsData (long startRow, long nElement)
 
- Public Member Functions inherited from gcp::util::PtSrcFitsReader
 PtSrcFitsReader (std::string catalogFile)
 
 PtSrcFitsReader ()
 
virtual ~PtSrcFitsReader ()
 
void initialize ()
 
void openCatalogFile ()
 
void closeCatalogFile ()
 
PtSrcReader::Source readNextEntry ()
 
void readNextChunk ()
 
bool eof ()
 
void setRaRange (HourAngle &ra, Declination &dec, Angle &radius)
 
void initRange ()
 
- Public Member Functions inherited from gcp::util::PtSrcReader
 PtSrcReader (std::string catalogFile)
 
 PtSrcReader ()
 
void initialize ()
 
virtual ~PtSrcReader ()
 
void setCatalogFile (std::string catalogFile)
 
std::vector< PtSrcReader::Source > findSources (HourAngle ra, Declination dec, Angle radius, Flux fMin=minFlux_, Flux fMax=maxFlux_, bool doPrint=true)
 
unsigned countSources (HourAngle ra, Declination dec, Angle radius, Flux fMin=minFlux_, Flux fMax=maxFlux_)
 
void indexSources ()
 
bool checkAngle (PtSrcReader::Source &src, HourAngle &ra, Declination &dec, Angle &radius)
 
void printHeader (std::ostream &os)
 

Static Public Member Functions

static void setGenericPositionErrors (PtSrcReader::Source &src, Angle &calRaErr, Angle &calDecErr)
 
static double effSnr2 (PtSrcReader::Source &src, double alphaMajor, double alphaMinor)
 

Additional Inherited Members

- Public Attributes inherited from gcp::util::PtSrcReader
HourAngle raMin_
 
HourAngle raMax_
 
- Static Public Attributes inherited from gcp::util::PtSrcReader
static Flux minFlux_
 
static Flux maxFlux_
 
- Protected Member Functions inherited from gcp::util::PtSrcFitsReader
void incrementRange ()
 
- Protected Member Functions inherited from gcp::util::PtSrcReader
void throwCfitsioError (int status)
 
- Protected Attributes inherited from gcp::util::PtSrcFitsReader
fitsfile * fitsFile_
 
int status_
 
long indices_ [25]
 
unsigned nChunk_
 
unsigned iChunk_
 
long nRow_
 
long nRowTotal_
 
unsigned iRow_
 
double ras_ [chunkSize_]
 
double decs_ [chunkSize_]
 
float peakFluxes_ [chunkSize_]
 
float rmsFluxes_ [chunkSize_]
 
float majorAxes_ [chunkSize_]
 
float minorAxes_ [chunkSize_]
 
float positionAngles_ [chunkSize_]
 
char * sourceNames_ [chunkSize_]
 
unsigned int nRange_
 
unsigned int iRange_
 
unsigned rowMin_
 
unsigned int rangeStartInd_ [2]
 
unsigned int rangeStopInd_ [2]
 
- Protected Attributes inherited from gcp::util::PtSrcReader
std::string catalogFile_
 
long sourceIndices_ [25]
 
unsigned nSrc_
 
- Static Protected Attributes inherited from gcp::util::PtSrcFitsReader
static const unsigned chunkSize_ = 10
 

Constructor & Destructor Documentation

NvssReader::NvssReader ( std::string  catalogFile)

Constructor.

....................................................................... Constructor.

NvssReader::NvssReader ( )

....................................................................... Constructor.

NvssReader::~NvssReader ( )
virtual

Destructor.

....................................................................... Destructor.

Member Function Documentation

void NvssReader::applyCorrections ( PtSrcReader::Source &  src)
virtual

....................................................................... NVSS main correction functions (translated from NVSSlist.f)

Reimplemented from gcp::util::PtSrcReader.

void NvssReader::correctForConfusionBias ( PtSrcReader::Source &  src)

....................................................................... NVSS notes say that to correct for confusion bias, the following correction should be applied:

void NvssReader::correctForPositionBias ( PtSrcReader::Source &  src)

....................................................................... Correct for positional biases

void NvssReader::deconvolve ( PtSrcReader::Source &  src)

This subroutine deconvolves a gaussian point-source response from a fitted elliptical gaussian to yield the gaussian source parameters.

This calculation is from AJ, 109, 2318 and Condon's notes of 010209. deconvolve() also determines whether each source axis is significantly (98% confidence = 2.33 sigma) resolved.

If so, the rms error in that axis is calculated. If not, the 98% confidence upper limit to that axis is calculated and its rms error is set to -1. These error calculations are based on the old DECONV subroutine and the NVSS paper (AJ, 115, 1693).

Input arguments:

ptsmaj = major-axis of pt src response (arcsec) ptsmin = minor-axis of pt src response (arcsec) ptspa = position angle of pt src response (DEG) (MUST BE IN RANGE -180 DEG TO +180 DEG) fitmaj = major-axis of raw fit (arcsec) ufitmaj = rms error in fitted major axis (arcsec) fitmin = minor-axis of raw fit (arcsec) ufitmin = rms error in fitted minor axis (arcsec) fitpa = position angle of raw fit (DEG) (MUST BE IN RANGE -180 DEG TO +180 DEG)

Output arguments:

srcmaj = major-axis of deconv source (arcsec) or 98% confidence upper limit usrcmaj = rms error in srcmaj if resolved, or -1 if srcmaj is an upper limit srcmin = minor-axis of deconv source (arcsec) or 98% confidence upper limit usrcmin = rms error in srcmin if resolved, or -1 if srcmaj in an upper limit srcpa = position angle of deconv source (deg) (RETURNED IN RANGE -90 TO +90 deg)

double NvssReader::effSnr2 ( PtSrcReader::Source &  src,
double  alphaMajor,
double  alphaMinor 
)
static

........................................................................ Equation (1) from catalog.ps, for the effective SNR in a given parameter

double NvssReader::getBiasErrInJy ( PtSrcReader::Source &  src)

....................................................................... Return the appropriate bias error for D or DnC configuration

PtSrcReader::Source NvssReader::parseData ( )
virtual

....................................................................... Read the next entry from the catalog file

Implements gcp::util::PtSrcFitsReader.

void NvssReader::readFitsData ( long  startRow,
long  nElement 
)
virtual

....................................................................... Read the next chunk of data from the FITS file

Data structure from the NVSS catalog.ps file:

1 RA(2000) J2000 Right Ascension (degrees) 2 DEC(2000) J2000 declination (degrees) 3 PEAK INT Peak Stokes I (Jy/beam) 4 MAJOR AX Fitted major axis (degrees) 5 MINOR AX Fitted minior axis (degrees) 6 POSANGLE Fitted position angle of the major axis (degrees) 7 Q CENTER Interpolated Q value at position of I peak (Jy/beam) 8 U CENTER Interpolated U value at position of I peak (Jy/beam) 9 P FLUX Integrated polarized (linear) flux (Jy) 10 I RMS RMS noise in Stokes I image (Jy/beam) 11 POL RMS RMS noise in Stokes Q,U image (Jy/beam) 12 RES RMS RMS Stokes I residual after fit 13 RES PEAK Peak Stokes I residual (Jy/beam) 14 RES FLUX Integrated Stokes I residual (Jy) 15 CENTER X X pixel in FIELD of center (pixels) 16 CENTER Y Y pixel in FIELD of center (pixels) 17 FIELD Name of the 4x4 degree field (string) 18 JD PROCESSED Julian date on which this entry was derived from an image (day)

Implements gcp::util::PtSrcFitsReader.

void NvssReader::setGenericPositionErrors ( PtSrcReader::Source &  src,
Angle &  calRaErr,
Angle &  calDecErr 
)
static

....................................................................... Set the error in the fitted RA and DEC (see catalog.ps, section 6.8 & 6.9)

void NvssReader::setIntegratedFlux ( PtSrcReader::Source &  src)

....................................................................... Calculate the integrated flux

void NvssReader::setMajAxisError ( PtSrcReader::Source &  src)

....................................................................... Set the error in the fitted major axis size (see catalog.ps, section 6.2)

void NvssReader::setMinAxisError ( PtSrcReader::Source &  src)

....................................................................... Set the error in the fitted minor axis size (see catalog.ps, section 6.2)

void NvssReader::setNvssPsFn ( PtSrcReader::Source &  src)

....................................................................... NVSS point-spread function calculation (translated from NVSSlist.f)

void NvssReader::setPartiallyResolvedIntegratedFlux ( PtSrcReader::Source &  src)

....................................................................... Set the integrated flux and error for a partially resolved source

void NvssReader::setPartiallyResolvedIntegratedFluxError ( PtSrcReader::Source &  src,
double  pseudoPeakJy 
)

....................................................................... Set the error in the integrated flux for an unresolved source (see catalog.ps, section 6.5)

void NvssReader::setPeakError ( PtSrcReader::Source &  src)

....................................................................... Set the error in the fitted peak flux value (see catalog.ps, section 6.1)

void NvssReader::setPositionAngleError ( PtSrcReader::Source &  src)

....................................................................... Set the error in the fitted minor axis size (see catalog.ps, section 6.4)

void NvssReader::setPositionErrors ( PtSrcReader::Source &  src)

....................................................................... Set the error in the fitted RA and DEC (see catalog.ps, section 6.8 & 6.9)

void NvssReader::setResolvedIntegratedFlux ( PtSrcReader::Source &  src)

....................................................................... Set the integrated flux and error for an unresolved source

void NvssReader::setResolvedIntegratedFluxError ( PtSrcReader::Source &  src)

....................................................................... Set the error in the integrated flux for an unresolved source (see catalog.ps, section 6.5)

void NvssReader::setUnresolvedIntegratedFlux ( PtSrcReader::Source &  src)

....................................................................... Set the integrated flux and error for an unresolved source

void NvssReader::setUnresolvedIntegratedFluxError ( PtSrcReader::Source &  src)

....................................................................... Set the error in the integrated flux for an unresolved source (see catalog.ps, section 6.6)

The documentation for this class was generated from the following files: