def METHOD_NAME(self, station_id, n=100):
    """
    gets set of n templates to allow for the calculation of average templates
    loops first over different viewing angles
    and then over azimuth angles of 0, 22.5 and 45 degree
    and then over zenith angles of 100, 120 and 140 degree
    """
    if station_id not in self.__ref_nu_templates.keys():
        self.__load_nu_template(station_id)
    if self.__nu_template_set == {}:
        self.logger.info("Getting set of templates for station ID {}".format(station_id))
        n_tmpl = 0
        zen_refs = np.deg2rad([100, 120, 140])
        az_refs = np.deg2rad([0, 22.5, 45])
        dCherenkovs = np.deg2rad([0, -0.5, -1, -1.5, -2, -3, -4, -5])
        for zen in zen_refs:
            for az in az_refs:
                for dCh in dCherenkovs:
                    self.__nu_template_set[n_tmpl] = self.__nu_templates[station_id][zen][az][dCh]
                    n_tmpl += 1
                    if n_tmpl >= n:
                        break
    return self.__nu_template_set