#include <TROOT.h>
#include <TH1.h>
Double_t Median(TH1F* h1)
{
  const Int_t nq = 100;
  Double_t xq[nq];  // position where to compute the quantiles in [0,1]
  Double_t yq[nq];  // array to contain the quantiles
  for (Int_t i=0;i<nq;i++) xq[i] = Float_t(i+1)/nq;
  
  h1->GetQuantiles(nq,yq,xq);
  Double_t med=yq[nq/2-1];

  return med;

}