[PIPE2D-1110] Make match filterName between pfsConfig and flux calibration Created: 07/Nov/22  Updated: 20/Dec/22  Resolved: 09/Dec/22

Status: Done
Project: DRP 2-D Pipeline
Component/s: None
Affects Version/s: None
Fix Version/s: None

Type: Task Priority: Normal
Reporter: Takuji Yamashita Assignee: sogo.mineo
Resolution: Done Votes: 0
Labels: flux-calibration
Remaining Estimate: Not Specified
Time Spent: Not Specified
Original Estimate: Not Specified

Attachments: File PIPE2D-1110-change-filter-names.py    
Issue Links:
Relates
relates to OBSPROC-33 Review the pfsDesign generation scrip... Done
Reviewers: price

 Description   

Flux calibration reads broad-band information from pfsConfig. But the current flux calibration assumes different names from targetDB/gaiaDB, and therefore we need to unify the names to run flux calibration. We change filterName in flux calibration and use ones of targetDB and gaiaDB.
 
{{filterName}}s passed down from targeted and gaiaDB to pfsDesign and pfsConfig are:
[band]_ps1 for PS1
[band]_gaia for Gaia
 
Similarly, we can use {{filterName}}s of [band]_hsc and [band]_sdss for HSC and SDSS, respectively.

EDIT

PIPE2D-1110-change-filter-names.py has to be run in order not to break the integration tests and weekly tests after this branch is merged to the main branch.



 Comments   
Comment by Takuji Yamashita [ 08/Nov/22 ]

In regard to HSC r2 and i2 filters:
I propose to use r2/i2 filters instead of the old r/i filters for flux calibration unless users explicitly set the filterName of the old ones.

The HSC filters are used to scale fluxes of science targets in the last part of flux calibration. Calculating flux calibration vectors does not use HSC photometries because FLUXSTD stars have PS1 photometries. Users do not distinguish which filters, r or r2(i2), are used in the HSC catalog. The new r2 and i2 filters have been in use since 2016. Therefore, I think that it is natural to use r2 and i2 filters for our flux calibration.

For example, if i_hsc is set in filterName of pfsConfig, the flux calibration reads the filter transmission of i2.

It might be useful to have an option in which the old i filter is read when i_old_hsc (or other nice names) is set. Does targetDB accept such names for the old filters, Kiyoto Yabe ?

Comment by price [ 08/Nov/22 ]

We should adopt a naming scheme for the photometric bands that makes the filter curve explicit, without having to assume anything.

Comment by Takuji Yamashita [ 08/Nov/22 ]

Do you mean that we should use the names of i2 and i, which corresponds with the filter curve names? In that case, targetDB needs to use a name including i2 to designate the i2 filter although the HSC catalog uses a name without "2" even for i2 band. Or do we choose the i2 filter curve internally in the pipeline?

Comment by price [ 08/Nov/22 ]

Oh, the problem is the HSC SSP catalog doesn't distinguish between i1 and i2 (though other reductions may). Yeah, that's annoying. Maybe we need to map i_HSC to i2, i2_HSC to i2 and i1_HSC to i1?

Comment by Takuji Yamashita [ 08/Nov/22 ]

Yes. That is because images taken with i1 and i2 are coadded in DUD, and i1 magnitudes in WIDE are corrected to i2. A user can fetch i1 magnitudes from the old data releases, but the band name is i not i1.
 

map i_HSC to i2, i2_HSC to i2 and i1_HSC to i1

I agree with this. I prefer i_old for i1. It is a bit more meaningful, but this is subjective.

Comment by sogo.mineo [ 18/Nov/22 ]

There is already a filter named "ENG-R1", so I don't think the name "r1" is good for the older "HSC-R". Yamashita-san and I think the filter names should be:

HSC: {g,r_old,r2,i_old,i2,z,y}_hsc
PS1: {g,r,i,z,y}_ps1
Gaia: {bp,rp,g}_gaia
SDSS: {u,g,r,i,z}_sdss

I am changing filter names in fitReference.py but I am afraid that the changes break the integration test (and weekly test). psfConfig in drp_stella_data must be modified. I'm not sure I am able to do it properly. 

Comment by sogo.mineo [ 18/Nov/22 ]

Also, I believe that i_hsc and r_hsc should result in an error. It would be better than having to memorize whether i means HSC-I or HSC-I2.

Comment by sogo.mineo [ 01/Dec/22 ]

Once this branch is merged to the main branch, pfsConfig*.fits and pfsDetection*.fits have to be modified (PIPE2D-1110-change-filter-names.py can be used to modify the files).

Comment by sogo.mineo [ 02/Dec/22 ]

Could you review this PR?

Comment by price [ 06/Dec/22 ]

We need to coordinate merging this ticket with the necessary changes to the integration test and weekly. I suggest you update the integration test inputs (in drp_stella_data), run the integration test to confirm it works, and then I'll update the weekly after you merge.

Comment by sogo.mineo [ 06/Dec/22 ]

> you update the integration test inputs (in drp_stella_data),
I have already run the modified integration test. Should I push the change of drp_stella_data?

Comment by price [ 06/Dec/22 ]

Yes, please.

Comment by sogo.mineo [ 06/Dec/22 ]

I pushed tickets/PIPE2D-1110 to drp_stella_data and created a pull request for it. The integration test ran successfully.

Comment by sogo.mineo [ 08/Dec/22 ]

price Can I merge the branches of drp_stella and drp_stella_data to the main branch?

Comment by price [ 08/Dec/22 ]

Please go ahead. Thanks!

Comment by price [ 09/Dec/22 ]

I ran the script on /projects/HSC/PFS/weekly-20221201:

(lsst-scipipe-3.0.0) pprice@tiger2-sumire:/projects/HSC/PFS/weekly-20221201 $ python ~/PIPE2D-1110-change-filter-names.py pfsConfig*.fits --goahead
Comment by sogo.mineo [ 09/Dec/22 ]

Merged. Thanks for the review and changing the weekly test.

Comment by price [ 20/Dec/22 ]

And now I've run the script on /projects/HSC/PFS/scienceSims/scienceSims-20221201:

(lsst-scipipe-3.0.0) pprice@tiger2-sumire:/projects/HSC/PFS/scienceSims/scienceSims-20221201 $ python ~/PIPE2D-1110-change-filter-names.py pfsConfig*.fits --goahead
Generated at Sat Feb 10 16:02:55 JST 2024 using Jira 8.3.4#803005-sha1:1f96e09b3c60279a408a2ae47be3c745f571388b.