Properties of Gutzwiller wave functions for multiband models