(March 2016, J. Skowron)
INPUT: (mag, error) AND COEFFICIENTS FOR THE GIVEN FIELD.CHIP: (GAMMA, EPSILON, UNIT, BACKG, MAG_01)
if error > 0.003:
new_error = sqrt( (GAMMA*error)^2 + EPSILON^2 )
else:
F = 10^(0.4*(UNIT - mag))
B = 10^(0.4*(UNIT - BACKG))
FLUX_NOISE = sqrt(F + B)
MAG_NOISE = 2.5/ln(10) * FLUX_NOISE/F
new_error = max( error, sqrt( (GAMMA*MAG_NOISE)^2 + EPSILON^2 ) , 0.01 * 10^( 0.6*(MAG_01 - mag) ) )
end
RETURN (mag, new_error)