#include "JTools/JCombinatorics.hh"
#include "TH1.h"
#include "TFitResult.h"

using namespace JPP;

namespace JNANOBEACON{
  
  typedef JCombinatorics::pair_type  pair_type;
  inline bool comparepair(const pair_type& A, const pair_type& B) {
    return( abs(A.first - A.second) > abs(B.first - B.second) ) ;
  }

  inline TFitResultPtr Fit(TH1D* h){
    TFitResultPtr fitresult;
    int MaxBin = h->GetMaximumBin();
    double maximum = h->GetBinCenter(MaxBin);
    new TF1("fit","gaus",maximum-15,maximum+15);
    fitresult = h->Fit("fit","QRS");
    return fitresult;
  }

}