void dpwMip() {
  gROOT->SetStyle("Plain");

  TCanvas *c1=new TCanvas("c1","Canvas",800,600);
  gStyle->SetOptStat(0);
  gStyle->SetPalette(1);
  gStyle->SetGridWidth(10);
  //gStyle->SetPadLeftMargin(0.12);
  //gStyle->SetPadRightMargin(0.12);
  //c1->UseCurrentStyle();

  TFile::Open("summaryDpw.root");

  hAmm->SetTitle(";MIP x position (#mum);MIP y position (#mum);Charge (electrons)");

  for(int x(1);x<151;x++) {
    for(int y(1);y<151;y++) {
      double b=1100.0*hAmm->GetBinContent(x,y);
      hAmm->SetBinContent(x,y,b);
    }
  }

  // For lego2
  //hAmm->GetXaxis()->SetTitleOffset(1.4);
  //hAmm->GetYaxis()->SetTitleOffset(1.5);
  //hAmm->GetZaxis()->SetTitleOffset(1.3);

  // For colz
  hAmm->GetXaxis()->SetTitleOffset(1.1);
  hAmm->GetYaxis()->SetTitleOffset(1.1);

  hAmm->GetZaxis()->SetRangeUser(0,670.0);
  //Trim0Sigmas2D->GetXaxis()->SetRangeUser(0,83);
  //Trim0Sigmas2D->GetZaxis()->SetRangeUser(0,15);

  //hAmm->Draw("lego2");
  //hAmm->Draw("surf2");
  hAmm->Draw("colz");

  c1->Update();
  c1->Print("../dpwMip.png");
}