1 #ifndef __JCOMPAREHISTOGRAMS__JTESTSIGNIFICANCE__
2 #define __JCOMPAREHISTOGRAMS__JTESTSIGNIFICANCE__
17 namespace JCOMPAREHISTOGRAMS {}
18 namespace JPP {
using namespace JCOMPAREHISTOGRAMS; }
20 namespace JCOMPAREHISTOGRAMS {
50 JTest_t(
"Significance",
"Significance")
65 const TH1* h1 =
dynamic_cast<const TH1*
>(o1);
66 const TH1* h2 =
dynamic_cast<const TH1*
>(o2);
68 if (h1 == NULL || h2 == NULL) {
69 THROW(JValueOutOfRange,
"JTestSignificance::test(): Could not cast given TObjects to TH1.");
72 if(h1->GetNbinsX() != h2->GetNbinsX() ||
73 h1->GetNbinsY() != h2->GetNbinsY() ||
74 h1->GetNbinsZ() != h2->GetNbinsZ()) {
75 THROW(JValueOutOfRange,
"JTestSignificance::test(): Histograms with different bining. The objects: " <<
76 h1->GetName() <<
" and " << h2->GetName() <<
" can not be compared." << endl);
80 K = h2->GetEntries() / h1->GetEntries();
83 TH1* h3 = (TH1*) h1->Clone(h1->GetName() == h2->GetName() ?
91 for (
int i=1 ;
i <= h1->GetNbinsX() ; ++
i){
92 for (
int j=1 ;
j <= h1->GetNbinsY() ; ++
j){
93 for (
int k=1 ;
k <= h1->GetNbinsZ() ; ++
k){
95 const double a = h1->GetBinContent(
i,
j,
k);
96 const double b = h2->GetBinContent(
i,
j,
k);
99 S += fabs((a -
K*b)/sqrt(a +
K*
K*b));
105 S /= (h1->GetNbinsX()*h1->GetNbinsY()*h1->GetNbinsZ());
111 h3->SetTitle(title.
getTitle().c_str());
113 const int Ndims = h3->GetDimension();
117 }
else if (Ndims == 2) {
136 std::istream&
read(std::istream&
in)
override
143 THROW(JValueOutOfRange,
"JTestSignificance::read(): Invalid threshold value " <<
threshold);
then fatal No hydrophone data file $HYDROPHONE_TXT fi sort gr k
Interface to read input and write output for TObject tests.
std::string getTitle() const
Returns a standard string to be used as title of a graphical root object.
Class dedicated to standardize the title of the graphical objects produced by the JTest_t() derived c...
#define THROW(JException_t, A)
Marco for throwing exception with std::ostream compatible message.
#define MAKE_CSTRING(A)
Make C-string.
Auxiliary class to handle file name, ROOT directory and object name.
void test(const TObject *o1, const TObject *o2) override
Applies Significance test for two ROOT TH1 histograms.
std::istream & read(std::istream &in) override
Read test parameters from input.
double K
normalization factor between histograms.
#define MAKE_STRING(A)
Make string.
const std::string resultType
test result type
const std::string testName
test name
double threshold
threshold p-value to decide if test is passed.
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
JTestSignificance()
Default constructor.
Structure containing the result of the test.