Hi Xin Soh,
Assuming your distribution is not affected by pixel noise (that is, you have traced enough rays), you can use NSDD operand in the merit function to return the standard deviation.
Look for “NSC Operands” page in the help. It contains a table of what values to use for the Pix# and Data#. In you case, Pix# should be -4, and Data# can be 0 or 1.
You can also access this data in ZPL scripts using NSDD function, and in extensions with NSCDetectorData command (use of extensions are deprecated).