void ZOffset3D( TH3* h, float z) {
//     This boosts the height for each bin of a histogram by a given constant value of z
   int nbinx = h->GetNbinsX();
   int nbiny = h->GetNbinsY();
   int nbinz = h->GetNbinsZ();
   for (int i=1; i<=nbinx; i++) {
      for (int j=1; j<=nbiny; j++) {
          for (int k=1; k<=nbinz; k++) {
              int h0=(h->GetBinContent(i,j,k));
            if (h0>0){
//             cout<<h->GetBinContent(i,j,k)<<"   ";
            h->SetBinContent(i,j,k,z+(h->GetBinContent(i,j,k)) );
//             cout<<h->GetBinContent(i,j,k)<<endl;
            }
          }
      }
   }
    h->SetMinimum(0);
}