{

  gStyle->SetOptStat("e");
  gStyle->SetOptTitle(0);
  TCanvas *c3 = new TCanvas("c3","c3",500,500);

  TChain *ch = new TChain("tree");
  ch->Add("/home/deepak/BsToPhiMuMu/MC/mumuphisel/cut1mcnogen/sel_BsToPhiMuMu_mumuphi_cut1_mc_nogen_mc.nogen_cut1_*.root");
  TTree *tr1 = ch;


  tr1->SetBranchStatus("*", 0);
  tr1->SetBranchStatus("Bvtxcl", 1);
  tr1->SetBranchStatus("Triggers", 1);
  tr1->SetBranchStatus("Bmass", 1);
  tr1->SetBranchStatus("Phimass", 1);
  tr1->SetBranchStatus("Mumumass",1);
  tr1->SetBranchStatus("Mumumasserr",1);

  TH1F *h1 = new TH1F("h1", "K^{+} Transverse momentum",80,0.0,1.0);
  tr1->Draw("Bvtxcl>>h1","Triggers==1 && Bmass > 5.1 && Bmass < 5.6 && Phimass > 1.01 && Phimass < 1.03 && (Mumumass < 3.096916-5.5*Mumumasserr || Mumumass > 3.096916+3.5*Mumumasserr) && (Mumumass < 3.686109-3.5*Mumumasserr || Mumumass > 3.686109+3.5*Mumumasserr) && (Bmass > 5.287 && Bmass < 5.447)");
  h1->SetDirectory(0);
  h1->Draw();
  gPad->Update();
  TPaveStats *tps1 = (TPaveStats*) h1->FindObject("stats");
  tps1->SetName("Hist1 Stats");
  double X1 = tps1->GetX1NDC();
  double Y1 = tps1->GetY1NDC();
  double X2 = tps1->GetX2NDC();
  double Y2 = tps1->GetY2NDC();

  TChain *ch = new TChain("tree");
  ch->Add("/home/deepak/BsToPhiMuMu/DATA_DoubleMuon/cut1/sel_BsToPhiMuMu_Double_2016*.root");
  ch->Add("/home/deepak/BsToPhiMuMu/DATA_Charmonium/Rereco/sel_BsToPhiMuMu_*.root");

  TTree *tr2 = ch;

  TH1F *h2 = new TH1F("h2", "K^{+} Transverse momentum",80,0.0,1.0);
  tr2->Draw("Bvtxcl>>h2","Triggers==1 && Bmass > 5.1 && Bmass < 5.6 && Phimass > 1.01 && Phimass < 1.03 && (Mumumass < 3.096916-5.5*Mumumasserr || Mumumass > 3.096916+3.5*Mumumasserr) && (Mumumass < 3.686109-3.5*Mumumasserr || Mumumass > 3.686109+3.5*Mumumasserr) && (Bmass > 5.191 && Bmass < 5.271 || Bmass > 5.463 && Bmass <5.543)");
  h2->SetDirectory(0);
  h2->Scale(1367.144);
  h2->Draw("E");
  gPad->Update();
  TPaveStats *tps2 = (TPaveStats*) h2->FindObject("stats");
  tps2->SetTextColor(kRed);
  tps2->SetLineColor(kRed);
  tps2->SetX1NDC(X1);
  tps2->SetX2NDC(X2);
  tps2->SetY1NDC(Y1-(Y2-Y1));
  tps2->SetY2NDC(Y1);

  h2->Scale(1.0 / h2->Integral());
  //  h2->Scale(1367.144);
  h2->Draw();
  h1->Scale(1.0 / h1->Integral());
  h2->SetLineColor(kRed);
  h1->Draw("same");
  tps1->Draw("same");
  tps2->Draw("same");
  h2->GetYaxis()->SetTitle("Events");
  h2->GetXaxis()->SetTitle("B_{s} Vertex CL");
  h2->GetYaxis()->SetTitleOffset(1.5);
  h2->GetXaxis()->SetLabelFont(63);
  h2->GetXaxis()->SetLabelSize(14);
  h2->GetYaxis()->SetLabelFont(63);
  h2->GetYaxis()->SetLabelSize(14);
  h2->GetXaxis()->SetTitleFont();
  h2->GetYaxis()->SetTitleFont();

  leg = new TLegend(0.1,0.8,0.2,0.9);
  leg->AddEntry("h1","MC");
  leg->AddEntry("h2","Data");
  leg->Draw();

  c3->SaveAs("bvtxcl.png");
}