/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_5/src/RecoMuon/DetLayers/src/GeneralBinFinderInR.h

Go to the documentation of this file.

#ifndef GeneralBinFinderInR_H
#define GeneralBinFinderInR_H

#include <Utilities/BinningTools/interface/BaseBinFinder.h>
#include "RBorderFinder.h"

#include <cmath>
#include <vector>

template <class T>
class GeneralBinFinderInR : public BaseBinFinder<T>{
public:
  
  typedef RBorderFinder::Det Det; //FIXME!!!

  GeneralBinFinderInR() : theNbins(0) {}

  GeneralBinFinderInR(const RBorderFinder& bf) {
    theBorders=bf.RBorders();
    theBins=bf.RBins();
    theNbins=theBins.size();
  }

  GeneralBinFinderInR(std::vector<Det*>::const_iterator first,
                      std::vector<Det*>::const_iterator last)
    : theNbins( last-first)
  {
    std::vector<const Det*> dets(first,last);
    RBorderFinder bf(dets);
    theBorders=bf.RBorders();
    theBins=bf.RBins();
    theNbins=theBins.size();
  }

  
  virtual int binIndex( T R) const {
    int i;
    for (i = 0; i<theNbins; i++) {
      if (R < theBorders[i]){
         break;
      }
    }
    return binIndex(i-1);
  }

  virtual int binIndex( int i) const {
    return std::min( std::max( i, 0), theNbins-1);
  }
   
  virtual T binPosition( int ind) const {
    return theBins[binIndex(ind)];
  }


private:
  int theNbins;
  std::vector<T> theBorders;
  std::vector<T> theBins;

};
#endif