I am not sure that I fully understand why you have the two goals in this case.
If you want to constrain the minimization to within your region of interest, then you can change the upper and lower bounds for the design variable in the inputs to simulannealbnd. If you want the objective function to be 0 outside the region of interest, then you can encode that in the objective function by testing if the x at that step is within the region. If it is not, you can return 0 or some factor of the distance from the x to the region.
Another option for global optimization for multiple object functions is gamultiobj which finds the Pareto front of multiple fitness functions using genetic algorithm. The documentation for that function is here.