Jpp
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Public Member Functions | Public Attributes | Private Attributes | List of all members
JTestKolmogorovSlice Class Reference

Implementation of the Kolmogorov test for 2D histograms. More...

#include <JTestKolmogorovSlice.hh>

Inheritance diagram for JTestKolmogorovSlice:
JTestKolmogorov_t JTest_t

Public Member Functions

 JTestKolmogorovSlice ()
 Default constructor. More...
 
std::istream & read (std::istream &in) override
 Read test parameters from input. More...
 
std::ostream & write (std::ostream &out, std::string delimiter=" ", bool onlyFailures=false) const override
 Write test result to output. More...
 
void test (TObject *o1, TObject *o2) override
 Applies Kolmogorov test for two ROOT TH2 histograms. More...
 
virtual void save (TFile *f, std::string path, bool onlyFailures=false) override
 Writes the test result to root file. More...
 
void clear () override
 Clear results. More...
 
std::string getName () override
 Get test name. More...
 
std::string getParameter () override
 Get parameter name. More...
 
JTestResult JKolmogorovTest (TH1 *h1, TH1 *h2, double threshold, std::string testName, std::string parameterName)
 Kolmogorov test for 1D histograms. More...
 
JTestResult JKolmogorovTestSlice (TH2 *h1, TH2 *h2, double threshold, double failuresThreshold, std::string testName, std::string parameterName, char slice)
 Kolmogorov test for sliced 2D histograms. More...
 
JTestResult JKolmogorovTest2D (TH2 *h1, TH2 *h2, double threshold, std::string testName, std::string parameterName)
 Kolmogorov test for 2D histograms. More...
 

Public Attributes

std::vector< JTestResult > results
 

Private Attributes

double threshold
 threshold p-value to decide if test is passed. More...
 
double failuresThreshold
 threshold p-value to decide if test is passed. More...
 
char slice
 axis to slice. x or X for x-axis, y or Y for y-axis, n or N for None. More...
 
const std::string name = "KS_Slice"
 test name. More...
 
const std::string parameter = "p-Value(KS)"
 parameter name. More...
 

Detailed Description

Implementation of the Kolmogorov test for 2D histograms.

Author
rgruiz
This class is derived from the abstract class JTest_t(). For a general description of the implementation of this and other tests derived from JTest_t(), see its documentation.
This test compares two 2D histograms. If the parameter slice() equals x, X y or Y, the histograms are sliced along the corresponding axis, and the JKolmogorovTest() test is applied to each slice.
If slice() equals n or N, the histograms are not sliced, and JKolmogorovTest2D() is applied.
The input parameter threshold(), is used to evaluate whether the test is passed or failed for each slice or for the full 2D distribution.
The evaluation is done by comparing the threshold() value with the result produced by JKolmogorovTest() or JKolmogorovTest2D(), which is a p-value.
The parameter threshold() should therefore be a real value between 0 and 1.

Definition at line 23 of file JTestKolmogorovSlice.hh.

Constructor & Destructor Documentation

JTestKolmogorovSlice::JTestKolmogorovSlice ( )
inline

Default constructor.

Definition at line 31 of file JTestKolmogorovSlice.hh.

31 {}

Member Function Documentation

std::istream& JTestKolmogorovSlice::read ( std::istream &  in)
inlineoverridevirtual

Read test parameters from input.

Parameters
ininput stream
Returns
input stream

Implements JTest_t.

Definition at line 39 of file JTestKolmogorovSlice.hh.

39  {
40  return in >> threshold >> slice >> failuresThreshold;
41  };
JTestKolmogorovSlice::slice
char slice
axis to slice. x or X for x-axis, y or Y for y-axis, n or N for None.
Definition: JTestKolmogorovSlice.hh:145
JTestKolmogorovSlice::failuresThreshold
double failuresThreshold
threshold p-value to decide if test is passed.
Definition: JTestKolmogorovSlice.hh:144
JTestKolmogorovSlice::threshold
double threshold
threshold p-value to decide if test is passed.
Definition: JTestKolmogorovSlice.hh:143
in
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable INPUT_FILE $argv[2] eval JPrintDetector a $DETECTOR O IDENTIFIER eval JPrintDetector a $DETECTOR O SUMMARY source JAcoustics sh $DETECTOR_ID typeset A TRIPODS get_tripods $WORKDIR tripod txt TRIPODS for EMITTER in
Definition: JCanberra.sh:36
std::ostream& JTestKolmogorovSlice::write ( std::ostream &  out,
std::string  delimiter = " ",
bool  onlyFailures = false 
) const
inlineoverridevirtual

Write test result to output.

Parameters
outoutput stream
delimiterfield delimiter
onlyFailuresIf true, write only failures.
Returns
output stream

Implements JTest_t.

Definition at line 51 of file JTestKolmogorovSlice.hh.

51  {
52  for (std::vector<JTestResult>::const_iterator r = results.begin() ; r != results.end() ; ++r) {
53  if (onlyFailures){
54  if (!r->passed)
55  out << std::scientific << std::setprecision(2) << (r->passed ? JPP::GREEN : JPP::RED) << r->print(delimiter) << std::endl;
56  }else{
57  out << std::scientific << std::setprecision(2) << (r->passed ? JPP::GREEN : JPP::RED) << r->print(delimiter) << std::endl;
58  }
59  }
60  return out;
61  };
JTest_t::results
std::vector< JTestResult > results
Definition: JTest_t.hh:251
JTOOLS::r
data_type r[M+1]
Definition: JPolint.hh:742
void JTestKolmogorovSlice::test ( TObject *  o1,
TObject *  o2 
)
inlineoverridevirtual

Applies Kolmogorov test for two ROOT TH2 histograms.

Parameters
o1First histogram
o2Second histogram

Implements JTest_t.

Definition at line 69 of file JTestKolmogorovSlice.hh.

69  {
70 
71  using namespace std;
72  using namespace JPP;
73 
74  if (!(dynamic_cast<TH2*>(o1) == NULL) && !(dynamic_cast<TH2*>(o2) == NULL)) {
75 
76  TH2D* h1 = dynamic_cast<TH2D*>(o1);
77  TH2D* h2 = dynamic_cast<TH2D*>(o2);
78 
79  if(slice == 'x' || slice == 'X' || slice == 'y' || slice == 'Y'){
80 
81  JTestResult r = JKolmogorovTestSlice(h1, h2, threshold, failuresThreshold, name, parameter, slice);
82 
83  results.push_back(r);
84 
85  } else {
86  ERROR ("Slice option should be: x, X, y or Y, and not " + slice);
87  }
88 
89  }else if (!(dynamic_cast<TH1*>(o1) == NULL) && !(dynamic_cast<TH1*>(o2) == NULL)) {
90  ERROR("For 1D histograms call JTestKolmogorov_1D: " << o1->GetName() << endl);
91  }
92  };
JTestKolmogorovSlice::parameter
const std::string parameter
parameter name.
Definition: JTestKolmogorovSlice.hh:147
JTestKolmogorovSlice::slice
char slice
axis to slice. x or X for x-axis, y or Y for y-axis, n or N for None.
Definition: JTestKolmogorovSlice.hh:145
JTestKolmogorovSlice::name
const std::string name
test name.
Definition: JTestKolmogorovSlice.hh:146
JTestKolmogorovSlice::failuresThreshold
double failuresThreshold
threshold p-value to decide if test is passed.
Definition: JTestKolmogorovSlice.hh:144
JTest_t::results
std::vector< JTestResult > results
Definition: JTest_t.hh:251
h1
then for HISTOGRAM in h0 h1
Definition: JMatrixNZ.sh:69
JTOOLS::r
data_type r[M+1]
Definition: JPolint.hh:742
JTestKolmogorov_t::JKolmogorovTestSlice
JTestResult JKolmogorovTestSlice(TH2 *h1, TH2 *h2, double threshold, double failuresThreshold, std::string testName, std::string parameterName, char slice)
Kolmogorov test for sliced 2D histograms.
Definition: JTestKolmogorov_t.hh:139
JTestResult
Structure containing the result of the test.
Definition: JTest_t.hh:164
JTestKolmogorovSlice::threshold
double threshold
threshold p-value to decide if test is passed.
Definition: JTestKolmogorovSlice.hh:143
ERROR
#define ERROR(A)
Definition: JMessage.hh:66
virtual void JTestKolmogorovSlice::save ( TFile *  f,
std::string  path,
bool  onlyFailures = false 
)
inlineoverridevirtual

Writes the test result to root file.

Parameters
fA ROOT file
pathPath in root file.
onlyFailuresIf true, write only failures.

Implements JTest_t.

Definition at line 100 of file JTestKolmogorovSlice.hh.

100  {
101 
102  using namespace std;
103 
104  for (vector<JTestResult>::const_iterator r = results.begin() ; r != results.end() ; ++r){
105 
106  if (onlyFailures){
107  if (!r->passed){
108  if (f -> GetDirectory(path.c_str())==0) f->mkdir(path.c_str());
109  f->cd(path.c_str());
110  r->obj->Write();
111  }
112  }else{
113  if (f -> GetDirectory(path.c_str())==0) f->mkdir(path.c_str());
114  f->cd(path.c_str());
115  r->obj->Write();
116  }
117  }
118  }
JTest_t::results
std::vector< JTestResult > results
Definition: JTest_t.hh:251
f
do set_array DAQHEADER JPrintDAQHeader f
Definition: JTuneHV.sh:79
JTOOLS::r
data_type r[M+1]
Definition: JPolint.hh:742
path
do cat driver txt<< EOFevent ev_configure{RC_EVT%< ev_configure.txt > RC_DWRT path
Definition: software/JDataWriter/JDataWriter.sh:143
void JTestKolmogorovSlice::clear ( )
inlineoverridevirtual

Clear results.

Implements JTest_t.

Definition at line 123 of file JTestKolmogorovSlice.hh.

123  {
124  results.clear();
125  }
JTest_t::results
std::vector< JTestResult > results
Definition: JTest_t.hh:251
std::string JTestKolmogorovSlice::getName ( )
inlineoverridevirtual

Get test name.

Implements JTest_t.

Definition at line 130 of file JTestKolmogorovSlice.hh.

130  {
131  return name;
132  }
JTestKolmogorovSlice::name
const std::string name
test name.
Definition: JTestKolmogorovSlice.hh:146
std::string JTestKolmogorovSlice::getParameter ( )
inlineoverridevirtual

Get parameter name.

Implements JTest_t.

Definition at line 137 of file JTestKolmogorovSlice.hh.

137  {
138  return parameter;
139  }
JTestKolmogorovSlice::parameter
const std::string parameter
parameter name.
Definition: JTestKolmogorovSlice.hh:147
JTestResult JTestKolmogorov_t::JKolmogorovTest ( TH1 *  h1,
TH1 *  h2,
double  threshold,
std::string  testName,
std::string  parameterName 
)
inlineinherited

Kolmogorov test for 1D histograms.


The input parameter threshold, is used to evaluate whether the test is passed or failed.
The evaluation is done by comparing the threshold value with the result of the Kolmogorov test. The output of a Kolmogorov test is a p-value.
The parameter threshold should therefore be a real value between 0 and 1.

Parameters
h1First histogram
h2Second historgram
thresholdp-value
parameterNameName of the parameter used to test the histograms
testNameName of the test used to compare the histograms
Returns
Test result

Definition at line 37 of file JTestKolmogorov_t.hh.

37  {
38 
39  using namespace std;
40  using namespace JPP;
41 
42  int n1 = h1 -> GetNbinsX();
43  int n2 = h2 -> GetNbinsX();
44 
45  if(n1 != n2)
46  ERROR("Histograms with different bining. The objects: " << h1 -> GetName() << " can not be compared." << endl);
47 
48  TH2D* h3 = (TH2D*)h1->Clone(h1->GetName()==h2->GetName() ?
49  MAKE_CSTRING(to_string(h1->GetName())) :
50  MAKE_CSTRING(to_string(h1->GetName()) + "_VS_" + to_string(h2->GetName())));
51 
52  bool afunc1 = false;
53  bool afunc2 = false;
54  double s1 = 1./h1->Integral();
55  double s2 = 1./h2->Integral();
56 
57  double ew1, ew2, w1 = 0, w2 = 0;
58 
59  for (int bin = 1; bin <= n1; ++bin){
60  ew1 = h1->GetBinError(bin);
61  ew2 = h2->GetBinError(bin);
62  w1 += ew1*ew1;
63  w2 += ew2*ew2;
64  }
65 
66  double esum1 = 0, esum2 = 0;
67  if (w1 > 0)
68  esum1 = 1./s1/s1/w1;
69  else
70  afunc1 = true;
71  if (w2 > 0)
72  esum2 = 1./s2/s2/w2;
73  else
74  afunc2 = true;
75  if (afunc2 && afunc1) {
76  ERROR("Errors are zero for both histograms\n");
77  }
78 
79  double c1 = 0, c2 = 0;
80 
81  double dmax = 0;
82  for (int bin=1 ; bin<=n1 ; ++bin){
83  c1 += s1*h1->GetBinContent(bin);
84  c2 += s2*h2->GetBinContent(bin);
85  double d = TMath::Abs(c1-c2)*TMath::Sqrt(esum1*esum2/(esum1+esum2));
86  double p = TMath::KolmogorovProb(d);
87  h3->SetBinContent(bin,p);
88  dmax = TMath::Max(dmax,TMath::Abs(c1-c2));
89  }
90 
91  double z;
92 
93  if (afunc1)
94  z = dmax*TMath::Sqrt(esum2);
95  else if (afunc2)
96  z = dmax*TMath::Sqrt(esum1);
97  else
98  z = dmax*TMath::Sqrt(esum1*esum2/(esum1+esum2));
99 
100  double pValue = TMath::KolmogorovProb(z);
101 
102  bool passed;
103 
104  (pValue < threshold ? passed = false : passed = true);
105 
106  JResultTitle title(testName, parameterName, passed , pValue);
107 
108  h3->SetTitle(title.getTitle().c_str());
109 
110  JTestResult r (testName,
111  string (h1->GetDirectory()->GetPath()).append(h1->GetName()),
112  string (h2->GetDirectory()->GetPath()).append(h2->GetName()),
113  h1->GetDirectory()->GetFile()->GetName(),
114  h2->GetDirectory()->GetFile()->GetName(),
115  parameterName, pValue, threshold, h3, passed);
116 
117  return r;
118  };
JResultTitle
Class dedicated to standardize the title of the graphical objects produced by the JTest_t() derived c...
Definition: JTest_t.hh:22
MAKE_CSTRING
#define MAKE_CSTRING(A)
Make C-string.
Definition: JPrint.hh:151
h1
then for HISTOGRAM in h0 h1
Definition: JMatrixNZ.sh:69
JTOOLS::r
data_type r[M+1]
Definition: JPolint.hh:742
JTestResult
Structure containing the result of the test.
Definition: JTest_t.hh:164
ERROR
#define ERROR(A)
Definition: JMessage.hh:66
c1
TCanvas * c1
Global variables to handle mouse events.
Definition: JDrawModule3D.cc:34
d
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:45
JLANG::to_string
std::string to_string(const T &value)
Convert value to string.
Definition: JLangToolkit.hh:192
JTestResult JTestKolmogorov_t::JKolmogorovTestSlice ( TH2 *  h1,
TH2 *  h2,
double  threshold,
double  failuresThreshold,
std::string  testName,
std::string  parameterName,
char  slice 
)
inlineinherited

Kolmogorov test for sliced 2D histograms.


The histograms are sliced along the axis specified by the slice parameter. A slice per bin is made.
For each of the slices, the input parameter threshold is used to evaluate whether the test is passed or failed.
The evaluation is done by comparing the threshold value with the result of the Kolmogorov test. The output of the Kolmogorov test is a p-value.
The parameter threshold should therefore be a real value between 0 and 1. The fraction of failed tests is compared to the input parameter failuresThreshold. If this fraction is larger than failuresThreshold, the test fails.

Parameters
h1First histogram
h2Second histogram
thresholdThreshold value for the test result
failuresThresholdThreshold value for the test result
parameterNameName of the parameter used to test the histograms
testNameName of the test used to compare the histograms
sliceThe axis along which the histogram is sliced.
Returns
Test result.

Definition at line 139 of file JTestKolmogorov_t.hh.

139  {
140 
141  using namespace std;
142  using namespace JPP;
143 
144  int nFailures = 0;
145 
146  JTestResult r;
147 
148  if(slice == 'x' || slice == 'X'){
149 
150  int nSlices1 = h1->GetNbinsX();
151  int nSlices2 = h2->GetNbinsX();
152 
153  TH1* h3 = h1->ProjectionX(h1->GetName()==h2->GetName() ?
154  MAKE_CSTRING(to_string(h1->GetName())) :
155  MAKE_CSTRING(to_string(h1->GetName()) + "_VS_" + to_string(h2->GetName())));
156 
157  if(nSlices1 != nSlices2)
158  ERROR("Histograms with different binning. The objects: " << h1->GetName() << " and " << h2->GetName() << " can not be compared." << endl);
159 
160  for (int i=1 ; i<=nSlices1 ; ++i){
161 
162  std::string sliceName = MAKE_STRING(h1->GetName() + to_string("_") + to_string(i));
163 
164  TH1D* s1 = h1->ProjectionY (sliceName.c_str(),i,i);
165  TH1D* s2 = h2->ProjectionY (sliceName.c_str(),i,i);
166 
167  double p = s1 -> KolmogorovTest (s2);
168 
169  bool passed;
170 
171  (p < threshold ? passed = false : passed = true);
172 
173  if (!passed) nFailures++;
174 
175  h3->SetBinContent(i,p);
176 
177  }
178 
179  bool passed;
180 
181  (nFailures/nSlices1 > failuresThreshold ? passed = false : passed = true);
182 
183  JResultTitle title(testName, parameterName, passed , nFailures);
184 
185  h3->SetTitle(title.getTitle().c_str());
186 
187  r = JTestResult (testName,
188  string (h1->GetDirectory()->GetPath()).append(h1->GetName()),
189  string (h2->GetDirectory()->GetPath()).append(h2->GetName()),
190  h1->GetDirectory()->GetFile()->GetName(),
191  h2->GetDirectory()->GetFile()->GetName(),
192  parameterName, nFailures, failuresThreshold, h3, passed);
193 
194  }else if (slice == 'y' || slice == 'Y'){
195 
196  int nSlices1 = h1->GetNbinsY();
197  int nSlices2 = h2->GetNbinsY();
198 
199  TH1* h3 = h1->ProjectionY(h1->GetName()==h2->GetName() ?
200  MAKE_CSTRING(to_string(h1->GetName())) :
201  MAKE_CSTRING(to_string(h1->GetName()) + "_VS_" + to_string(h2->GetName())));
202 
203  if(nSlices1 != nSlices2)
204  ERROR("Histograms with different binning. The objects: " << h1->GetName() << " can not be compared." << endl);
205 
206  for (int i=1 ; i<=nSlices1 ; ++i){
207 
208  std::string sliceName = MAKE_STRING(h1->GetName() + to_string("_") + to_string(i));
209 
210  TH1D* s1 = h1->ProjectionX (sliceName.c_str(),i,i);
211  TH1D* s2 = h2->ProjectionX (sliceName.c_str(),i,i);
212 
213  double p = s1 -> KolmogorovTest (s2);
214 
215  bool passed;
216 
217  (p < threshold ? passed = false : passed = true);
218 
219  if (!passed) nFailures++;
220 
221  h3->SetBinContent(i,p);
222  }
223 
224  bool passed;
225 
226  (nFailures/nSlices1 > failuresThreshold ? passed = false : passed = true);
227 
228  JResultTitle title(testName, parameterName, passed , nFailures);
229 
230  h3->SetTitle(title.getTitle().c_str());
231 
232  r = JTestResult (testName,
233  string (h1->GetDirectory()->GetPath()).append(h1->GetName()),
234  string (h2->GetDirectory()->GetPath()).append(h2->GetName()),
235  h1->GetDirectory()->GetFile()->GetName(),
236  h2->GetDirectory()->GetFile()->GetName(),
237  parameterName, nFailures, failuresThreshold, h3, passed);
238 
239  }
240 
241  return r;
242  };
JResultTitle
Class dedicated to standardize the title of the graphical objects produced by the JTest_t() derived c...
Definition: JTest_t.hh:22
MAKE_CSTRING
#define MAKE_CSTRING(A)
Make C-string.
Definition: JPrint.hh:151
h1
then for HISTOGRAM in h0 h1
Definition: JMatrixNZ.sh:69
JTOOLS::r
data_type r[M+1]
Definition: JPolint.hh:742
JTestResult
Structure containing the result of the test.
Definition: JTest_t.hh:164
MAKE_STRING
#define MAKE_STRING(A)
Make string.
Definition: JPrint.hh:142
ERROR
#define ERROR(A)
Definition: JMessage.hh:66
JLANG::to_string
std::string to_string(const T &value)
Convert value to string.
Definition: JLangToolkit.hh:192
JTestResult JTestKolmogorov_t::JKolmogorovTest2D ( TH2 *  h1,
TH2 *  h2,
double  threshold,
std::string  testName,
std::string  parameterName 
)
inlineinherited

Kolmogorov test for 2D histograms.


The input parameter threshold, is used to evaluate whether the test is passed or failed.
The evaluation is done by comparing the threshold value with the result of the Kolmogorov test. The output of a Kolmogorov test is a p-value.
The parameter threshold should therefore be a real value between 0 and 1.

Parameters
h1First histogram
h2Second historgram
thresholdp-value
parameterNameName of the parameter used to test the histograms
testNameName of the test used to compare the histograms
Returns
Test result

Definition at line 260 of file JTestKolmogorov_t.hh.

260  {
261 
262  using namespace std;
263  using namespace JPP;
264 
265  int n1x = h1 -> GetNbinsX();
266  int n2x = h2 -> GetNbinsX();
267  int n1y = h1 -> GetNbinsY();
268  int n2y = h2 -> GetNbinsY();
269 
270  if(n1x != n2x || n1y != n2y)
271  ERROR("Histograms with different bining. The objects: " << h1 -> GetName() << " can not be compared." << endl);
272 
273  if(h1->Integral()==0 || h2->Integral()==0)
274  ERROR("Empty histogram: " << h1 -> GetName() << " can not be compared." << endl);
275 
276  double s1 = 1./h1->Integral();
277  double s2 = 1./h2->Integral();
278 
279  TH2D* h3 = (TH2D*)h1->Clone(h1->GetName()==h2->GetName() ?
280  MAKE_CSTRING(to_string(h1->GetName())) :
281  MAKE_CSTRING(to_string(h1->GetName()) + "_VS_" + to_string(h2->GetName())));
282 
283  bool afunc1 = false;
284  bool afunc2 = false;
285 
286  double ew1, ew2, w1 = 0, w2 = 0;
287 
288  for (int i = 1; i <= n1x; ++i){
289  for (int j = 1; j <= n1y; ++j){
290  ew1 = h1->GetBinError(i,j);
291  ew2 = h2->GetBinError(i,j);
292  w1 += ew1*ew1;
293  w2 += ew2*ew2;
294  }
295  }
296 
297  double esum1 = 0, esum2 = 0;
298  if (w1 > 0)
299  esum1 = 1./s1/s1/w1;
300  else
301  afunc1 = true;
302  if (w2 > 0)
303  esum2 = 1./s2/s2/w2;
304  else
305  afunc2 = true;
306  if (afunc2 && afunc1) {
307  ERROR("Errors are zero for both histograms\n");
308  }
309 
310  double c1 = 0, c2 = 0;
311 
312  double dmax1 = 0;
313  for (int i=1 ; i<=n1x ; ++i){
314  for (int j=1 ; j<=n1y ; ++j){
315  c1 += s1*h1->GetBinContent(i,j);
316  c2 += s2*h2->GetBinContent(i,j);
317  double d = TMath::Abs(c1-c2)*TMath::Sqrt(esum1*esum2/(esum1+esum2));
318  h3->Fill(i,j,d);
319  dmax1 = TMath::Max(dmax1,TMath::Abs(c1-c2));
320  }
321  }
322 
323  c1 = 0, c2 = 0;
324 
325  double dmax2 = 0;
326  for (int j=1 ; j<=n1y ; ++j){
327  for (int i=1 ; i<=n1x ; ++i){
328  c1 += s1*h1->GetBinContent(i,j);
329  c2 += s2*h2->GetBinContent(i,j);
330  double d = TMath::Abs(c1-c2)*TMath::Sqrt(esum1*esum2/(esum1+esum2));
331  h3->Fill(i,j,d);
332  dmax1 = TMath::Max(dmax2,TMath::Abs(c1-c2));
333  }
334  }
335 
336  double dmax = 0.5*(dmax1+dmax2);
337  double z;
338 
339  if (afunc1)
340  z = dmax*TMath::Sqrt(esum2);
341  else if (afunc2)
342  z = dmax*TMath::Sqrt(esum1);
343  else
344  z = dmax*TMath::Sqrt(esum1*esum2/(esum1+esum2));
345 
346  double pValue = TMath::KolmogorovProb(z);
347 
348  for (int i=1 ; i<=n1x ; ++i){
349  for (int j=1 ; j<=n1y ; ++j){
350  h3->SetBinContent(i,j,TMath::KolmogorovProb(0.5 * h3->GetBinContent(i,j)));
351  }
352  }
353 
354  bool passed;
355 
356  (pValue < threshold ? passed = false : passed = true);
357 
358  JResultTitle title(testName, parameterName, passed , pValue);
359 
360  h3->SetTitle(title.getTitle().c_str());
361 
362  JTestResult r (testName,
363  string (h1->GetDirectory()->GetPath()).append(h1->GetName()),
364  string (h2->GetDirectory()->GetPath()).append(h2->GetName()),
365  h1->GetDirectory()->GetFile()->GetName(),
366  h2->GetDirectory()->GetFile()->GetName(),
367  parameterName, pValue, threshold, h3, passed);
368 
369  return r;
370  };
JResultTitle
Class dedicated to standardize the title of the graphical objects produced by the JTest_t() derived c...
Definition: JTest_t.hh:22
MAKE_CSTRING
#define MAKE_CSTRING(A)
Make C-string.
Definition: JPrint.hh:151
h1
then for HISTOGRAM in h0 h1
Definition: JMatrixNZ.sh:69
JTOOLS::r
data_type r[M+1]
Definition: JPolint.hh:742
JTestResult
Structure containing the result of the test.
Definition: JTest_t.hh:164
ERROR
#define ERROR(A)
Definition: JMessage.hh:66
c1
TCanvas * c1
Global variables to handle mouse events.
Definition: JDrawModule3D.cc:34
d
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
Definition: JMuonPath.sh:45
JLANG::to_string
std::string to_string(const T &value)
Convert value to string.
Definition: JLangToolkit.hh:192
JTOOLS::j
int j
Definition: JPolint.hh:666

Member Data Documentation

double JTestKolmogorovSlice::threshold
private

threshold p-value to decide if test is passed.

Definition at line 143 of file JTestKolmogorovSlice.hh.

double JTestKolmogorovSlice::failuresThreshold
private

threshold p-value to decide if test is passed.

Definition at line 144 of file JTestKolmogorovSlice.hh.

char JTestKolmogorovSlice::slice
private

axis to slice. x or X for x-axis, y or Y for y-axis, n or N for None.

Definition at line 145 of file JTestKolmogorovSlice.hh.

const std::string JTestKolmogorovSlice::name = "KS_Slice"
private

test name.

Definition at line 146 of file JTestKolmogorovSlice.hh.

const std::string JTestKolmogorovSlice::parameter = "p-Value(KS)"
private

parameter name.

Definition at line 147 of file JTestKolmogorovSlice.hh.

std::vector<JTestResult> JTest_t::results
inherited

Buffer to store results of multiple tests.

Definition at line 251 of file JTest_t.hh.

The documentation for this class was generated from the following file:
Generated by   doxygen 1.8.5