//
// Note: one MUST use this code as (pre)compiled (e.g. by ACLiC) because
//       the "storage" class inherits from a compiled class ("TObject"):
//       root [0] .L breakingExample.C++
//       root [1] breakingExample();
//

#include <iostream>

#include <TFile.h>
#include <TTree.h>
#include <TClonesArray.h>
#include <TVector3.h>
#include <TRandom.h>

class storage : public TObject
{
public:
  TTree* tree;
  
  storage() { tree = 0; }
  
  ClassDef(storage, 1); // my storage class
};

void breakingExample()
{

	TFile* file = TFile::Open("breaking.root", "RECREATE");
	storage stor;

	TClonesArray* datum = new TClonesArray("TVector3", 10);

	stor.tree = new TTree("tree", "tree");
	stor.tree->SetDirectory(0); // make it "memory resident, unregistered"
	stor.tree->Branch("brName", "TClonesArray", &datum, 256000, 99);

	std::cout<<"point_0"<<std::endl;

	const int nEvents = 10000;
	for(long i = 0; i < nEvents; ++i) {
		for(unsigned int j = 0; j < 10; ++j) {
			new ((*datum)[j]) TVector3(gRandom->Uniform(-10., 10.), gRandom->Uniform(-10., 10.), gRandom->Uniform(-10., 10.));
		}
		stor.tree->Fill();
	}
	stor.Write("my_storage_object");
	file->Close();

	std::cout<<"point_1"<<std::endl;

	file = TFile::Open("breaking.root", "READ");
	std::cout<<"point_2"<<std::endl;
	storage* storPtr = (storage*)file->Get("my_storage_object");
	std::cout<<"point_3"<<std::endl;
	datum = 0;
	storPtr->tree->SetBranchAddress("brName", &datum);
	std::cout<<"point_4"<<std::endl;
	for(long i = 0; i < storPtr->tree->GetEntries(); ++i) {
		std::cout<<"point_5 ... "<<i<<std::endl;
		storPtr->tree->GetEntry(i);
		std::cout<<"point_6"<<std::endl;
		for(unsigned j = 0; j < 10; ++j) {
			((TVector3*)(*datum)[j])->Print();
		}
	}
	file->Close();

}