{
  int nptbins  = 11; 
  float ptbin[nptbins+1]={0,1,2,3,4,5,6,7,8,9,10,20};
  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("hmapselalleff","hmapselalleff",netabins,etabin,nptbins,ptbin);        
  Den_EFmu4_map->Sumw2();
  Num_EFmu4_map->Sumw2();
  hmapden->Divide(Num_EFmu4_map,Den_EFmu4_map,1.,1.,"b") ;
  hmapden->SaveAs("Efmu4_trkjpsi_2011_periodG.root");
}