  We present first results using the stochastic Laplace-Heaviside (sLapH)
method in the twisted mass formulation. The calculations are performed on gauge
configurations provided by the European Twisted Mass (ETM) collaboration with
2+1+1 dynamical quark flavours at a single value of the lattice spacing. In
particular, we compute disconnected contributions to flavour singlet
pseudo-scalar mesons and compare sLapH to standard volume noise methods.
