RectangleFilter.h

Go to the documentation of this file.

#ifndef _utl_RectangleFilter_h_
#define _utl_RectangleFilter_h_

#include <complex>
#include <utl/AnalyticFilter.h>
#include <utl/AugerException.h>
#include <utl/AugerUnits.h>
#include <utl/MathConstants.h>

namespace utl {

  class RectangleFilter : public AnalyticFilter {
    public: 
      RectangleFilter(double lowerCutoff, double upperCutoff, double delay);
      virtual ~RectangleFilter();
      virtual std::complex<double> TransmissivityAtFrequency(double parFrequency) const;
      virtual double EffectiveBandwidth() const { return (fUpperCutoff-fLowerCutoff); }
  };

  inline
  std::complex<double> RectangleFilter::TransmissivityAtFrequency(double parFrequency) const
  {
    if ((parFrequency >= fLowerCutoff) && (parFrequency <= fUpperCutoff))
    {
      return std::complex<double>(1.0, 0.0) * std::exp(-2 * utl::kPi * std::complex<double>(0.0,1.0) * parFrequency/hertz * fDelay/second);
    }
    else
    {
      return std::complex<double>(0.0, 0.0);
    }
  }


} // utl

#endif