1 #ifndef __JTESTSIGNIFICANCE_T__
2 #define __JTESTSIGNIFICANCE_T__
56 if(h1 -> GetNbinsX() != h2 -> GetNbinsX())
57 ERROR(
"Histograms with different bining. The objects: " << h1 -> GetName() <<
" can not be compared." << endl);
60 T* h3 = (
T*)h1->Clone(h1->GetName()==h2->GetName() ?
68 for (
int i=0 ; i < h1->GetNbinsX() ; ++i){
69 for (
int j=0 ;
j< h1->GetNbinsY() ; ++
j){
70 for (
int k=0 ;
k< h1->GetNbinsY() ; ++
k){
72 double a = h1->GetBinContent(i+1,
j+1,
k+1);
73 double b = h2->GetBinContent(i+1,
j+1,
k+1);
76 S += fabs((a - K*b)/sqrt(a + K*K*b));
82 S /= (h1->GetNbinsX()*h1->GetNbinsY()*h1->GetNbinsZ());
86 (
S < threshold ? passed =
false : passed =
true);
90 h3->SetTitle(title.
getTitle().c_str());
93 string (h1->GetDirectory()->GetPath()).append(h1->GetName()),
94 string (h2->GetDirectory()->GetPath()).append(h2->GetName()),
95 h1->GetDirectory()->GetFile()->GetName(),
96 h2->GetDirectory()->GetFile()->GetName(),
97 parameterName,
S , threshold, h3, passed);
static const uint32_t K[64]
Class dedicated to standardize the title of the graphical objects produced by the JTest_t() derived c...
JTestResult JSignificanceTest(T *h1, T *h2, double threshold, double K, std::string testName, std::string parameterName)
Significance test for histograms.
then fatal No sound hydrophone file $HYDROPHONE_TXT fi JGraph f $HYDROPHONE_TXT o $HYDROPHONE_ROOT sort gr k
JTestSignificance_t()
Default constructor.
#define MAKE_CSTRING(A)
Make C-string.
then for HISTOGRAM in h0 h1
Structure containing the result of the test.
do set_variable OUTPUT_DIRECTORY $WORKDIR T
std::string getTitle()
Returns a standard string to be used as title of a graphical root object.
do set_variable SIGMA_NS set_variable OUTLIERS set_variable OUTPUT_FILE matrix[${ALPHA_DEG}\deg\] root $JPP JMatrixNZ a $DETECTOR f $INPUT_FILE o $OUTPUT_FILE S
std::string to_string(const T &value)
Convert value to string.