Frequently asked questions¶
This page contains answers to common questions about the data access and pipeline middleware, such as the Butler and PipelineTask classes.
The lsst.daf.butler package documentation includes a number of overview documentation pages (especially Organizing and identifying datasets) that provide an introduction to many of the concepts referenced here.
When should I use each of the query methods/commands?¶
The Butler Python class and butler command-line tool support five major query operations that can be used to inspect a data repository:
Command-line tool |
Python API |
|---|---|
|
|
|
|
|
|
|
|
|
These operations share many optional arguments that constrain what is returned, but their return types each reflect a different aspect of how datasets are organized).
Collections¶
From the command-line, butler query-collections generally provides the best high-level view of the contents of a data repository, with the --chains=tree format.
For all but the smallest repos, it’s best to start with some kind of guess at what you’re looking for, or the results will still be overwhelming large.
$ butler query-collections /repo/main HSC/runs/RC2/* --chains=tree
Name Type
-------------------------------------- -----------
HSC/runs/RC2/w_2021_02/DM-28282/sfm RUN
HSC/runs/RC2/w_2021_02/DM-28282/rest RUN
HSC/runs/RC2/w_2021_06/DM-28654 CHAINED
HSC/runs/RC2/w_2021_06/DM-28654/rest RUN
HSC/runs/RC2/w_2021_06/DM-28654/sfm RUN
HSC/raw/RC2/9615 TAGGED
HSC/raw/RC2/9697 TAGGED
HSC/raw/RC2/9813 TAGGED
HSC/calib/gen2/20180117 CALIBRATION
HSC/calib/DM-28636 CALIBRATION
HSC/calib/gen2/20180117/unbounded RUN
HSC/calib/DM-28636/unbounded RUN
HSC/masks/s18a RUN
skymaps RUN
refcats/DM-28636 RUN
HSC/runs/RC2/w_2021_02/DM-28282 CHAINED
HSC/runs/RC2/w_2021_02/DM-28282/rest RUN
HSC/runs/RC2/w_2021_02/DM-28282/sfm RUN
HSC/raw/RC2/9615 TAGGED
HSC/raw/RC2/9697 TAGGED
HSC/raw/RC2/9813 TAGGED
HSC/calib/gen2/20180117 CALIBRATION
HSC/calib/DM-28636 CALIBRATION
HSC/calib/gen2/20180117/unbounded RUN
HSC/calib/DM-28636/unbounded RUN
HSC/masks/s18a RUN
skymaps RUN
refcats/DM-28636 RUN
HSC/runs/RC2/w_2021_06/DM-28654/sfm RUN
HSC/runs/RC2/w_2021_06/DM-28654/rest RUN
Note that some collections appear multiple times here - once as a top-level collection, and again later as some child of a CHAINED collection (that’s what the indentation means here).
In the future we may be able to remove some of this duplication.
Similar functionality can be accessed from Python using Butler.collections.query_info.
Dataset Types¶
Registry.queryDatasetTypes reports the dataset types that have been registered with a data repository, even if there aren’t any datasets of that type actually present.
That makes it less useful for exploring a data repository generically, but it’s an important tool when you know the name of the dataset type already and want to see how it’s defined.
Dimension Records¶
Butler.query_dimension_records is the best way to inspect the metadata records associated with data ID keys (“dimensions”).
Those metadata tables include observations (the exposure and visit dimensions), instruments (instrument, physical_filter, detector), and regions on the sky (skymap, tract, patch, htm7).
That isn’t an exhaustive list of dimension tables (actually pseudo-tables in some cases), but you can get one in Python with:
>>> print(butler.dimensions.getStaticDimensions())
And while query_dimension_records shows you the schema of those tables with each record it returns, you can also get it without querying for any data with (e.g.)
>>> print(butler.dimensions["exposure"].schema)
exposure:
instrument: string
id: int
physical_filter: string
obs_id: string
exposure_time: float
Duration of the exposure with shutter open (seconds).
dark_time: float
Duration of the exposure with shutter closed (seconds).
observation_type: string
The observation type of this exposure (e.g. dark, bias, science).
observation_reason: string
The reason this observation was taken. (e.g. science, filter scan,
unknown).
day_obs: int
Day of observation as defined by the observatory (YYYYMMDD
format).
seq_num: int
Counter for the observation within a larger sequence. Context of
the sequence number is observatory specific. Can be a global
counter or counter within day_obs.
group_name: string
String group identifier associated with this exposure by the
acquisition system.
group_id: int
Integer group identifier associated with this exposure by the
acquisition system.
target_name: string
Object of interest for this observation or survey field name.
science_program: string
Observing program (survey, proposal, engineering project)
identifier.
tracking_ra: float
Tracking ICRS Right Ascension of boresight in degrees. Can be NULL
for observations that are not on sky.
tracking_dec: float
Tracking ICRS Declination of boresight in degrees. Can be NULL for
observations that are not on sky.
sky_angle: float
Angle of the instrument focal plane on the sky in degrees. Can be
NULL for observations that are not on sky, or for observations
where the sky angle changes during the observation.
zenith_angle: float
Angle in degrees from the zenith at the start of the exposure.
timespan: timespan
For most dimensions and most data repositories, the number of records is quite large, so you’ll almost always want a very constraining where argument to control what’s returned, e.g.:
$ butler query-dimension-records /repo/main detector \
--where "instrument='HSC' AND detector.id IN (6..8)"
instrument id full_name name_in_raft raft purpose
---------- --- --------- ------------ ---- -------
HSC 6 1_44 44 1 SCIENCE
HSC 7 1_45 45 1 SCIENCE
HSC 8 1_46 46 1 SCIENCE
When working with repositories of transient, cached datasets, note that dimension values may be retained in the registry for datasets that no longer exist (e.g. for provenance purposes) and may sometimes be present for datasets that do not yet exist.
Datasets¶
Butler.query_datasets is used to query for DatasetRef objects - handles that point directly to something at least approximately like a file on disk.
These correspond directly to what can be retrieved with Butler.get.
Because there are usually many datasets in a data repository (even in a single collection), this also isn’t a great tool for general exploration; it’s perhaps most useful as a way to explore things like the thing you’re looking for (perhaps because a call to Butler.get unexpectedly failed), by looking with similar collections, dataset types, or data IDs.
query_datasets usually isn’t what you want if you’re looking for raw-image metadata (use query_dimension_records instead); it’s easy to confuse the dimensions that represent observations with instances of the raw dataset type, because they are always ingested into the data repository together.
Data Ids¶
Butler.query_data_ids is used to query for combinations of dimension values that could be used to identify datasets.
The most important thing to know about query_data_ids is when not to use it:
It’s usually not what you want if you’re looking for datasets that already exist (use
query_datasetsinstead).It’s usually not what you want if you’re looking for metadata associated with those data ID values (use
query_dimension_records). Whilequery_data_idscan do that, too (via thewith_dimension_recordsparameter), it’s overkill if you’re looking for metadata that corresponds to a single dimension rather than all of them.
query_data_ids is most useful when you want to query for future datasets that could exist, such as when debugging empty QuantumGraphs.
Where can I find documentation for command-line butler queries?¶
The butler command line tool uses a plugin system to allow packages downstream of daf_butler to define their own butler subcommands.
Unfortunately, this means there’s no single documentation page that lists all subcommands; each package has its own page documenting the subcommands it provides.
The daf_butler and obs_base pages contain most subcommands, but the best way to find them all is to use --help on the command-line.
The pipetask tool is implemented entirely within ctrl_mpexec.
Why do queries using Registry methods return duplicate results?¶
Note
Modern Butler query methods (Butler.query_datasets, Butler.query_data_ids, and Butler.query_dimension_records) no longer return duplicate results.
The information in this section only applies when using the old Registry interface.
The Registry.queryDataIds, queryDatasets, and queryDimensionRecords methods can sometimes return true duplicate values, simply because the SQL queries used to implement them do.
You can always remove those duplicates by wrapping the calls in set(); the DataCoordinate, DatasetRef, and DimensionRecord objects in the returned iterables are all hashable.
This is a conscious design choice; these methods return lazy iterables in order to handle large results efficiently, and that rules out removing duplicates inside the methods themselves.
We similarly don’t want to always remove duplicates in SQL via SELECT DISTINCT, because that can be much less efficient than deduplication in Python, but in the future we may have a way to turn this on explicitly (and may even make it the default).
We do already remove these duplicates automatically in the butler command-line interface.
It is also possible for queryDatasets (and the butler query-datasets command) to return datasets that have the same dataset type and data ID from different collections,
This can happen even if the users passes only collection to search, if that collection is a CHAINED collection (because this evaluates to searching one or more child collections).
These results are not true duplicates, and will not be removed by wrapping the results in set().
They are best removed by passing findFirst=True (or --find-first), which will return - for each data ID and dataset type - the dataset from the first collection with a match.
For example, from the command-line, this command returns one calexp from each of the given collections:
$ butler query-datasets /repo/main calexp \
--collections HSC/runs/RC2/w_2021_06/DM-28654 \
--collections HSC/runs/RC2/w_2021_02/DM-28282 \
--where "instrument='HSC' AND visit=1228 AND detector=40"
type run id band instrument detector physical_filter visit_system visit
------ ----------------------------------- ------- ---- ---------- -------- --------------- ------------ -----
calexp HSC/runs/RC2/w_2021_02/DM-28282/sfm 5928697 i HSC 40 HSC-I 0 1228
calexp HSC/runs/RC2/w_2021_06/DM-28654/sfm 5329565 i HSC 40 HSC-I 0 1228
(with no guaranteed order!) while adding --find-first yields only the calexp found in the first collection:
$ butler query-datasets /repo/main calexp --find-first \
--collections HSC/runs/RC2/w_2021_06/DM-28654 \
--collections HSC/runs/RC2/w_2021_02/DM-28282 \
--where "instrument='HSC' AND visit=1228 AND detector=40"
type run id band instrument detector physical_filter visit_system visit
------ ----------------------------------- ------- ---- ---------- -------- --------------- ------------ -----
calexp HSC/runs/RC2/w_2021_06/DM-28654/sfm 5329565 i HSC 40 HSC-I 0 1228
Passing findFirst=True or --find-first requires the list of collections to be clearly ordered, however, ruling out wildcards like ... (“all collections”), globs, and regular expressions.
Single-dataset search methods like Butler.get and Butler.find_dataset always use the find-first logic (and hence always require ordered collections).
Why are some keys (usually filters) sometimes missing from data IDs?¶
While most butler methods accept regular dictionaries as data IDs, internally we standardize them into instances of the DataCoordinate class, and that’s also what will be returned by Butler and Registry methods.
Printing a DataCoordinate can sometimes yield results with a confusing ... in it:
>>> dataId = butler.registry.expandDataId(instrument="HSC", exposure=903334)
>>> print(dataId)
{instrument: 'HSC', exposure: 903334, ...}
And similarly asking for its keys doesn’t show everything you’d expect (same for values or items); in particular, there are no physical_filter or band keys here, either:
>>> print(dataId.keys())
{instrument, exposure}
The quick solution to these problems is to use DataCoordinate.full, which is another more straightforward Mapping that contains all of those keys:
>>> print(dataId.full)
{band: 'r', instrument: 'HSC', physical_filter: 'HSC-R', exposure: 903334}
You can also still use expressions like dataId["band"], even though those keys seem to be missing:
>>> print(dataId["band"])
r
The catch is these solutions only work if DataCoordinate.hasFull returns True; when it doesn’t, accessing DataCoordinate.full will raise AttributeError, essentially saying that the DataCoordinate doesn’t know what the filter values are, even though it knows other values (i.e. the exposure ID) that could be used to fetch them.
The terminology we use for this is that {instrument, exposure} are the required dimensions for this data ID and {physical_filter, band} are implied dimensions:
>>> dataId.graph.required
{instrument, exposure}
>>> dataId.graph.implied
{band, physical_filter}
The good news is that any DataCoordinate returned by the Registry query methods will always have hasFull return True, and you can use Registry.expandDataId to transform any other DataCoordinate or dict data ID into one that contains everything the database knows about those values.
The obvious follow-up question is why DataCoordinate.keys and stringification don’t just report all of they key-value pairs the object actually knows, instead of hiding them.
The answer is that DataCoordinate is trying to satisfy a conflicting set of demands on it:
We want it to be a
collections.abc.Mapping, so it behaves much like thedictobjects often used informally for data IDs.We want a
DataCoordinatethat only knows the value for required dimensions to compare as equal to any data ID with the same values for those dimensions, regardless of whether those other data IDs also have values for implied dimensions.collections.abc.Mappingdefines equality to be equivalent to equality overitems(), so if one mapping includes more keys than the other, they can’t be equal.
Our solution was to make it so DataCoordinate is always a Mapping over just its required keys, with full available sometimes as a Mapping over all of them.
And because the Mapping interface doesn’t prohibit us from allowing __getitem__ to succeed even when the given value isn’t in keys, we support that for implied dimensions as well.
It’s possible it would have been better to just not make it a Mapping at all (i.e. remove keys, values, and items in favor of other ways to access those things).
DataCoordinate has already been through a number of revisions, though, and it’s not clear it’s worth yet another try.
How do I avoid errors involving queries for calibration datasets?¶
Note
The modern Butler.query_datasets method is able to search in calibration collections.
The information in this section only applies when using the old Registry interface.
Registry.queryDatasets currently has a major limitation in that it can’t query for datasets within a CALIBRATION collection; the error message looks like this:
NotImplementedError: Query for dataset type 'flat' in CALIBRATION-type collection 'HSC/calib' is not yet supported.
We do expect to fix this limitation in the future, but it may take a while. In the meantime, there are a few ways to work around this problem.
First, if you don’t actually want to search for calibrations at all, but this exception is still getting in your way, you can make your query more specific.
If you use a dataset type list or pattern (a shell-style glob on the command line, or re.compile in the Python interface) that doesn’t match any calibration datasets, this error should not occur.
Similarly, if you can use a list of collections or a collection pattern that doesn’t include any CALIBRATION collections, that will avoid the problem as well - but this is harder, because CHAINED collections that include CALIBRATION collections are quite common.
For example, both processing-output collections with names like “HSC/runs/w_2025_06/DM-50000” and per-instrument default collections like “HSC/defaults” include a CALIBRATION child collection.
You can recursively expand a collection list and filter out any child CALIBRATION collections from it with this snippet:
expanded = list(
butler.registry.queryCollections(
original,
flattenChains=True,
collectionTypes=(CollectionType.all - {CollectionType.CALIBRATION}),
)
)
where original is the original, unexpanded list of collections to search.
The equivalent command-line invocation is:
$ butler query-collections /repo/main --chains=flatten \
--collection-type RUN \
--collection-type CHAINED \
--collection-type TAGGED \
HSC/defaults
Name Type
--------------------------------- ----
HSC/raw/all RUN
HSC/calib/gen2/20180117/unbounded RUN
HSC/calib/DM-28636/unbounded RUN
HSC/masks/s18a RUN
refcats/DM-28636 RUN
skymaps RUN
Another possible workaround is to make the query much more general - passing collections=... to search all collections in the repository will avoid this limitation even for calibration datasets, because it will take advantage of the fact that all datasets are in exactly one RUN collection (even if they can also be in one or more other kinds of collection) by searching only all of the RUN collections.
That same feature of RUN collections can also be used with Registry.queryCollections (and our naming conventions) to find calibration datasets that might belong to particular CALIBRATION collections.
For example, if “HSC/calib” is a CALIBRATION collection (or a pointer to one), the datasets in it will usually also be present in RUN collections that start with “HSC/calib/”, so logic like this might be useful:
run_collections = list(
butler.registry.queryCollections(
re.compile("HSC/calib/.+"),
collectionTypes={CollectionTypes.RUN},
)
)
Or, from the command-line,
$ butler query-collections /repo/main --collection-type RUN \
HSC/calib/gen2/20200115/*
Name Type
---------------------------------------- ----
HSC/calib/gen2/20200115/20170821T000000Z RUN
HSC/calib/gen2/20200115/20160518T000000Z RUN
HSC/calib/gen2/20200115/20170625T000000Z RUN
HSC/calib/gen2/20200115/20150417T000000Z RUN
HSC/calib/gen2/20200115/20181207T000000Z RUN
HSC/calib/gen2/20200115/20190407T000000Z RUN
HSC/calib/gen2/20200115/20150407T000000Z RUN
HSC/calib/gen2/20200115/20160114T000000Z RUN
HSC/calib/gen2/20200115/20170326T000000Z RUN
...
The problem with this approach is that it may return many datasets that aren’t in “HSC/calib”, including datasets that were not certified, and (like all of the previous workarounds) it doesn’t tell you anything about the validity ranges of the datasets that it returns.
If you just want to load the calibration dataset appropriate for a particular raw (and you have the data ID for that raw in hand), the right solution is to use Butler.get with that raw data ID, which takes care of everything for you:
flat = butler.get(
"flat",
instrument="HSC", exposure=903334, detector=0,
collections="HSC/calib"
)
The lower-level Butler.find_dataset method can also perform this search without actually reading the dataset, but you’ll need to be explicit about how to do the temporal lookup:
raw_data_id = butler.registry.expandDataId(
instrument="HSC",
exposure=903334,
detector=0,
)
ref = butler.find_dataset(
"flat",
raw_data_id,
timespan=raw_data_id.timespan,
)
It’s worth noting that find_dataset doesn’t need or use the exposure key in the raw_data_id argument that is passed to it - a master flat isn’t associated with an exposure - but it’s happy to ignore it, and we do need it (or something else temporal) in order to get a data ID with a timespan for the last argument.
Finally, if you need to query for calibration datasets and their validity ranges, and don’t have a point in time you’re starting from, the only option is Registry.queryDatasetAssociations.
That’s a bit less user-friendly - it only accepts one dataset type at a time, and doesn’t let you restrict the data IDs at all - but it can query CALIBRATION collections and it returns the associated validity ranges as well.
It actually only exists as a workaround for the fact that queryDatasets can’t do those things, and it will probably be removed sometime after those limitations are lifted.
How do I fix an empty QuantumGraph?¶
The pipetask tool attempts to predict all of the processing a pipeline will perform in advance, representing the results as a QuantumGraph object that can be saved or directly executed.
When that graph is empty, it means it thinks there’s no work to be done, and unfortunately this is both a common and hard-to-diagnose problem.
The QuantumGraph generation algorithm begins with a large SQL query (a complicated invocation of Registry.queryDataIds, actually), where the result rows are essentially data IDs and the result columns are all of the dimensions referenced by any task or dataset type in the pipeline.
Queries for all "regular input" datasets (i.e. not PrerequisiteInputs”) are included as subqueries, spatial and temporal joins are automatically included, and the user-provided query expression is translated into an equivalent SQL WHERE clause.
That means there are many ways to get no result rows - and hence an empty graph.
Sometimes we can tell what will go wrong even before the query is executed - the butler maintains a summary of which dataset types are present each each collection, so if the input collections don’t have any datasets of a needed type at all, a warning log message will be generated stating the problem. This will also catch most cases where a pipeline is misconfigured such that what should be an intermediate dataset isn’t actually being produced in the pipeline, because it will appear instead as an overall input that (usually) won’t be present in those input collections.
We also perform some follow-up queries after generating an empty QuantumGraph, to see if any needed dimensions are lacking records entirely (the most common example of this case is forgetting to define visits after ingesting raws in a new data repository).
If you get an empty QuantumGraph without any clear explanations in the warning logs, it means something more complicated went wrong in that initial query, such as the input datasets, available dimensions, and boolean expression being mutually inconsistent (e.g. not having any bands in common, or tracts and visits not overlapping spatially).
In this case, the arguments to queryDataIds will be logged again as warnings, and the next step in debugging is to try that call manually with slight adjustments.
To guide this process, it can be very helpful to first use pipetask build --show pipeline-graph to create a diagram of the pipeline graph - a simpler directed acyclic graph that relates tasks to dataset types, without any data IDs:
$ pipetask build ... --show pipeline-graph
○ camera
│
○ │ raw
│ │
◍ │ │ yBackground, transmission_sensor, transmi...[1]
├─┼─┤
■ │ │ isr
│ │ │
○ │ │ postISRCCD
│ │ │
■ │ │ characterizeImage
│ │ │
◍ │ │ icSrc, icExpBackground, icExp
│ │ │
○ │ │ │ ps1_pv3_3pi_20170110
├─┤ │ │
│ ■ │ │ calibrate
(...)
The --pipeline-dot argument can also be used to create a version of this graph in the GraphViz dot language, and you can use the ubiquitous dot command-line tool to transform that into a PNG, SVG, or other graphical format file:
$ pipetask build ... --pipeline-dot pipeline.dot
$ dot pipeline.dot -Tsvg > pipeline.svg
That ... should be replaced by most of the arguments you’d pass to pipetask that describe what to run (which tasks, pipelines, configuration, etc.), but not the ones that describe how, or what to use as inputs (no collection options).
See pipetask build --help for details.
This graph will often reveal some unexpected input dataset types, tasks, or relationships between the two that make it obvious what’s wrong.
Another useful approach is to try to simplify the pipeline, ideally removing all but the first task; if that works, you can generally rule it out as the cause of the problem, add the next task in, and repeat.
Because the big initial query only involves regular inputs, it can also be helpful to change regular Input connections into PrerequisiteInput connections - when a prerequisite input is missing, pipetask should provide more useful diagnostics.
This is only possible when the dataset type is already in your input collections, rather than something to be produced by another task within the same pipeline.
But if you work through your pipeline task-by-task, and run each single-task pipeline as well as produce a QuantumGraph for it, this should be true each step of the way as well.
How can I make QuantumGraph generation faster?¶
QuantumGraph generation can be slow in several different ways for different pipelines and datasets, and the first step in speeding it up is to look at the logs to see where it’s spending its time.
We strongly recommend passing --long-log to include timestamps in all logging, and passing --log-level lsst.pipe.base.quantum_graph_builder=VERBOSE can provide more information about QuantumGraph generation in particular.
If you’re running BPS, logs for this step are written to quantumGraphGeneration.out in the submit directory.
Here’s what’s going on after a few important log messages (all INFO level):
Processing pipeline subgraph X of Y with N task(s).: we’re running the “big initial query” for all of the data IDs that might appear in the graph. This step is usually quite fast (seconds or minutes for large graphs, not hours), but occasionally catastrophically slow (days) when the database’s query optimizer chooses a bad plan, so it’s the step most amenable to big speedups (more on this below).Iterating over query results to associate quanta with datasets.: we’re processing the result rows of that big query, each of which will correspond to an edge or a set of similar edges in the graph. This step is pure Python (no database queries), and the only way to make it faster is to shrink the size of the problem by splitting it up. Splitting the task into steps may help more than splitting up data when this is the bottleneck, but only slightly.Initial bipartite graph has 290189 quanta, 1073224 dataset nodes, and 3591767 edges from 234155 query row(s).: the preliminary graph is built, and now we’re performing many smaller database queries to look for input datasets (or outputs that may be in the way, in some cases), and asking each task if each of its quanta should be kept or pruned out. This step is usually very close to linear in the number of quanta and is typically dominated by Python logic, but it does involve some database queries. TheVERBOSElogging can provide information about exactly which dataset it’s querying for, and if any of these seem to be unusually slow, please report it to the middleware team (with logs and a link to the pipeline you’re running). There’s not much a user can do about slowdowns here (aside from splitting up the problem).When the graph has been built,
pipetaskwill print a table with the number of quanta for each task. If it pauses a long time after this, it’s probably spending a long time writing the graph to disk. This takes longer than it should (this is a known issue we have plans to fix, but it’ll require some deep changes), but it should be linear in the number of quanta in the graph.
When the “big initial query” is catastrophically slow, it’s almost always because the query is complex enough that the database’s query optimizer chose to execute it in a way that didn’t take advantage of the right index, and our goal is to give it an equivalent or nearly-equivalent query that’s simpler.
By default, the query includes both the --data-query expression provided by the user and joins to a subqueries for each regular input dataset in the pipeline (but not prerequisites).
The best way to simplify the query is to eliminate as many of those dataset subqueries as you can via the --dataset-query-constraint option, which provides direct control over the dataset types to join against.
If you can easily write a --data-query argument that includes all of the data IDs you want to process and almost no data IDs you don’t want to process (like an explicit tract or visit range), pass --dataset-query-constraint off to get rid of all of the dataset subqueries.
When that’s not easy, try to identify one input dataset type whose existence strongly implies the others (perhaps because they’re all produced together by some previous processing), and pass that as the argument to --dataset-query-constraint.
Visualizing the pipeline as a graph (see e.g. pipetask build --show pipeline-graph, as described in How do I fix an empty QuantumGraph?) is the best way to do this.
Dataset types with data IDs that are more similar to the data IDs of the quanta are probably best, and dataset types with coarser data IDs are probably better choices than those with finer data IDs (e.g. prefer tract over patch, visit over {visit, detector}), but this is based on intuition, not experience, and the most important thing is to reduce the number of dataset types down to zero or one.
In most cases, a complex --data-query argument is preferable to even one input dataset constraint, but there are exceptions:
If the
--data-queryreferences a dimension that is completely irrelevant to the graph (e.g. putting anexposureconstraint into a graph that only uses{tract, patch}data IDs), it can really slow things down, because it still gets included in the query and the number of result rows is multiplied by the number of matching irrelevant-dimension values (e.g. the numberexposures). The fact that theexposuredimension is not spatial (butvisitis) interacts with this in a particularly dramatic way: while it’s fine to add a constraint ontractorpatchto spatially control the avisit-based pipeline, if you do this on a pipeline that only referencesexposure, notvisit(like ISR alone), the query system will not recognize that it needs to usevisitto mediate betweenexposureandtract/patch, and a disastrously huge query will be the result.If the
--data-queryreferences dimension metadata fields rather than primary key values (e.g.visit.exposure_timerather than justvisit), we may not have indexes in place to make those selections fast. Note that this includes theseqnumfield ofvisitandexposure, and - until the repositories are migrated to the latest dimension universe -day_obsas well. We haven’t actually observed this ever leading to catastrophic query performance, so it’s not worth worrying about unless you’re trying to fix a graph-generation problem that you know is slow, and if you do think this is a problem for you, please report it so we can add indexes in the future.
Finally, while we haven’t seen this problem in the wild (perhaps because --dataset-query-constraint is underused), if the combination of the --data-query and --dataset-query-constraint arguments leave the query underconstrained, it might run quickly but return many more result rows than we need.
For example, if one passes --dataset-query-constraint off and the --data-query matches 1000 visits while only 10 of those have inputs, the initial query will return a factor of 100 more result rows than it might need - and while the initial query may still be fast enough to avoid being the bottleneck, this will result in a preliminary graph that is too big and needs to be pruned considerably by the follow-up queries for input datasets, making later steps of the process 100x slower.
What do I do if a query method/command is slow?¶
Adding the --log-level sqlalchemy.engine=DEBUG option to the butler or pipetask command will allow the SQL queries issued by the command to be inspected.
Similarly, for a slow query method, adding logging.getLogger("sqlalchemy.engine").setLevel(logging.DEBUG) can help.
The resulting query logs can be useful for developers and database administrators to determine what, if anything, is going wrong.
How can I get a report on all failures and missing datasets in a run?¶
The pipetask report tool can be used to analyze executed quantum graphs, troubleshoot, diagnose failures, and confirm that fail-and-recovery attempts (such as when using --skip-existing-in) are effective.
When analyzing multiple graphs with pipetask report, all graphs should be attempts to execute the same pipeline with the same dataquery.
The recommended usage is
pipetask report --full-output-filename <path/to/output_file>.json --force-v2 REPO QGRAPHS
The
--full-output-filename <path/to/output_file>.jsonoption provides the path to a file where the output of the full summary information can be storedThe
--force-v2option makes sure that the most recent version of the tool is used even when the user passes only one graphThe
REPOargument is theButlerrepo where the output from the processing is storedThe
QGRAPHSargument is aSequenceof `QuantumGraph`s to be analyzed, separated by spaces and passed in order of first to last executed
Note
If the --full-output-filename argument is not provided, information on the error messages and data IDs for all failed quanta will be output to the command line, which can be overwhelming in the case of many failures.
If you do not need error information, consider the --brief option, which prints the same tables to the command line as --full-output-filename, but does not save additional information to a file.
This will print two bps report-style tables, one for quanta and one for output datasets.
In the output JSON file will be
A summary under every task with:
Every failed data ID and corresponding error message
Every run containing failing data IDs, and their status
A list of data IDs which have been “recovered”; i.e., successes from fail-and-recovery attempts
A list of the data IDs associated with every missing dataset
A field called
producerconnecting eachdatasetTypefor each data ID to the task which produced itCounts of quanta and datasets in all possible states.
Currently, the --force-v2 option is the suggested usage until version 1 of pipetask report (using the QuantumGraphExecutionReport instead of the QuantumProvenanceGraph) is deprecated.
With that, we’ll go into how to read the output of pipetask report.
The Quanta Table¶
The table for Quanta from the w_2024_28 DC2 test-med-1 reprocessing run, before recoveries, looked like:
which indicates that there are 6 failed templateGen quanta, 14 failed assembleCoadd quanta, 14 blocked detection quanta, and so on.
These totals are all determined by examining the status of each task run on each data ID from run to run.
Status Definitions for Task Quanta¶
Unknown¶
The Unknown category could mean one of two things.
Being the default value for the status of a task run on a particular data ID, it could technically mean that no attempt has been made to execute this part of the graph.
However, in a table that does not look empty, it much more likely denotes that the task’s metadata is missing for that data ID.
This makes it impossible for us to tell whether the task succeeded or failed, but does point to possible infrastructure problems.
Successful¶
This status is the trademark of a successful quantum.
The specific path to being marked as successful by the QuantumProvenanceGraph is that metadata and log datasets exist for the task for the data ID in question.
Blocked¶
Blocked quanta are the successors of unsuccesful quanta which cannot be executed because their inputs are not present.
They do not represent failure per se, but since they could not execute, they also do not produce output data products.
The way that the QuantumProvenanceGraph identifies a Blocked quantum is that it has no metadata and no logs (so it did not start or finish execution) and that the task run on that particular data ID is a successor (according to the quantum graph) of a Failed quantum.
Failed¶
There are log datasets but no metadata datasets for this task and data ID, indicating that the task started but did not finish execution.
Wonky¶
A Wonky quantum is the result of infrastructure problems or concerning middleware mismatches.
The category is intended to halt processing and require human intervention to proceed. As such, when a task is marked as Wonky for a particular data ID, no further successes (or any other statuses, for that matter) will change the overall status out of Wonky.
A quantum can only exit a Wonky state via human intervention. Currently, there are three paths to the Wonky state for a quantum:
A quantum which was marked as
Successfulon a previous processing attempt (run) has a more recent attempt (run) which theQuantumProvenanceGraphidentifies as unsuccessful (i.e., Graph 1 says task a ran successfully on data ID x, but Graph 2 says task a’s attempt at data ID x wasFailed,Blocked, orUnknown)Logs are missing for at least one of the attempts to run this task on this data ID
Registry.queryDatasetsfor the output datasets of this quantum return outputs from multiple different processing attempts (runs)
The outputs of Registry.queryDatasets are important because they are the datasets which will be used as inputs to downstream tasks. If the inputs to a downstream task are from different processing attempts, the Butler cannot ensure that they have been processed in the same way, with the same inputs, dataqueries, etc.
Total¶
The sum of all the previous categories.
Expected¶
The expected number of quanta, according to the quantum graphs.
The Dataset Table¶
The (abbreviated because there are many datasets) table for the output datasets from the same w_2024_38 DC2 test-med-1 run, before recoveries, looked like
Dataset Visible Shadowed Predicted Only Unsuccessful Cursed TOTAL EXPECTED
-------------------------------------------- ------- -------- -------------- ------------ ------ ----- --------
deepCoadd_directWarp 6384 0 457 0 0 6841 6841
deepCoadd_psfMatchedWarp 6369 0 472 0 0 6841 6841
makeWarp_metadata 6841 0 0 0 0 6841 6841
makeWarp_log 6841 0 0 0 0 6841 6841
deepCoaddVisits 294 0 0 0 0 294 294
selectDeepCoaddVisits_metadata 294 0 0 0 0 294 294
selectDeepCoaddVisits_log 294 0 0 0 0 294 294
goodSeeingVisits 294 0 0 0 0 294 294
selectGoodSeeingVisits_metadata 294 0 0 0 0 294 294
selectGoodSeeingVisits_log 294 0 0 0 0 294 294
goodSeeingCoadd 288 0 0 6 0 294 294
goodSeeingCoadd_nImage 288 0 0 6 0 294 294
templateGen_metadata 288 0 0 6 0 294 294
templateGen_log 288 0 0 6 0 294 294
deepCoadd 280 0 0 14 0 294 294
deepCoadd_inputMap 280 0 0 14 0 294 294
deepCoadd_nImage 280 0 0 14 0 294 294
assembleCoadd_metadata 280 0 0 14 0 294 294
assembleCoadd_log 280 0 0 14 0 294 294
deepCoadd_det 280 0 0 14 0 294 294
deepCoadd_calexp 280 0 0 14 0 294 294
deepCoadd_calexp_background 280 0 0 14 0 294 294
...
which shows Unsuccessful datasets for all the outputs of the failed getTemplate and assembleCoadd tasks, as well as the outputs of their (Blocked) successors. It also shows Predicted Only for some deepCoadd_directWarp and deepCoadd_psfMatchedWarp datasets, which is a result of these outputs being predicted when the QuantumGraph was built, but found to be unnecessary by the Science Pipelines during execution.
Status Definitions for Datasets¶
The statuses for datasets all pertain to whether a datasetType associated with a certain data ID exists in the Butler, and whether it could be used as an input to downstream tasks.
Visible¶
The dataset exists, and comes up in a find-first search (Registry.queryDatasets). This means it will be used as an input for downstream tasks.
Shadowed¶
The dataset exists in the Butler, but does not come up in a find-first search (Registry.queryDatasets). Therefore, we don’t have to worry about it being used as an input.
Predicted Only¶
This status occurs when the graph predicts an output dataset that ultimately is not produced by the Science Pipelines because it was deemed unnecessary when the graph was executed. These are commonly referred to as “NoWorkFound cases” because their output messages say “No work found for task” (and hence said task does not produce output datasets).
Unsuccessful¶
An unsuccessful status for a dataset just means it does not exist in the Butler. Unsuccessful datasets are the results of Failed and Blocked quanta.
Cursed¶
A Cursed dataset is the result of an unsuccessful quantum which would otherwise be marked as Visible. This means that if any datasets which were produced by non-Successful quanta (this includes Unknown, Failed, Blocked and Wonky quanta) come up in a Registry.queryDatasets find-first search, we flag them as Cursed in order to halt processing, start investigation, and prevent them from being used as inputs later. Ideally Cursed datasets should be removed before processing continues.
Total¶
The sum of all the previous categories.
Expected¶
The expected number of datasets, according to the quantum graphs.
Note
The number in Total should always match the number in Expected.
If it does not, this is a problem with this algorithm and should be reported to the developers.
Combining group-level pipetask report summaries into a step-level rollup¶
pipetask report works on the group-level (on attempts to execute the same pipeline with the same data-query). However, it is possible to collate the JSON output from multiple groups into one file. This is intended for combining group-level summaries into summaries over processing steps, and answering questions like “What are all the errors that occurred in step 1?” or “How many quanta were blocked in step 5?”
Recommended usage is
pipetask aggregate-reports --full-output-filename <path/to/combined/output_file>.json <path/to/group/file/1.json> <path/to/group/file/2.json>...
where the argument to --full-output-filename is a filepath to store the combined QuantumProvenanceGraph summary.
If no argument is passed, the combined summary will be printed to the screen.
How do I clean up processing runs I don’t need anymore?¶
Because a data repository stores information on both a filesystem or object store and a SQL database, deleting datasets completely requires using butler commands, even if you know where the associated files are stored on disk.
For processing runs that follow our usual conventions (following them is automatic if you use --output and don’t override --output-run when running pipetask), two different collections are created:
a
RUNcollection that directly holds your outputsa
CHAINEDcollection that points to that RUN collection as well as all of your input collections.
If you perform multiple processing runs with the same --output, you’ll get multiple RUN collections in the same CHAINED collection.
The CHAINED collection will have the name you passed to --output, and the RUN collections will start with that and end with a timestamp.
You can see this structure for your own collections with a command like this one:
$ butler query-collections /repo/main --chains=tree u/jbosch/*
u/jbosch/DM-30649 CHAINED
u/jbosch/DM-30649/20210614T191615Z RUN
HSC/raw/RC2/9813 TAGGED
HSC/calib/gen2/20180117 CALIBRATION
HSC/calib/DM-28636 CALIBRATION
HSC/calib/gen2/20180117/unbounded RUN
HSC/calib/DM-28636/unbounded RUN
HSC/masks/s18a RUN
HSC/fgcmcal/lut/RC2/DM-28636 RUN
refcats/DM-28636 RUN
skymaps RUN
u/jbosch/DM-30649/20210614T191615Z RUN
The RUN collections that directly hold the datasets are what we want to remove in order to free up space, but we have to start by deleting the CHAINED collections that hold them first:
$ butler remove-collections /repo/main u/jbosch/DM-30649
You can add the --no-confirm option to skip the confirmation prompt if you like.
If you’re only deleting one collection at a time, it doesn’t tell you anything new.
Not deleting the CHAINED collection¶
If you don’t want to remove the CHAINED collection - you just want to remove the RUN collection from it - you can instead do
$ butler collection-chain /repo/main –remove u/jbosch/DM-30649 u/jbosch/DM-20210614T191615Z
Or, if you know the RUN is the first one in the chain,
$ butler collection-chain /repo/main –pop u/jbosch/DM-30649
In any case, once the CHAINED collection is out of the way, we can delete the RUN collections that start with the same prefix using a glob pattern:
$ butler remove-runs /repo/main u/jbosch/DM-30649/*
The following RUN collections will be removed:
u/jbosch/DM-30649/20210614T191615Z
The following datasets will be removed:
calexp(18222), calexpBackground(18222), calexp_camera(168), calibrate_config(1), calibrate_metadata(18222), characterizeImage_config(1), characterizeImage_metadata(18231), consolidateSourceTable_config(1), consolidateVisitSummary_config(1), consolidateVisitSummary_metadata(168), fgcmBuildStarsTable_config(1), fgcmFitCycle_config(1), fgcmOutputProducts_config(1), icExp(18231), icExpBackground(18231), icSrc(18231), icSrc_schema(1), isr_config(1), isr_metadata(18232), postISRCCD(18232), skyCorr(17304), skyCorr_config(1), skyCorr_metadata(168), source(18222), src(18222), srcMatch(18222), srcMatchFull(18222), src_schema(1), transformSourceTable_config(1), visitSummary(168), writeSourceTable_config(1), writeSourceTable_metadata(18222)
Continue? [y/N]: y
Removed collections
Here we’ve left the default confirmation behavior on because we used a glob, just to be safe. You can write one or more full RUN collection names explicitly (separated by commas), too, and that’s what you’ll need to do if you didn’t follow the naming convention well enough for a glob to work.
Removing RUN collections always removes the files within them, but it does not remove the directory structure, because in the presence of arbitrary path templates (including any that may have been used in the past) and possible concurrent writes, it’s difficult for the butler to recognize efficiently when a directory will end up empty.
You’re welcome to delete empty directories on your own after using remove-runs; they’re typically in subdirectories of the main repository directory named after the collection (it’s possible to configure the butler such that this isn’t the case, but rare).
It’s also completely fine to just leave them.
Note
If you delete files from the filesystem before using butler commands to remove entries from the database, the commands for cleaning up the database are actually exactly the same. The butler won’t know that the files are gone until you try to use or delete them, but when you try to delete them, it will just log this at debug level.
Deleting only some datasets¶
If you don’t want to delete the full RUN collection, just some datasets within it, you can generally use the prune-datasets subcommand:
$ butler prune-datasets /repo/main --purge u/jbosch/DM-29776/singleFrame/20210426T161854Z --datasets postISRCCD u/jbosch/DM-29776/singleFrame/20210426T161854Z
The following datasets will be removed:
type run id band instrument detector physical_filter exposure
---------- ---------------------------------------------- ------------------------------------ ---- ---------- -------- --------------- --------
postISRCCD u/jbosch/DM-29776/singleFrame/20210426T161854Z c45a177f-24e8-4dc9-9268-5895decb7989 y HSC 0 HSC-Y 318
postISRCCD u/jbosch/DM-29776/singleFrame/20210426T161854Z 461d0293-3c80-45ea-9f06-21a90525c185 y HSC 1 HSC-Y 318
postISRCCD u/jbosch/DM-29776/singleFrame/20210426T161854Z 1572dd02-c959-4d23-ba03-91cf235e1291 y HSC 2 HSC-Y 318
postISRCCD u/jbosch/DM-29776/singleFrame/20210426T161854Z b38afec9-1970-478d-80d8-4f61c5a992d0 y HSC 3 HSC-Y 318
postISRCCD u/jbosch/DM-29776/singleFrame/20210426T161854Z 769bb9ce-9267-4e57-812f-82fee3fd0afa y HSC 4 HSC-Y 318
(...)
Continue? [y/N]: y
The datasets were removed.
Note that here you have to know the exact RUN collection that holds the datasets, and specify it twice (the argument to --purge is the collection to delete from, while the positional argument is the collection to query within - the latter could be some other kind of collection, but it’s rare for that to be useful).
The Python Butler.pruneDatasets method can be used for even greater control of what you want to delete, as it accepts an arbitrary DatasetRef iterable indicating what to delete.