{
  gStyle->SetPadGridX(1);
  gStyle->SetPadGridY(1);

  TH2F* b1_p1 = new TH2F("b1_p1","recoil proton; momentum (GeV/c); angle (degree)",100,0,1.5,100,0,180);
  TH2F* b1_p2 = new TH2F("b1_p2","proton; momentum (GeV/c); angle (degree)",100,0,1.5,100,0,180);
  TH2F* b1_pi = new TH2F("b1_pi","pion; momentum (GeV/c); angle (degree)",100,0,1.5,100,0,180);
  TH2F* b1_k = new TH2F("b1_k","kaon; momentum (GeV/c); angle (degree)",100,0,1.5,100,0,180);
  TH2F* b1_k1 = new TH2F("b1_k1","kaon+; momentum (GeV/c); angle (degree)",100,0,1.5,100,0,180);
  TH2F* b1_k2 = new TH2F("b1_k2","kaon-; momentum (GeV/c); angle (degree)",100,0,1.5,100,0,180);

  //  TTree* T = (TTree*)gROOT->FindObject("T_bl"); 
  //  TTree* T = (TTree*)gROOT->FindObject("T_bs");
  //  TTree* T = (TTree*)gROOT->FindObject("T_dl");
  //  TTree* T = (TTree*)gROOT->FindObject("T_ds");
  //  TTree* T = (TTree*)gROOT->FindObject("T_b");
  //  TTree* T = (TTree*)gROOT->FindObject("T_d");
  //  TTree* T = (TTree*)gROOT->FindObject("T_nb");
  //  TTree* T = (TTree*)gROOT->FindObject("T_l");
  TTree* T = (TTree*)gROOT->FindObject("T_af");

  //  char* weight = "w.pf*w.lk*w.ppi*w.nf";
  //  char* weight = "w.pf*w.sk*w.lg*w.ppi*w.nf";
  //  char* weight = "w.lk*w.ppi";
  //  char* weight = "w.sk*w.lg*w.ppi";
  //  char* weight = "w.pf*w.pkk*w.nf";
  //  char* weight = "w.pkk";
  //  char* weight = "w.pf*w.kk";
  //  char* weight = "w.lk*w.pk";
  char* weight = "w.paf*w.kk";


  c1 = new TCanvas("c1","",800,600);
  c1->Divide(2,2);
  c1->cd(1);
  T->Draw("p1.af.Th*180/3.1416:p1.af.P>>b1_p1",weight,"cont");
  c1->cd(2);
  //  T->Draw("p2.af.Th*180/3.1416:p2.af.P>>b1_p2",weight,"cont");
  c1->cd(3);
  //  T->Draw("pi.af.Th*180/3.1416:pi.af.P>>b1_pi",weight,"cont");
  T->Draw("k1.af.Th*180/3.1416:k1.af.P>>b1_k1",weight,"cont");
  c1->cd(4);
  //  T->Draw("k.af.Th*180/3.1416:k.af.P>>b1_k",weight,"cont");
  T->Draw("k2.af.Th*180/3.1416:k2.af.P>>b1_k2",weight,"cont");
}