#include <iostream>

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

class storage : public TObject
{

  public:

	TTree* tree;

	ClassDef(storage, 1);

};



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->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("storage");
	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("storage");
	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"<<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();

}