Dev

pyproject.toml

@@ -0,0 +1,30 @@
+[build-system]
+requires = ["setuptools>=61.0.0"]
+build-backend = "setuptools.build_meta"
+
+[project]
+name = "sn_tools"
+version = "1.0.2"
+requires-python = ">=3.5"
+dependencies = ['astropy_healpix>=0.6','descartes>=1.0','sncosmo>=1.0','aplpy']
+authors = [
+  {name = "Philippe Gris", email = "philippe.gris@clermont.in2p3.fr"}
+]
+maintainers = [
+  {name = "Philippe Gris", email = "philippe.gris@clermont.in2p3.fr"}
+]
+description = "Set of tools to run sn_pipe"
+readme = "README.md"
+license = {file = "LICENSE"}
+classifiers = ["Programming Language :: Python",]
+
+[tool.setuptools.packages.find]
+where = [".","sn_script_input"]
+include = ['sn_tools', 'sn_analysis','sn_script_input']
+
+[tool.setuptools.package-data]
+"*" = ["*.txt"]
+
+[project.urls]
+Homepage = "http://github.com/lsstdesc/sn_tools"
+Repository = "https://github.com/lsstdesc/sn_tools"

requirements.txt

@@ -0,0 +1 @@
+###### Requirements without Version Specifiers ######

setup.py → setup_deprecated.py

@@ -1,7 +1,7 @@
 from setuptools import setup
 
 # get the version here
-pkg_vars  = {}
+pkg_vars = {}
 
 with open("version.py") as fp:
     exec(fp.read(), pkg_vars)
@@ -14,7 +14,10 @@
     author='Philippe Gris',
     author_email='philippe.gris@clermont.in2p3.fr',
     license='BSD',
-    packages= ['sn_tools'],
+    packages= ['sn_tools','sn_script_input'],
+    # All files from folder sn_script_input
+    package_data={'sn_script_input': ['*.txt']},
     python_requires='>=3.5',
-    zip_safe=False
+    zip_safe=False,
+    install_requires=['astropy_healpix>=1.0','descartes>=1.0','sncosmo>=1.0'],
 )

sn_script_input/__init__.py

@@ -0,0 +1 @@
+from .version import __version__

sn_script_input/config_cadence.txt

@@ -0,0 +1,115 @@
+## RA and Dec
+RA 0.0 float # RA
+Dec 0.0 float #Dec
+
+## fieldName (to be set as 'note' in the file)
+field DD:COSMOS str # field
+
+## seasons
+seasons 1 str # seasons
+seasonLength 180 str #season_length (days)
+MJDmin 60321. float #MJD min for observations - 1/12/2024
+
+## bands to consider
+bands grizy str #bands to consider
+
+## cadence[day] per band
+cadence g 1. float # cadence g band
+cadence r 1. float # cadence r band
+cadence i 1. float # cadence i band
+cadence z 1. float # cadence z band
+cadence y 1. float # cadence y band
+
+## obs file
+obsCondFile m5_night.csv str # Observing condition file
+obsFromSimu 0 int # to take obs. cond. from obsCondFile
+
+## fiveSigmaDepth per band
+fiveSigmaDepth u 23.48 str #u 5-sigma depth values
+fiveSigmaDepth g 24.26 str #g 5-sigma depth values
+fiveSigmaDepth r 23.87 str #r 5-sigma depth values
+fiveSigmaDepth i 23.46 str #i 5-sigma depth values
+fiveSigmaDepth z 22.76 str #z 5-sigma depth values
+fiveSigmaDepth y 22.02 str #y 5-sigma depth values
+
+## number of visits per band
+Nvisits u 10 int #Nvisits u
+Nvisits g 10 int #Nvisits g
+Nvisits r 20 int #Nvisits r
+Nvisits i 20  int #Nvisits i
+Nvisits z 26 int #Nvisits z
+Nvisits y 20 int #Nvisits y
+
+## Single visit exposure time
+visitExposureTime u  30. float #single visit exp time u-band
+visitExposureTime g  30. float #single visit exp time g-band
+visitExposureTime r  30. float #single visit exp time r-band
+visitExposureTime i  30. float #single visit exp time i-band
+visitExposureTime z  30. float #single visit exp time z-band
+visitExposureTime y  30. float #single visit exp time y-band
+
+## numExposures
+numExposures u  2 int #single visit exp time u-band
+numExposures g  2 int #single visit exp time g-band
+numExposures r  2 int #single visit exp time r-band
+numExposures i  2 int #single visit exp time i-band
+numExposures z  2 int #single visit exp time z-band
+numExposures y  2 int #single visit exp time y-band
+
+##shift between visits per band per day
+shiftDays 0.0069444 float #in days: 10./(24.*60.)
+
+##seeing
+seeingFwhmEff u 0.90 float # seeingFwhmEff u-band
+seeingFwhmEff g 0.87 float # seeingFwhmEff g-band
+seeingFwhmEff r 0.83 float # seeingFwhmEff r-band
+seeingFwhmEff i 0.80  float # seeingFwhmEff i-band
+seeingFwhmEff z 0.78 float # seeingFwhmEff z-band
+seeingFwhmEff y 0.76 float # seeingFwhmEff y-band
+
+seeingFwhmGeom u 0.90 float # seeingFwhmGeom u-band
+seeingFwhmGeom g 0.87 float # seeingFwhmGeom g-band
+seeingFwhmGeom r 0.83 float # seeingFwhmGeom r-band
+seeingFwhmGeom i 0.80  float # seeingFwhmGeom i-band
+seeingFwhmGeom z 0.78 float # seeingFwhmGeom z-band
+seeingFwhmGeom y 0.76 float # seeingFwhmGeom y-band
+
+##airmass
+airmass u 1.2 float #airmass u-band
+airmass g 1.2 float #airmass g-band
+airmass r 1.2 float #airmass r-band
+airmass i 1.2 float #airmass i-band
+airmass z 1.2 float #airmass z-band
+airmass y 1.2 float #airmass y-band
+
+##sky
+sky u 1.2 float #sky g-band
+sky g 1.2 float #sky g-band
+sky r 1.2 float #sky r-band
+sky i 1.2 float #sky i-band
+sky z 1.2 float #sky z-band
+sky y 1.2 float #sky y-band
+
+
+## moonPhase
+moonPhase u 0.5 float #moonPhase u-band
+moonPhase g 0.5 float #moonPhase g-band
+moonPhase r 0.5 float #moonPhase r-band
+moonPhase i 0.5 float #moonPhase i-band
+moonPhase z 0.5 float #moonPhase z-band
+moonPhase y 0.5 float #moonPhase y-band
+moonPhaseu 20. float #max moon phase for u-obs
+moonswapFilter y str #filter to be replaced by when no moon
+mooncompensate 0 int # to compensate for moon effects
+
+## rotTelPos
+rotTelPos 0 int #tel rotation dithering
+
+## m5file from realistic simus
+# m5File NoData str #m5 file from realistic simus
+# healpixID -1 int #healpixID for m5 values
+
+## Save data infos
+saveData 0 int #flag to save data
+outDir ../Fake_Observations str #output directory
+outName Fakes str #output file name

sn_script_input/config_process.txt

@@ -0,0 +1,35 @@
+## generic metric parameters
+dbName alt_sched str #db name
+dbExtens npy str #db extension
+dbDir /sps/lsst/cadence/LSST_SN_CADENCE/cadence_db str #db dir 
+templateDir /sps/lsst/data/dev/pgris/Templates_final_new str #template dir
+nproc 8 int #number of proc
+nproc_pixels 8 int #number of proc to process pixels
+remove_dithering 0 int #remove dithering for DDF
+simuType 0 int #flag for new simulations
+saveData 0 int #flag to dump data on disk
+dirRefs reference_files str #dir of reference files for the metric
+dirFake input/Fake_cadence str #dir of fake files for the metric
+pixelmap_dir None str #dir where to find pixel maps
+nclusters 0 int #number of clusters in data (DD only)
+radius 4.0 float #radius around clusters (DD and Fakes)
+pixelList None str #list of healpixIds to process
+RAmin 0.0 float #RA min for obs area - for WDF only
+RAmax 360. float #RA max for obs area - for WDF only
+Decmin -90.0 float #Dec min for obs area - for WDF only
+Decmax 40.0 float #Dec max for obs area - for WDF only
+npixels -1 int #number of pixels to process
+FoV 9.6 float #telescope field of view
+# VRO_FP circular str #VRO Focal Plane
+# project_FP gnomonic str #type of projection in FP
+telrot 0 int #telescope rotation angle
+fp_level ccd str #fp level(raft,ccd,sensor)
+display 0 int #display results
+fieldType DD str #field type DD or WFD
+nside 64 int #healpix nside
+DD_list COSMOS,CDFS,EDFS,ELAISS1,XMM-LSS str #list of DDFs
+fieldName COSMOS str #fieldName - for DD only
+ebvofMW_pixel -1.0 float #E(B-V) for pixel processing
+lookup_ddf input/simulation/lookup_ddf.csv str #DDF lookup table
+noteCol target_name str #obs col to grab field name
+code new str#code to use (old/new)

sn_script_input/config_sel.txt

@@ -0,0 +1,2 @@
+selection_params input/lc_selection/light_curve_selection.csv str #csv file of LC selection
+nproc_sel 8 int # number of procs for multiprocessing

sn_script_input/config_test_simtofit_obs.txt

@@ -0,0 +1,112 @@
+## config file for observation
+config_obs input/sim_to_simu_test/combi_obs.csv str #config file for obs
+
+## fieldName (to be set as 'note' in the file)
+field DD:COSMOS str # field
+
+## config file for simu+fit scheme
+config_simu input/sim_to_simu_test/combi_simufit.csv str #config file for simu
+
+## RA and Dec
+RA 0.0 float # RA
+Dec 0.0 float #Dec
+
+## seasons
+seasons 1 str # seasons
+seasonLength 180 str #season_length (days)
+MJDmin 60321. str #MJD min for observations - 1/12/2024
+
+## bands to consider
+bands grizy str #bands to consider
+
+## obs file
+obsCondFile m5_night.csv str # Observing condition file
+obsFromSimu 0 int # to take obs. cond. from obsCondFile
+
+## cadence[day] per band - see config_obs
+## cadence g 1. float # cadence g band
+## cadence r 1. float # cadence r band
+## cadence i 1. float # cadence i band
+## cadence z 1. float # cadence z band
+## cadence y 1. float # cadence y band
+
+## m5 per band
+## m5 per band
+m5 u 23.48 str #u 5-sigma depth values
+m5 g 24.26 str #g 5-sigma depth values
+m5 r 23.87 str #r 5-sigma depth values
+m5 i 23.46 str #i 5-sigma depth values
+m5 z 22.76 str #z 5-sigma depth values
+m5 y 22.02 str #y 5-sigma depth values
+
+## number of visits per band - see config_obs
+#Nvisits g 10,2 str #Nvisits g
+#Nvisits r  20,9 str #Nvisits r
+#Nvisits i 20,37  str #Nvisits i
+#Nvisits z 26,52 str #Nvisits z
+#Nvisits y 20,21 str #Nvisits y
+
+## Single visit exposure time
+ExposureTime u  30. float #single visit exp time g-band
+ExposureTime g  30. float #single visit exp time g-band
+ExposureTime r  30. float #single visit exp time r-band
+ExposureTime i  30. float #single visit exp time i-band
+ExposureTime z  30. float #single visit exp time z-band
+ExposureTime y  30. float #single visit exp time y-band
+
+##shift between visits per band per day
+shiftDays 0.0069444 float #in days: 10./(24.*60.)
+
+##seeing
+seeingEff u 0.90 float # seeingEff u-band
+seeingEff g 0.87 float # seeingEff g-band
+seeingEff r 0.83 float # seeingEff r-band
+seeingEff i 0.80  float # seeingEff i-band
+seeingEff z 0.78 float # seeingEff z-band
+seeingEff y 0.76 float # seeingEff y-band
+
+seeingGeom g 0.90 float # seeingGeom u-band
+seeingGeom g 0.87 float # seeingGeom g-band
+seeingGeom r 0.83 float # seeingGeom r-band
+seeingGeom i 0.80  float # seeingGeom i-band
+seeingGeom z 0.78 float # seeingGeom z-band
+seeingGeom y 0.76 float # seeingGeom y-band
+
+##airmass
+airmass u 1.2 float #airmass u-band
+airmass g 1.2 float #airmass g-band
+airmass r 1.2 float #airmass r-band
+airmass i 1.2 float #airmass i-band
+airmass z 1.2 float #airmass z-band
+airmass y 1.2 float #airmass y-band
+
+##sky
+sky u 1.2 float #sky u-band
+sky g 1.2 float #sky g-band
+sky r 1.2 float #sky r-band
+sky i 1.2 float #sky i-band
+sky z 1.2 float #sky z-band
+sky y 1.2 float #sky y-band
+
+
+## moonPhase
+moonphase g 0.5 float #moonphase g-band
+moonphase r 0.5 float #moonphase r-band
+moonphase i 0.5 float #moonphase i-band
+moonphase z 0.5 float #moonphase z-band
+moonphase y 0.5 float #moonphase y-band
+moonPhaseu 20. float #max moon phase for u-obs
+moonswapFilter z str #filter to be replaced by when no moon
+mooncompensate 0 int # to compensate for moon effects
+
+## rotTelPos
+rotTelPos 0 int #tel rotation dithering
+
+## m5file from realistic simus
+# m5File NoData str #m5 file from realistic simus
+# healpixID -1 int #healpixID for m5 values
+
+## Save data infos
+saveData 0 int #flag to save data
+outDir ../Fake_Observations str #output directory
+outName Fakes str #output file name

sn_script_input/config_test_simtofit_simu.txt

@@ -0,0 +1,24 @@
+#x1
+SN x1 type unique str #SN x1 type
+SN x1 min -2.0 float #SN x1 min
+SN x1 max 0.2 float #SN x1 max
+SN x1 step 0.1 float #SN x1 step
+# color
+SN color type unique str #SN color type
+SN color min 0.2 float #SN color min
+SN color max 1.0 float #SN color max
+SN color step 0.1 float #SN color step
+# redshift
+SN z type uniform str #SN z type
+SN z min 0.2 float #SN z min
+SN z max 1.1 float #SN z max
+SN z step 0.01 float #SN z step
+# daymax
+SN daymax type unique str #SN daymax type
+SN daymax step 1. float #SN daymax step
+# other parameters
+SN NSNfactor 1 int #NSN*factor for simulation
+SN NSNabsolute 1 int #absolute number of SN to produce
+#Multiprocessing
+MultiprocessingSimu nproc 1 int #multiprocessing number of procs
+MultiprocessingFit nproc 1 int #multiprocessing number of procs

sn_script_input/version.py

@@ -0,0 +1 @@
+__version__='v1.0.0'

sn_tools/__init__.py

@@ -1 +1,18 @@
 from .version import __version__
+
+"""
+import os
+
+# requested infos for throughputs
+throughputs_dir = os.path.join(os.getenv("PWD"), "throughputs")
+if not os.path.isdir(throughputs_dir):
+    cmd = 'git clone https://github.com/lsst/throughputs'
+    os.system(cmd)
+
+os.environ['LSST_THROUGHPUTS_BASELINE'] = '{}/{}'.format(
+    throughputs_dir, 'baseline')
+os.environ['THROUGHPUTS_DIR'] = throughputs_dir
+
+print('Reading throughputs from',
+      os.environ['LSST_THROUGHPUTS_BASELINE'], os.environ['THROUGHPUTS_DIR'])
+"""