def setup_options(options): if options is None: parser= get_options() options, args= parser.parse_args([]) #Set up to default fit options.potential= 'mpdiskplhalofixbulgeflat' options.aAmethod= 'staeckelg' options.dfeh=0.1 options.dafe=0.05 options.ngl= 20 options.singles= True return options
r'$[\alpha/\mathrm{Fe}] = %.3f$' % (afes[bin]), size=16.,top_right=True) left, bottom, width, height= 0.1, 0.1, 0.8, 0.25 ax2= pyplot.axes([left,bottom,width,height]) fig= pyplot.gcf() fig.sca(ax2) bovy_plot.bovy_plot(rs,derivProps[:,2],'k-',lw=2., xrange=[4.,10.], overplot=True) indx= numpy.argmin(numpy.fabs(derivProps[:,2])) bovy_plot.bovy_plot([4.,10.],[0.,0.],'-',color='0.5',overplot=True) bovy_plot.bovy_plot([rs[indx],rs[indx]],[derivProps[indx,2],1000.], 'k--',overplot=True) thisax= pyplot.gca() pyplot.ylabel(r'$\mathrm{Correlation\ between}$'+'\n'+r'$\!\!\!\!\!\!\!\!R_d\ \&\ \Sigma_{1.1}(R)$') pyplot.xlabel(r'$R\ (\mathrm{kpc})$') pyplot.xlim(4.,10.) pyplot.ylim(-.5,.5) bovy_plot._add_ticks() pyplot.sca(axTop) bovy_plot.bovy_plot([rs[indx],rs[indx]],[0.01,sigs[3,indx]], 'k--',overplot=True) bovy_plot.bovy_end_print(options.outfilename) if __name__ == '__main__': parser= get_options() options,args= parser.parse_args() illustrateBestR(options,args)
#First do all dens comparisons if not options.justvel: options.type= 'z' plotDensComparisonDFMulti(options,args) #options.type= 'R' #plotDensComparisonDFMulti(options,args) options.type= 'vt' plotVelComparisonDFMulti(options,args) options.type= 'vz' plotVelComparisonDFMulti(options,args) options.type= 'vr' plotVelComparisonDFMulti(options,args) return None if __name__ == '__main__': (options,args)= get_options().parse_args() try: os.mkdir(os.path.join('..',args[0][0:-1])) except OSError: pass args[0]= os.path.join('..',args[0][0:-1],args[0]) if options.allgroups: options.group= 'aenhanced' plotAllComparisonDFMulti(options,args) options.group= 'apoor' plotAllComparisonDFMulti(options,args) options.group= 'apoorfpoor' plotAllComparisonDFMulti(options,args) options.group= 'aintermediate' plotAllComparisonDFMulti(options,args) else: