Uncertainty.cxx

Go to the documentation of this file.

//***************************************************************************************
//***************************************************************************************
 
#include "NAGASH/Uncertainty.h"
 
using namespace NAGASH;
 
double Uncertainty::AplusB(double A_err, double B_err)
{
    return sqrt(A_err * A_err + B_err * B_err);
}
 
double Uncertainty::AminusB(double A_err, double B_err)
{
    return sqrt(A_err * A_err + B_err * B_err);
}
 
double Uncertainty::AoverAplusB(double A, double AplusB, double A_err, double AplusB_err)
{
    double B = AplusB - A;
    double B_err = A_err < AplusB_err ? sqrt(pow(AplusB_err, 2) - pow(A_err, 2)) : 0;
 
    if (AplusB != 0)
        return sqrt(B * B * A_err * A_err + A * A * B_err * B_err) / pow(A + B, 2);
    else
        return 0;
}
 
double Uncertainty::AoverB(double A, double B, double A_err, double B_err)
{
    if (B != 0)
        return sqrt(B * B * A_err * A_err + A * A * B_err * B_err) / pow(B, 2);
    else
        return 0;
}
 
double Uncertainty::AminusBoverAplusB(double A, double B, double A_err, double B_err)
{
    if (A != 0 || B != 0)
        return 2 * sqrt(B * B * A_err * A_err + A * A * B_err * B_err) / pow(A + B, 2);
    else
        return 0;
}