  Thermal states of light are widely used in quantum optics for various quantum
phenomena testing. Particularly, they can be utilized for characterization of
photon creation and photon annihilation operations. During the last decade the
problem of photon subtraction from multimode quantum states become of much
significance. Therefore, in this work we present a technique for statistical
parameter estimation of multimode multiphoton subtracted thermal states of
light, which can be used for multimode photon annihilation test.
