{

  gStyle->SetOptStat("emr");

  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/cut1/sel_BsToPhiMuMu_*.root");
  TTree *tr = ch;

  TH2D *h_kl   = new TH2D("h_kl"  ,"cos#theta_{K} vs cos#theta_{L}; cos#theta_{K} ; cos#theta_{L}",75, -1.0, 1.0, 75, -1.0, 1.0);

  //  h_kl->SetMarkerColor(kBlack);
  h_kl->SetMarkerStyle(20);
  h_kl->GetYaxis()->SetTitleOffset(1.6);

  TCut trig   = "Triggers==1";
  TCut bm = "Bmass > 5.1 && Bmass < 5.6";
  TCut phm = "Phimass > 1.01 && Phimass < 1.03";

  TCut jpR = "(Mumumass < 3.096916-5.5*Mumumasserr || Mumumass > 3.096916+3.5*Mumumasserr)";
  TCut psR = "(Mumumass < 3.686109-3.5*Mumumasserr || Mumumass > 3.686109+3.5*Mumumasserr)";

  
  TCut AR = "abs(Bmass - Mumumass - 2.270) > 0.170 && abs(Bmass - Mumumass - 1.681) > 0.080";

  TCut selbm = trig && bm && phm && jpR && psR  && AR;
  
  TCanvas *cv1 = new TCanvas("cv1","",500,500);
  tr->Draw("CosThetaL : CosThetaK >> h_kl",selbm);

  
  cv1->SaveAs("coskl_scatr_jpsr_ar.png");

  //  h_kl::GetCovariance(Double_t axis1 = 1, Double_t axis2 = 2 );const
  
  h_kl.GetCovariance(Double_t axis=1 , Double_t axis=2);
  
  
}