Figure 4: Single-photon Wigner function.

The reconstructed density matrix of our measurement is given below. 

dm = {{0.292986, 0.010693 + 0.010368*I,
    0.042302 - 0.016842*I, -0.000181 - 0.027813*I,
    0.004984 - 0.008119*I,
    0.003445 - 0.002869*I, -0.012434 - 0.000518*I, -0.002187 +
     0.000711*I, 0.002833 + 0.003162*I}, {0.010693 - 0.010369*I,
    0.678311, -0.045697 - 0.051679*I, 0.043459 + 0.017763*I,
    0.010088 - 0.006043*I, 0.011128 - 0.025010*I,
    0.007945 + 0.004867*I, -0.015916 + 0.005969*I, -0.005254 +
     0.002235*I}, {0.042302 + 0.016842*I, -0.045697 + 0.051679*I,
    0.014932, -0.003577 - 0.001030*I, 0.000479 - 0.000908*I,
    0.002749 + 0.002735*I, -0.003315 - 0.000340*I,
    0.000532 - 0.001805*I,
    0.000312 + 0.000123*I}, {-0.000181 + 0.027813*I,
    0.043459 - 0.017763*I, -0.003577 + 0.001030*I, 0.007636,
    0.000401 + 0.001061*I, 0.000023 - 0.002966*I,
    0.001328 - 0.000563*I, -0.001325 + 0.000539*I, -0.000451 +
     0.000656*I}, {0.004984 + 0.008119*I, 0.010088 + 0.006043*I,
    0.000479 + 0.000908*I, 0.000401 - 0.001061*I,
    0.001977, -0.000401 + 0.000762*I, -0.000004 -
     0.001003*I, -0.000189 + 0.000253*I, -0.000151 -
     0.000048*I}, {0.003445 + 0.002869*I, 0.011127 + 0.025010*I,
    0.002749 - 0.002735*I,
    0.000023 + 0.002966*I, -0.000401 - 0.000762*I,
    0.002398, -0.000763 + 0.000669*I, -0.000370 -
     0.000815*I, -0.000276 + 0.000014*I}, {-0.012434 + 0.000518*I,
    0.007945 - 0.004867*I, -0.003315 + 0.000340*I,
    0.001328 + 0.000563*I, -0.000004 + 0.001003*I, -0.000763 -
     0.000669*I,
    0.001090, -0.000240 + 0.000250*I, -0.000090 -
     0.000084*I}, {-0.002187 - 0.000711*I, -0.015916 - 0.005969*I,
    0.000532 + 0.001805*I, -0.001325 - 0.000539*I, -0.000189 -
     0.000253*I, -0.000370 + 0.000815*I, -0.000240 - 0.000250*I,
    0.000537,
    0.000093 - 0.000050*I}, {0.002833 - 0.003162*I, -0.005254 -
     0.002235*I,
    0.000312 - 0.000123*I, -0.000451 - 0.000656*I, -0.000151 +
     0.000048*I, -0.000276 - 0.000014*I, -0.000090 + 0.000084*I,
    0.000093 + 0.000050*I, 0.000132}};