> How can we see the estimated irf function for the contrast between two conditions?
Since you've already had the response function for the two conditions, you can use 3dcalc to create the time series of the difference between the two. Another approach would be to set up a symbolic contrast coding (assuming 6 basis functions)
-gltsym 'SYM: +first_condition[[0..5]] -second_condition[[0..5]]'
and then you can extract the time series of the HRF difference.
Gang