PeriodicBinFinderInZ.h

Go to the documentation of this file.

#ifndef DetLayers_PeriodicBinFinderInZ_H
#define DetLayers_PeriodicBinFinderInZ_H

#include "Utilities/BinningTools/interface/BaseBinFinder.h"
#include "Geometry/CommonDetUnit/interface/GeomDet.h"
#include <cmath>

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

  PeriodicBinFinderInZ() : theNbins(0), theZStep(0), theZOffset(0) {}

  PeriodicBinFinderInZ(std::vector<const GeomDet*>::const_iterator first,
               std::vector<const GeomDet*>::const_iterator last) :
    theNbins( last-first) 
  {
    float zFirst = (**first).surface().position().z();
    theZStep = ((**(last-1)).surface().position().z() - zFirst) / (theNbins-1);
    theZOffset = zFirst - 0.5*theZStep;
  }

  virtual int binIndex( T z) const {
    return binIndex( int((z-theZOffset)/theZStep));
  }

  virtual int binIndex( int i) const {
    return std::min( std::max( i, 0), theNbins-1);
  }
   
  virtual T binPosition( int ind) const {
    return theZOffset + theZStep * ( ind + 0.5);
  }

private:

  int theNbins;
  T theZStep;
  T theZOffset;

};
#endif