{
gStyle->SetStatBorderSize(1);
gStyle->SetStatColor(0);
gStyle->SetTitleFillColor(0);
gStyle->SetCanvasColor(0);
gStyle->SetDrawBorder(0);
gStyle->SetPadColor(0);
gStyle->SetPadBorderMode(0);
gStyle->SetFrameBorderMode(0);
gStyle->SetPalette(1,0);
gStyle->SetOptStat(0);
gStyle->SetOptTitle(0);

TFile* ff = new TFile("tree_cu63.root","r");
TTree* T_b = (TTree*)gROOT->FindObject("T_b");
TTree* T_d = (TTree*)gROOT->FindObject("T_d");
TTree* T_nb = (TTree*)gROOT->FindObject("T_nb");
TTree* T_l = (TTree*)gROOT->FindObject("T_l");
TTree* T_af = (TTree*)gROOT->FindObject("T_af");

TH2F* hb = new TH2F("hb","Momentum Separation",600,0,1.2,600,0,1.2);
TH2F* hb1 = new TH2F("hb1","Momentum Separation",600,0,1.2,600,0,1.2);
TH2F* hd = new TH2F("hd","Momentum Separation",600,0,1.2,600,0,1.2);
TH2F* hnb = new TH2F("hnb","Momentum Separation",600,0,1.2,600,0,1.2);
TH2F* hl = new TH2F("hl","Momentum Separation",600,0,1.2,600,0,1.2);
TH2F* haf = new TH2F("haf","Momentum Separation",600,0,1.2,600,0,1.2);

TH1F* hhb = new TH1F("hhb","Momentum Sum",500,0,2);
TH1F* hhb1 = new TH1F("hhb1","Momentum Sum",500,0,2);
TH1F* hhd = new TH1F("hhd","Momentum Sum",500,0,2);
TH1F* hhnb = new TH1F("hhnb","Momentum Sum",500,0,2);
TH1F* hhl = new TH1F("hhl","Momentum Sum",500,0,2);
TH1F* hhaf = new TH1F("hhaf","Momentum Sum",500,0,2);

c1 = new TCanvas("c1","",600,600);
T_b->Draw("p2.af.P:k1.af.P+k2.af.P>>hb","w.pf*w.pkk*w.nf");
T_b->Draw("p1.af.P:k1.af.P+k2.af.P>>hb1","w.pf*w.pkk*w.nf");
T_d->Draw("p1.af.P:k1.af.P+k2.af.P>>hd","w.pkk");
T_nb->Draw("p1.af.P:k1.af.P+k2.af.P>>hnb","w.pf*w.kk");
T_l->Draw("p1.af.P:k1.af.P+k2.af.P>>hl","w.lk*w.pk");
T_af->Draw("p1.af.P:k1.af.P+k2.af.P>>haf","w.paf*w.kk");

T_b->Draw("p2.af.P+k1.af.P+k2.af.P>>hhb","w.pf*w.pkk*w.nf/50.");
T_b->Draw("p1.af.P+k1.af.P+k2.af.P>>hhb1","w.pf*w.pkk*w.nf/50.");
T_d->Draw("p1.af.P+k1.af.P+k2.af.P>>hhd","w.pkk/100.");
T_nb->Draw("p1.af.P+k1.af.P+k2.af.P>>hhnb","w.pf*w.kk/100.");
T_l->Draw("p1.af.P+k1.af.P+k2.af.P>>hhl","w.lk*w.pk/100.");
T_af->Draw("p1.af.P+k1.af.P+k2.af.P>>hhaf","w.paf*w.kk/100.");

hnb->SetXTitle("|P_{K^{+}}|+|P_{K^{-}}| (GeV/c)");
hnb->SetYTitle("|P_{proton}| (GeV/c)");
hnb->SetMarkerColor(4);
hd->SetMarkerColor(8);
hl->SetMarkerColor(1);
haf->SetMarkerColor(7);
hb->SetMarkerColor(2);
hb1->SetMarkerColor(6);

hnb->SetFillColor(4);
hd->SetFillColor(8);
hl->SetFillColor(1);
haf->SetFillColor(7);
hb->SetFillColor(2);
hb1->SetFillColor(6);

hhnb->SetXTitle("|P_{proton}|+|P_{K^{+}}|+|P_{K^{-}}| (GeV/c)");
hhnb->SetYTitle("Arbitary Scale");
hhnb->SetLineColor(4);
hhd->SetLineColor(8);
hhl->SetLineColor(1);
hhaf->SetLineColor(7);
hhb->SetLineColor(2);
hhb1->SetLineColor(6);

hhnb->SetLineWidth(1);
hhd->SetLineWidth(1);
hhl->SetLineWidth(1);
hhaf->SetLineWidth(1);
hhb->SetLineWidth(1);
hhb1->SetLineWidth(1);

hnb->Draw();
hl->Draw("same");
haf->Draw("same");
hd->Draw("same");
hb1->Draw("same");
hb->Draw("same");

TLegend* leg = new TLegend(0.62,0.6,0.92,0.92,"|P_{proton}|:|P_{K^{+}}|+|P_{K^{-}}|");
leg->SetFillColor(0);
leg->AddEntry(hb,"Bound state","f");
leg->AddEntry(hb1,"B.S. w/ recoil proton","f");
leg->AddEntry(hnb,"No bound state","f");
leg->AddEntry(hd,"Direct production","f");
leg->AddEntry(hl,"#Lambda_{1520} production","f");
leg->AddEntry(haf,"a0/f0 production","f");
leg->Draw();

c2 = new TCanvas("c2","",600,600);
c2->cd();
hhnb->Draw();
hhl->Draw("same");
hhaf->Draw("same");
hhd->Draw("same");
hhb1->Draw("same");
hhb->Draw("same");

TLegend* leg2 = new TLegend(0.15,0.56,0.45,0.88,"|P_{proton}|+|P_{K1}|+|P_{K2}|");
leg2->SetFillColor(0);
leg2->AddEntry(hhb,"Bound state","l");
leg2->AddEntry(hhb1,"B.S. w/ recoil proton","l");
leg2->AddEntry(hhnb,"No bound state","l");
leg2->AddEntry(hhd,"Direct production","l");
leg2->AddEntry(hhl,"#Lambda_{1520} production","l");
leg2->AddEntry(hhaf,"a0/f0 production","l");
leg2->Draw();

}