ATMWVRMeasurement.h

Go to the documentation of this file.

#ifndef _ATM_WVRMEASUREMENT_H
#define _ATM_WVRMEASUREMENT_H
/*******************************************************************************
 * ALMA - Atacama Large Millimiter Array
 * (c) Instituto de Estructura de la Materia, 2009
 *
 * 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: ATMWVRMeasurement.h Exp $"
 *
 * who       when      what
 * --------  --------  ----------------------------------------------
 * pardo     24/03/09  created
 */

#ifndef __cplusplus
#error "This is a C++ include file and cannot be used from plain C"
#endif

#include "ATMAngle.h"
#include "ATMCommon.h"
#include "ATMLength.h"
#include "ATMTemperature.h"
#include <vector>
#include <math.h>

using std::vector;

ATM_NAMESPACE_BEGIN

class WVRMeasurement
{
public:
  WVRMeasurement();
  WVRMeasurement(const Angle &elevation,
                 const vector<Temperature> &measuredSkyBrightness);
  WVRMeasurement(const Angle &elevation,
                 const vector<Temperature> &measuredSkyBrightness,
                 const vector<Temperature> &fittedSkyBrightness,
                 const Length &retrievedWaterVaporColumn,
                 const Temperature &sigma_fittedSkyBrightness);

  virtual ~WVRMeasurement();

  Angle getElevation() const { return elevation_; }
   double getAirMass() const { return 1.0 / sin(elevation_.get()); }
   vector<Temperature> getmeasuredSkyBrightness() const { return v_measuredSkyBrightness_; }
   vector<Temperature> getfittedSkyBrightness() const { return v_fittedSkyBrightness_; }
   void setfittedSkyBrightness(const vector<Temperature> &a) { v_fittedSkyBrightness_ = a; }
   Length getretrievedWaterVaporColumn() const { return retrievedWaterVaporColumn_; }
   void setretrievedWaterVaporColumn(const Length &a) { retrievedWaterVaporColumn_ = a; }
   Temperature getSigmaFit() const { return sigma_fittedSkyBrightness_; }
   void setSigmaFit(const Temperature &a) { sigma_fittedSkyBrightness_ = a; }

protected:
  Angle elevation_;
  vector<Temperature> v_measuredSkyBrightness_;
  vector<Temperature> v_fittedSkyBrightness_; // !< Fitted sky brightness temperatures over all WVR channels for each event
  Length retrievedWaterVaporColumn_; // !< Retrieved zenith water vapor column for each event
  Temperature sigma_fittedSkyBrightness_; // !< Sigma on the fitted sky brightness temperatures (average sigma over the WVR channels for each event).
}; // class WVRMeasurement

ATM_NAMESPACE_END

#endif