####################################### # Simulate GRID with # script: /home/s1/rkessler/snana/util/sim_SNgrid.pl # sim-job: /home/s1/rkessler/snana/bin/snlc_sim.exe # ####################################### ************************************************** Begin execution of snlc_sim.exe SNDATA_ROOT = /data/dp62.b/data/SNDATA_ROOT SNANA_DIR = /home/s1/rkessler/snana ######################################################## INIT_SNVAR: Init variables. ######################################################## HOST MACHINE = sdssdp62.fnal.gov (i386) SNDATA_ROOT = /data/dp62.b/data/SNDATA_ROOT SNANA_DIR = /home/s1/rkessler/snana (v10_31c ) sizeof(INPUTS) = 0.613 MB sizeof(GENLC) = 6.808 MB -------------------------------------------------------- Read default input from: /data/dp62.b/data/SNDATA_ROOT/models/snlc_sim.defaults -------------------------------------------------------- Read user input file: SDSS_SIMGEN_S11DM15.INPUT PROCESS COMMAND LINE ARG: GRID_FORMAT SIMLIB file : SDSS_3year.SIMLIB (start LIBID=0) Generation Version : GRID_SDSS_S11DM15 Generation source : GRID Generation model : S11DM15.v01 Number of LC to Generate: 0 OMEGA_(MATTER,LAMBDA)= 0.270, 0.730, W= -1.00 H0= 70.0 KCOR file : SDSS_Doi2010/kcor_SDSS_Bessell90_BD17.fits Observer Gen-FILTERS : ugriz Random number seed: 1 Gen-Range for RA(deg) : -360.000 to 360.000 Gen-Range for DECL(deg): -360.000 to 360.000 Gen-Range for REDSHIFT : 0.010 to 0.700 (sigma= 0.0000) Gen-Range for PEAKMJD : 0.0 to 0.0 Gen-Range for Trest : -20.0 to 90.0 days Gen-Range for SHAPEPAR : 0.6 to 1.7 Gen-Range for AV : 0.00 to 1.00 (dN/dAv = exp(-AV/0.33) Gen-Mean for RV : 2.18 Gen-sigma for RV : 0.00 , 0.00 (lower , upper ) Gen-Range for RV : 3.10 to 3.10 AB mag offsets (ugriz) : MODEL mag offsets (ugriz) : 0.00 0.00 0.00 0.00 0.00 ************************************************** SIMLIB_open(): Open conditions library for simulation. Opened /data/dp62.b/data/SNDATA_ROOT/simlib/SDSS_3year.SIMLIB SIMLIB Survey : SDSS SIMLIB Filters : ugriz SIMLIB telescope : SDSS SIMLIB created by: rkessler SIMLIB MJD-SKIP RANGE: 53705 - 53975 SIMLIB MJD-SKIP RANGE: 54060 - 54345 SIMLIB Found FLUXERR_ADD: for 'ugriz' SIMLIB Found flux-error correction map. SIMLIB PSF unit: PIXEL_SIGMA SIMLIB SKYSIG unit: ADU_PER_SQPIXEL ************************************************** init_GRIDsource for SURVEY=SDSS and GENMODEL=S11DM15.v01 : 160 GRID-BINS for LOGZ : -2.00 -1.99 -1.98 -1.97 -1.95 -1.94 -1.93 -1.92 -1.91 -1.90 -1.88 -1.87 -1.86 -1.85 -1.84 -1.83 -1.81 -1.80 -1.79 -1.78 -1.77 -1.76 -1.74 -1.73 -1.72 -1.71 -1.70 -1.69 -1.68 -1.66 -1.65 -1.64 -1.63 -1.62 -1.61 -1.59 -1.58 -1.57 -1.56 -1.55 -1.54 -1.52 -1.51 -1.50 -1.49 -1.48 -1.47 -1.45 -1.44 -1.43 -1.42 -1.41 -1.40 -1.38 -1.37 -1.36 -1.35 -1.34 -1.33 -1.32 -1.30 -1.29 -1.28 -1.27 -1.26 -1.25 -1.23 -1.22 -1.21 -1.20 -1.19 -1.18 -1.16 -1.15 -1.14 -1.13 -1.12 -1.11 -1.09 -1.08 -1.07 -1.06 -1.05 -1.04 -1.03 -1.01 -1.00 -0.99 -0.98 -0.97 -0.96 -0.94 -0.93 -0.92 -0.91 -0.90 -0.89 -0.87 -0.86 -0.85 -0.84 -0.83 -0.82 -0.80 -0.79 -0.78 -0.77 -0.76 -0.75 -0.74 -0.72 -0.71 -0.70 -0.69 -0.68 -0.67 -0.65 -0.64 -0.63 -0.62 -0.61 -0.60 -0.58 -0.57 -0.56 -0.55 -0.54 -0.53 -0.51 -0.50 -0.49 -0.48 -0.47 -0.46 -0.45 -0.43 -0.42 -0.41 -0.40 -0.39 -0.38 -0.36 -0.35 -0.34 -0.33 -0.32 -0.31 -0.29 -0.28 -0.27 -0.26 -0.25 -0.24 -0.22 -0.21 -0.20 -0.19 -0.18 -0.17 -0.15 2 GRID-BINS for AV : 0.00 1.00 1 GRID-BINS for RV : 3.10 13 GRID-BINS for DM15 : 0.60 0.69 0.78 0.88 0.97 1.06 1.15 1.24 1.33 1.43 1.52 1.61 1.70 5 GRID-BINS for FILTER : 0.00 1.00 2.00 3.00 4.00 56 GRID-BINS for TREST : -20.00 -18.00 -16.00 -14.00 -12.00 -10.00 -8.00 -6.00 -4.00 -2.00 0.00 2.00 4.00 6.00 8.00 10.00 12.00 14.00 16.00 18.00 20.00 22.00 24.00 26.00 28.00 30.00 32.00 34.00 36.00 38.00 40.00 42.00 44.00 46.00 48.00 50.00 52.00 54.00 56.00 58.00 60.00 62.00 64.00 66.00 68.00 70.00 72.00 74.00 76.00 78.00 80.00 82.00 84.00 86.00 88.00 90.00 Generate 4160 light curves, each with ugriz and 56 epochs. GRID storage (mag+err, excluding header): 4.726 MB. ============================================= RDKCOR_FITS: Read Calib/Kcor tables. Opened /data/dp62.b/data/SNDATA_ROOT/kcor/SDSS_Doi2010/kcor_SDSS_Bessell90_BD17.fits with Blocksize= 1 and ISTAT= 0 RDKCOR_FITS_HEAD Read VERSION = 3 (kcor.exe version) Read NPRIM = 2 primary refs Read NFILTERS = 16 filters Read NKCOR = 22 K-COR tables RDKCOR_FITS_ZPT RDKCOR_FITS_SNSED RDKCOR_FITS_KCOR Added Bessell-BX to rest-frame list. RDKCOR_FITS_MAGS RDKCOR_FITS_FILTERS RDKCOR_FITS_PRIMARYSED FTCLOSE ---------------------- !!! RDKCOR_SUMMARY !!! ---------------------- Primary spectrum: AB FILTER SUMMARY (RDKCOR) : internal LAM LAM rest-lam primary index name AVG RMS kcor range mag ZP ------------------------------------------------------------------ 01 SDSS-u 3594 250 0- 0 0.000 0.000 02 SDSS-g 4718 386 0- 0 0.000 0.000 03 SDSS-r 6187 343 0- 0 0.000 0.000 04 SDSS-i 7506 382 0- 0 0.000 0.000 05 SDSS-z 8918 510 0- 0 0.000 0.000 11 Bessell-U 3592 220 1000- 3900 9.724 -0.760 12 Bessell-B 4384 353 3900- 4850 9.907 0.130 13 Bessell-V 5489 357 4850- 5850 9.464 0.012 14 Bessell-R 6524 628 5850- 7050 9.166 -0.160 15 Bessell-I 8033 462 7050-12000 8.846 -0.398 16 Bessell-BX 4405 355 3615- 5595 9.907 0.140 ------------------------------------------------------------------ NFILTDEF[SURVEY,REST] = 5 6 KCOR SUMMARY ( 11 tables) 1 Found K_Uu table (11-> 1 : Bessell-U -> SDSS-u) 2 Found K_Ug table (11-> 2 : Bessell-U -> SDSS-g) 3 Found K_Bg table (12-> 2 : Bessell-B -> SDSS-g) 4 Found K_Br table (12-> 3 : Bessell-B -> SDSS-r) 5 Found K_Vr table (13-> 3 : Bessell-V -> SDSS-r) 6 Found K_Vi table (13-> 4 : Bessell-V -> SDSS-i) 7 Found K_Rr table (14-> 3 : Bessell-R -> SDSS-r) 8 Found K_Ri table (14-> 4 : Bessell-R -> SDSS-i) 9 Found K_Rz table (14-> 5 : Bessell-R -> SDSS-z) 10 Found K_Ii table (15-> 4 : Bessell-I -> SDSS-i) 11 Found K_Iz table (15-> 5 : Bessell-I -> SDSS-z) Rest-frame filters : UBVRI (ifiltdef= 11 12 13 14 15 16) K-CORR TABLE uses 23.85% of 8000000 available bins ( 32.0 MB) AVWARP TABLE uses 38.16% of 1000000 available bins ( 4.0 MB) LCMAG TABLE uses 19.08% of 2000000 available bins ( 8.0 MB) MWXT TABLE uses 15.90% of 2000000 available bins ( 8.0 MB) MWXT TABLE params: RV= 3.10 and OPT_MWCOLORLAW= 94 Read 921 Lambda bins, 2100. - 11300. binSize= 10.0 Read 106 Trest bins, -20.00 - 85.00 binSize=1.000 Read 24 Redshift bins, 0.00 - 0.46 binSize=0.020 Read 25 AVwarp bins, -6.00 - 6.00 binSize=0.500 RDKCOR: Done. ---------------------------------------------- Zeropoint offsets (ZPOFF) from : KCOR file Will apply net MAGOBS_SHIFT_ZP(u) = -0.0679 Will apply net MAGOBS_SHIFT_ZP(g) = 0.0203 Will apply net MAGOBS_SHIFT_ZP(r) = 0.0049 Will apply net MAGOBS_SHIFT_ZP(i) = 0.0178 Will apply net MAGOBS_SHIFT_ZP(z) = 0.0102 ************************************************** init_genmag_S11DM15 : Initialize S11DM15.v01 Read primary ref 'AB' with 921 lambda bins Defined Lambda AB AB Filter Range(step) mag ZP ---------------------------------------------------------------------- id=01 ( SDSS-u) ( 3594.3 + 0) 3000- 7880( 10) 0.000 12.185 id=02 ( SDSS-g) ( 4717.6 + 0) 3610- 5630( 10) 0.000 14.341 id=03 ( SDSS-r) ( 6186.8 + 0) 5370- 7050( 10) 0.000 14.230 id=04 ( SDSS-i) ( 7506.2 + 0) 6590- 8390( 10) 0.000 13.886 id=05 ( SDSS-z) ( 8918.3 + 0) 7700-11160( 10) 0.000 11.912 ---------------------------------------------------------------------- Read S11DM15 model parameters from /data/dp62.b/data/SNDATA_ROOT/models/S11DM15/S11DM15.v01/S11DM15.INFO S11DM15.INFO RESTLAMBDA_RANGE: 1000 - 25000 A X = DM15 - 1.10 a(lam) = 1.248000*LAM^0 -0.000105*LAM^1 0.000000*LAM^2 b_PAR = 0.633000 TAUPOLY_STRETCH = 3.455000e+00 1.372000e+01 -3.601000e+00 9.460000e-01 Read S11DM15 SED from : /data/dp62.b/data/SNDATA_ROOT/models/S11DM15/S11DM15.v01/SED.DAT LAM(MIN,MAX,STEP)= 1000, 25000, 10 DAY(MIN,MAX,STEP)= -20, 85, 1.0 Max B-flux at Trest = 0.000000 days stretch(dm15=0.6) = 1.416 stretch(dm15=0.8) = 1.189 stretch(dm15=1.0) = 1.033 stretch(dm15=1.2) = 0.916 stretch(dm15=1.4) = 0.824 stretch(dm15=1.6) = 0.748 Host Transmission with RV=3.1 and AV=1.0 : Host-Trans(Lam= 1000.0) = 0.00803 Host-Trans(Lam= 2000.0) = 0.07294 Host-Trans(Lam= 3000.0) = 0.18739 Host-Trans(Lam= 4000.0) = 0.25952 Host-Trans(Lam= 5000.0) = 0.35571 Host-Trans(Lam= 6000.0) = 0.43382 Host-Trans(Lam= 7000.0) = 0.50111 Host-Trans(Lam= 8000.0) = 0.57677 Host-Trans(Lam= 9000.0) = 0.64339 Host-Trans(Lam=10000.0) = 0.68929 ************************************************** init_modelSmear: init intrinsic SN smearing with model= Smear-mode: interpolate from central filter wavelengths. ************************************************** Fill comments for README doc-file ************************************************** clr_VERSION LCMERGE Version GRID_SDSS_S11DM15 does not exist. SIM Version GRID_SDSS_S11DM15 exists. PHOTOMETRY Version GRID_SDSS_S11DM15 does not exist. Removing version GRID_SDSS_S11DM15 files ... execute command: rm -rf /data/dp62.b/data/SNDATA_ROOT/SIM/GRID_SDSS_S11DM15/ Done. Opened GRIDGEN fits file (istat=0): /data/dp62.b/data/SNDATA_ROOT/SIM/GRID_SDSS_S11DM15/GRID_SDSS_S11DM15.GRID create GRIDGEN header keys. create GRIDGEN info-table for each parameter : create GRIDGEN table of values for : LOGZ AV RV DM15 FILTER TREST create GRIDGEN table of pointers to 4160 light curves create GRIDGEN table for model MAGS and ERRORS ************************************************** Write Filter Transmision files mkdir /data/dp62.b/data/SNDATA_ROOT/SIM/GRID_SDSS_S11DM15/GRID_SDSS_S11DM15.FILTERS ************************************************** Begin Generating Lightcurves. Malloc S11DM15 TABLES for FLUXSCALE & XTHOST Finished writing CID = 1 (NEP=280) Finished writing CID = 100 (NEP=280) Finished writing CID = 200 (NEP=280) Finished writing CID = 300 (NEP=280) Finished writing CID = 400 (NEP=280) Finished writing CID = 500 (NEP=280) Finished writing CID = 600 (NEP=280) Finished writing CID = 700 (NEP=280) Finished writing CID = 800 (NEP=280) Finished writing CID = 900 (NEP=280) Finished writing CID = 1000 (NEP=280) Finished writing CID = 1100 (NEP=280) Finished writing CID = 1200 (NEP=280) Finished writing CID = 1300 (NEP=280) Finished writing CID = 1400 (NEP=280) Finished writing CID = 1500 (NEP=280) Finished writing CID = 1600 (NEP=280) Finished writing CID = 1700 (NEP=280) Finished writing CID = 1800 (NEP=280) Finished writing CID = 1900 (NEP=280) Finished writing CID = 2000 (NEP=280) Finished writing CID = 2100 (NEP=280) Finished writing CID = 2200 (NEP=280) Finished writing CID = 2300 (NEP=280) Finished writing CID = 2400 (NEP=280) Finished writing CID = 2500 (NEP=280) Finished writing CID = 2600 (NEP=280) Finished writing CID = 2700 (NEP=280) Finished writing CID = 2800 (NEP=280) Finished writing CID = 2900 (NEP=280) Finished writing CID = 3000 (NEP=280) Finished writing CID = 3100 (NEP=280) Finished writing CID = 3200 (NEP=280) Finished writing CID = 3300 (NEP=280) Finished writing CID = 3400 (NEP=280) Finished writing CID = 3500 (NEP=280) Finished writing CID = 3600 (NEP=280) Finished writing CID = 3700 (NEP=280) Finished writing CID = 3800 (NEP=280) Finished writing CID = 3900 (NEP=280) Finished writing CID = 4000 (NEP=280) Finished writing CID = 4100 (NEP=280) ************************************************** Done generating 4160 SN lightcurves from GRID source. (4160 lightcurves requested => 4160 were written) ************************************************** Fill comments for README doc-file ************************************************** DUMP README CONTENTS TO SCREEN USER COMMENT: BLANK BRIEF_DESCRIPTION: simulate GRID SNe with GENMODEL = S11DM15.v01 HOST MACHINE: sdssdp62.fnal.gov USERNAME: rkessler SNDATA_ROOT: /data/dp62.b/data/SNDATA_ROOT SNANA_DIR: /home/s1/rkessler/snana Current Dir: /data/dp62.b/data/dp47.b/data/www/html/sdsssn_data/Papers/DataRelease/WebPage/SDSS_dataRelease-snana/SNgrid FULL_DESCRIPTION: SIMLIB filename : SDSS_3year.SIMLIB (start at LIBID=0) SIMLIB SURVEY : SDSS (TELESCOPE=SDSS) SIMLIB UNITS : PIXEL_SIGMA for PSF, ADU_PER_SQPIXEL for SKYSIG NEWMJD_DIF : 10.08 minutes Generation VERSION: GRID_SDSS_S11DM15 Generation source : GRID Generation FILTERS: ugriz KCOR lookup tables: SDSS_Doi2010/kcor_SDSS_Bessell90_BD17.fits Flux-smearing is ON Zeropt-smearing is ON Host-galaxy noise is OFF MilkyWay extinction is ON Color law: CCM89+ODonell94 (OPT_MWCOLORLAW=94) E(B-V): FILE value (SIMLIB or data header) (OPT_MWEBV=1) sigma(MWEBV) = 0.16*MWEBV + 0.00 shift(MWEBV) = 0.000 mag GENERATION RANGES: Generate Redshift : 0.010 to 0.700 using POWERLAW distribution. REDSHIFT_FINAL is ZGEN smeared by : 0.00000 ugriz AB offsets : -0.068 0.020 0.005 0.018 0.010 ugriz MODEL offsets : 0.000 0.000 0.000 0.000 0.000 ugriz exposure times: 1.0000 1.0000 1.0000 1.0000 1.0000 RA : -360.00 to 360.00 deg PEAKMJD : 0.0 to 0.0 sigma-smear(SEARCH_PEAKMJD) : 0.00 days Trest : -20.00 to 90.00 days RISETIME-SHIFT(days) SIGMA(lo,hi) : 0.0 , 0.0 (Mean= 0.0) FALLTIME-SHIFT(days) SIGMA(lo,hi) : 0.0 , 0.0 (Mean= 0.0) Shape-par(DM15): Peak=0.000 SIG-+= 0.000,0.000 BND=0.60,1.70 Z-dependent SN-parameter shifts: None. Intrinsic MAG-smearing models (sigma clip -3.0 to 3.0) : Model 1: Coherent MAG-smearing (GENMAG_SMEAR) : 0.000 Model 2: passband MAG-smearing is OFF Model 3: model-smear is OFF Model 4: intrinsic scatter matrix is OFF Model 5: GENMAG_SMEAR_USRFUN is OFF AVWARP_OVERFLOWS: NONE. Cosmology Parameters: H0 = 70.00 km/s per MPc Omega_{M,L} = 0.27, 0.73 w = -1.00 --------------------------------------------------- Software-Pipeline Search Efficiency (MINOBS=2) from Spectroscopic Efficiency : APPLY_SEARCHEFF_OPT: 0 => Do NOT Apply trigger efficiency SNTYPE(Ia) = 1(SPEC) and 101(PHOT) --------------------------------------------------- HOSTLIB Summary: None. Fudges on SIMLIB Seeing Conditions: Fudge-factor for SIMLIB NOISE(SKY,READ) : 1.00, 1.00 Fudge-factor for SIMLIB PSF : 1.00 Fudge-factor for SIMLIB ZPT : 0.00 mag Fudge-factor for FLUX-ERROR : 1.00 ============ END OF SIMULATION SUMMARY ============== Random Number Sync: RANDOM SEED: 1 (RANLIST_START_GENSMEAR: 1) Generation Statistics: Generated 4160 simulated light curves. Wrote 4160 simulated light curves to SNDATA files. Spectroscopic-type: 4160 -> 4160 (before -> after cuts) Photometric-type: 0 -> 0 (before -> after cuts) Rejection Statistics: 0 rejected by GEN-RANGE 0 rejected by SEARCH-TRIGGER 0 rejected by CUTWIN-SELECTION SEARCH+CUTS Efficiency: 1.0000 +- 0.0002 Number of SNe per season AFTER CUTS : 0 +- 0 ===== END OF README FILE ====== AUXILIARY FILES: /data/dp62.b/data/SNDATA_ROOT/SIM/GRID_SDSS_S11DM15/GRID_SDSS_S11DM15.LIST /data/dp62.b/data/SNDATA_ROOT/SIM/GRID_SDSS_S11DM15/GRID_SDSS_S11DM15.README /data/dp62.b/data/SNDATA_ROOT/SIM/GRID_SDSS_S11DM15/GRID_SDSS_S11DM15.FILTERS /data/dp62.b/data/SNDATA_ROOT/SIM/GRID_SDSS_S11DM15/GRID_SDSS_S11DM15.GRID ************************************************** end_GRIDfile Total number of generated points: 1181440 Actual PTR_GRIDGEN_LC[last ilc= 4160] = 1181157 Expected PTR_GRIDGEN_LC[last ilc] = 1181157 DONE with snlc_sim.