56 const TH2* h1 =
dynamic_cast<const TH2*
>(o1);
57 const TH2* h2 =
dynamic_cast<const TH2*
>(o2);
59 if (h1 == NULL || h2 == NULL) {
65 const char*
const h3name = (h1->GetName() == h2->GetName() ?
74 const int nSlices1 = h1->GetNbinsX();
75 const int nSlices2 = h2->GetNbinsX();
77 if (nSlices1 != nSlices2) {
78 THROW(
JValueOutOfRange,
"JTestRuns_Slice2D::test(): Histograms with different binning. The objects: " << h1->GetName() <<
" and " << h2->GetName() <<
" can not be compared." << endl);
83 h3 = h1->ProjectionX(h3name);
85 for (
int i=1 ; i<=nSlices1 ; ++i){
87 const std::string sliceName =
MAKE_STRING(h3->GetName() <<
"_" << to_string(i));
89 const TH1* s1 = h1->ProjectionY (sliceName.c_str(),i,i);
90 const TH1* s2 = h2->ProjectionY (sliceName.c_str(),i,i);
92 if (!(s1->GetSumOfWeights() > 0 && s2->GetSumOfWeights() > 0)) {
continue; }
94 const double R = s1->Integral();
95 const double T = s2->Integral();
97 bool a = ((T/R)*s1->GetBinContent(1) - s2->GetBinContent(1)) < 0;
100 double p = (a ? 1 : 0);
101 double q = (a ? 0 : 1);
103 for (
int i = 2 ; i<s1->GetNbinsX() ; ++i){
105 const bool b = ((T/R)*s1->GetBinContent(i) - s2->GetBinContent(i)) < 0;
115 const double N = 1 + 2*p*q/(p+q) ;
116 const double s = sqrt( 2*p*q*(2*p*q-p-q)/(p+q)/(p+q)/(p+q-1) );
117 const double d = (n-N)/s;
121 h3->SetBinContent(i,fabs(d));
126 const int nSlices1 = h1->GetNbinsX();
127 const int nSlices2 = h2->GetNbinsX();
129 if (nSlices1 != nSlices2) {
130 THROW(
JValueOutOfRange,
"JTestRuns_Slice2D::test(): Histograms with different binning. The objects: " << h1->GetName() <<
" and " << h2->GetName() <<
" can not be compared." << endl);
135 h3 = h1->ProjectionY(h3name);
137 for (
int i=1 ; i<=nSlices1 ; ++i){
139 const std::string sliceName =
MAKE_STRING(h3->GetName() <<
"_" << to_string(i));
141 const TH1* s1 = h1->ProjectionX (sliceName.c_str(),i,i);
142 const TH1* s2 = h2->ProjectionX (sliceName.c_str(),i,i);
144 if (!(s1->GetSumOfWeights() > 0 && s2->GetSumOfWeights() > 0)) {
continue; }
146 const double R = s1->Integral();
147 const double T = s2->Integral();
149 bool a = ((T/R)*s1->GetBinContent(1) - s2->GetBinContent(1)) < 0;
152 double p = (a ? 1 : 0);
153 double q = (a ? 0 : 1);
155 for (
int i = 2 ; i<s1->GetNbinsX() ; ++i){
157 const bool b = ((T/R)*s1->GetBinContent(i) - s2->GetBinContent(i)) < 0;
167 const double N = 1 + 2*p*q/(p+q) ;
168 const double s = sqrt( 2*p*q*(2*p*q-p-q)/(p+q)/(p+q)/(p+q-1) );
169 const double d = (n-N)/s;
173 h3->SetBinContent(i,fabs(d));
181 h3->SetTitle(title.
getTitle().c_str());