code: def create_avg_edisp(skydir, ltc, event_class, event_types, erec,
                       egy, cth_bins, npts=None):
      
      return create_avg_rsp(create_edisp, skydir, ltc,
                            event_class, event_types,
                            erec, egy,  cth_bins, npts)
docs: Generate model for exposure-weighted DRM averaged over incidence
      angle.
  
      Parameters
      ----------
      egy : `~numpy.ndarray`
          True energies in MeV.
  
      cth_bins : `~numpy.ndarray`
          Bin edges in cosine of the incidence angle.