Hsc Survey Data

2. Set Credentials

For Internal Data Release (DR4/S23B):

export SSP_IDR_USR="your_username"
export SSP_IDR_PWD="your_password"

For Public Data Release (PDR3):

export SSP_PDR_USR="your_username"
export SSP_PDR_PWD="your_password"

Tool 1: SQL Catalog Query (hscSspQuery3.py)

Location in repo: dr4/catalogQuery/hscSspQuery3.py

Source: https://hsc-gitlab.mtk.nao.ac.jp/ssp-software/data-access-tools/blob/master/dr4/catalogQuery/hscSspQuery3.py

Basic Usage

python hscSspQuery3.py \
    --user $SSP_IDR_USR \
    --release-version dr4 \
    --format csv \
    query.sql > output.csv

Command-Line Options

API Endpoints (Default)

• Internal: https://hscdata.mtk.nao.ac.jp/datasearch/api/catalog_jobs/
• Public: https://hsc-release.mtk.nao.ac.jp/datasearch/api/catalog_jobs/

SQL Query Examples

Cone search around coordinates:

SELECT object_id, ra, dec, i_cmodel_flux, i_cmodel_mag
FROM s23b_wide.forced
WHERE coneSearch(coord, 136.47, -0.05, 60)
AND i_cmodel_flux > 0
ORDER BY i_cmodel_flux DESC

Box search:

SELECT *
FROM s23b_wide.forced
WHERE boxSearch(coord, 136.0, 137.0, -0.5, 0.5)
LIMIT 1000

Join with photoz:

SELECT f.object_id, f.ra, f.dec, f.i_cmodel_mag, p.photoz_best
FROM s23b_wide.forced AS f
JOIN s23b_wide.photoz_mizuki AS p ON f.object_id = p.object_id
WHERE coneSearch(f.coord, 150.0, 2.0, 120)

Python Module Usage

You can also import the script as a module:

import sys
sys.path.insert(0, '/path/to/data-access-tools/dr4/catalogQuery')
import hscSspQuery3

# Set arguments programmatically
hscSspQuery3.args = type('Args', (), {
    'user': os.environ['SSP_IDR_USR'],
    'release_version': 'dr4',
    'out_format': 'csv',
    'delete_job': True,
    'nomail': True,
    'skip_syntax_check': False,
    'api_url': 'https://hscdata.mtk.nao.ac.jp/datasearch/api/catalog_jobs/',
    'password_env': 'SSP_IDR_PWD'
})()

# Read SQL
with open('query.sql', 'r') as f:
    sql = f.read()

# Get credentials
credential = {
    'account_name': hscSspQuery3.args.user,
    'password': hscSspQuery3.getPassword()
}

# Submit job
job = hscSspQuery3.submitJob(credential, sql, 'csv')
print(f"Job ID: {job['id']}")

# Wait for completion
hscSspQuery3.blockUntilJobFinishes(credential, job['id'])

# Download to file
with open('output.csv', 'wb') as f:
    hscSspQuery3.download(credential, job['id'], f)

# Clean up
hscSspQuery3.deleteJob(credential, job['id'])

Tool 2: Cutout Images (downloadCutout.py)

Location in repo: dr4/downloadCutout/downloadCutout.py

Basic Usage

# Download single cutout
python downloadCutout.py \
    --ra 136.471165 \
    --dec -0.046470 \
    --sw 10arcsec \
    --sh 10arcsec \
    --filter HSC-I \
    --rerun s23b_wide \
    --type coadd/bg \
    --image true \
    --mask true \
    --variance true \
    --user $SSP_IDR_USR

Batch Download via List File

Create a coordinate list file coords.txt:

#? ra dec sw sh filter
136.471165 -0.046470 10arcsec 10arcsec HSC-G
136.471165 -0.046470 10arcsec 10arcsec HSC-R
136.471165 -0.046470 10arcsec 10arcsec HSC-I
136.471165 -0.046470 10arcsec 10arcsec HSC-Z
136.471165 -0.046470 10arcsec 10arcsec HSC-Y

Then run:

python downloadCutout.py \
    --list coords.txt \
    --rerun s23b_wide \
    --type coadd/bg \
    --image true \
    --mask true \
    --variance true \
    --user $SSP_IDR_USR

Key Parameters

Tool 3: PSF Download (downloadPsf.py)

Location in repo: dr4/downloadPsf/downloadPsf.py

⚠️ Important: S23B data requires the psf/9 endpoint, while older releases use psf/8. Edit the script's api_url if needed.

Basic Usage

python downloadPsf.py \
    --ra 136.471165 \
    --dec -0.046470 \
    --filter HSC-I \
    --rerun s23b_wide \
    --type coadd \
    --centered true \
    --user $SSP_IDR_USR

Batch PSF Download

Create list file psf_coords.txt:

#? ra dec filter
136.471165 -0.046470 HSC-G
136.471165 -0.046470 HSC-R
136.471165 -0.046470 HSC-I
136.471165 -0.046470 HSC-Z
136.471165 -0.046470 HSC-Y

Run:

python downloadPsf.py \
    --list psf_coords.txt \
    --rerun s23b_wide \
    --type coadd \
    --centered true \
    --user $SSP_IDR_USR

Database Schema Reference

Main Catalog Tables (S23B/DR4)

Note: The forced table is split into 6 parts due to PostgreSQL column limits. To get both CModel and PSF fluxes, JOIN on object_id:

SELECT f.object_id, f.ra, f.dec, 
       f.i_cmodel_flux, f2.i_psfflux_flux
FROM s23b_wide.forced AS f
LEFT JOIN s23b_wide.forced2 AS f2 
    ON f.object_id = f2.object_id
WHERE coneSearch(f.coord, 136.47, -0.05, 60)

Useful SQL Functions

Flux Units

• Fluxes: nano-Janskys (nJy)
• Magnitudes: AB mag (stored as _mag columns)
• Coordinates: degrees (J2000)

Complete Workflow Example

#!/bin/bash

# Set credentials
export SSP_IDR_USR="your_username"
export SSP_IDR_PWD="your_password"

# Target coordinates
RA=136.471165
DEC=-0.046470

# 1. Query catalog
cat > query.sql << EOF
SELECT object_id, ra, dec, 
    g_cmodel_mag, r_cmodel_mag, i_cmodel_mag, z_cmodel_mag, y_cmodel_mag
FROM s23b_wide.forced
WHERE coneSearch(coord, $RA, $DEC, 10)
AND i_cmodel_flux > 0
ORDER BY i_cmodel_flux DESC
LIMIT 10
EOF

python hscSspQuery3.py \
    --user $SSP_IDR_USR \
    --release-version dr4 \
    --format csv \
    query.sql > catalog.csv

# 2. Download cutouts
cat > cutout_list.txt << EOF
#? ra dec sw sh filter type
$RA $DEC 10arcsec 10arcsec HSC-G coadd/bg
$RA $DEC 10arcsec 10arcsec HSC-R coadd/bg
$RA $DEC 10arcsec 10arcsec HSC-I coadd/bg
$RA $DEC 10arcsec 10arcsec HSC-Z coadd/bg
$RA $DEC 10arcsec 10arcsec HSC-Y coadd/bg
EOF

python downloadCutout.py \
    --list cutout_list.txt \
    --rerun s23b_wide \
    --image true --mask true --variance true \
    --user $SSP_IDR_USR

# 3. Download PSFs
cat > psf_list.txt << EOF
#? ra dec filter
$RA $DEC HSC-G
$RA $DEC HSC-R
$RA $DEC HSC-I
$RA $DEC HSC-Z
$RA $DEC HSC-Y
EOF

python downloadPsf.py \
    --list psf_list.txt \
    --rerun s23b_wide \
    --type coadd \
    --centered true \
    --user $SSP_IDR_USR

Troubleshooting

PSF Download 404 Error

If PSF download fails with 404 for S23B data, edit downloadPsf.py and change:

api_url = "https://hscdata.mtk.nao.ac.jp/psf/8"