spectral_bound.h

Go to the documentation of this file.

#ifndef SPECTRAL_BOUND_H
#define SPECTRAL_BOUND_H
#include "printable.h"
#include <vector>
#include "fp_exception.h"
#include "double_with_unit.h"
#include "array_with_unit.h"

namespace FullPhysics {
/****************************************************************/
class SpectralBound : public Printable<SpectralBound> {
public:
//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  SpectralBound() { }

  SpectralBound(const std::vector<DoubleWithUnit>& Lower_bound,
                const std::vector<DoubleWithUnit>& Upper_bound);
  SpectralBound(const ArrayWithUnit<double, 2>& Bound);
  virtual ~SpectralBound() {}

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  int number_spectrometer() const { return (int) lower_b.size(); }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  DoubleWithUnit center(int Spec_index, const Unit& U) const
  {
    return ((lower_bound(Spec_index) + upper_bound(Spec_index)) / 2.0).
      convert_wave(U);
  }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  DoubleWithUnit lower_bound(int Spec_index) const 
  {
    range_check(Spec_index, 0, number_spectrometer());
    return lower_b[Spec_index];
  }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  DoubleWithUnit lower_bound(int Spec_index, const Unit& U) const
  {
    range_check(Spec_index, 0, number_spectrometer());
    return DoubleWithUnit(std::min(upper_b[Spec_index].convert_wave(U).value,
                                   lower_b[Spec_index].convert_wave(U).value),
                          U);
  }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  DoubleWithUnit upper_bound(int Spec_index) const 
  {
    range_check(Spec_index, 0, number_spectrometer());
    return upper_b[Spec_index];
  }

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

  DoubleWithUnit upper_bound(int Spec_index, const Unit& U) const
  {
    range_check(Spec_index, 0, number_spectrometer());
    return DoubleWithUnit(std::max(upper_b[Spec_index].convert_wave(U).value,
                                   lower_b[Spec_index].convert_wave(U).value),
                          U);
  }

  int spectral_index(const DoubleWithUnit& W) const;
  virtual void print(std::ostream& Os) const;
private:
  std::vector<DoubleWithUnit> lower_b, upper_b;
};
}
#endif