{
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");

TH1F* hb = new TH1F("hb","Invariant Mass of pKK",300,1.9,2.2);
TH1F* hb1 = new TH1F("hb1","Invariant Mass of pKK",300,1.9,2.2);
TH1F* hd = new TH1F("hd","Invariant Mass of pKK",300,1.9,2.2);
TH1F* hnb = new TH1F("hnb","Invariant Mass of pKK",300,1.9,2.2);
TH1F* hl = new TH1F("hl","Invariant Mass of pKK",300,1.9,2.2);
TH1F* haf = new TH1F("haf","Invariant Mass of pKK",300,1.9,2.2);

c1 = new TCanvas("c1","",800,600);
T_b->Draw("M2>>hb","w.pf*w.pkk*w.nf/100.");
T_b->Draw("M1>>hb1","w.pf*w.pkk*w.nf/100.");
T_d->Draw("M1>>hd","w.pkk/100.");
T_nb->Draw("M1>>hnb","w.pf*w.kk/100.");
T_l->Draw("M1>>hl","w.lk*w.pk/100.");
T_af->Draw("M1>>haf","w.paf*w.kk/100.");

hnb->SetXTitle("Invariant Mass (GeV)");
hnb->SetYTitle("Arbitary Scale");
hnb->SetMaximum(40);
hnb->SetLineColor(4);
hd->SetLineColor(8);
hl->SetLineColor(1);
haf->SetLineColor(7);
hb->SetLineColor(2);
hb1->SetLineColor(6);

hnb->SetLineWidth(1);
hd->SetLineWidth(1);
hl->SetLineWidth(1);
haf->SetLineWidth(1);
hb->SetLineWidth(1);
hb1->SetLineWidth(1);

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

TLegend* leg = new TLegend(0.55,0.5,0.85,0.85,"Invariant Mass of PK^{+}K^{-}");
leg->SetFillColor(0);
leg->AddEntry(hb,"Bound State","l");
leg->AddEntry(hb1,"B.S. w/ recoil proton","l");
leg->AddEntry(hnb,"No bound state","l");
leg->AddEntry(hd,"Direct  production","l");
leg->AddEntry(hl,"#Lambda_{1520} production","l");
leg->AddEntry(haf,"a0/f0 production","l");
leg->Draw();
}