ATMWaterVaporRadiometer.h

Go to the documentation of this file.

#ifndef _ATM_WATERVAPORRADIOMETER_H
#define _ATM_WATERVAPORRADIOMETER_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: ATMWaterVaporRadiometer.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 "ATMCommon.h"
#include "ATMPercent.h"
#include "ATMTemperature.h"

#include <vector>

using std::vector;

ATM_NAMESPACE_BEGIN

class WaterVaporRadiometer
{
public:

  WaterVaporRadiometer() {}

  WaterVaporRadiometer(const vector<unsigned int> &IdChannels);

  WaterVaporRadiometer(const vector<unsigned int> &IdChannels,
                       const vector<double> &skyCoupling);

  inline WaterVaporRadiometer(const vector<unsigned int> &IdChannels,
                              double skyCoupling)
  {
    WaterVaporRadiometer(IdChannels, vector<double> (IdChannels.size(),
                                                     skyCoupling));
  }

  WaterVaporRadiometer(const vector<unsigned int> &IdChannels,
                       const vector<Percent> &signalGain);

  /* WaterVaporRadiometer(vector<unsigned int> IdChannels, Percent signalGain); */

  WaterVaporRadiometer(const vector<unsigned int> &IdChannels,
                       const vector<double> &skyCoupling,
                       const vector<Percent> &signalGain);

  WaterVaporRadiometer(const vector<unsigned int> &IdChannels,
                       const Temperature &spilloverTemperature);

  WaterVaporRadiometer(const vector<unsigned int> &IdChannels,
                       const vector<double> &skyCoupling,
                       const Temperature &spilloverTemperature);

  WaterVaporRadiometer(const vector<unsigned int> &IdChannels,
                       const vector<Percent> &signalGain,
                       const Temperature &spilloverTemperature);

  WaterVaporRadiometer(const vector<unsigned int> &IdChannels,
                       const vector<double> &skyCoupling,
                       const vector<Percent> &signalGain,
                       const Temperature &spilloverTemperature);

  virtual ~WaterVaporRadiometer();

  inline vector<unsigned int> getIdChannels() const { return IdChannels_; }

  inline vector<double> getSkyCoupling() const { return skyCoupling_; }

  inline void multiplySkyCoupling(double factor)
  {
    for(unsigned int i = 0; i < skyCoupling_.size(); i++) {
      skyCoupling_[i] = skyCoupling_[i] * factor;
    }
  }
  inline void multiplySkyCouplingChannel(unsigned int ichan, double factor)
  {
    skyCoupling_[ichan] = skyCoupling_[ichan] * factor;
  }

  inline void setSkyCoupling(double factor)
  {
    /* for(unsigned int i=0; i<skyCoupling_.size(); i++){ skyCoupling_[i] = factor;} */
    for(vector<double>::iterator it(skyCoupling_.begin()),
        itmax(skyCoupling_.end()); it != itmax; ++it) {
      (*it) = factor;
    }
  }
  inline void setSkyCoupling(unsigned int i, double factor)
  {
    if(i < skyCoupling_.size()) {
      skyCoupling_[i] = factor;
    } else { /* throw(Error()); */
    }
  }

  inline vector<Percent> getsignalGain() const { return signalGain_; }
  inline Temperature getSpilloverTemperature() const { return spilloverTemperature_; }
  inline void setSpilloverTemperature(Temperature spilloverTemperature) { spilloverTemperature_ = spilloverTemperature; }

protected:
  vector<unsigned int> IdChannels_;
  vector<double> skyCoupling_;
  vector<Percent> signalGain_;
  Temperature spilloverTemperature_;

private:
  //     /** Inaccessible Copy Constructor */
  //     WaterVaporRadiometer(const WaterVaporRadiometer&);

  //      /** Inaccessible assignment operator */
  //      WaterVaporRadiometer& operator=(const WaterVaporRadiometer&);


}; // class WaterVaporRadiometer

ATM_NAMESPACE_END

#endif