void gdsNodpwMip() {
  bool leg(false);

  gROOT->SetStyle("Plain");

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

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

  hAmm->SetTitle(";MIP x position (#mum);MIP y position (#mum);Charge (electrons)");
  for(int x(1);x<31;x++) {
    for(int y(1);y<31;y++) {
      double b=1100.0*hAmm->GetBinContent(x,y);
      hAmm->SetBinContent(x,y,b);
    }
  }

  if(leg) {
    hAmm->GetXaxis()->SetTitleOffset(1.4);
    hAmm->GetYaxis()->SetTitleOffset(1.5);
    hAmm->GetZaxis()->SetTitleOffset(1.3);
  } else {
    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);

  if(leg) hAmm->Draw("lego2");
  else    hAmm->Draw("colz");

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