rubin-sim 2.6.1.dev2__py3-none-any.whl

rubin_sim/maf/batches/col_map_dict.py

@@ -0,0 +1,274 @@
+__all__ = ("col_map_dict",)
+
+
+def col_map_dict(dict_name=None):
+    if dict_name is None:
+        dict_name = "FBS"
+    dict_name = dict_name.lower()
+
+    if dict_name == "fbs":
+        col_map = {}
+        col_map["ra"] = "fieldRA"
+        col_map["dec"] = "fieldDec"
+        col_map["raDecDeg"] = True
+        col_map["mjd"] = "observationStartMJD"
+        col_map["exptime"] = "visitExposureTime"
+        col_map["visittime"] = "visitTime"
+        col_map["alt"] = "altitude"
+        col_map["az"] = "azimuth"
+        col_map["lst"] = "observationStartLST"
+        col_map["filter"] = "band"
+        col_map["band"] = "band"
+        col_map["fiveSigmaDepth"] = "fiveSigmaDepth"
+        col_map["night"] = "night"
+        col_map["slewtime"] = "slewTime"
+        col_map["slewdist"] = "slewDistance"
+        col_map["seeingEff"] = "seeingFwhmEff"
+        col_map["seeingGeom"] = "seeingFwhmGeom"
+        col_map["skyBrightness"] = "skyBrightness"
+        col_map["moonDistance"] = "moonDistance"
+        col_map["slewactivities"] = {}
+        col_map["metadataList"] = [
+            "airmass",
+            "normairmass",
+            "seeingEff",
+            "skyBrightness",
+            "fiveSigmaDepth",
+            "HA",
+            "moonDistance",
+            "solarElong",
+            "saturation_mag",
+        ]
+        col_map["metadataAngleList"] = ["rotSkyPos"]
+        col_map["scheduler_note"] = "scheduler_note"
+        col_map["scheduler_note_root"] = "scheduler_note_root"
+
+    elif dict_name == "consdb":
+        col_map = {}
+        col_map["ra"] = "s_ra"
+        col_map["dec"] = "s_dec"
+        col_map["raDecDeg"] = True
+        col_map["mjd"] = "obs_start_mjd"
+        col_map["exptime"] = "exp_time"
+        col_map["visittime"] = "dark_time"
+        col_map["alt"] = "altitude"
+        col_map["az"] = "azimuth"
+        col_map["filter"] = "band"
+        col_map["band"] = "band"
+        col_map["fiveSigmaDepth"] = "stats_mag_lim_median"
+        col_map["night"] = "day_obs"
+        col_map["slewtime"] = "slew_time"
+        col_map["slewdist"] = "slew_distance"
+        col_map["seeingEff"] = "fwhm_eff"
+        col_map["seeingGeom"] = "fwhm_geom"
+        col_map["skyBrightness"] = "sky_bg_median_mag"
+        col_map["moonDistance"] = "moon_distance"
+        col_map["slewactivities"] = {}
+        col_map["metadataList"] = [
+            "airmass",
+            "normairmass",
+            "fwhm_eff",
+            "sky_bg_median_mag",
+            "stats_mag_lim_median",
+            "HA",
+            "moon_distance",
+            "solarElong",
+            "saturation_mag",
+        ]
+        col_map["metadataAngleList"] = ["sky_rotation"]
+        col_map["scheduler_note"] = "scheduler_note"
+        col_map["scheduler_note_root"] = "scheduler_note_root"
+
+    elif dict_name == "fbs_sim" or dict_name == "opsimfbs":
+        col_map = {}
+        col_map["ra"] = "fieldRA"
+        col_map["dec"] = "fieldDec"
+        col_map["raDecDeg"] = True
+        col_map["mjd"] = "observationStartMJD"
+        col_map["exptime"] = "visitExposureTime"
+        col_map["visittime"] = "visitTime"
+        col_map["alt"] = "altitude"
+        col_map["az"] = "azimuth"
+        col_map["lst"] = "observationStartLST"
+        col_map["filter"] = "filter"
+        col_map["fiveSigmaDepth"] = "fiveSigmaDepth"
+        col_map["night"] = "night"
+        col_map["slewtime"] = "slewTime"
+        col_map["slewdist"] = "slewDistance"
+        col_map["seeingEff"] = "seeingFwhmEff"
+        col_map["seeingGeom"] = "seeingFwhmGeom"
+        col_map["skyBrightness"] = "skyBrightness"
+        col_map["moonDistance"] = "moonDistance"
+        col_map["slewactivities"] = {}
+        col_map["metadataList"] = [
+            "airmass",
+            "normairmass",
+            "seeingEff",
+            "skyBrightness",
+            "fiveSigmaDepth",
+            "HA",
+            "moonDistance",
+            "solarElong",
+            "saturation_mag",
+        ]
+        col_map["metadataAngleList"] = ["rotSkyPos"]
+        col_map["scheduler_note"] = "scheduler_note"
+        col_map["scheduler_note_root"] = "scheduler_note_root"
+
+    elif dict_name == "opsimv4":
+        col_map = {}
+        col_map["ra"] = "fieldRA"
+        col_map["dec"] = "fieldDec"
+        col_map["raDecDeg"] = True
+        col_map["mjd"] = "observationStartMJD"
+        col_map["exptime"] = "visitExposureTime"
+        col_map["visittime"] = "visitTime"
+        col_map["alt"] = "altitude"
+        col_map["az"] = "azimuth"
+        col_map["lst"] = "observationStartLST"
+        col_map["filter"] = "filter"
+        col_map["fiveSigmaDepth"] = "fiveSigmaDepth"
+        col_map["night"] = "night"
+        col_map["slewtime"] = "slewTime"
+        col_map["slewdist"] = "slewDistance"
+        col_map["seeingEff"] = "seeingFwhmEff"
+        col_map["seeingGeom"] = "seeingFwhmGeom"
+        col_map["skyBrightness"] = "skyBrightness"
+        col_map["moonDistance"] = "moonDistance"
+        col_map["fieldId"] = "fieldId"
+        col_map["proposalId"] = "proposalId"
+        # slew speeds table
+        col_map["slewSpeedsTable"] = "SlewMaxSpeeds"
+        # slew states table
+        col_map["slewStatesTable"] = "SlewFinalState"
+        # slew activities list
+        col_map["slewActivitiesTable"] = "SlewActivities"
+        # Slew columns
+        col_map["Dome Alt Speed"] = "domeAltSpeed"
+        col_map["Dome Az Speed"] = "domeAzSpeed"
+        col_map["Tel Alt Speed"] = "telAltSpeed"
+        col_map["Tel Az Speed"] = "telAzSpeed"
+        col_map["Rotator Speed"] = "rotatorSpeed"
+        col_map["Tel Alt"] = "telAlt"
+        col_map["Tel Az"] = "telAz"
+        col_map["Rot Tel Pos"] = "rotTelPos"
+        col_map["Dome Alt"] = "domeAlt"
+        col_map["Dome Az"] = "domeAz"
+        col_map["slewactivities"] = {
+            "Dome Alt": "domalt",
+            "Dome Az": "domaz",
+            "Dome Settle": "domazsettle",
+            "Tel Alt": "telalt",
+            "Tel Az": "telaz",
+            "Tel Rot": "telrot",
+            "Tel Settle": "telsettle",
+            "TelOptics CL": "telopticsclosedloop",
+            "TelOptics OL": "telopticsopenloop",
+            "Readout": "readout",
+            "Filter": "filter",
+        }
+        col_map["metadataList"] = [
+            "airmass",
+            "normairmass",
+            "seeingEff",
+            "skyBrightness",
+            "fiveSigmaDepth",
+            "HA",
+            "moonDistance",
+            "solarElong",
+        ]
+        col_map["metadataAngleList"] = ["rotSkyPos"]
+
+    elif dict_name == "opsimv3":
+        col_map = {}
+        col_map["ra"] = "fieldRA"
+        col_map["dec"] = "fieldDec"
+        col_map["raDecDeg"] = False
+        col_map["mjd"] = "expMJD"
+        col_map["exptime"] = "visitExpTime"
+        col_map["visittime"] = "visitTime"
+        col_map["alt"] = "altitude"
+        col_map["az"] = "azimuth"
+        col_map["lst"] = "lst"
+        col_map["filter"] = "filter"
+        col_map["fiveSigmaDepth"] = "fiveSigmaDepth"
+        col_map["night"] = "night"
+        col_map["slewtime"] = "slewTime"
+        col_map["slewdist"] = "slewDist"
+        col_map["seeingEff"] = "FWHMeff"
+        col_map["seeingGeom"] = "FWHMgeom"
+        col_map["skyBrightness"] = "filtSkyBrightness"
+        col_map["moonDistance"] = "dist2Moon"
+        col_map["fieldId"] = "fieldID"
+        col_map["proposalId"] = "propID"
+        # slew speeds table
+        col_map["slewSpeedsTable"] = "SlewMaxSpeeds"
+        # slew states table
+        col_map["slewStatesTable"] = "SlewStates"
+        # Slew activities list
+        col_map["slewActivitiesTable"] = "SlewActivities"
+        col_map["Dome Alt Speed"] = "domeAltSpeed"
+        col_map["Dome Az Speed"] = "domeAzSpeed"
+        col_map["Tel Alt Speed"] = "telAltSpeed"
+        col_map["Tel Az Speed"] = "telAzSpeed"
+        col_map["Rotator Speed"] = "rotatorSpeed"
+        col_map["Tel Alt"] = "telAlt"
+        col_map["Tel Az"] = "telAz"
+        col_map["Rot Tel Pos"] = "rotTelPos"
+        col_map["Dome Alt"] = "domAlt"
+        col_map["Dome Az"] = "domAz"
+        col_map["slewactivities"] = {
+            "Dome Alt": "DomAlt",
+            "Dome Az": "DomAz",
+            "Tel Alt": "TelAlt",
+            "Tel Az": "TelAz",
+            "Tel Rot": "Rotator",
+            "Settle": "Settle",
+            "TelOptics CL": "TelOpticsCL",
+            "TelOptics OL": "TelOpticsOL",
+            "Readout": "Readout",
+            "Filter": "Filter",
+        }
+        col_map["metadataList"] = [
+            "airmass",
+            "normairmass",
+            "seeingEff",
+            "skyBrightness",
+            "fiveSigmaDepth",
+            "HA",
+            "moonDistance",
+            "solarElong",
+        ]
+        col_map["metadataAngleList"] = ["rotSkyPos"]
+
+    elif dict_name == "barebones":
+        col_map = {}
+        col_map["ra"] = "ra"
+        col_map["dec"] = "dec"
+        col_map["raDecDeg"] = True
+        col_map["mjd"] = "mjd"
+        col_map["exptime"] = "exptime"
+        col_map["visittime"] = "exptime"
+        col_map["alt"] = "alt"
+        col_map["az"] = "az"
+        col_map["filter"] = "filter"
+        col_map["fiveSigmaDepth"] = "fivesigmadepth"
+        col_map["night"] = "night"
+        col_map["slewtime"] = "slewtime"
+        col_map["slewdist"] = None
+        col_map["seeingGeom"] = "seeing"
+        col_map["seeingEff"] = "seeing"
+        col_map["metadataList"] = [
+            "airmass",
+            "normairmass",
+            "seeingEff",
+            "skyBrightness",
+            "fiveSigmaDepth",
+            "HA",
+        ]
+        col_map["metadataAngleList"] = ["rotSkyPos"]
+
+    else:
+        raise ValueError(f"No built in column dict with name {dict_name}")
+
+    return col_map

rubin_sim/maf/batches/common.py

@@ -0,0 +1,388 @@
+__all__ = (
+    "combine_info_labels",
+    "filter_list",
+    "standard_summary",
+    "extended_summary",
+    "standard_metrics",
+    "extended_metrics",
+    "lightcurve_summary",
+    "standard_angle_metrics",
+    "summary_completeness_at_time",
+    "summary_completeness_over_h",
+    "fraction_population_at_threshold",
+    "microlensing_summary",
+)
+
+
+import rubin_sim.maf.metrics as metrics
+
+
+def combine_info_labels(info1, info2):
+    if info1 is not None and info2 is not None:
+        info = info1 + " " + info2
+    elif info1 is not None:
+        info = info1
+    elif info2 is not None:
+        info = info2
+    else:
+        info = None
+    return info
+
+
+def filter_list(all=True, extra_sql=None, extra_info_label=None, band_col="band"):
+    """Return a list of filters, plot colors and orders.
+
+    Parameters
+    ----------
+    all : `bool`, optional
+        Include 'all' in the list of filters and as part of the
+        colors/orders dictionaries.
+    extra_sql : `str`, optional
+        Additional sql constraint to add to constraints returned per filter.
+    extra_info_label : `str`, optional
+        Substitute info_label to add to info_label strings composed per band.
+    band_col : `str`, optional
+        The column name to use for the band (or filter).
+
+    Returns
+    -------
+    filterlist : `list` [`str]
+    colors : `dict` {`str`: `str`}
+    orders : `dict` {`str`: int}
+    sqls : `dict` {`str`: `str`}
+    info_labels : `dict` {`str`: `str}
+    """
+    if all:
+        bandlist = ("all", "u", "g", "r", "i", "z", "y")
+    else:
+        bandlist = ("u", "g", "r", "i", "z", "y")
+    colors = {"u": "cyan", "g": "g", "r": "orange", "i": "r", "z": "m", "y": "b"}
+    orders = {"u": 1, "g": 2, "r": 3, "i": 4, "z": 5, "y": 6}
+    if all:
+        colors["all"] = "k"
+        orders["all"] = 0
+    sqls = {}
+    info_labels = {}
+    if extra_info_label is None:
+        if extra_sql is None or len(extra_sql) == 0:
+            md = ""
+        else:
+            md = "%s " % extra_sql
+    else:
+        md = "%s " % extra_info_label
+    for b in bandlist:
+        if b == "all":
+            sqls[b] = ""
+            info_labels[b] = md + "all bands"
+        else:
+            sqls[b] = f"{band_col} = '{b}'"
+            info_labels[b] = md + f"{b} band"
+    if extra_sql is not None and len(extra_sql) > 0:
+        for s in sqls:
+            if s == "all":
+                sqls[s] = extra_sql
+            else:
+                sqls[s] = f"({extra_sql}) and ({sqls[s]})"
+    return bandlist, colors, orders, sqls, info_labels
+
+
+def standard_summary(with_count=True):
+    """A set of standard summary metrics, to calculate
+    Mean, RMS, Median, #, Max/Min, and # 3-sigma outliers.
+
+    Parameters
+    ----------
+    with_count : `bool`, optional
+        Include the "Count" metric in the set of summary metrics or not.
+
+    Returns
+    -------
+    standardSummary : `list` [`maf.BaseMetric`]
+        List of metrics appropriate to use to summarize the results from
+        another metric.
+    """
+    standardSummary = [
+        metrics.MeanMetric(),
+        metrics.RmsMetric(),
+        metrics.MedianMetric(),
+        metrics.MaxMetric(),
+        metrics.MinMetric(),
+        metrics.NoutliersNsigmaMetric(metric_name="N(+3Sigma)", n_sigma=3),
+        metrics.NoutliersNsigmaMetric(metric_name="N(-3Sigma)", n_sigma=-3.0),
+    ]
+    if with_count:
+        standardSummary += [metrics.CountMetric()]
+    return standardSummary
+
+
+def extended_summary():
+    """An extended set of summary metrics, to calculate all that is in
+    the standard summary stats, plus 25/75 percentiles.
+
+    Returns
+    --------
+    extendedSummary : `list` [`maf.BaseMetric`]
+        List of metrics appropriate to use to summarize the results
+        from another metric.
+    """
+    extendedStats = standard_summary()
+    extendedStats += [
+        metrics.PercentileMetric(metric_name="25th%ile", percentile=25),
+        metrics.PercentileMetric(metric_name="75th%ile", percentile=75),
+    ]
+    return extendedStats
+
+
+def lightcurve_summary():
+    lightcurveSummary = [
+        metrics.SumMetric(metric_name="Total detected"),
+        metrics.CountMetric(metric_name="Total lightcurves in footprint"),
+        metrics.CountMetric(metric_name="Total lightcurves on sky", mask_val=0),
+        metrics.MeanMetric(metric_name="Fraction detected in footprint (mean)"),
+        metrics.MeanMetric(mask_val=0, metric_name="Fraction detected of total (mean)"),
+    ]
+    return lightcurveSummary
+
+
+def standard_metrics(colname, replace_colname=None):
+    """A set of standard simple metrics for some quantity.
+    Typically would be applied with unislicer.
+
+    Parameters
+    ----------
+    colname : `str`
+        The column name to apply the metrics to.
+    replace_colname: `str` or None, optional
+        Value to replace colname with in the metric_name.
+        i.e. if replace_colname='' then metric name is Mean,
+        instead of Mean Airmass, or
+        if replace_colname='seeingGeom', then metric name is
+        Mean seeingGeom instead of Mean seeingFwhmGeom.
+
+    Returns
+    -------
+    standardMetrics : `list` [`maf.BaseMetric`]
+        List of appropriate MAF metrics to evaluate a distribution.
+    """
+    standardMetrics = [
+        metrics.MeanMetric(colname),
+        metrics.MedianMetric(colname),
+        metrics.MinMetric(colname),
+        metrics.MaxMetric(colname),
+    ]
+    if replace_colname is not None:
+        for m in standardMetrics:
+            if len(replace_colname) > 0:
+                m.name = m.name.replace("%s" % colname, "%s" % replace_colname)
+            else:
+                m.name = m.name.rstrip(" %s" % colname)
+    return standardMetrics
+
+
+def extended_metrics(colname, replace_colname=None):
+    """An extended set of simple metrics for some quantity.
+    Typically applied with unislicer.
+
+    Parameters
+    ----------
+    colname : `str`
+        The column name to apply the metrics to.
+    replace_colname: `str` or None, optional
+        Value to replace colname with in the metric_name.
+        i.e. if replace_colname='' then metric name is Mean,
+        instead of Mean Airmass, or
+        if replace_colname='seeingGeom', then metric name is
+        Mean seeingGeom instead of Mean seeingFwhmGeom.
+
+    Returns
+    -------
+    extendedMetrics : `list` [`maf.BaseMetric`]
+        List of appropriate MAF metrics to evaluate a distribution.
+    """
+    extendedMetrics = standard_metrics(colname, replace_colname=None)
+    extendedMetrics += [
+        metrics.RmsMetric(colname),
+        metrics.NoutliersNsigmaMetric(colname, metric_name="N(+3Sigma) " + colname, n_sigma=3),
+        metrics.NoutliersNsigmaMetric(colname, metric_name="N(-3Sigma) " + colname, n_sigma=-3),
+        metrics.PercentileMetric(colname, percentile=25),
+        metrics.PercentileMetric(colname, percentile=75),
+        metrics.CountMetric(colname),
+    ]
+    if replace_colname is not None:
+        for m in extendedMetrics:
+            if len(replace_colname) > 0:
+                m.name = m.name.replace("%s" % colname, "%s" % replace_colname)
+            else:
+                m.name = m.name.rstrip(" %s" % colname)
+    return extendedMetrics
+
+
+def standard_angle_metrics(colname, replace_colname=None):
+    """A set of standard simple metrics for a wrap-around angle quantity.
+
+    Parameters
+    ----------
+    colname : `str`
+        The column name to apply the metrics to.
+    replace_colname: `str` or None, optional
+        Value to replace colname with in the metric_name.
+        i.e. if replace_colname='' then metric name is Mean,
+        instead of Mean Airmass, or
+        if replace_colname='seeingGeom', then metric name is
+        Mean seeingGeom instead of Mean seeingFwhmGeom.
+
+    Returns
+    -------
+    standardAngleMetrics : `list` [`maf.BaseMetric`]
+        List of appropriate MAF metrics for angle distributions.
+    """
+    standardAngleMetrics = [
+        metrics.MeanAngleMetric(colname),
+        metrics.RmsAngleMetric(colname),
+        metrics.FullRangeAngleMetric(colname),
+        metrics.MinMetric(colname),
+        metrics.MaxMetric(colname),
+    ]
+    if replace_colname is not None:
+        for m in standardAngleMetrics:
+            if len(replace_colname) > 0:
+                m.name = m.name.replace("%s" % colname, "%s" % replace_colname)
+            else:
+                m.name = m.name.rstrip(" %s" % colname)
+    return standardAngleMetrics
+
+
+def summary_completeness_at_time(times, h_val, h_index=0.33):
+    """A simple list of summary metrics to be applied to the Discovery_Time
+    or PreviouslyKnown metrics.
+    (can be used with any moving object metric which returns the
+    time of discovery).
+
+    Parameters
+    ----------
+    times : `np.ndarray `or list` [`float`]
+        The times at which to evaluate the completeness @ Hval.
+    h_val : `float`
+        The H value at which to evaluate the completeness
+        (cumulative and differential).
+    h_index : `float`, optional
+        The index of the power law to integrate H over
+        (for cumulative completeness).
+
+    Returns
+    -------
+    summaryMetrics : `list` [`maf.MoCompletenessAtTimeMetric`]
+        List of completeness metrics to be evaluated at the specified times.
+    """
+    summaryMetrics = [
+        metrics.MoCompletenessAtTimeMetric(times=times, hval=h_val, hindex=h_index, cumulative=False),
+        metrics.MoCompletenessAtTimeMetric(times=times, hval=h_val, hindex=h_index, cumulative=True),
+    ]
+    return summaryMetrics
+
+
+def summary_completeness_over_h(requiredChances=1, Hindex=0.33):
+    """A simple list of summary metrics to be applied to the
+    Discovery_N_Chances metric.
+
+    Parameters
+    ----------
+    requiredChances : `int`, optional
+        Number of discovery opportunities required to consider an
+        object 'discovered'.
+    Hindex : `float`, optional
+        The index of the power law to integrate H over
+        (for cumulative completeness).
+
+    Returns
+    -------
+    summaryMetrics : `list` [`maf.MoCompletenessMetric`]
+        List of moving object MoCompleteness metrics
+        (cumulative and differential)
+    """
+    summaryMetrics = [
+        metrics.MoCompletenessMetric(threshold=requiredChances, cumulative=False, hindex=Hindex),
+        metrics.MoCompletenessMetric(threshold=requiredChances, cumulative=True, hindex=Hindex),
+    ]
+    return summaryMetrics
+
+
+def fraction_population_at_threshold(thresholds, optnames=None):
+    """Creates a list of summary metrics to be applied to any moving object
+    metric which reports a float value, calculating the fraction of the
+    population above X.
+
+    Parameters
+    ----------
+    thresholds : `list` [`float`]
+        The thresholds at which to calculate what fraction of the population
+        exceeds these values.
+    optnames : `list` [`str`], optional
+        If provided, these names will be used instead of the threshold values
+        when constructing the metric names.
+        This allows more descriptive summary statistic names.
+
+    Returns
+    -------
+    fracMetrics : `list` [`maf.MoCompletenessMetric`]
+        List of moving object MoCompleteness metrics
+        (differential fractions of the population).
+    """
+    fracMetrics = []
+    for i, threshold in enumerate(thresholds):
+        if optnames is not None:
+            o = optnames[i]
+        else:
+            o = threshold
+        m = metrics.MoCompletenessMetric(
+            threshold=threshold, cumulative=False, metric_name=f"FractionPop {o}"
+        )
+        fracMetrics.append(m)
+    return fracMetrics
+
+
+def microlensing_summary(metric_type, npts_required=10, Fisher_sigmatE_tE_cutoff=0.1):
+    """Calculate summary metrics for the microlensing population metrics.
+
+    Parameters
+    -----------
+    metric_type : `str`
+        Identify whether the metric is "Npts" or Fisher"
+    npts_required : `int`, optional
+        Count the fraction of microlensing events with more than npts_required
+        observations
+    Fisher_sigmatE_tE_cutoff : `float`, optional
+        Count the fraction of microlensing events with characterization
+        uncertainty less than this
+
+    Returns
+    -------
+    microlensingSummary : `list` [`maf.BaseMetric`]
+        List of appropriate MAF metrics for this type of microlensing
+        metric with the specified threshold values.
+    """
+    if metric_type != "Npts" and metric_type != "Fisher":
+        raise Exception('metric_type must be "Npts" or "Fisher"')
+    if metric_type == "Npts":
+        microlensingSummary = [
+            metrics.FracAboveMetric(
+                cutoff=npts_required,
+                metric_name=f"Fraction w/ at least {npts_required} points",
+            ),
+            metrics.CountMetric(metric_name="Total lightcurves in footprint"),
+            metrics.CountMetric(metric_name="Total lightcurves on sky", mask_val=0),
+            metrics.MeanMetric(metric_name="Mean number of points per lightcurves in footprint"),
+            metrics.MeanMetric(mask_val=0, metric_name="Mean number of points per lightcurves in total"),
+        ]
+    elif metric_type == "Fisher":
+        microlensingSummary = [
+            metrics.FracBelowMetric(
+                cutoff=Fisher_sigmatE_tE_cutoff,
+                metric_name=f"Fraction w/ sigma_tE/tE < {Fisher_sigmatE_tE_cutoff}",
+            ),
+            metrics.CountMetric(metric_name="Total lightcurves in footprint"),
+            metrics.CountMetric(metric_name="Total lightcurves on sky", mask_val=0),
+            metrics.RealMeanMetric(metric_name="Mean sigma_tE/tE in footprint (mean)"),
+            metrics.RealMeanMetric(mask_val=0, metric_name="Mean sigma_tE/tE of total (mean)"),
+        ]
+    return microlensingSummary