At each grid point (τ1, τ3) at a given waiting time Tw, the photodissociation yield of 27 phase-cycling steps were weighted by the appropriate coefficient for rephasing, nonrephasing or absorptive pathways and summed up. After obtaining full 2D time domain data, it was Fourier transformed along τ3 axis first and then τ1 axis to be converted into frequency domain. Zero padding with five times the length of signal (220 zeroes) was employed along both τ1 and τ3 axes to interpolate the data. A Hanning window was used apodize the time-domain data. Since all four pulses were generated with the same pulse shaper and their phases were controlled, phase correction was not required in data processing.