Esempio n. 1
0
func tree0(f *hio.File) {
	t, err := f.CreateTuple("events")
	if err != nil {
		panic(err)
	}

	// fill some events with random numbers
	nevents := *evtmax
	for iev := 0; iev != nevents; iev++ {
		if iev%1000 == 0 {
			fmt.Printf(":: processing event %d...\n", iev)
		}

		e := Event{}
		// the two energies follow a gaussian distribution
		ea, eb := croot.GRandom.Rannord()
		e.A.E = ea
		e.B.E = eb

		e.A.T = croot.GRandom.Rndm(1)
		e.B.T = e.A.T * croot.GRandom.Gaus(0., 1.)

		if iev%1000 == 0 {
			fmt.Printf("evt.a.e= %8.3f\n", e.A.E)
			fmt.Printf("evt.a.t= %8.3f\n", e.A.T)
			fmt.Printf("evt.b.e= %8.3f\n", e.B.E)
			fmt.Printf("evt.b.t= %8.3f\n", e.B.T)
		}
		err = t.Write(e)
		if err != nil {
			panic(err)
		}
	}
}
Esempio n. 2
0
func tree0(f *hio.File) {
	t, err := f.OpenTuple("events")
	if err != nil {
		panic(err)
	}

	// fill some events with random numbers
	nevents := int64(*evtmax)
	if nevents < 0 {
		nevents = t.Entries()
	}
	if nevents > t.Entries() {
		nevents = t.Entries()
	}
	for iev := int64(0); iev != nevents; iev++ {
		if iev%1000 == 0 || true {
			fmt.Printf(":: processing event %d...\n", iev)
		}
		e := Event{}
		err = t.Read(&e)
		if err != nil {
			panic(err)
		}
		if iev%1000 == 0 || true {
			fmt.Printf("ievt= %d\n", iev)
			fmt.Printf("evt.a.e= %8.3f\n", e.A.E)
			fmt.Printf("evt.a.t= %8.3f\n", e.A.T)
			fmt.Printf("evt.b.e= %8.3f\n", e.B.E)
			fmt.Printf("evt.b.t= %8.3f\n", e.B.T)
		}
	}
}