TH2D *hEGammaVsKKpi = new TH2D("hEGamma","",100,1.2,1.5,60,5.4,11.4)
 mass_pi0kpkm__B4_M7_Tree->Draw("beam_energy:m_kpkmpi0>>hEGamma","(m_kppi0 < 0.81 || m_kppi0 > 0.97)&&(m_kmpi0 < 0.81 || m_kmpi0 > 0.97)")
TFile *hFile = new TFile("hFile.root","recreate")
hEGamma->Write()
hFile->Close()