detectorfluxanalysis2d_tests.cpp

Go to the documentation of this file.

/**
 * @file detectorfluxanalysis2d_tests.cpp
 * @brief Implements regression and invariant tests for detectorfluxanalysis2d tests.
 * @ingroup qpe_tests
 *
 * Source: @c tests/detectorfluxanalysis2d_tests.cpp. This file is part of Quantum Phenomena Explorer (QPE),
 * maintained by Officina Turini. See LICENSE_IT.txt / LICENSE.txt for the
 * QPE freeware source-available terms and THIRD_PARTY_NOTICES.txt for
 * third-party components.
 */
 
#include "detectorfluxanalysis2d_tests.h"
 
#include "core/quantum/detectorfluxanalysis2d.h"
 
#include <cmath>
#include <vector>
 
namespace qpe { namespace tests {
 
/**
 * @brief Runs detector flux analysis2 d tests.
 * @param error Output diagnostic text on failure.
 */
bool runDetectorFluxAnalysis2DTests(std::string &error)
{
    std::vector<double> profile(9u, 0.0);
    profile[3] = 1.0;
    profile[5] = 3.0;
 
    const std::vector<double> rebinned =
        DetectorFluxAnalysis2D::rebinNormalizedProfile(profile, 9u);
    if (rebinned.size() != 9u ||
        std::fabs(rebinned[3] - 0.25) > 1.0e-12 ||
        std::fabs(rebinned[5] - 0.75) > 1.0e-12) {
        error = "detector flux profile rebin/normalization failed";
        return false;
    }
 
    std::vector<PositionMeasurement2DPoint> hits;
    for (int i = 0; i < 25; ++i) {
        PositionMeasurement2DPoint p;
        p.x = 8.0f;
        p.y = -1.0f; // index 3 for dy=1 on 9 samples
        hits.push_back(p);
    }
    for (int i = 0; i < 75; ++i) {
        PositionMeasurement2DPoint p;
        p.x = 8.0f;
        p.y = 1.0f; // index 5
        hits.push_back(p);
    }
 
    const DetectorFluxAnalysis2DResult match =
        DetectorFluxAnalysis2D::analyze(profile, 1.0f, hits, 9u);
    if (!match.expectedValid || !match.observedValid ||
        std::fabs(match.totalVariationDistance) > 1.0e-12) {
        error = "TVD should be zero for a matching arrival ensemble";
        return false;
    }
 
    hits.clear();
    for (int i = 0; i < 100; ++i) {
        PositionMeasurement2DPoint p;
        p.x = 8.0f;
        p.y = -1.0f;
        hits.push_back(p);
    }
    const DetectorFluxAnalysis2DResult mismatch =
        DetectorFluxAnalysis2D::analyze(profile, 1.0f, hits, 9u);
    if (!(mismatch.totalVariationDistance > 0.70 &&
          mismatch.totalVariationDistance <= 1.0)) {
        error = "TVD did not detect a strongly mismatched arrival ensemble";
        return false;
    }
 
    error.clear();
    return true;
}
 
} }