Ratio of two Tprofile histograms

See the possible options here: ROOT: TRatioPlot Class Reference
(diff, diffsig …)

Yes now it looks fine I guess. Yes, I will check. Thank you very much!

The question is … which options make sense for TProfile histograms (e.g., which options will return “properly” calculated errors)?

I tried these.
Ratio2.C (848 Bytes)

A->ProjectionX();
What this does?
Also, I want to know what A->Sumw2(); does?

Not sure about it. @moneta ?

https://root.cern/doc/master/classTProfile.html#a4d3fa02873a069643c90918d6c5a29bd

https://root.cern/doc/master/classTProfile.html#a9591898784d69e20b76a1841a58bac07

Hi,
You cannot perform a direct ratio of two TProfile’s histogram using TRatioPlot or TProfile::Divide. The resulting ratio of two TProfile’s is not anymore a TProfile, so the correct way of computing the ratio is to convert first the TProfile in a TH1 object using TProfile::ProjectionX()and then use TH1::Divide or TRatioPlot.
I will add an error message in TProfile::Divide explaining this.
However, note also that computing the correct error on a ratio of means is tricky, using a normal approximation as it is done in TH1::Divide can give also incorrect results.
If one is really interest in the correct confidence interval in the ratio, it should be estimated using simulation (i.e with bootstrapping)

Thank you for raising this issue.

Lorenzo

1 Like

This topic was automatically closed 14 days after the last reply. New replies are no longer allowed.