sum( W_iRQ_iRU_i )/(Nsum( W_i ))
RQ_I is the ith residual from the Q map, and RU_I is the ith residual from the U map. A residual is the bolometer Q or U sample minus the Q or U map pixel value (the pixel value is the weighted mean of the Q or U sample values that fall in the pixel). These residuals can be dumped by makemap when the Q and U maps are created.
W_i is the weight for the ith sample. It is the geometric mean of the weights associated with RQ_i and RU_i. These can also be dumped by makemap, but are created initially by calcuq as the residual of the fitting process for each sample.
The input NDFs should be created by running makemap twice on the CALCQU output for a single
POL2 observation - once to create a Q map and once to create a U map. On each invocation of
makemap, the config should include
" . This will causes two extra NDFs
to be created by each invocation of makemap with suffixes
" . These should be supplied as input to this command.