Hi,

I am using the honeycomb lattice and boson hubbard model from the alps library on vispython. I wish to measure the staggered structure factor given by the equation

                                               S = 1/N * Summation [ e(i)e(j) <n(i)n(j)> ]

where e(i) = 1 or -1 for i on sublattice A(B) in terms of density-density correlations. I know that alps measure the density structure factor using wave vectors. How can i make alps measure the structure factor given above ?

Regards,

Mike