{
  int nptbins  = 7; 
  float ptbin[8]={0,3,4,5,7,9,15,25};
  int netabins = 18; 
  float etabin[19]={-2.4,-2.0,-1.7,-1.5,-1.3,-1.05,-0.8,-0.5,-0.1,0,0.1,0.5,0.8,1.05,1.3,1.5,1.7,2.0,2.4};
  TH2F *hmapden = new TH2F("hmapden","hmapden",netabins,etabin,nptbins,ptbin);        
  mapden_all->Sumw2();
  mapnum_sel->Sumw2();
  hmapden->Divide(mapnum_sel,mapden_all,1.,1.,"b") ;
  hmapden->SaveAs("TMapDenData2010.root");
}