ugali.utils.skymap.writeSparseHealpixMap(self.roi.pixels_target, data_dict, self.config['coords']['nside_pixel'], outfile, distance_modulus_array=self.distance_modulus_array, coordsys='NULL', ordering='NULL', header_dict=header_dict) ############################################################ if __name__ == "__main__": import ugali.utils.parser description = "Script for executing the likelihood scan." parser = ugali.utils.parser.Parser(description=description) parser.add_config() parser.add_argument('outfile',metavar='outfile.fits',help='Output fits file.') parser.add_debug() parser.add_verbose() parser.add_coords(required=True,radius=False) opts = parser.parse_args() if len(opts.coords) != 1: raise Exception('Must specify exactly one coordinate.') lon,lat,radius = opts.coords[0] grid = createGridSearch(opts.config,lon,lat) if not opts.debug: result = grid.search() grid.write(opts.outfile)
logger.info("Fit parameter values:") for d in dtype: logger.info('\t%s: %s'%(d[0], results[i][d[0]])) if (i%self.config['simulate']['save'])==0 and outfile: ugali.utils.fileio.write(outfile,results) if outfile: ugali.utils.fileio.write(outfile,results,clobber=True) return results if __name__ == "__main__": import ugali.utils.parser parser = ugali.utils.parser.Parser(description=__doc__) parser.add_config() parser.add_argument('-p','--popfile',default=None, help='simulated population input file') parser.add_argument('-c','--catfile',default=None, help='simulated catalog input file') parser.add_argument('-o','--outfile',default=None, help='output results file') parser.add_argument('-i','--mc-source-id',default=None,type=int,action='append', help='specific source id to run') parser.add_force() parser.add_debug() parser.add_verbose() args = parser.parse_args() analyzer = Analyzer(args.config,args.catfile,args.popfile) if args.mc_source_id is None: basename = os.path.splitext(args.catfile)[0]
if samples is not None: mcmc.load_samples(samples) return mcmc if __name__ == "__main__": import ugali.utils.parser description = "Script for running MCMC followup" parser = ugali.utils.parser.Parser(description=description) parser.add_config() parser.add_debug() parser.add_verbose() #parser.add_coords(required=True) parser.add_coords() parser.add_name() parser.add_argument('--srcmdl',help='Source model file') parser.add_argument('--grid',action='store_true',help='Grid search for intial parameters') parser.add_argument('outfile',default=None,help="Output file name") opts = parser.parse_args() if opts.coords is not None and len(opts.coords) != 1: raise Exception('Must specify exactly one coordinate.') #try: # # Initialize the MPI-based pool used for parallelization. # from emcee.utils import MPIPool # pool = MPIPool(loadbalance=True) #except (ImportError,ValueError) as e: # logger.warn(e.message) # pool = None
clobber=True) fitsio.write(outfile, dict(DISTANCE_MODULUS=self.distance_modulus_array.astype( 'f4', copy=False)), extname='DISTANCE_MODULUS', clobber=False) if __name__ == "__main__": import ugali.utils.parser description = "Script for executing the likelihood scan." parser = ugali.utils.parser.Parser(description=description) parser.add_config() parser.add_argument('outfile', metavar='outfile.fits', help='Output fits file.') parser.add_argument('-m', '--mlimit', nargs='?', default=-1, type=int, help='Memory limit (KB)') parser.add_debug() parser.add_verbose() parser.add_coords(required=True, radius=False) opts = parser.parse_args() if len(opts.coords) != 1: raise Exception('Must specify exactly one coordinate.') lon, lat, radius = opts.coords[0]
class Marigo2017(Padova): defaults = dict(defaults_30) defaults['isoc_kind'] = 'parsec_CAF09_v1.2S_NOV13' def factory(name, **kwargs): from ugali.utils.factory import factory return factory(name, module=__name__, **kwargs) if __name__ == "__main__": import ugali.utils.parser description = "Download isochrones" parser = ugali.utils.parser.Parser(description=description) parser.add_verbose() parser.add_force() parser.add_argument('-a','--age',default=None,type=float) parser.add_argument('-z','--metallicity',default=None,type=float) parser.add_argument('-k','--kind',default='Bressan2012') parser.add_argument('-s','--survey',default='des') parser.add_argument('-o','--outdir',default=None) parser.add_argument('-n','--njobs',default=10,type=int) args = parser.parse_args() if args.verbose: try: from http.client import HTTPConnection except ImportError: from httplib import HTTPConnection HTTPConnection.debuglevel = 1 if args.outdir is None: args.outdir = os.path.join(args.survey.lower(),args.kind.lower())
outfile.write('s.WAVG_SPREAD_MODEL_R, s.CLASS_STAR_R \n') outfile.write('FROM %s s \n' % (table)) outfile.write('WHERE %s; \n' % (select)) outfile.write(' > %s \n' % (db)) outfile.close() if __name__ == "__main__": import ugali.utils.parser description = "Download data set." parser = ugali.utils.parser.Parser(description=description) parser.add_config() parser.add_debug() parser.add_verbose() parser.add_argument('pixfile', metavar='pixels.dat', default=None, nargs='?', help='Input pixel file.') opts = parser.parse_args() from ugali.utils.config import Config config = Config(opts.config) survey = config['data']['survey'].lower() release = config['data']['release'].lower() outdir = config['data']['dirname'] db = databaseFactory(survey, release) db.run(pixfile=opts.pixfile, outdir=outdir)
results.load_samples(samples) return results def write_results(filename, config, srcfile, samples): """ Package everything nicely """ results = createResults(config, srcfile, samples=samples) results.write(filename) if __name__ == "__main__": import ugali.utils.parser parser = ugali.utils.parser.Parser(description=__doc__) parser.add_config() parser.add_verbose() parser.add_argument('--srcmdl', required=True, help='Source model file') parser.add_argument('--section', default='source', help='Section of source file') parser.add_argument('--samples', required=True, help='Posterior samples file') parser.add_argument('outfile', default=None, help="Output file name") args = parser.parse_args() #write_results(args.outfile,args.config,args.srcmdl,args.samples) results = createResults(args.config, args.srcmdl, samples=args.samples) results.write(args.outfile)
import ugali.utils.parser import numpy as np from ugali.utils.projector import gal2cel, ang2const, ang2iau #CATALOGS = ['McConnachie12','Rykoff14', 'Harris96', 'Corwen04', 'Nilson73', 'Webbink85', 'Kharchenko13', 'WEBDA14'] CATALOGS = [ 'McConnachie12', 'Harris96', 'Corwen04', 'Nilson73', 'Webbink85', 'Kharchenko13', 'Bica08', 'WEBDA14', 'ExtraDwarfs', 'ExtraClusters' ] if __name__ == "__main__": import argparse description = "python script" parser = ugali.utils.parser.Parser(description=description) parser.add_coords(required=True, radius=True, targets=True) parser.add_argument('-n', '--nnearest', default=1, type=int) opts = parser.parse_args() catalog = ugali.candidate.associate.SourceCatalog() for i in CATALOGS: catalog += ugali.candidate.associate.catalogFactory(i) for name, (glon, glat, radius) in zip(opts.names, opts.coords): ra, dec = gal2cel(glon, glat) iau = ang2iau(glon, glat) const = ang2const(glon, glat)[0] if radius <= 0: radius = None idx1, idx2, sep = catalog.match([glon], [glat], tol=radius, nnearest=opts.nnearest)
results.load_samples(samples) return results def write_results(filename,config,srcfile,samples): """ Package everything nicely """ results = createResults(config,srcfile,samples=samples) results.write(filename) if __name__ == "__main__": import ugali.utils.parser parser = ugali.utils.parser.Parser(description=__doc__) parser.add_config() parser.add_verbose() parser.add_argument('--srcmdl',required=True, help='Source model file') parser.add_argument('--section',default='source', help='Section of source file') parser.add_argument('--samples',required=True, help='Posterior samples file') parser.add_argument('outfile',default=None, help="Output file name") args = parser.parse_args() #write_results(args.outfile,args.config,args.srcmdl,args.samples) results = createResults(args.config,args.srcmdl,samples=args.samples) results.write(args.outfile)
defaults = dict(defaults_30) defaults['isoc_kind'] = 'parsec_CAF09_v1.2S_NOV13' def factory(name, **kwargs): from ugali.utils.factory import factory return factory(name, module=__name__, **kwargs) if __name__ == "__main__": import ugali.utils.parser description = "Download isochrones" parser = ugali.utils.parser.Parser(description=description) parser.add_verbose() parser.add_force() parser.add_argument('-a', '--age', default=None, type=float) parser.add_argument('-z', '--metallicity', default=None, type=float) parser.add_argument('-k', '--kind', default='Bressan2012') parser.add_argument('-s', '--survey', default='des') parser.add_argument('-o', '--outdir', default=None) parser.add_argument('-n', '--njobs', default=10, type=int) args = parser.parse_args() if args.verbose: try: from http.client import HTTPConnection except ImportError: from httplib import HTTPConnection HTTPConnection.debuglevel = 1 if args.outdir is None:
assert np.allclose(age,float(a),atol=1e-5) except: msg = "Age does not match:\n"+lines[13] raise Exception(msg) def factory(name, **kwargs): from ugali.utils.factory import factory return factory(name, module=__name__, **kwargs) if __name__ == "__main__": import ugali.utils.parser description = "Download Dotter isochrones" parser = ugali.utils.parser.Parser(description=description) parser.add_verbose() parser.add_force() parser.add_argument('-a','--age',default=None,type=float,action='append') parser.add_argument('-z','--metallicity',default=None,type=float,action='append') parser.add_argument('-k','--kind',default='Dotter2008') parser.add_argument('-s','--survey',default='des') parser.add_argument('-o','--outdir',default=None) parser.add_argument('-n','--njobs',default=1,type=int) args = parser.parse_args() if args.verbose: try: from http.client import HTTPConnection except ImportError: from httplib import HTTPConnection HTTPConnection.debuglevel = 1 if args.outdir is None:
outfile.write('s.WAVG_MAG_PSF_G, s.WAVG_MAGERR_PSF_G, \n') outfile.write('s.MAG_PSF_R, s.MAGERR_PSF_R, \n') outfile.write('s.WAVG_MAG_PSF_R, s.WAVG_MAGERR_PSF_R, \n') outfile.write('s.SPREAD_MODEL_R, s.SPREADERR_MODEL_R, \n') outfile.write('s.WAVG_SPREAD_MODEL_R, s.CLASS_STAR_R \n') outfile.write('FROM %s s \n'%(table)) outfile.write('WHERE %s; \n'%(select)) outfile.write(' > %s \n'%(db)) outfile.close() if __name__ == "__main__": import ugali.utils.parser description = "Download data set." parser = ugali.utils.parser.Parser(description=description) parser.add_config() parser.add_debug() parser.add_verbose() parser.add_argument('pixfile',metavar='pixels.dat',default=None, nargs='?',help='Input pixel file.') opts = parser.parse_args() from ugali.utils.config import Config config = Config(opts.config) survey = config['data']['survey'].lower() release = config['data']['release'].lower() outdir = config['data']['dirname'] db = databaseFactory(survey,release) db.run(pixfile=opts.pixfile,outdir=outdir)
mcmc.load_samples(samples) return mcmc if __name__ == "__main__": import ugali.utils.parser description = "Script for running MCMC followup" parser = ugali.utils.parser.Parser(description=description) parser.add_config() parser.add_debug() parser.add_verbose() #parser.add_coords(required=True) parser.add_coords() parser.add_name() parser.add_argument('--srcmdl', help='Source model file') parser.add_argument('--grid', action='store_true', help='Grid search for intial parameters') parser.add_argument('outfile', default=None, help="Output file name") opts = parser.parse_args() if opts.coords is not None and len(opts.coords) != 1: raise Exception('Must specify exactly one coordinate.') #try: # # Initialize the MPI-based pool used for parallelization. # from emcee.utils import MPIPool # pool = MPIPool(loadbalance=True) #except (ImportError,ValueError) as e:
clobber=True) fitsio.write(outfile, dict(DISTANCE_MODULUS=self.distance_modulus_array.astype( 'f4', copy=False)), extname='DISTANCE_MODULUS', clobber=False) if __name__ == "__main__": import ugali.utils.parser description = "Script for executing the likelihood scan." parser = ugali.utils.parser.Parser(description=description) parser.add_config() parser.add_argument('outfile', metavar='outfile.fits', help='Output fits file.') parser.add_debug() parser.add_verbose() parser.add_coords(required=True, radius=False) opts = parser.parse_args() if len(opts.coords) != 1: raise Exception('Must specify exactly one coordinate.') lon, lat, radius = opts.coords[0] grid = createGridSearch(opts.config, lon, lat) if not opts.debug: result = grid.search() grid.write(opts.outfile)
logger.info('\t%s: %s' % (d[0], results[i][d[0]])) if (i % self.config['simulate']['save']) == 0 and outfile: ugali.utils.fileio.write(outfile, results) if outfile: ugali.utils.fileio.write(outfile, results, clobber=True) return results if __name__ == "__main__": import ugali.utils.parser parser = ugali.utils.parser.Parser(description=__doc__) parser.add_config() parser.add_argument('-p', '--popfile', default=None, help='simulated population input file') parser.add_argument('-c', '--catfile', default=None, help='simulated catalog input file') parser.add_argument('-o', '--outfile', default=None, help='output results file') parser.add_argument('-i', '--mc-source-id', default=None, type=int, action='append', help='specific source id to run')
results[i]['FIT_MASS_ERR'] = grid.stellar_mass_conversion*richness_err distance_modulus_err = (err['distance_modulus'][1]-err['distance_modulus'][0])/2. results[i]['FIT_DISTANCE_ERR'] = distance_modulus_err """ return results run = runall if __name__ == "__main__": import ugali.utils.parser parser = ugali.utils.parser.Parser(description=__doc__) parser.add_config() parser.add_argument('-p', '--popfile', default=None, help='simulated population input file') parser.add_argument('-c', '--catfile', default=None, help='simulated catalog input file') parser.add_argument('-o', '--outfile', default=None, help='output results file') parser.add_argument('-i', '--mc-source-id', default=None, type=int, action='append', help='specific source id to run')
def factory(name, **kwargs): from ugali.utils.factory import factory return factory(name, module=__name__, **kwargs) if __name__ == "__main__": import ugali.utils.parser description = "Download Dotter isochrones" parser = ugali.utils.parser.Parser(description=description) parser.add_verbose() parser.add_force() parser.add_argument('-a', '--age', default=None, type=float, action='append') parser.add_argument('-z', '--metallicity', default=None, type=float, action='append') parser.add_argument('-k', '--kind', default='Dotter2008') parser.add_argument('-s', '--survey', default='des') parser.add_argument('-o', '--outdir', default=None) parser.add_argument('-n', '--njobs', default=1, type=int) args = parser.parse_args() if args.verbose: try:
#!/usr/bin/env python import ugali.candidate.associate import ugali.utils.parser import numpy as np from ugali.utils.projector import gal2cel,ang2const,ang2iau #CATALOGS = ['McConnachie12','Rykoff14', 'Harris96', 'Corwen04', 'Nilson73', 'Webbink85', 'Kharchenko13', 'WEBDA14'] CATALOGS = ['McConnachie12', 'Harris96', 'Corwen04', 'Nilson73', 'Webbink85', 'Kharchenko13', 'Bica08', 'WEBDA14', 'ExtraDwarfs','ExtraClusters'] if __name__ == "__main__": import argparse description = "python script" parser = ugali.utils.parser.Parser(description=description) parser.add_coords(required=True,radius=True,targets=True) parser.add_argument('-n','--nnearest',default=1,type=int) opts = parser.parse_args() catalog = ugali.candidate.associate.SourceCatalog() for i in CATALOGS: catalog += ugali.candidate.associate.catalogFactory(i) for name,(glon,glat,radius) in zip(opts.names, opts.coords): ra,dec = gal2cel(glon,glat) iau = ang2iau(glon,glat) const = ang2const(glon,glat)[0] if radius <= 0: radius = None idx1,idx2,sep = catalog.match([glon],[glat],tol=radius,nnearest=opts.nnearest) match = catalog[idx2] if len(match) > 0:
defaults['isoc_kind'] = 'gi10b' class Bressan2012(Padova): defaults = dict(defaults_27) defaults['isoc_kind'] = 'parsec_CAF09_v1.2S' if __name__ == "__main__": import ugali.utils.parser description = "Download isochrones" parser = ugali.utils.parser.Parser(description=description) parser.add_config() parser.add_debug() parser.add_verbose() parser.add_force() parser.add_argument('-a','--age',default=None,type=float) parser.add_argument('-z','--metallicity',default=None,type=float) opts = parser.parse_args() from ugali.utils.config import Config config = Config(opts.config) #outdir = '/u/ki/kadrlica/des/isochrones/v1' #outdir = '/u/ki/kadrlica/sdss/isochrones/v2' #outdir = '/u/ki/kadrlica/des/isochrones/v2' #outdir = '/u/ki/kadrlica/des/isochrones/v3' #outdir = '/u/ki/kadrlica/des/isochrones/v4' outdir = '/u/ki/kadrlica/des/isochrones/v5' #outdir = '/u/ki/kadrlica/des/isochrones/v6' survey = config['data']['survey']