HannWindow.h

Go to the documentation of this file.

#ifndef _utl_HannWindow_h_
#define _utl_HannWindow_h_ 

#include <utl/AnalyticWindow.h>
#include <utl/AugerException.h>
#include <utl/PhysicalConstants.h>
#include <cmath>

namespace utl {

  class HannWindow : public AnalyticWindow {
    public: 
      HannWindow(double relativeWidthOnEachSide, long traceLength = 0);
      virtual ~HannWindow();
      virtual double GetWeightAtBin(long parBin) const;
      virtual void SetTraceLength(long parTraceLength);
    protected:
      double fRise;
      double fFall;
      long fWindowSize;
  };

  inline
  double HannWindow::GetWeightAtBin(long parBin) const
  {
    if ((parBin < 0) || (parBin >= fTraceLength))
      throw(OutOfBoundException("Bin is outside of trace."));
    if (parBin < fRise)
      return 0.5 * (1 - std::cos(2 * kPi * parBin / (fWindowSize-1)));
    else if (parBin > fFall)
      return 0.5 * (1 - std::cos(2 * kPi * (long(parBin - fFall + fRise)) / (fWindowSize-1)));
    else
      return 1.0;
  }

  inline
  void HannWindow::SetTraceLength(long parTraceLength)
  {
    fTraceLength = parTraceLength;
    fRise = (fTraceLength * fRelativeWidthOnEachSide);
    fFall = (fTraceLength * (1. - fRelativeWidthOnEachSide));
    fWindowSize = 2 * long(fRise);
    fRenormalizationFactor = 0.0; // delete cached renormalization factor (will usually only change slightly with trace length, but anyway ...)
  }

} // utl

#endif