void test_5_cound_z_to_x_to_c (void)
{
  TStopwatch watcher;
  watcher.Start();

  cout << "hi" << endl;

  TFile *f = new TFile ("caf_zcc.root", "READ");

  TTree *tree = (TTree*) f->Get("TMBTree");
  tree->SetBranchStatus("*",1);

  //tree->Draw("MCpart.getDMCvtx()!=0>>hev", "MCpart.pdgid()==4");
  int entries_found = tree->Draw(">>hev", "(MCpart.pdgid()==4)&&(MCpart.getPMCvtx()!=0)&&MCpart.getPMCvtx()!=0&&MCpart.getPMCvtx().getParent(0)!=0&&MCpart.getPMCvtx().getParent(0).getPMCvtx()!=0&&MCpart.getPMCvtx().getParent(0).getPMCvtx().getParent(0).pdgid()==23");
  
  cout << "Entries is " << entries_found << endl;
  cout << "Events is: " << hev->GetN() << endl;

  watcher.Stop();
  watcher.Print ("m");
  
}