66 const TH1* h1 =
dynamic_cast<const TH1*
>(o1);
67 const TH1* h2 =
dynamic_cast<const TH1*
>(o2);
69 if (h1 == NULL || h2 == NULL) {
73 if(h1->GetNbinsX() != h2->GetNbinsX() ||
74 h1->GetNbinsY() != h2->GetNbinsY() ||
75 h1->GetNbinsZ() != h2->GetNbinsZ()) {
76 THROW(
JValueOutOfRange,
"JTestSignificance::test(): Histograms with different bining. The objects: " <<
77 h1->GetName() <<
" and " << h2->GetName() <<
" can not be compared." << endl);
81 K = h2->GetEntries() / h1->GetEntries();
84 TH1* h3 = (TH1*) h1->Clone(h1->GetName() == h2->GetName() ?
92 for (
int i=1 ; i <= h1->GetNbinsX() ; ++i){
93 for (
int j=1 ; j <= h1->GetNbinsY() ; ++j){
94 for (
int k=1 ; k <= h1->GetNbinsZ() ; ++k){
96 const double a = h1->GetBinContent(i,j,k);
97 const double b = h2->GetBinContent(i,j,k);
100 S += fabs((a -
K*b)/sqrt(a +
K*
K*b));
106 S /= (h1->GetNbinsX()*h1->GetNbinsY()*h1->GetNbinsZ());
112 h3->SetTitle(title.
getTitle().c_str());
114 const int Ndims = h3->GetDimension();
118 }
else if (Ndims == 2) {