void realDpwMip() {
  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.12);
  //c1->UseCurrentStyle();

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

  hAmm->SetTitle(";MIP x position (#mum);MIP y position (#mum);Laser signal (DTU)");
  for(int x(1);x<31;x++) {
    for(int y(1);y<31;y++) {
      double b=500.0*hAmm->GetBinContent(x,y);
      if(b<0.01) b=0.01;
      hAmm->SetBinContent(x,y,b);
    }
  }

  if(leg) {
    hAmm->GetXaxis()->SetTitleOffset(1.4);
    hAmm->GetYaxis()->SetTitleOffset(1.5);
    hAmm->GetZaxis()->SetTitleOffset(1.3);
    hAmm->GetZaxis()->SetRangeUser(0,287.0);
  } else {
    hAmm->GetXaxis()->SetTitleOffset(1.1);
    hAmm->GetYaxis()->SetTitleOffset(1.1);
    hAmm->GetZaxis()->SetRangeUser(0,299.0);
  }
  //Trim0Sigmas2D->GetXaxis()->SetRangeUser(0,83);
  //Trim0Sigmas2D->GetZaxis()->SetRangeUser(0,15);

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

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