How can define the subdomain

First of all, please use markdown syntax to highlight your code, see: Problem with bool in solve() - #3 by kamensky

Also note that you should make your code complete, i.e. include all definitions so that people can copy the code and reproduce the behavior.

Finally, you should visualize domains in Paraview for verification, by for instance calling File("domains.pvd")<<domains