def getmax(dset,cols):
ldata = l.getcols(dset,cols)
maxlogL = max(np.exp(np.sum(ldata, axis=1))) #sum rows (likelihood contributions for each point) then find maximum
return maxlogL
line = ax.plot(data[:,0],data[:,1],'.') #ax.set_ylim(logL[100000],max(logL)) """ """ t0 = time.time() data = l.getcols(txtdset,['M0','neg2LogL']) print 't:', time.time()-t0 fig = plt.figure() ax = fig.add_subplot(111) line = ax.plot(data[:,0],data[:,1],'.') #ax.set_ylim(logL[100000],max(logL)) """ t0 = time.time() pars=['M0','TanB','neg2LogL'] data = l.getcols(txtdset,pars) print 't:', time.time()-t0 #l.bin2d(data,binop='min') fig2 = plt.figure() ax = fig2.add_subplot(111) plot = l.profplot(ax,data,labels=pars[0:2]) plt.show() """ fig2 = plt.figure() ax = fig2.add_subplot(111) plot = l.chi2scatplot(ax,data,labels=pars[0:2]) cbar = fig2.colorbar(plot, ticks=[0, 1, 2, 3, 4, 5]) cbar.ax.set_yticklabels(['0','1','4','9','16',r'$\geq$25'])# vertically oriented colorbar
#root='/media/backup1/temp_CMSSM_backup/pysusyAnalysis/CMSSMpoints/CMSSMproject_tmp/m2/' #outdirroot='/home/farmer/Dropbox/Projects/pysusy_pytools/tmp_plots/m2/' #dirs = ['5SUSY','6SUSYCCR','7higgs','8higgsCCR'] #dsets = ['SUSY/log/mu+/txt','SUSY/CCR/mu+/txt','higgs/log/mu+/txt','higgs/CCR/mu+/txt'] #dirs = ['5SUSY','6SUSYCCR','8higgsCCR'] #dsets = ['SUSY/log/mu+/txt','SUSY/CCR/mu+/txt','higgs/CCR/mu+/txt'] #dirs = ['SUSYCCR'] #dsets = ['SUSY/CCR/mu+/txt'] #------------------------------------------------- # TASKS #-------------------------------------------------- #----------FIND LOWEST HIGGS MASS OBTAINED IN PRE-LEP DATASET----------- f = h5py.File(root+'CMSSMev.hdf5','r') dset = f['preLEP/log/mu+/txt'] data = l.getcols(dset,['higgs','neg2LogL']) print 'Minimum Higgs mass in dataset: ', min(data[:,0]) fig= plt.figure() plt.plot(data[:,0],np.exp(-0.5*(data[:,1] - min(data[:,1]))),'.') plt.show()