#include "TFile.h"
#include "TTree.h"
#include "TH2F.h"

Int_t PlotStable() {
  TFile *f1 = new TFile("StableNuclei.root","READ");
  TTree *tree1 = (TTree*)f1->Get("StableNuclei");
  TH2F  *hist = new TH2F("hist", "Stable Nuclei;Neutron number N; Atomic number Z;",141,-0.5,140.5,91,-0.5,90.5);
  tree1->Draw("Z:N>>hist");
  
  gStyle->SetPalette(0);

  hist->GetXaxis()->CenterTitle();
  hist->GetYaxis()->CenterTitle();
  hist->Draw("COL");
  
  //  delete hist;
  //  delete f;
  return 0;
}