Shifts in BCFW method for QED

We study the application of BCFW recursion relations to the QED processes $0\to e^- e^+ n \gamma$. Based on 6-point amplitudes (both MHVA and NMHVA) computed from Feynman diagrams in the Berends-Giele gauge, we conduct a comprehensive study on all different shifts. Then we propose a new shift (LLYZ shift) which can lead to the full amplitudes for these processes and can have some realistic computation advantages. We compare the number of terms and the independent amplitudes of this novel shift with a few typical shifts.