Hi Adam,
I wrote a series of robust estimators, including Median, Mode
(discrete and peak fit), and the corresponding measures of spread:
ADev (sum of absolute deviations), and FWHM (with adjustable height
instead of just half max); and the corresponding routines to turn a 2-
d histogram into a 1-d histograms of the center and spread, like
TH2::ProfileX(), and TH2::FitSlicesY(). There are also normalization
functions. They can be found in 'blast/exp/analysis/utils/hist.C',
which is in the root path, and the names should be self-explanatory.
So instead of filling a TProfile, fill a regular TH2F, and then
call ProfileMedianY(), ProfileModeY(), or FitPeaksY() on the
histogram, similar to the TH2 methods.
--Chris
TA-53/MPF-1/D111 P-23 MS H803
LANL, Los Alamos, NM 87545
505-665-9804(o) 665-4121(f) 662-0639(h)
On May 4, 2006, at 05:50:23, Adam DeGrush wrote:
>
> Hey Chris,
>
> You said that you had written something that does what TProfile
> does except that it uses the maximum in the bin instead of the bin
> average. Can you send it to me? It sounds like it could be useful
> in the future.
>
> Hope all is well,
> Adam
This archive was generated by hypermail 2.1.2 : Mon Feb 24 2014 - 14:07:33 EST