Jpp  18.0.0-rc.3
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Public Member Functions | Public Attributes | Protected Attributes | Private Attributes | List of all members
JCOMPAREHISTOGRAMS::JTestKolmogorov_1D Class Reference

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

#include <JTestKolmogorov_1D.hh>

Inheritance diagram for JCOMPAREHISTOGRAMS::JTestKolmogorov_1D:
JCOMPAREHISTOGRAMS::JTest_t JCOMPAREHISTOGRAMS::JTestKolmogorov_t

Public Member Functions

 JTestKolmogorov_1D ()
 Default constructor. More...
 
std::istream & read (std::istream &in) override
 Read test parameters from input. More...
 
void test (TObject *o1, TObject *o2) override
 Applies Kolmogorov test for two ROOT TH1 histograms. More...
 
std::ostream & write (std::ostream &out, const char delimiter= ' ', const bool onlyFailures=false) const
 Write test result to output. More...
 
virtual void save (TFile *f, const std::string &path, const bool onlyFailures=false)
 Writes the test result to root file. More...
 
virtual void clear ()
 Clear results. More...
 
const std::string & getTestName ()
 Get test name. More...
 
const std::string & getResultType ()
 Get result type. 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
 Buffer to store results of multiple tests.*/. More...
 

Protected Attributes

const std::string testName
 test name More...
 
const std::string resultType
 result type More...
 

Private Attributes

double threshold
 threshold p-value to decide if test is passed. More...
 

Detailed Description

Implementation of the Kolmogorov test for 1D histograms.


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 uses the input parameter threshold() to evaluate whether the test is passed or failed.
The evaluation is done by comparing the threshold value with the result of the JKolmogorovTest() test. This is a p-value.
The parameter threshold should therefore be a real value between 0 and 1.

Definition at line 24 of file JTestKolmogorov_1D.hh.

Constructor & Destructor Documentation

JCOMPAREHISTOGRAMS::JTestKolmogorov_1D::JTestKolmogorov_1D ( )
inline

Default constructor.

Definition at line 33 of file JTestKolmogorov_1D.hh.

33  :
34  JTest_t("KS_1D", "p-Value(KS)"),
35  JTestKolmogorov_t()
36  {}
JCOMPAREHISTOGRAMS::JTestKolmogorov_t::JTestKolmogorov_t
JTestKolmogorov_t()
Default constructor.
Definition: JTestKolmogorov_t.hh:27
JCOMPAREHISTOGRAMS::JTest_t::JTest_t
JTest_t(const std::string &name, const std::string &type)
Constructor.
Definition: JTest_t.hh:57

Member Function Documentation

std::istream& JCOMPAREHISTOGRAMS::JTestKolmogorov_1D::read ( std::istream &  in)
inlineoverridevirtual

Read test parameters from input.

Parameters
ininput stream
Returns
input stream

Implements JCOMPAREHISTOGRAMS::JTest_t.

Definition at line 44 of file JTestKolmogorov_1D.hh.

44  {
45  return in >> threshold;
46  };
JCOMPAREHISTOGRAMS::JTestKolmogorov_1D::threshold
double threshold
threshold p-value to decide if test is passed.
Definition: JTestKolmogorov_1D.hh:70
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 JAcoustics sh $DETECTOR_ID source JAcousticsToolkit sh CHECK_EXIT_CODE typeset A EMITTERS get_tripods $WORKDIR tripod txt EMITTERS get_transmitters $WORKDIR transmitter txt EMITTERS for EMITTER in
Definition: JCanberra.sh:46
void JCOMPAREHISTOGRAMS::JTestKolmogorov_1D::test ( TObject *  o1,
TObject *  o2 
)
inlineoverridevirtual

Applies Kolmogorov test for two ROOT TH1 histograms.

Parameters
o1First histogram
o2Second histogram

Implements JCOMPAREHISTOGRAMS::JTest_t.

Definition at line 54 of file JTestKolmogorov_1D.hh.

54  {
55 
56  using namespace std;
57 
58  if (!(dynamic_cast<TH2*>(o1) == NULL) || !(dynamic_cast<TH2*>(o2) == NULL)) {
59  ERROR("For 2D histograms call JTestKolmogorov_2D: " << o1->GetName() << endl);
60 
61  } else if (!(dynamic_cast<TH1*>(o1) == NULL) && !(dynamic_cast<TH1*>(o2) == NULL)) {
62 
63  TH1D* h1 = dynamic_cast<TH1D*>(o1);
64  TH1D* h2 = dynamic_cast<TH1D*>(o2);
65 
66  JTestResult r = JKolmogorovTest(h1, h2, threshold, testName, resultType);
67 
68  results.push_back(r);
69  }
70  };
JTOOLS::r
data_type r[M+1]
Definition: JPolint.hh:779
JCOMPAREHISTOGRAMS::JTestKolmogorov_t::JKolmogorovTest
JTestResult JKolmogorovTest(TH1 *h1, TH1 *h2, double threshold, std::string testName, std::string parameterName)
Kolmogorov test for 1D histograms.
Definition: JTestKolmogorov_t.hh:43
JCOMPAREHISTOGRAMS::JTest_t::results
std::vector< JTestResult > results
Buffer to store results of multiple tests.*/.
Definition: JTest_t.hh:162
JCOMPAREHISTOGRAMS::JTestKolmogorov_1D::threshold
double threshold
threshold p-value to decide if test is passed.
Definition: JTestKolmogorov_1D.hh:70
ERROR
#define ERROR(A)
Definition: JMessage.hh:66
JCOMPAREHISTOGRAMS::JTest_t::resultType
const std::string resultType
result type
Definition: JTest_t.hh:167
JCOMPAREHISTOGRAMS::JTest_t::testName
const std::string testName
test name
Definition: JTest_t.hh:166
std::ostream& JCOMPAREHISTOGRAMS::JTest_t::write ( std::ostream &  out,
const char  delimiter = ' ',
const bool  onlyFailures = false 
) const
inlineinherited

Write test result to output.

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

Definition at line 81 of file JTest_t.hh.

83  {
84 
85  for (std::vector<JTestResult>::const_iterator r = results.begin() ; r != results.end() ; ++r) {
86 
87  if (onlyFailures && r->passed) {
88  continue;
89  }
90 
91  print(out, *r, delimiter, true);
92  }
93 
94  return out;
95  }
JCOMPAREHISTOGRAMS::print
std::ostream & print(std::ostream &out, const JTestSummary &summary, const char delimiter= ' ', const bool useColors=true)
Print test summary.
Definition: JTestSummary.hh:255
JTOOLS::r
data_type r[M+1]
Definition: JPolint.hh:779
JCOMPAREHISTOGRAMS::JTest_t::results
std::vector< JTestResult > results
Buffer to store results of multiple tests.*/.
Definition: JTest_t.hh:162
virtual void JCOMPAREHISTOGRAMS::JTest_t::save ( TFile *  f,
const std::string &  path,
const bool  onlyFailures = false 
)
inlinevirtualinherited

Writes the test result to root file.

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

Definition at line 113 of file JTest_t.hh.

115  {
116 
117  using namespace std;
118 
119  if (f -> GetDirectory(path.c_str())==0) {
120  f->mkdir(path.c_str());
121  }
122 
123  f->cd(path.c_str());
124 
125  for (vector<JTestResult>::const_iterator r = results.begin() ; r != results.end() ; ++r){
126 
127  if (!onlyFailures || !r->passed){
128  r->obj->Write();
129  }
130  }
131  }
f
o $QUALITY_ROOT d $DEBUG!CHECK_EXIT_CODE JPlot1D f
Definition: JDataQuality.sh:76
JTOOLS::r
data_type r[M+1]
Definition: JPolint.hh:779
JCOMPAREHISTOGRAMS::JTest_t::results
std::vector< JTestResult > results
Buffer to store results of multiple tests.*/.
Definition: JTest_t.hh:162
path
$WORKDIR driver txt done cat $WORKDIR driver txt<< EOFprocess ${DATAFILTER}$FILTER_HOST csh-c '(setenv ROOTSYS $ROOTSYS &&source $JPP_DIR/setenv.csh $JPP_DIR &&(JDataFilter-H\$SERVER\$-M\$LOGGER\$-d $DEBUG-u ${DATAFILTER}-P $PORT</dev/null > &/dev/null &))';process ${DATAWRITER}$WRITER_HOST csh-c '(setenv ROOTSYS $ROOTSYS &&source $JPP_DIR/setenv.csh $JPP_DIR &&(JDataWriter-H\$SERVER\$-M\$LOGGER\$-d $DEBUG-u ${DATAWRITER}</dev/null > &/dev/null &))';print enterevent ev_init{RC_CMD}event ev_reset{RC_CMD}event ev_init{RC_CMD}event ev_configure{RC_DFLTR%<$WORKDIR/ev_configure_datafilter.txt > RC_DQSIM<$WORKDIR/ev_configure_dqsimulator.txt > RC_DWRT path
Definition: software/JDataFilter/JDataFilter.sh:199
virtual void JCOMPAREHISTOGRAMS::JTest_t::clear ( )
inlinevirtualinherited

Clear results.

Definition at line 137 of file JTest_t.hh.

137  {
138  results.clear();
139  }
JCOMPAREHISTOGRAMS::JTest_t::results
std::vector< JTestResult > results
Buffer to store results of multiple tests.*/.
Definition: JTest_t.hh:162
const std::string& JCOMPAREHISTOGRAMS::JTest_t::getTestName ( )
inlineinherited

Get test name.

Returns
test name

Definition at line 147 of file JTest_t.hh.

147  {
148  return testName;
149  }
JCOMPAREHISTOGRAMS::JTest_t::testName
const std::string testName
test name
Definition: JTest_t.hh:166
const std::string& JCOMPAREHISTOGRAMS::JTest_t::getResultType ( )
inlineinherited

Get result type.

Returns
result type

Definition at line 157 of file JTest_t.hh.

157  {
158  return resultType;
159  }
JCOMPAREHISTOGRAMS::JTest_t::resultType
const std::string resultType
result type
Definition: JTest_t.hh:167
JTestResult JCOMPAREHISTOGRAMS::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 43 of file JTestKolmogorov_t.hh.

43  {
44 
45  using namespace std;
46  using namespace JPP;
47 
48  int n1 = h1 -> GetNbinsX();
49  int n2 = h2 -> GetNbinsX();
50 
51  if(n1 != n2)
52  ERROR("Histograms with different bining. The objects: " << h1 -> GetName() << " can not be compared." << endl);
53 
54  TH2D* h3 = (TH2D*)h1->Clone(h1->GetName()==h2->GetName() ?
55  MAKE_CSTRING(to_string(h1->GetName())) :
56  MAKE_CSTRING(to_string(h1->GetName()) + "_VS_" + to_string(h2->GetName())));
57 
58  bool afunc1 = false;
59  bool afunc2 = false;
60  double s1 = 1./h1->Integral();
61  double s2 = 1./h2->Integral();
62 
63  double ew1, ew2, w1 = 0, w2 = 0;
64 
65  for (int bin = 1; bin <= n1; ++bin){
66  ew1 = h1->GetBinError(bin);
67  ew2 = h2->GetBinError(bin);
68  w1 += ew1*ew1;
69  w2 += ew2*ew2;
70  }
71 
72  double esum1 = 0, esum2 = 0;
73  if (w1 > 0)
74  esum1 = 1./s1/s1/w1;
75  else
76  afunc1 = true;
77  if (w2 > 0)
78  esum2 = 1./s2/s2/w2;
79  else
80  afunc2 = true;
81  if (afunc2 && afunc1) {
82  ERROR("Errors are zero for both histograms\n");
83  }
84 
85  double c1 = 0, c2 = 0;
86 
87  double dmax = 0;
88  for (int bin=1 ; bin<=n1 ; ++bin){
89  c1 += s1*h1->GetBinContent(bin);
90  c2 += s2*h2->GetBinContent(bin);
91  double d = TMath::Abs(c1-c2)*TMath::Sqrt(esum1*esum2/(esum1+esum2));
92  double p = TMath::KolmogorovProb(d);
93  h3->SetBinContent(bin,p);
94  dmax = TMath::Max(dmax,TMath::Abs(c1-c2));
95  }
96 
97  double z;
98 
99  if (afunc1)
100  z = dmax*TMath::Sqrt(esum2);
101  else if (afunc2)
102  z = dmax*TMath::Sqrt(esum1);
103  else
104  z = dmax*TMath::Sqrt(esum1*esum2/(esum1+esum2));
105 
106  double pValue = TMath::KolmogorovProb(z);
107 
108  bool passed;
109 
110  (pValue < threshold ? passed = false : passed = true);
111 
112  JResultTitle title(testName, parameterName, passed , pValue);
113 
114  h3->SetTitle(title.getTitle().c_str());
115 
116  JTestResult r (testName,
117  JRootObjectID(MAKE_STRING(h1->GetDirectory()->GetPath() << h1->GetName())),
118  JRootObjectID(MAKE_STRING(h2->GetDirectory()->GetPath() << h1->GetName())),
119  parameterName, pValue, threshold, h3, passed);
120 
121  return r;
122  };
MAKE_CSTRING
#define MAKE_CSTRING(A)
Make C-string.
Definition: JPrint.hh:136
JTOOLS::r
data_type r[M+1]
Definition: JPolint.hh:779
MAKE_STRING
#define MAKE_STRING(A)
Make string.
Definition: JPrint.hh:127
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:47
JLANG::to_string
std::string to_string(const T &value)
Convert value to string.
Definition: JLangToolkit.hh:207
JTestResult JCOMPAREHISTOGRAMS::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 143 of file JTestKolmogorov_t.hh.

143  {
144 
145  using namespace std;
146  using namespace JPP;
147 
148  int nFailures = 0;
149 
150  JTestResult r;
151 
152  if(slice == 'x' || slice == 'X'){
153 
154  int nSlices1 = h1->GetNbinsX();
155  int nSlices2 = h2->GetNbinsX();
156 
157  TH1* h3 = h1->ProjectionX(MAKE_CSTRING(h1->GetName() << "_VS_" <<
158  h2->GetName() << "_KolmogorovTestSliceX"));
159 
160  if(nSlices1 != nSlices2) {
161  ERROR("Histograms with different binning. The objects: " << h1->GetName() << " and " << h2->GetName() << " can not be compared." << endl);
162  }
163 
164  for (int i=1 ; i<=nSlices1 ; ++i){
165 
166  std::string sliceName = MAKE_STRING(h3->GetName() << "_" << to_string(i));
167 
168  TH1D* s1 = h1->ProjectionY (sliceName.c_str(),i,i);
169  TH1D* s2 = h2->ProjectionY (sliceName.c_str(),i,i);
170 
171  if (!(s1->GetSumOfWeights() > 0 && s2->GetSumOfWeights() > 0)) { continue; }
172 
173  double p = s1 -> KolmogorovTest (s2);
174 
175  bool passed = (p < threshold ? false : true);
176 
177  if (!passed) nFailures++;
178 
179  h3->SetBinContent(i,p);
180 
181  }
182 
183  bool passed = (nFailures/nSlices1 > failuresThreshold ? false : true);
184 
185  JResultTitle title(testName, parameterName, passed, nFailures);
186 
187  h3->SetTitle(title.getTitle().c_str());
188 
189  r = JTestResult (testName,
190  JRootObjectID(MAKE_STRING(h1->GetDirectory()->GetPath() << h1->GetName())),
191  JRootObjectID(MAKE_STRING(h2->GetDirectory()->GetPath() << h1->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(MAKE_CSTRING(h1->GetName() << "_VS_" <<
200  h2->GetName() << "_KolmogorovTestSliceY"));
201 
202  if(nSlices1 != nSlices2) {
203  ERROR("Histograms with different binning. The objects: " << h1->GetName() << " can not be compared." << endl);
204  }
205 
206  for (int i=1 ; i<=nSlices1 ; ++i){
207 
208  std::string sliceName = MAKE_STRING(h3->GetName() << "_" << 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  if (!(s1->GetSumOfWeights() > 0 && s2->GetSumOfWeights() > 0)) { continue; }
214 
215  double p = s1 -> KolmogorovTest (s2);
216 
217  bool passed = (p < threshold ? false : true);
218 
219  if (!passed) nFailures++;
220 
221  h3->SetBinContent(i,p);
222  }
223 
224  bool passed = (nFailures/nSlices1 > failuresThreshold ? false : true);
225 
226  JResultTitle title(testName, parameterName, passed, nFailures);
227 
228  h3->SetTitle(title.getTitle().c_str());
229 
230  r = JTestResult (testName,
231  JRootObjectID(MAKE_STRING(h1->GetDirectory()->GetPath() << h1->GetName())),
232  JRootObjectID(MAKE_STRING(h2->GetDirectory()->GetPath() << h1->GetName())),
233  parameterName, nFailures, failuresThreshold, h3, passed);
234 
235  }
236 
237  return r;
238  };
MAKE_CSTRING
#define MAKE_CSTRING(A)
Make C-string.
Definition: JPrint.hh:136
JTOOLS::r
data_type r[M+1]
Definition: JPolint.hh:779
MAKE_STRING
#define MAKE_STRING(A)
Make string.
Definition: JPrint.hh:127
ERROR
#define ERROR(A)
Definition: JMessage.hh:66
string
then awk string
Definition: pre-calibration.sh:45
JLANG::to_string
std::string to_string(const T &value)
Convert value to string.
Definition: JLangToolkit.hh:207
JTestResult JCOMPAREHISTOGRAMS::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 256 of file JTestKolmogorov_t.hh.

256  {
257 
258  using namespace std;
259  using namespace JPP;
260 
261  int n1x = h1 -> GetNbinsX();
262  int n2x = h2 -> GetNbinsX();
263  int n1y = h1 -> GetNbinsY();
264  int n2y = h2 -> GetNbinsY();
265 
266  if(n1x != n2x || n1y != n2y)
267  ERROR("Histograms with different bining. The objects: " << h1 -> GetName() << " can not be compared." << endl);
268 
269  if(h1->Integral()==0 || h2->Integral()==0)
270  ERROR("Empty histogram: " << h1 -> GetName() << " can not be compared." << endl);
271 
272  double s1 = 1./h1->Integral();
273  double s2 = 1./h2->Integral();
274 
275  TH2D* h3 = (TH2D*)h1->Clone(h1->GetName()==h2->GetName() ?
276  MAKE_CSTRING(to_string(h1->GetName())) :
277  MAKE_CSTRING(to_string(h1->GetName()) + "_VS_" + to_string(h2->GetName())));
278 
279  bool afunc1 = false;
280  bool afunc2 = false;
281 
282  double ew1, ew2, w1 = 0, w2 = 0;
283 
284  for (int i = 1; i <= n1x; ++i){
285  for (int j = 1; j <= n1y; ++j){
286  ew1 = h1->GetBinError(i,j);
287  ew2 = h2->GetBinError(i,j);
288  w1 += ew1*ew1;
289  w2 += ew2*ew2;
290  }
291  }
292 
293  double esum1 = 0, esum2 = 0;
294  if (w1 > 0)
295  esum1 = 1./s1/s1/w1;
296  else
297  afunc1 = true;
298  if (w2 > 0)
299  esum2 = 1./s2/s2/w2;
300  else
301  afunc2 = true;
302  if (afunc2 && afunc1) {
303  ERROR("Errors are zero for both histograms\n");
304  }
305 
306  double c1 = 0, c2 = 0;
307 
308  double dmax1 = 0;
309  for (int i=1 ; i<=n1x ; ++i){
310  for (int j=1 ; j<=n1y ; ++j){
311  c1 += s1*h1->GetBinContent(i,j);
312  c2 += s2*h2->GetBinContent(i,j);
313  double d = TMath::Abs(c1-c2)*TMath::Sqrt(esum1*esum2/(esum1+esum2));
314  h3->Fill(i,j,d);
315  dmax1 = TMath::Max(dmax1,TMath::Abs(c1-c2));
316  }
317  }
318 
319  c1 = 0, c2 = 0;
320 
321  double dmax2 = 0;
322  for (int j=1 ; j<=n1y ; ++j){
323  for (int i=1 ; i<=n1x ; ++i){
324  c1 += s1*h1->GetBinContent(i,j);
325  c2 += s2*h2->GetBinContent(i,j);
326  double d = TMath::Abs(c1-c2)*TMath::Sqrt(esum1*esum2/(esum1+esum2));
327  h3->Fill(i,j,d);
328  dmax1 = TMath::Max(dmax2,TMath::Abs(c1-c2));
329  }
330  }
331 
332  double dmax = 0.5*(dmax1+dmax2);
333  double z;
334 
335  if (afunc1)
336  z = dmax*TMath::Sqrt(esum2);
337  else if (afunc2)
338  z = dmax*TMath::Sqrt(esum1);
339  else
340  z = dmax*TMath::Sqrt(esum1*esum2/(esum1+esum2));
341 
342  double pValue = TMath::KolmogorovProb(z);
343 
344  for (int i=1 ; i<=n1x ; ++i){
345  for (int j=1 ; j<=n1y ; ++j){
346  h3->SetBinContent(i,j,TMath::KolmogorovProb(0.5 * h3->GetBinContent(i,j)));
347  }
348  }
349 
350  bool passed;
351 
352  (pValue < threshold ? passed = false : passed = true);
353 
354  JResultTitle title(testName, parameterName, passed , pValue);
355 
356  h3->SetTitle(title.getTitle().c_str());
357 
358  JTestResult r (testName,
359  JRootObjectID(MAKE_STRING(h1->GetDirectory()->GetPath() << h1->GetName())),
360  JRootObjectID(MAKE_STRING(h2->GetDirectory()->GetPath() << h1->GetName())),
361  parameterName, pValue, threshold, h3, passed);
362 
363  return r;
364  };
MAKE_CSTRING
#define MAKE_CSTRING(A)
Make C-string.
Definition: JPrint.hh:136
JTOOLS::r
data_type r[M+1]
Definition: JPolint.hh:779
MAKE_STRING
#define MAKE_STRING(A)
Make string.
Definition: JPrint.hh:127
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:47
JLANG::to_string
std::string to_string(const T &value)
Convert value to string.
Definition: JLangToolkit.hh:207
JTOOLS::j
int j
Definition: JPolint.hh:703

Member Data Documentation

double JCOMPAREHISTOGRAMS::JTestKolmogorov_1D::threshold
private

threshold p-value to decide if test is passed.

Definition at line 70 of file JTestKolmogorov_1D.hh.

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

Buffer to store results of multiple tests.*/.

Definition at line 162 of file JTest_t.hh.

const std::string JCOMPAREHISTOGRAMS::JTest_t::testName
protectedinherited

test name

Definition at line 166 of file JTest_t.hh.

const std::string JCOMPAREHISTOGRAMS::JTest_t::resultType
protectedinherited

result type

Definition at line 167 of file JTest_t.hh.

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