PeriodicBinFinderInPhi.h

Go to the documentation of this file.

#ifndef PeriodicBinFinderInPhi_H
#define PeriodicBinFinderInPhi_H

#include "Utilities/BinningTools/interface/BaseBinFinder.h"

#include <cmath>

template <class T>
class PeriodicBinFinderInPhi : public BaseBinFinder<T> {
public:

  PeriodicBinFinderInPhi() : theNbins(0), thePhiStep(0), thePhiOffset(0) {}

  PeriodicBinFinderInPhi( T firstPhi, int nbins) :
    theNbins( nbins), thePhiStep( twoPi() / nbins),
    thePhiOffset( firstPhi - thePhiStep/2.) {}

  virtual  int binIndex( T phi) const {
    T tmp = fmod((phi - thePhiOffset), twoPi()) / thePhiStep;
    if ( tmp < 0) tmp += theNbins;
    return std::min( int(tmp), theNbins-1);
  }

  virtual int binIndex( int i) const {
    int ind = i % theNbins;
    return ind < 0 ? ind+theNbins : ind;
  }
   
  virtual T binPosition( int ind) const {
    return thePhiOffset + thePhiStep * ( ind + 0.5);
  }

  static T pi() { return 3.141592653589793238;}
  static T twoPi() { return 2.*pi();}

private:

  int theNbins;
  T thePhiStep;
  T thePhiOffset;

};
#endif