// sntools_genSmear.h // // Mar 30 2018: MXLAM_GENSMEAR_SALT2 --> 4000 (was 1000) // Oct 18 2019: add COVSED model #define MASK_GENSMEAR_APPLY 1 // apply genSmear, old or new #define MASK_GENSMEAR_NEW 2 // re-compute genSmear void init_genSmear_FLAGS(int MSKOPT, char *SCALE_STRING); int istat_genSmear(void) ; void init_genSmear_SCALE(char *SCALE_STRING); void init_genSmear_USRFUN(int NPAR, double *parList, double *LAMRANGE ) ; void init_genSmear_SALT2(char *version, char *dispFile, double SIGCOH, double *GENRANGE_REDSHIFT); void init_genSmear_randoms(int NRANGauss, int NRANFlat) ; void read_genSmear_SALT2disp(char *smearFile, double COLOR_DISP_MAX) ; void read_genSmear_SALT2INFO(char *versionSALT2, GRIDMAP1D_DEF *SIGCOH_LAM, double *COLOR_DISP_MAX ) ; void parse_SIGCOH_SALT2(char *KEYNAME, char *KEYARG, GRIDMAP1D_DEF *SIGCOH_LAM); void getFileName_SALT2colorDisp(char *fileName) ; // added Jan 2017 void init_genSmear_Chotard11(int OPT_farUV) ; void init_genSmear_VCR(char *VCR_version, int index_SNmodel); void init_genSmear_CCM89(double *LAMRANGE) ; void init_genSmear_COH(char *stringArg) ; void init_genSmear_biModalUV(void) ; void init_genSmear_OIR(char *OIR_version); void init_genSmear_COVSED(char *COVSED_version, int OPTMASK); void parse_COVSED_INFO_FILE(void); void readFits_genSmear_COVSED(char *covFileName); void init_genSmear_phaseCor(double magSmear, double exptau); void get_genSmear_phaseCor(int CID, double phase, double *magSmear ); void check_genSmear_phaseCor(void); void init_genSmear_private(void) ; int nval_genSmear_override(char *inputKey, char *parName); void store_genSmear_override(char *parName, int NVAL, double *tmpList); int exec_genSmear_override(int ipar, char *keyName, double *val) ; void read_OIR_INFO(void) ; // read misc. info/parameters void sort_OIR_BANDS(void); void prep_OIR_COVAR(void) ; void read_VCR_VSI(void) ; // read v_Si distribution void read_VCR_INFO(void) ; // read misc. info/parameters void prep_VCR_forSALT2(void); // inits specific to SALT2 void sort_VCR_BANDS(void); void prep_VCR_COVAR(void) ; void get_VCR_colorOffset(int ic, double *OFFSET, double *RMS ); void parse_VCR_colorString(int ic); void get_NRAN_genSmear(int *NRANGauss, int *NRANFlat); // returns NRANxxx void SETSNPAR_genSmear(double shape, double color, double redshift) ; void get_genSmear(double *parList, int NLam, double *Lam, double *magSmear) ; int repeat_genSmear(double Trest, int NLam, double *Lam); void load_genSmear_randoms(int CID, double rmin, double rmax, double RANFIX); void init_genSmear_COVLAM_debug(double *lam, double COVMAT[2][2] ); void update_genSmear_COVLAM_debug(double *magSmear); double get_genSmear_SCALE(double *parList); void get_genSmear_USRFUN(double Trest, int NLam, double *Lam, double *magSmear ) ; void get_genSmear_SALT2(double Trest, int NLam, double *Lam, double *magSmear ) ; void get_genSmear_Chotard11(double Trest, int NLam, double *Lam, double *magSmear ) ; void get_genSmear_VCR(double Trest, int NLam, double *Lam, double *magSmear ) ; double get_CCM89_RV(void) ; void get_genSmear_CCM89(double Trest, int NLam, double *Lam, double *magSmear ) ; void get_genSmear_COH(double Trest, int NLam, double *Lam, double *magSmear) ; void get_genSmear_biModalUV(double Trest, int NLam, double *Lam, double *magSmear ) ; void get_genSmear_OIR(double Trest, int NLam, double *Lam, double *magSmear ) ; void get_genSmear_COVSED(double Trest, int NLam, double *Lam, double *magSmear ) ; void get_genSmear_private(double Trest, int NLam, double *Lam, double *magSmear ) ; int INODE_LAMBDA(double LAM, int NNODE, double *LAM_NODES); void extraFilters_4genSmear(char *modelName, int *NFILT_extra, int *IFILTOBS_extra) ; #define MXBAND_GENSMEAR 20 void genSmear_nodes(int NBAND, double *LAMCEN, double CholMatrix[MXBAND_GENSMEAR][MXBAND_GENSMEAR], int NLam, double *Lam, double *magSmear ) ; // ============= GLOBALS ================ #define MXRAN_GENSMEAR 500 // 100->500 Oct 18 2019 // define central wavelength for commonly used Bessell filters #define LAMCEN_FILT_v 2500.0 #define LAMCEN_FILT_U 3560.0 #define LAMCEN_FILT_B 4390.0 #define LAMCEN_FILT_V 5490.0 #define LAMCEN_FILT_R 6545.0 #define LAMCEN_FILT_I 8045.0 struct GENSMEAR { int NUSE; // number of used models with init_genSmear_XXX call int NCALL; // number of calls to get_genSmear int NGEN_RANGauss ; int NGEN_RANFlat ; int NSET_RANGauss ; int NSET_RANFlat ; double *RANGauss_LIST; // [MXRAN_GENSMEAR] ; double *RANFlat_LIST; // [MXRAN_GENSMEAR] ; int CID; // CID for each set of randoms double SHAPE, COLOR ; // allows more complex magSmear models double REDSHIFT ; // allows for redshift evolution (Jan 2014) double ZPOW[10]; // store powers of redshift for faster calculations. int MSKOPT ; // generic opt-mask // debug option to check scatter at 'CHECK' wavelength double SUMSMEAR_CHECK[2], SQSUMSMEAR_CHECK[2], SUMCROSS; int NCHECK; // Nov 2019: allow up to two independent COH scatter terms double MAGSMEAR_COH[2]; // keep track of last values to check of magSmear needs // to be re-computed int CID_LAST, NLAM_LAST; double TREST_LAST, LAMMIN_LAST, LAMMAX_LAST; // define global mag-vs-lanbda to allow repeat function to re-use (Jan 2020) double *MAGSMEAR_LIST ; } GENSMEAR ; // Mar 22 2020: // define scaling of mag-smearing as global scale or ploynomial func struct { int USE ; double GLOBAL ; // global scale ... or ... GENPOLY_DEF POLY ; // polynominal function of ... char VARNAME[20]; // this variable name } GENSMEAR_SCALE ; // ---------- USRFUN struct ----------------- #define NPAR_GENSMEAR_USRFUN_REQUIRED 8 struct GENSMEAR_USRFUN { int USE ; double SIGCOH ; // sigma for coherent scatter double A5500 ; double LAMINTERVAL; double LAMPHASE ; double TAU_LAM ; double TAU_EPOCH ; double CORR_LAM ; // wavelength correlation, Angstroms double CORR_EPOCH ; // epoch correlation, days double LAMREF ; // Ampl defined here double LAMRANGE[2]; // define temp variables in case we need to debug int NNODE; double LAM_NODE[MXRAN_GENSMEAR] ; // wavelength at each node double AMPL_NODE[MXRAN_GENSMEAR] ; // Amplitude at each node } GENSMEAR_USRFUN ; // ---------- G10 struct ----------------- #define MXLAM_GENSMEAR_SALT2 4000 struct GENSMEAR_SALT2 { int USE ; char FILE[200] ; // mag smear vs. wavelength int NLAM ; // no. lambda bins defining color smearing double LAM[MXLAM_GENSMEAR_SALT2]; // lambda at each bin double MINLAM, MAXLAM; // min,max wavelenght double SIGMA[MXLAM_GENSMEAR_SALT2]; // sigma-smear value, mags double SIGMA_SCALE ; // SIGMA_INT from SALT2.INFO file or from user input GRIDMAP1D_DEF SIGCOH_LAM ; // May 30 2018 : sigma_coh vs. wavelenth double SIGCOH ; // traditional coherent term double LAMSEP_NODE ; double LAM_NODE[MXRAN_GENSMEAR]; double SIG_NODE[MXRAN_GENSMEAR]; int NNODE ; double FUDGE_dSmear_dLAM[2] ; // Fudge G10 curve to match broadband smear double FUDGE_LAMSWITCH ; } GENSMEAR_SALT2 ; // ---------- Chotard11 struct ----------------- #define NBAND_C11 6 // UBVRI correlations struct GENSMEAR_C11 { int USE ; CHOLESKY_DECOMP_DEF DECOMP ; int OPT_farUV; // see sub-models C11_0, C11_1, C11_2 } GENSMEAR_C11 ; // ------------ OIR struct ---------------------- #define NBAND_OIR 7 // BgriYJH correlations struct GENSMEAR_OIR { int USE ; double Cholesky[NBAND_OIR][NBAND_OIR]; // path and filenames char MODELPATH[MXPATHLEN] ; char INFO_FILE[MXPATHLEN] ; // inputs from OIR.INFO file int NCOLOR ; char COLOR_STRING[NBAND_OIR][8]; double COLOR_SLOPE[NBAND_OIR] ; double SIGMACOH ; double COLOR_SIGMA_SCALE ; // multiplies all the COLOR_SIGMA values. double COLOR_SIGMA[NBAND_OIR] ; double COLOR_CORMAT[NBAND_OIR][NBAND_OIR] ; double LAMCEN_BAND[NBAND_OIR*2]; // define warp nodes double SPECSHIFT_SCALE; // fudge-scale for spectral shifts // define list of filters ordered by wavelength int NBAND_OIRDEF ; // number of unique filters = size of ordered list int ORDERED_IFILTDEF[NBAND_OIR*2]; // sparse/ordered list double ORDERED_LAMCEN[NBAND_OIR*2]; // sparse/ordered list int LAMCEN[MXFILTINDX] ; // vs. IFILTDEF int IFILT_B ; } GENSMEAR_OIR; struct GENSMEAR_COVSED { int USE ; int OPTMASK ; char MODEL_PATH[MXPATHLEN], VERSION[80] ; char COVSED_FILE[MXPATHLEN]; int REBIN_LAM; double LAMPAIR_DEBUG[2]; // COVSED and binning read from FITS file int NBIN_WAVE, NBIN_EPOCH, NBIN_WAVExEPOCH, NBIN_COVMAT ; double *WAVE, *EPOCH; // matrix used to generate correlated randoms CHOLESKY_DECOMP_DEF DECOMP; // xxx marl delete double **Cholesky; // magSmear scatter values for each event. double *SCATTER_VALUES; } GENSMEAR_COVSED; // ------------- CCM89 struct --------------- struct GENSMEAR_CCM89 { int USE ; double SIGCOH ; // coherent term double RV_REF; // reference RV; RV=RV_REF -> no fluctuation double RV_sigma ; // sigma for random RV double RV_range[2] ; // generate withing this range int NLAM_MAP ; double *LAM_MAP ; double *XTDIF_MAP; // XT-XT0 vs. LAM at current RV double RV ; // randome RV for each SN } GENSMEAR_CCM89; #define MXCOLOR_VCR 10 struct GENSMEAR_VCR { int USE ; // path and filenames char MODELPATH[MXPATHLEN] ; char INFO_FILE[MXPATHLEN] ; char VSI_FILE[MXPATHLEN]; // inputs from VCR.INFO file int NCOLOR ; char COLOR_STRING[MXCOLOR_VCR][8]; double COLOR_SLOPE[MXCOLOR_VCR] ; double SIGMACOH_MB ; double COLOR_SIGMA_SCALE ; // multiplies all the COLOR_SIGMA values. double COLOR_SIGMA[MXCOLOR_VCR] ; double COLOR_CORMAT[MXCOLOR_VCR][MXCOLOR_VCR] ; double LAMCEN_BAND[MXCOLOR_VCR*2]; // define warp nodes double SPECSHIFT_SCALE; // fudge-scale for spectral shifts int USE_ZVAR ; // logical flag double ZVARIATION_POLY_VSI[4]; // VSI += sum [z^i * coef[i] ] // VSI distribution read from file int NBIN_VSI; // number of bins for v_Si distribution double *VSI, *PROB ; // v_Si prob distribution double *SUMPROB ; // cumulative prov (for efficiency MC generation) double VSI_MEAN ; // misc. double COLOR_OFFSET[MXCOLOR_VCR] ; double COLOR_RMS[MXCOLOR_VCR] ; char COLOR_BAND[MXCOLOR_VCR][2][2]; // each band int COLOR_IFILTDEF[MXCOLOR_VCR][2] ; // absolute filter indices double COLOR_COVAR2[MXCOLOR_VCR][MXCOLOR_VCR] ; // covariances double COLOR_COVAR1[MXCOLOR_VCR*MXCOLOR_VCR] ; // covariances double Cholesky[MXCOLOR_VCR][MXCOLOR_VCR] ; int COLOR_IFILT[MXCOLOR_VCR] ; // points to ORDERED_IFILTDEF // define list of filters ordered by wavelength int NFILTDEF ; // number of unique filters = size of ordered list int ORDERED_IFILTDEF[MXCOLOR_VCR*2]; // sparse/ordered list double ORDERED_LAMCEN[MXCOLOR_VCR*2]; // sparse/ordered list int LAMCEN[MXFILTINDX] ; // vs. IFILTDEF int IFILT_B ; // generated quantieis for SIMGEN-DUMP file double GENRAN_VSI ; // generated VSI double GENRAN_COLORSHIFT[MXCOLOR_VCR]; // broadband color shift } GENSMEAR_VCR ; struct GENSMEAR_COH { int USE ; int NSIGMA; double MAGSIGMA[2] ; } GENSMEAR_COH ; struct GENSMEAR_BIMODAL_UV { int USE ; double SIGMACOH ; double LAMU_MIN, LAMU_CEN, LAMU_MAX; double MAGU_SPLIT, MAGU_SIGMA ; } GENSMEAR_BIMODAL_UV ; // see sim input GENMAG_SMEAR_ADDPHASECOR: struct { int USE; double INPUT_MAGSMEAR, INPUT_EXPTAU; // copied from user input int NBIN ; double BINSIZE ; double RANGE[2] ; // min & max phase to store COVMAT double *GRID_PHASE ; // phase value at each grid point double *GRID_MAGSMEAR ; // for each event double *RANGauss_LIST ; // new randoms for each event int CID_LAST ; // used to identify new event CHOLESKY_DECOMP_DEF DECOMP; int NCHECK, NSUM; double SUMCHECK[10]; // for internal COV check double sumCHECK ; } GENSMEAR_PHASECOR ; // --------- private struct for testing -------------- struct GENSMEAR_PRIVATE { int USE ; } GENSMEAR_PRIVATE ; // -------------------------------------------------------------------- // Override params passed from simulation to allow command-line overrides // See sim-input key GENMAG_SMEARPAR_OVERRIDE (Mar 22 2014) #define MXSMEARPAR_OVERRIDE 20 int NSMEARPAR_OVERRIDE ; struct { int NVAL ; // number of values passed per SMEARPAR (default=1) char NAME[60]; // name of parameter to override double VALUE[MXSMEARPAR_OVERRIDE] ; // list of value(s) } GENMAG_SMEARPAR_OVERRIDE[MXSMEARPAR_OVERRIDE] ;