Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

DM-20382: Add description of RC2 ap_pipe repos #12

Closed
wants to merge 1 commit into from

Conversation

yalsayyad
Copy link
Contributor

No description provided.

@yalsayyad yalsayyad changed the title Add description of RC2 ap_pipe repos DM-20382: Add description of RC2 ap_pipe repos Jun 28, 2019
@yalsayyad
Copy link
Contributor Author

Includes all and "best seeing" visit lists for tract 9813.

Copy link
Member

@kfindeisen kfindeisen left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I have some things I'd like clarified, but looks basically good.

@@ -0,0 +1,84 @@
# Running ap_pipe on HSC RC2 data from PDR1

When PDR2 hits lsst-dev we'll recreate an AP-style repo using PDR1 templates to difference new visits in PDR2. Until then, this repo will demonstrate running ap_pipe on a dataset more like DRP mode, where the templates are the best seeing images within a data release and all visits in that data release are differenced.
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I don't know what a PDR is, and other readers might not, either.

* `/datasets/hsc/repo/rerun/private/yusra/RC2/ap_pipe_PDR1_bestcase` (HSC-I band)
* `/datasets/hsc/repo/rerun/private/yusra/RC2/ap_pipe_PDR1_bestcase_noSkyCorr` (HSC-G, HSC-R, HSC-I, HSC-Z)

These instructions go into detail on the second was generated.
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

missing "how"?


These instructions go into detail on the second was generated.

To remember the name, the dataset only contails "PDR1", is "bestcase" in that the templates are chosen as the best seeing images, the best case scenario. The only config difference was that the HSC-specific skyCorrection was turned off.
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

"contains", plus the grammar makes it hard to follow. Consider rewording.

false_positives/ap_pipe_HSC_RC2.md Outdated Show resolved Hide resolved
export COSMOS_Z_SEEING=17928^17908^17906^17926^17904^17902^17930^17962^17946^17944
export COSMOS_Y_SEEING=22602^22604^22644^22632^22628^22642^22626^22630^22664^11730^22662^22606^11736^22646^22660^11718^11732
```
LSST has no plans for a narrowband filter, so let's put that data aside.
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I'm confused. What narrowband filter? Why bring it up in the first place?




Choose the 1/3 of visits with the best seeing. These visit lists are sorted by ascending visit FHWM, so reducing the number of visits can be accomplished by removing visits from the end:
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

At first I was confused because I assumed this would be done with a BestSeeingSelector. Maybe say "Manually choose" to make it more explicit?


```
mkdir /datasets/hsc/repo/rerun/private/yusra/RC2/ap_pipe_PDR1_bestcase_noSkyCorr
make_ppdb.py --config ppdb.isolation_level=READ_UNCOMMITTED ppdb.db_url="sqlite:////datasets/hsc/repo/rerun/private/yusra/RC2/ap_pipe_PDR1_bestcase_noSkyCorr/association.db
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I suggest using the same syntax as for the ap_pipe.py call:

Suggested change
make_ppdb.py --config ppdb.isolation_level=READ_UNCOMMITTED ppdb.db_url="sqlite:////datasets/hsc/repo/rerun/private/yusra/RC2/ap_pipe_PDR1_bestcase_noSkyCorr/association.db
make_ppdb.py --configfile ./configfile.py

Writing the calls in the same way helps ensure the database setup matches what ap_pipe.py expects. The config is safe; make_ppdb.py will simply ignore any non-PPDB fields.

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

(Though, looking closer, I do find it a bit odd that the output repo is specified both in the config file and on the ap_pipe.py command line...)

false_positives/ap_pipe_HSC_RC2.md Outdated Show resolved Hide resolved
config.ppdb.db_url="sqlite:////datasets/hsc/repo/rerun/private/yusra/RC2/ap_pipe_PDR1_bestcase_noSkyCorr/association.db"
config.ppdb.isolation_level="READ_UNCOMMITTED"
config.ppdb.connection_timeout=240
config.differencer.doWriteMatchedExp=True
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I'm a bit confused by this config. I think you told me at one point not to use doSelectSources with coadd templates, but you seem to be doing so here. What's the difference?

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

You're right that the filterMap configs are a no-op if doSelectSources is False. I'll clarify.

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Now I'm more confused. Is doSelectSources False from some other config?

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

D'oh! 🤦‍♂

config.differencer.refObjLoader.retarget(LoadIndexedReferenceObjectsTask)
config.differencer.refObjLoader.load(os.path.join(getPackageDir("obs_subaru"), "config", "filterMap.py"))

config.differencer.refObjLoader.ref_dataset_name = "ps1_pv3_3pi_20170110"
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

After DM-18394 I don't think you need differencer.refObjLoader anymore, because kernelSourcesFromRef is False.

@yalsayyad
Copy link
Contributor Author

Didn't force push after rebasing master, so GitHub is not automatically registering the merge.

@yalsayyad yalsayyad closed this Jul 8, 2019
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

2 participants