For a non-interacting many particle Bose system whose energy operator is diagonal in the number of occupation operatorsn j upper bounds on the thermal averages n j are obtained. These bounds lead to the proof of Bose-Einstein condensation for finite values of the inverse temperature β and for chemical potential µ = 0. Finally for µ < 0, in the case of a generalization of the studied model system, the property of Local Gaussian Domination for the grand canonical partition function is proved.