void ReverseXAxis(TH1 *h)
{
   // Redraw the new axis
   gPad->Update();
   TGaxis *newaxis = new TGaxis(gPad->GetUxmax(),
                                gPad->GetUymin(),
                                gPad->GetUxmin(),
                                gPad->GetUymin(),
                                h->GetXaxis()->GetXmin(),
                                h->GetXaxis()->GetXmax(),
                                510,"-");
   newaxis->SetLabelOffset(-0.03);
   newaxis->Draw();
}
void ReverseYAxis(TH1 *h)
{
   // Redraw the new axis
   gPad->Update();
   TGaxis *newaxis = new TGaxis(gPad->GetUxmin(),
                                gPad->GetUymax(),
                                gPad->GetUxmin()-0.001,
                                gPad->GetUymin(),
                                h->GetYaxis()->GetXmin(),
                                h->GetYaxis()->GetXmax(),
                                510,"+");
   newaxis->SetLabelOffset(-0.05);
   newaxis->Draw();
}

void yepes2()
{
   TCanvas *C = new TCanvas("C","C",0.,0,800,800);

   TPad *Overlay = new TPad("Overlay","",0,0,1,1);
   Overlay->SetFillColor(0);
   Overlay->SetFillStyle(4000);
   Overlay->SetBorderMode(0);
   Overlay->SetFrameBorderMode(0);
   Overlay->SetFrameFillStyle(0);
   Overlay->SetFrameBorderMode(0);

   TFile *f = new TFile("example-figures.root");

   TH2F *h1 = (TH2F *)f->Get("patientRho_yx");
   TH2F *h2 = (TH2F *)f->Get("DoseW-fdc_yx");

   h1->SetStats(0);
   h2->SetStats(0);
   h1->GetXaxis()->SetLabelSize(0);
   h2->GetXaxis()->SetLabelSize(0);
   h1->GetYaxis()->SetLabelSize(0);
   h2->GetYaxis()->SetLabelSize(0);

   TExec *ex1 = new TExec("ex1","gStyle->SetPalette(kGreyScale,0,1.0);");
   TExec *ex2 = new TExec("ex2","gStyle->SetPalette(kRainBow,0,0.5);");

   h2->SetMaximum(65);

   h1->Draw("col");
   ReverseXAxis(h1);
   ReverseYAxis(h1);

   ex1->Draw();
   h1->Draw("col same");

   ex2->Draw();
   h2->SetTitle("");
   Overlay->Draw();
   Overlay->cd();
   h2->Draw("colz same");
}