Added selection.cpp. TRying to produce fig1 from AN.