{
  gStyle->SetOptStat("emr");

  //  TChain *tr = new TChain("tree");
  TChain *ch = new TChain("tree");
  //  ch->Add("/home/deepak/Desktop/Xcheck_bin3/macro_aod/sel_BsToPhiMuMu_*.root");
    ch->Add("/home/deepak/Desktop/Xcheck_bin3/macro_aod/sel_BsToPhiMuMu_phimumu_2016_mc_cutopt_s16000000.root");
  TTree *tr = ch;
  // tr->Print();

  tr->SetBranchStatus("*", 0);
  tr->SetBranchStatus("Bmass", 1);
  tr->SetBranchStatus("Mumumass", 1);
  tr->SetBranchStatus("Phimass", 1);
  tr->SetBranchStatus("Q2", 1);

  //const char *varmass[] = {"Bmass", "Mumumass", "Phimass"};
  const char *varmass[] = {"Bmass", "Mumumass"};
  Int_t n = sizeof(varmass) / sizeof(char*);

  TCut Q[] = {"Q2 > 1.0 && Q2 < 4.30",
	      "Q2 > 4.30 && Q2 < 8.68",
	      "Q2 > 8.68 && Q2 < 10.09",
	      "Q2 > 10.09 && Q2 < 12.86",
	      "Q2 > 12.86 && Q2 < 14.18",
	      "Q2 > 14.18 && Q2 < 16.0",
	      "Q2 > 16.0 && Q2 < 19.0"}
  Int_t m = sizeof(Q) / sizeof(TCut);

  //TCanvas *c = new TCanvas("c", "c");
  //  c->Divide(m, n);
  // c->DivideSquare(n * m);

  Int_t p = 1;

  TH1D *h_bs   = new TH1D("h_bs","B_{s} mass ; B_{s} Mass; Events/0.00625", 80, 5.1,5.6);
  
  for (Int_t i = 0; i < n; i++) {
    for (Int_t j = 0; j < m; j++) {
      //      c->cd(p++);
      /*   TCanvas *c[j] = new TCanvas;
      c[j]->cd(p++);
      c[j] = new TCanvas(Form("c%d",j));*/
      TCanvas* c=new TCanvas(Form("c_%d_%d",i,j),"",600,600);
      //      tr->Draw(varmass[i], Q[j]);
      tr->Project("h_bs",varmass[i], Q[j]);
      c->SaveAs(Form("multvar_%s_%d.pdf",varmass[i],j));
    }
  }
  //  c->cd(0);
  //c->SaveAs("multvar.pdf");
}