{
gSystem->Load("./libMyClass.so");
gSystem->Load("libPhysics.so");
TFile f("framework.root", "RECREATE");


TTree *pTree = new TTree("MyTree","My analysis tree");
TMyEvent *pEvent = new TMyEvent(0);
pTree->Branch("TMyEvent", "TMyEvent", &pEvent, 256000/8, 1);
TMyCand *pCand;

int n(0), nCand(0);
// -- Run over 100 'events' 
for (Int_t i=0; i<100; i++) {
  pEvent->Clear();
  n = 20*gRandom->Rndm();       // each event has a random number of candidates
  for (Int_t j=0; j<n; j++) {
    pCand = pEvent->AddCand();
    pCand->setEraw(3.1415); 
    pCand->setCharge(1);
    pCand->setPlab(10.*gRandom->Gaus(2., 0.10), 10.*gRandom->Gaus(2., 0.12), 10.*gRandom->Gaus(2., 0.20));
  }
  pEvent->SetNCand(n);
  pTree->Fill();
}

pTree->Write();
f.Close();
}