code: def model_counts_map(self, name=None, exclude=None, use_mask=False):
                  
          maps = [c.model_counts_map(name, exclude, use_mask=use_mask)
                  for c in self.components]
          return skymap.coadd_maps(self.geom, maps)
docs: Return the model counts map for a single source, a list of
          sources, or for the sum of all sources in the ROI.  The
          exclude parameter can be used to exclude one or more
          components when generating the model map.
  
          Parameters
          ----------
          name : str or list of str
  
             Parameter controlling the set of sources for which the
             model counts map will be calculated.  If name=None the
             model map will be generated for all sources in the ROI.
  
          exclude : str or list of str
  
             List of sources that will be excluded when calculating the
             model map.
  
          use_mask : bool
  
             Parameter that specifies in the model counts map should include
             mask pixels (i.e., ones whose weights are <= 0)
  
          Returns
          -------
          map : `~gammapy.maps.Map`