go-hep.org/x/hep@v0.38.1/groot/rtree/branch.go

// Copyright ©2017 The go-hep Authors. All rights reserved.
// Use of this source code is governed by a BSD-style
// license that can be found in the LICENSE file.

package rtree

import (
	"fmt"
	"reflect"
	"strings"

	"go-hep.org/x/hep/groot/rbase"
	"go-hep.org/x/hep/groot/rbytes"
	"go-hep.org/x/hep/groot/rcont"
	"go-hep.org/x/hep/groot/rdict"
	"go-hep.org/x/hep/groot/riofs"
	"go-hep.org/x/hep/groot/root"
	"go-hep.org/x/hep/groot/rtypes"
	"go-hep.org/x/hep/groot/rvers"
)

const (
	defaultBasketSize = 32 * 1024 // default basket size in bytes
	defaultSplitLevel = 99        // default split-level for branches
	defaultMaxBaskets = 10        // default number of baskets
)

type tbranch struct {
	named          rbase.Named
	attfill        rbase.AttFill
	compress       int         // compression level and algorithm
	basketSize     int         // initial size of Basket buffer
	entryOffsetLen int         // initial length of entryOffset table in the basket buffers
	writeBasket    int         // last basket number written
	entryNumber    int64       // current entry number (last one filled in this branch)
	iobits         tioFeatures // IO features for newly-created baskets
	offset         int         // offset of this branch
	maxBaskets     int         // maximum number of baskets so far
	splitLevel     int         // branch split level
	entries        int64       // number of entries
	firstEntry     int64       // number of the first entry in this branch
	totBytes       int64       // total number of bytes in all leaves before compression
	zipBytes       int64       // total number of bytes in all leaves after compression
	branches       []Branch    // list of branches of this branch
	leaves         []Leaf      // list of leaves of this branch
	baskets        []Basket    // list of baskets of this branch

	basketBytes []int32 // length of baskets on file
	basketEntry []int64 // table of first entry in each basket
	basketSeek  []int64 // addresses of baskets on file

	fname string // named of file where buffers are stored (empty if in same file as Tree header)

	ctx basketCtx // basket context for the current basket

	tree *ttree          // tree header
	btop Branch          // top-level parent branch in the tree
	bup  Branch          // parent branch
	dir  riofs.Directory // directory where this branch's buffers are stored
}

func newBranchFromWVar(w *wtree, name string, wvar WriteVar, parent Branch, lvl int, cfg wopt) (Branch, error) {
	base := &tbranch{
		named:    *rbase.NewNamed(name, ""),
		attfill:  *rbase.NewAttFill(),
		compress: int(cfg.compress),

		iobits:      w.ttree.iobits,
		basketSize:  int(cfg.bufsize),
		maxBaskets:  defaultMaxBaskets,
		basketBytes: make([]int32, 0, defaultMaxBaskets),
		basketEntry: make([]int64, 1, defaultMaxBaskets),
		basketSeek:  make([]int64, 0, defaultMaxBaskets),

		tree: &w.ttree,
		btop: btopOf(parent),
		bup:  parent,
		dir:  w.dir,
	}

	var (
		b Branch = base

		title = new(strings.Builder)
		rt    = reflect.TypeOf(wvar.Value).Elem()
	)

	title.WriteString(wvar.Name)
	switch k := rt.Kind(); k {
	case reflect.Array:
		et, shape := flattenArrayType(rt)
		for _, dim := range shape {
			fmt.Fprintf(title, "[%d]", dim)
		}
		rt = et

	case reflect.Slice:
		switch wvar.Count {
		case "":
			// write as a std::vector<T>
			return newBranchElementFromWVar(w, base, wvar, parent, lvl, cfg)
		default:
			fmt.Fprintf(title, "[%s]", wvar.Count)
			rt = rt.Elem()
		}
		base.entryOffsetLen = 1000 // slice, so we need an offset array

	case reflect.String:
		base.entryOffsetLen = 1000 // string, so we need an offset array

	case reflect.Struct:
		return newBranchElementFromWVar(w, base, wvar, parent, lvl, cfg)
	}

	isBranchElem := false
	if parent != nil {
		if parent, ok := parent.(*tbranchElement); ok {
			isBranchElem = true
			be := &tbranchElement{
				tbranch: *base,
				class:   parent.class,
				parent:  parent.class,
			}
			b = be
			base = &be.tbranch
		}
	}

	if !isBranchElem && rt.Kind() != reflect.Struct {
		code := gotypeToROOTTypeCode(rt)
		fmt.Fprintf(title, "/%s", code)
	}

	_, err := newLeafFromWVar(w, b, wvar, lvl, cfg)
	if err != nil {
		return nil, err
	}

	base.named.SetTitle(title.String())
	base.createNewBasket()

	return b, nil
}

func (b *tbranch) RVersion() int16 {
	return rvers.Branch
}

func (b *tbranch) Name() string {
	return b.named.Name()
}

func (b *tbranch) Title() string {
	return b.named.Title()
}

func (b *tbranch) Class() string {
	return "TBranch"
}

func (b *tbranch) getTree() *ttree {
	return b.tree
}

func (b *tbranch) setTree(t *ttree) {
	b.tree = t
	for _, sub := range b.branches {
		sub.setTree(t)
	}
}

func (b *tbranch) Branches() []Branch {
	return b.branches
}

func (b *tbranch) Leaves() []Leaf {
	return b.leaves
}

func (b *tbranch) Branch(name string) Branch {
	for _, bb := range b.Branches() {
		if bb.Name() == name {
			return bb
		}
	}
	return nil
}

func (b *tbranch) Leaf(name string) Leaf {
	for _, lf := range b.Leaves() {
		if lf.Name() == name {
			return lf
		}
	}
	switch {
	case b.bup != nil:
		return b.bup.Leaf(name)
	case b.btop != nil:
		return b.btop.Leaf(name)
	case b.tree != nil:
		return b.tree.Leaf(name)
	}
	return nil
}

func (b *tbranch) GoType() reflect.Type {
	if len(b.Leaves()) == 1 {
		return b.leaves[0].Type()
	}
	fields := make([]reflect.StructField, len(b.leaves))
	for i, leaf := range b.leaves {
		ft := &fields[i]
		ft.Name = "ROOT_" + leaf.Name()
		etype := leaf.Type()
		switch {
		case leaf.LeafCount() != nil:
			etype = reflect.SliceOf(etype)
		case leaf.Len() > 1 && leaf.Kind() != reflect.String:
			etype = reflect.ArrayOf(leaf.Len(), etype)
		}
		ft.Type = etype
	}
	return reflect.StructOf(fields)
}

func (b *tbranch) getReadEntry() int64 {
	return b.ctx.entry
}

func (b *tbranch) getEntry(i int64) {
	if b.ctx.entry == i {
		return
	}
	err := b.loadEntry(i)
	if err != nil {
		panic(fmt.Errorf("rtree: branch [%s] failed to load entry %d: %w", b.Name(), i, err))
	}
}

func (b *tbranch) MarshalROOT(w *rbytes.WBuffer) (int, error) {
	if w.Err() != nil {
		return 0, w.Err()
	}

	maxBaskets := b.maxBaskets
	defer func() { b.maxBaskets = maxBaskets }()
	b.maxBaskets = max(b.writeBasket+1, defaultMaxBaskets)

	hdr := w.WriteHeader(b.Class(), b.RVersion())
	w.WriteObject(&b.named)
	w.WriteObject(&b.attfill)
	w.WriteI32(int32(b.compress))
	w.WriteI32(int32(b.basketSize))
	w.WriteI32(int32(b.entryOffsetLen))
	w.WriteI32(int32(b.writeBasket))
	w.WriteI64(b.entryNumber)
	w.WriteObject(&b.iobits)
	w.WriteI32(int32(b.offset))
	w.WriteI32(int32(b.maxBaskets))
	w.WriteI32(int32(b.splitLevel))
	w.WriteI64(b.entries)
	w.WriteI64(b.firstEntry)
	w.WriteI64(b.totBytes)
	w.WriteI64(b.zipBytes)

	{
		branches := rcont.NewObjArray()
		if len(b.branches) > 0 {
			elems := make([]root.Object, len(b.branches))
			for i, v := range b.branches {
				elems[i] = v
			}
			branches.SetElems(elems)
		}
		if n, err := branches.MarshalROOT(w); err != nil {
			return n, err
		}
	}
	{
		leaves := rcont.NewObjArray()
		if len(b.leaves) > 0 {
			elems := make([]root.Object, len(b.leaves))
			for i, v := range b.leaves {
				elems[i] = v
			}
			leaves.SetElems(elems)
		}
		if n, err := leaves.MarshalROOT(w); err != nil {
			return n, err
		}
	}
	{
		baskets := rcont.NewObjArray()
		if len(b.baskets) > 0 {
			elems := make([]root.Object, len(b.baskets))
			for i := range b.baskets {
				elems[i] = &b.baskets[i]
			}
			baskets.SetElems(elems)
		}
		if n, err := baskets.MarshalROOT(w); err != nil {
			return n, err
		}
		baskets.SetElems(nil)
	}

	{
		sli := b.basketBytes[:b.writeBasket]
		w.WriteI8(1)
		w.WriteArrayI32(sli)
		if n := b.maxBaskets - len(sli); n > 0 {
			// fill up with zeros.
			w.WriteArrayI32(make([]int32, n))
		}
	}

	{
		sli := b.basketEntry[:b.writeBasket+1]
		w.WriteI8(1)
		w.WriteArrayI64(sli)
		if n := b.maxBaskets - len(sli); n > 0 {
			// fill up with zeros.
			w.WriteArrayI64(make([]int64, n))
		}
	}

	{
		sli := b.basketSeek[:b.writeBasket]
		w.WriteI8(1)
		w.WriteArrayI64(sli)
		if n := b.maxBaskets - len(sli); n > 0 {
			// fill up with zeros.
			w.WriteArrayI64(make([]int64, n))
		}
	}

	w.WriteString(b.fname)

	return w.SetHeader(hdr)
}

// ROOTUnmarshaler is the interface implemented by an object that can
// unmarshal itself from a ROOT buffer
func (b *tbranch) UnmarshalROOT(r *rbytes.RBuffer) error {
	if r.Err() != nil {
		return r.Err()
	}

	hdr := r.ReadHeader(b.Class(), b.RVersion())

	b.tree = nil
	b.ctx.bk = nil
	b.ctx.entry = -1
	b.ctx.first = -1
	b.ctx.next = -1

	const minVers = 6
	switch {
	case hdr.Vers >= 10:

		r.ReadObject(&b.named)
		r.ReadObject(&b.attfill)

		b.compress = int(r.ReadI32())
		b.basketSize = int(r.ReadI32())
		b.entryOffsetLen = int(r.ReadI32())
		b.writeBasket = int(r.ReadI32())
		b.entryNumber = r.ReadI64()
		if hdr.Vers >= 13 {
			if err := b.iobits.UnmarshalROOT(r); err != nil {
				return err
			}
		}
		b.offset = int(r.ReadI32())
		b.maxBaskets = int(r.ReadI32())
		b.splitLevel = int(r.ReadI32())
		b.entries = r.ReadI64()
		if hdr.Vers >= 11 {
			b.firstEntry = r.ReadI64()
		}
		b.totBytes = r.ReadI64()
		b.zipBytes = r.ReadI64()

		{
			var branches rcont.ObjArray
			if err := branches.UnmarshalROOT(r); err != nil {
				return err
			}
			b.branches = make([]Branch, branches.Last()+1)
			for i := range b.branches {
				br := branches.At(i).(Branch)
				b.branches[i] = br
			}
		}

		{
			var leaves rcont.ObjArray
			if err := leaves.UnmarshalROOT(r); err != nil {
				return err
			}
			b.leaves = make([]Leaf, leaves.Last()+1)
			for i := range b.leaves {
				leaf := leaves.At(i).(Leaf)
				leaf.setBranch(b)
				b.leaves[i] = leaf
			}
		}
		{
			var baskets rcont.ObjArray
			if err := baskets.UnmarshalROOT(r); err != nil {
				return err
			}

			b.baskets = make([]Basket, baskets.Last()+1)
			for i := range b.baskets {
				bkt := baskets.At(i)
				// FIXME(sbinet) check why some are nil
				if bkt == nil {
					continue
				}
				bk := bkt.(*Basket)
				b.baskets[i] = *bk
			}
		}

		b.basketBytes = nil
		b.basketEntry = nil
		b.basketSeek = nil

		/*isArray*/
		_ = r.ReadI8()
		b.basketBytes = make([]int32, b.maxBaskets)
		r.ReadArrayI32(b.basketBytes)
		b.basketBytes = b.basketBytes[:b.writeBasket:b.writeBasket]

		/*isArray*/
		_ = r.ReadI8()
		b.basketEntry = make([]int64, b.maxBaskets)
		r.ReadArrayI64(b.basketEntry)
		b.basketEntry = b.basketEntry[: b.writeBasket+1 : b.writeBasket+1]

		/*isArray*/
		_ = r.ReadI8()
		b.basketSeek = make([]int64, b.maxBaskets)
		r.ReadArrayI64(b.basketSeek)
		b.basketSeek = b.basketSeek[:b.writeBasket:b.writeBasket]

		b.fname = r.ReadString()

	case hdr.Vers >= 6:
		r.ReadObject(&b.named)

		if hdr.Vers > 7 {
			r.ReadObject(&b.attfill)
		}

		b.compress = int(r.ReadI32())
		b.basketSize = int(r.ReadI32())
		b.entryOffsetLen = int(r.ReadI32())
		b.writeBasket = int(r.ReadI32())
		b.entryNumber = int64(r.ReadI32())
		b.offset = int(r.ReadI32())
		b.maxBaskets = int(r.ReadI32())
		if hdr.Vers > 6 {
			b.splitLevel = int(r.ReadI32())
		}
		b.entries = int64(r.ReadF64())
		b.totBytes = int64(r.ReadF64())
		b.zipBytes = int64(r.ReadF64())

		{
			var branches rcont.ObjArray
			if err := branches.UnmarshalROOT(r); err != nil {
				return err
			}
			b.branches = make([]Branch, branches.Last()+1)
			for i := range b.branches {
				br := branches.At(i).(Branch)
				b.branches[i] = br
			}
		}

		{
			var leaves rcont.ObjArray
			if err := leaves.UnmarshalROOT(r); err != nil {
				return err
			}
			b.leaves = make([]Leaf, leaves.Last()+1)
			for i := range b.leaves {
				leaf := leaves.At(i).(Leaf)
				leaf.setBranch(b)
				b.leaves[i] = leaf
			}
		}
		{
			var baskets rcont.ObjArray
			if err := baskets.UnmarshalROOT(r); err != nil {
				return err
			}
			b.baskets = make([]Basket, baskets.Last()+1)
			for i := range b.baskets {
				bkt := baskets.At(i)
				// FIXME(sbinet) check why some are nil
				if bkt == nil {
					continue
				}
				bk := bkt.(*Basket)
				b.baskets[i] = *bk
			}
		}

		b.basketBytes = nil
		b.basketEntry = nil
		b.basketSeek = nil

		/*isArray*/
		_ = r.ReadI8()
		b.basketBytes = make([]int32, b.maxBaskets)
		r.ReadArrayI32(b.basketBytes)

		/*isArray*/
		_ = r.ReadI8()
		{
			slice := make([]int32, b.maxBaskets)
			r.ReadArrayI32(slice)
			b.basketEntry = make([]int64, len(slice))
			for i, v := range slice {
				b.basketEntry[i] = int64(v)
			}
		}

		switch r.ReadI8() {
		case 2:
			b.basketSeek = make([]int64, b.maxBaskets)
			r.ReadArrayI64(b.basketSeek)
		default:
			slice := make([]int32, b.maxBaskets)
			r.ReadArrayI32(slice)
			b.basketSeek = make([]int64, len(slice))
			for i, v := range slice {
				b.basketSeek[i] = int64(v)
			}
		}

		b.fname = r.ReadString()

	default:
		panic(fmt.Errorf("rtree: too old TBranch version (%d<%d)", hdr.Vers, minVers))
	}

	if b.splitLevel == 0 && len(b.branches) > 0 {
		b.splitLevel = 1
	}

	r.CheckHeader(hdr)
	return r.Err()
}

func (b *tbranch) loadEntry(ientry int64) error {
	var err error

	if len(b.basketBytes) == 0 {
		return nil
	}

	err = b.loadBasket(ientry)
	if err != nil {
		return err
	}
	b.firstEntry = b.ctx.first

	jentry := ientry - b.ctx.first
	switch len(b.leaves) {
	case 1:
		err = b.ctx.bk.loadLeaf(jentry, b.leaves[0])
		if err != nil {
			return err
		}
	default:
		for _, leaf := range b.leaves {
			err = b.ctx.bk.loadLeaf(jentry, leaf)
			if err != nil {
				return err
			}
		}
	}

	return nil
}

func (b *tbranch) loadBasket(entry int64) error {
	ib := b.findBasketIndex(entry)
	if ib < 0 {
		return fmt.Errorf("rtree: no basket for entry %d", entry)
	}
	var (
		err   error
		bufsz = int(b.basketBytes[ib])
		seek  = b.basketSeek[ib]
		f     = b.tree.getFile()
	)
	b.ctx.id = ib
	b.ctx.entry = entry
	b.ctx.next = b.basketEntry[ib+1]
	b.ctx.first = b.basketEntry[ib]
	if ib < len(b.baskets) {
		b.ctx.bk = &b.baskets[ib]
		if b.ctx.bk.rbuf == nil {
			err = b.ctx.inflate(bufsz, seek, f)
			if err != nil {
				return fmt.Errorf("rtree: could not inflate basket: %w", err)
			}

			return b.setupBasket(&b.ctx, entry)
		}
		return nil
	}

	b.baskets = append(b.baskets, Basket{})
	b.ctx.bk = &b.baskets[len(b.baskets)-1]
	err = b.ctx.inflate(bufsz, seek, f)
	if err != nil {
		return fmt.Errorf("rtree: could not inflate basket: %w", err)
	}
	return b.setupBasket(&b.ctx, entry)
}

func (b *tbranch) findBasketIndex(entry int64) int {
	switch {
	case entry == 0:
		return 0
	case b.ctx.first <= entry && entry < b.ctx.next:
		return b.ctx.id
	}
	/*
		    // binary search is not efficient for small slices (like basketEntry)
			// TODO(sbinet): test at which length of basketEntry it starts to be efficient.
			entries := b.basketEntry[1:]
			i := sort.Search(len(entries), func(i int) bool { return entries[i] >= entry })
			if b.basketEntry[i+1] == entry {
				return i + 1
			}
			return i
	*/

	for i := b.ctx.id; i < len(b.basketEntry); i++ {
		v := b.basketEntry[i]
		if v > entry && v > 0 {
			return i - 1
		}
	}
	if entry == b.basketEntry[len(b.basketEntry)-1] {
		return -2 // len(b.basketEntry) - 1
	}
	return -1
}

func (b *tbranch) setupBasket(ctx *basketCtx, entry int64) error {
	ib := ctx.id
	switch ctx.keylen {
	case 0: // FIXME(sbinet): from trial and error. check this is ok for all cases
		b.basketEntry[ib] = 0
		b.basketEntry[ib+1] = int64(ctx.bk.nevbuf)

	default:
		if b.entryOffsetLen <= 0 {
			return nil
		}

		last := int64(ctx.bk.last)
		ctx.bk.rbuf.SetPos(last)
		n := int(ctx.bk.rbuf.ReadI32())
		ctx.bk.offsets = rbytes.ResizeI32(ctx.bk.offsets, n)
		ctx.bk.rbuf.ReadArrayI32(ctx.bk.offsets)
		if err := ctx.bk.rbuf.Err(); err != nil {
			return err
		}
	}

	return nil
}

func (b *tbranch) setStreamer(s rbytes.StreamerInfo, ctx rbytes.StreamerInfoContext) {
	// no op
}

func (b *tbranch) setStreamerElement(s rbytes.StreamerElement, ctx rbytes.StreamerInfoContext) {
	// no op
}

func (b *tbranch) createNewBasket() {
	cycle := int16(b.writeBasket) + 1
	bk := newBasketFrom(b.tree, b, cycle, b.basketSize, b.entryOffsetLen)
	b.ctx.bk = &bk
	if n := b.writeBasket; n > b.maxBaskets {
		b.maxBaskets = n
	}
}

func (b *tbranch) write() (int, error) {
	b.entries++
	b.entryNumber++

	szOld := b.ctx.bk.wbuf.Len()
	b.ctx.bk.update(szOld)
	_, err := b.writeToBuffer(b.ctx.bk.wbuf)
	szNew := b.ctx.bk.wbuf.Len()
	n := int(szNew - szOld)
	if err != nil {
		return n, fmt.Errorf("could not write to buffer (branch=%q): %w", b.Name(), err)
	}
	if n > b.ctx.bk.nevsize {
		b.ctx.bk.nevsize = n
	}

	// FIXME(sbinet): harmonize or drive via "auto-flush" ?
	if szNew+int64(n) >= int64(b.basketSize) {
		err = b.flush()
		if err != nil {
			return n, fmt.Errorf("could not flush branch (auto-flush): %w", err)
		}

		b.createNewBasket()
	}
	return n, nil
}

func (b *tbranch) writeToBuffer(w *rbytes.WBuffer) (int, error) {
	var tot int
	for i, leaf := range b.leaves {
		n, err := leaf.writeToBuffer(w)
		if err != nil {
			return tot, fmt.Errorf("could not write leaf[%d] name=%q of branch %q: %w", i, leaf.Name(), b.Name(), err)
		}
		tot += n
	}
	return tot, nil
}

func (b *tbranch) flush() error {
	for i, sub := range b.branches {
		err := sub.flush()
		if err != nil {
			return fmt.Errorf("could not flush subbranch[%d]=%q of branch %q: %w", i, sub.Name(), b.Name(), err)
		}
	}

	f := b.tree.getFile()
	totBytes, zipBytes, err := b.ctx.bk.writeFile(f, int32(b.compress))
	if err != nil {
		return fmt.Errorf("could not marshal basket[%d] (branch=%q): %w", b.writeBasket, b.Name(), err)
	}
	b.totBytes += totBytes
	b.zipBytes += zipBytes

	b.basketBytes = append(b.basketBytes, b.ctx.bk.key.Nbytes())
	b.basketEntry = append(b.basketEntry, b.entryNumber)
	b.basketSeek = append(b.basketSeek, b.ctx.bk.key.SeekKey())
	b.writeBasket++
	b.ctx.bk = nil

	return nil
}

// tbranchObject is a Branch for objects.
type tbranchObject struct {
	tbranch
	class string // class name of referenced object
}

func (b *tbranchObject) RVersion() int16 {
	return rvers.BranchObject
}

func (b *tbranchObject) Class() string {
	return "TBranchObject"
}

func (b *tbranchObject) MarshalROOT(w *rbytes.WBuffer) (int, error) {
	if w.Err() != nil {
		return 0, w.Err()
	}

	hdr := w.WriteHeader(b.Class(), b.RVersion())
	w.WriteObject(&b.tbranch)
	w.WriteString(b.class)

	return w.SetHeader(hdr)
}

// ROOTUnmarshaler is the interface implemented by an object that can
// unmarshal itself from a ROOT buffer
func (b *tbranchObject) UnmarshalROOT(r *rbytes.RBuffer) error {
	if r.Err() != nil {
		return r.Err()
	}

	hdr := r.ReadHeader(b.Class(), b.RVersion())

	if hdr.Vers < 1 {
		r.SetErr(fmt.Errorf("rtree: TBranchObject version too old (%d < 8)", hdr.Vers))
		return r.Err()
	}

	r.ReadObject(&b.tbranch)

	for _, leaf := range b.leaves {
		switch leaf := leaf.(type) {
		case *tleaf:
			leaf.branch = b
		case *tleafObject:
			leaf.branch = b
		case *tleafElement:
			leaf.branch = b
		}
	}

	b.class = r.ReadString()

	r.CheckHeader(hdr)
	return r.Err()
}

// tbranchElement is a Branch for objects.
type tbranchElement struct {
	tbranch
	class   string          // class name of referenced object
	parent  string          // name of parent class
	clones  string          // named of class in TClonesArray (if any)
	chksum  uint32          // checksum of class
	clsver  uint16          // version number of class
	id      int32           // element serial number in fInfo
	btype   int32           // branch type
	stype   int32           // branch streamer type
	max     int32           // maximum entries for a TClonesArray or variable array
	stltyp  int32           // STL container type
	bcount1 *tbranchElement // pointer to primary branchcount branch
	bcount2 *tbranchElement // pointer to secondary branchcount branch

	streamer  rbytes.StreamerInfo
	estreamer rbytes.StreamerElement
}

func newBranchElementFromWVar(w *wtree, base *tbranch, wvar WriteVar, parent Branch, lvl int, cfg wopt) (Branch, error) {
	var (
		rv       = reflect.ValueOf(wvar.Value)
		streamer = rdict.StreamerOf(w.ttree.f, reflect.Indirect(rv).Type())
		pclass   = ""
	)

	if parent, ok := parent.(*tbranchElement); ok {
		pclass = parent.class
	}

	b := &tbranchElement{
		tbranch:  *base,
		class:    streamer.Name(),
		parent:   pclass,
		chksum:   uint32(streamer.CheckSum()),
		clsver:   uint16(streamer.ClassVersion()),
		id:       -1,
		btype:    0,
		stype:    -1,
		streamer: streamer,
	}
	switch reflect.TypeOf(wvar.Value).Elem().Kind() {
	case reflect.Struct:
		b.tbranch.entryOffsetLen = 20
	case reflect.Slice:
		b.tbranch.entryOffsetLen = 400
	}

	w.ttree.f.RegisterStreamer(b.streamer)

	_, err := newLeafFromWVar(w, b, wvar, lvl, cfg)
	if err != nil {
		return nil, err
	}
	b.named.SetTitle(wvar.Name)
	b.createNewBasket()
	return b, nil
}

func (b *tbranchElement) RVersion() int16 {
	return rvers.BranchElement
}

func (b *tbranchElement) Class() string {
	return "TBranchElement"
}

func (b *tbranchElement) MarshalROOT(w *rbytes.WBuffer) (int, error) {
	if w.Err() != nil {
		return 0, w.Err()
	}

	hdr := w.WriteHeader(b.Class(), b.RVersion())
	w.WriteObject(&b.tbranch)
	w.WriteString(b.class)
	w.WriteString(b.parent)
	w.WriteString(b.clones)
	w.WriteU32(b.chksum)
	w.WriteU16(b.clsver)
	w.WriteI32(b.id)
	w.WriteI32(b.btype)
	w.WriteI32(b.stype)
	w.WriteI32(b.max)

	{
		var obj root.Object
		if b.bcount1 != nil {
			obj = b.bcount1
		}
		w.WriteObjectAny(obj)
	}
	{
		var obj root.Object
		if b.bcount2 != nil {
			obj = b.bcount2
		}
		w.WriteObjectAny(obj)
	}

	return w.SetHeader(hdr)
}

// ROOTUnmarshaler is the interface implemented by an object that can
// unmarshal itself from a ROOT buffer
func (b *tbranchElement) UnmarshalROOT(r *rbytes.RBuffer) error {
	if r.Err() != nil {
		return r.Err()
	}

	hdr := r.ReadHeader(b.Class(), b.RVersion())

	if hdr.Vers < 1 {
		r.SetErr(fmt.Errorf("rtree: TBranchElement version too old (%d < 8)", hdr.Vers))
		return r.Err()
	}

	r.ReadObject(&b.tbranch)

	for _, leaf := range b.leaves {
		switch leaf := leaf.(type) {
		case *tleaf:
			leaf.branch = b
		case *tleafObject:
			leaf.branch = b
		case *tleafElement:
			leaf.branch = b
		}
	}

	b.class = r.ReadString()
	if hdr.Vers > 1 {
		b.parent = r.ReadString()
		b.clones = r.ReadString()
		b.chksum = r.ReadU32()
	}
	if hdr.Vers >= 10 {
		b.clsver = r.ReadU16()
	} else {
		b.clsver = uint16(r.ReadU32())
	}
	b.id = r.ReadI32()
	b.btype = r.ReadI32()
	b.stype = r.ReadI32()
	if hdr.Vers > 1 {
		b.max = r.ReadI32()

		bcount1 := r.ReadObjectAny()
		if bcount1 != nil {
			b.bcount1 = bcount1.(*tbranchElement)
		}

		bcount2 := r.ReadObjectAny()
		if bcount2 != nil {
			b.bcount2 = bcount2.(*tbranchElement)
		}
	}

	r.CheckHeader(hdr)
	return r.Err()
}

func (b *tbranchElement) loadEntry(ientry int64) error {
	if len(b.branches) > 0 {
		for _, sub := range b.branches {
			err := sub.loadEntry(ientry)
			if err != nil {
				return err
			}
		}
	}
	return b.tbranch.loadEntry(ientry)
}

func (b *tbranchElement) setupReadStreamer(sictx rbytes.StreamerInfoContext) error {
	streamer, err := sictx.StreamerInfo(b.class, int(b.clsver))
	if err != nil {
		streamer, err = sictx.StreamerInfo(b.class, -1)
		if err != nil {
			return fmt.Errorf("rtree: no StreamerInfo for class=%q version=%d checksum=%d", b.class, b.clsver, b.chksum)
		}
	}
	b.streamer = streamer

	for _, leaf := range b.tbranch.leaves {
		leaf, ok := leaf.(*tleafElement)
		if !ok {
			continue
		}
		leaf.streamers = b.streamer.Elements()
	}

	for _, sub := range b.branches {
		sub, ok := sub.(*tbranchElement)
		if !ok {
			continue
		}
		err := sub.setupReadStreamer(sictx)
		if err != nil {
			return err
		}
	}

	return nil
}

func (b *tbranchElement) GoType() reflect.Type {
	typ, err := rdict.TypeFromSI(b.tree.getFile(), b.streamer)
	if err != nil {
		panic(err)
	}
	return typ
}

func (b *tbranchElement) setStreamer(s rbytes.StreamerInfo, ctx rbytes.StreamerInfoContext) {
	b.streamer = s
	if len(b.tbranch.leaves) == 1 {
		typ, err := rdict.TypeFromSI(ctx, s)
		if err != nil {
			panic(err)
		}
		tle := b.tbranch.leaves[0].(*tleafElement)
		tle.streamers = s.Elements()
		tle.src = reflect.New(typ).Elem()
	}
	err := b.setupReadStreamer(ctx)
	if err != nil {
		panic(err)
	}
}

func (b *tbranchElement) setStreamerElement(se rbytes.StreamerElement, ctx rbytes.StreamerInfoContext) {
	b.estreamer = se
	if len(b.Leaves()) == 1 {
		tle := b.Leaves()[0].(*tleafElement)
		tle.streamers = []rbytes.StreamerElement{se}
		typ, err := rdict.TypeFromSE(ctx, se)
		if err != nil {
			panic(err)
		}
		tle.src = reflect.New(typ).Elem()
	}
	err := b.setupReadStreamer(ctx)
	if err != nil {
		panic(err)
	}
}

func (b *tbranchElement) write() (int, error) {
	b.entries++
	b.entryNumber++

	szOld := b.ctx.bk.wbuf.Len()
	b.ctx.bk.update(szOld)
	_, err := b.writeToBuffer(b.ctx.bk.wbuf)
	szNew := b.ctx.bk.wbuf.Len()
	n := int(szNew - szOld)
	if err != nil {
		return n, fmt.Errorf("could not write to buffer (branch=%q): %w", b.Name(), err)
	}
	if n > b.ctx.bk.nevsize {
		b.ctx.bk.grow(n)
	}

	// FIXME(sbinet): harmonize or drive via "auto-flush" ?
	if szNew+int64(n) >= int64(b.basketSize) {
		err = b.flush()
		if err != nil {
			return n, fmt.Errorf("could not flush branch (auto-flush): %w", err)
		}

		b.createNewBasket()
	}
	return n, nil
}

func (b *tbranchElement) writeToBuffer(w *rbytes.WBuffer) (int, error) {
	var tot int
	for i, leaf := range b.leaves {
		n, err := leaf.writeToBuffer(w)
		if err != nil {
			return tot, fmt.Errorf("could not write leaf[%d] name=%q of branch %q: %w", i, leaf.Name(), b.Name(), err)
		}
		tot += n
	}
	return tot, nil
}

func btopOf(b Branch) Branch {
	if b == nil {
		return nil
	}
	const max = 1<<31 - 1
	for range max {
		switch bb := b.(type) {
		case *tbranch:
			if bb.bup == nil {
				return bb
			}
			b = bb.bup
		case *tbranchElement:
			if bb.bup == nil {
				return bb
			}
			b = bb.bup
		default:
			panic(fmt.Errorf("rtree: unknown branch type %T", b))
		}
	}
	panic("impossible")
}

func init() {
	{
		f := func() reflect.Value {
			o := &tbranch{}
			return reflect.ValueOf(o)
		}
		rtypes.Factory.Add("TBranch", f)
	}
	{
		f := func() reflect.Value {
			o := &tbranchObject{}
			return reflect.ValueOf(o)
		}
		rtypes.Factory.Add("TBranchObject", f)
	}
	{
		f := func() reflect.Value {
			o := &tbranchElement{}
			return reflect.ValueOf(o)
		}
		rtypes.Factory.Add("TBranchElement", f)
	}
}

type basketCtx struct {
	id    int     // current basket number when reading
	entry int64   // current entry number when reading
	first int64   // first entry in the current basket
	next  int64   // next entry that will require us to go to TBranchElement next basket
	bk    *Basket // pointer to the current basket
	buf   []byte  // scratch space for the current basket

	keylen uint32
}

func (ctx *basketCtx) inflate(bufsz int, seek int64, f *riofs.File) error {
	ctx.buf = rbytes.ResizeU8(ctx.buf, bufsz)
	ctx.bk.key.SetFile(f)
	ctx.keylen = 0

	var (
		bk    = ctx.bk
		buf   = ctx.buf[:bufsz]
		sictx = f
		err   error
	)

	switch {
	case len(buf) == 0 && ctx.bk != nil: // FIXME(sbinet): from trial and error. check this is ok for all cases

	default:
		_, err = f.ReadAt(buf, seek)
		if err != nil {
			return fmt.Errorf("rtree: could not read basket buffer from file: %w", err)
		}

		err = bk.UnmarshalROOT(rbytes.NewRBuffer(buf, nil, 0, sictx))
		if err != nil {
			return fmt.Errorf("rtree: could not unmarshal basket buffer from file: %w", err)
		}

		ctx.keylen = uint32(bk.key.KeyLen())
	}

	buf = rbytes.ResizeU8(buf, int(bk.key.ObjLen()))
	_, err = bk.key.Load(buf)
	if err != nil {
		return err
	}

	bk.rbuf = rbytes.NewRBuffer(buf, nil, ctx.keylen, sictx)
	return nil
}

var (
	_ root.Object        = (*tbranch)(nil)
	_ root.Named         = (*tbranch)(nil)
	_ Branch             = (*tbranch)(nil)
	_ rbytes.Marshaler   = (*tbranch)(nil)
	_ rbytes.Unmarshaler = (*tbranch)(nil)

	_ root.Object        = (*tbranchObject)(nil)
	_ root.Named         = (*tbranchObject)(nil)
	_ Branch             = (*tbranchObject)(nil)
	_ rbytes.Marshaler   = (*tbranchObject)(nil)
	_ rbytes.Unmarshaler = (*tbranchObject)(nil)

	_ root.Object        = (*tbranchElement)(nil)
	_ root.Named         = (*tbranchElement)(nil)
	_ Branch             = (*tbranchElement)(nil)
	_ rbytes.Marshaler   = (*tbranchElement)(nil)
	_ rbytes.Unmarshaler = (*tbranchElement)(nil)
)