pyROOT
gwpyのFrequencySeriesをROOTのTH1にする
1 mean = FrequencySeries(...)
2 error = FrequencySeries(...)
3
4 nbin = mean.size
5 value = mean.value
6 freqs = mean.frequencies.value
7 df = mean.df.value
8 fmin = freqs[0] - df*0.5
9 fmax = freqs[-1]+ df*0.5
10
11 h = TH1F('h', ';Frequency [Hz]', nbin, fmin, fmax)
12 h.Sumw2()
13 h.SetContent(array('d', np.insert( value, 0,0).tolist() ) )
14 h.SetError( array('d', np.insert(error.value, 0,0).tolist() ) )
SetContent()とSetError()はゼロbin目(アンダーフロー)から詰めるため、ndarrayの先頭にゼロを追加しておく必要がある。