databaseFile = '/home/jeff/work/rotNSruns/mass-shed-models.db' connection = sqlite3.connect(databaseFile) c = connection.cursor() c.execute("SELECT DISTINCT a FROM MODELS") ahs = c.fetchall() databaseFile2 = '/home/jeff/work/rotNSruns/ls-cold-manual2.db' connection2 = sqlite3.connect(databaseFile2) c2 = connection2.cursor() for a in ahs: sequencePlot(["edMax", "baryMass"], c2, ("RedMax=0", "eos='LS220' ", " a=" + str(a[0])), grid=True, title="", suppressShow=True, c=matplotlib.pyplot.cm.jet(a[0])) sequencePlot(["edMax", "baryMass"], c2, ("RedMax>0", "eos='LS220' ", " a=" + str(a[0])), grid=True, title="", suppressShow=True, c=matplotlib.pyplot.cm.jet(a[0]), linestyle='--') sequencePlot(["edMax", "baryMass"], c, ("RedMax=0", "eos='HS'", "T=0.5"), "a", grid=True, title="", suppressShow=True,
#print edMaxes #c.execute("SELECT DISTINCT a FROM MODELS ORDER BY a") #aes=c.fetchall() #print aes # Here we try and do Ansorg plots SHED vs RPOE with a curves aes = arange(0.0, 1.1, 0.1) print aes filtersT = nearValueFilter("T", 0.5) filtersEd = nearValueFilter("edMax", 528787878700000.0) filters = filtersT + filtersEd for a in aes: thisFilter = filters + nearValueFilter("a", a) sequencePlot(["rpoe", "shed"], c, thisFilter, suppressShow=True, c=(a, a / 2.0, 0.5)) sequencePlot(["rpoe", "shed"], c, thisFilter, c=(a, a / 2.0, 0.5)) exit() T = 0.5 filtersT = nearValueFilter("T", T, 0.001) #("T>%s" % (T -.01), "T<%s" % (T +.01)) a = 0.7 filtersA = nearValueFilter("a", a, 0.001) #("a>%s" % (a -.001), "a<%s" % (a +.001)) filtersEd = ("edMax>1.5e14", ) # nearValueFilter("edMax",528787878700000.0,0.001 ) filters = filtersA + filtersT + filtersEd
resultTov=[] for T in cTov.execute("SELECT DISTINCT T FROM models"): T=T[0] print T answer = edOfMaxMass(equalsFiltersFromDict({'T': T}) ,connectionTov) resultTov.append((-.10, T, answer[plotVar])) databaseFile = '/home/jeff/work/rotNSruns/allRuns3-20-13.db' connection=sqlite3.connect(databaseFile) c = connection.cursor() sequencePlot(["edMax", "baryMass"], c, equalsFiltersFromDict({'a': 0.8, 'RedMax': 0.0}),"rpoe",grid=True) c.execute("SELECT DISTINCT a,T FROM models") grid = c.fetchall() c.execute("SELECT DISTINCT T FROM MODELS") Ts=c.fetchall() c.execute("SELECT DISTINCT a FROM MODELS") ahs=c.fetchall() resultS =[] resultT =[] for a,T in grid: print a,T
from sqlPlotRoutines import sequencePlot databaseFile = '/home/jeff/work/rotNSruns/mass-shed-models.db' connection=sqlite3.connect(databaseFile) c=connection.cursor() c.execute("SELECT DISTINCT a FROM MODELS") ahs=c.fetchall() databaseFile2 = '/home/jeff/work/rotNSruns/ls-cold-manual2.db' connection2=sqlite3.connect(databaseFile2) c2=connection2.cursor() for a in ahs: sequencePlot(["edMax","baryMass"],c2,("RedMax=0","eos='LS220' "," a="+str(a[0])), grid=True,title="",suppressShow=True,c=matplotlib.pyplot.cm.jet(a[0])) sequencePlot(["edMax","baryMass"],c2,("RedMax>0","eos='LS220' "," a="+str(a[0])), grid=True,title="",suppressShow=True,c=matplotlib.pyplot.cm.jet(a[0]),linestyle='--') sequencePlot(["edMax","baryMass"],c,("RedMax=0","eos='HS'","T=0.5"),"a", grid=True,title="",suppressShow=True,marker='o',s=40,vmin=0.0,vmax=1.0) sequencePlot(["edMax","baryMass"],c,("RedMax>0","eos='HS'","T=0.5"),"a", grid=True,title="$M_b$ vs $Ed_{max}$.\n Lines: $Cold LS220$ /w Dashed => toroid \n Points: $Cold HS$ /w Triangle => toroid ", suppressShow=False,marker='^',s=40,vmin=0.0,vmax=1.0) connection.close() del connection databaseFile = '/home/jeff/work/rotNSruns/hotTov.db' connectionTov=sqlite3.connect(databaseFile) cTov=connectionTov.cursor()
chs=connectionhs.cursor() databaseFile = '/home/jeff/work/rotNSruns/ls-TOV-new.db' connectionls=sqlite3.connect(databaseFile) cls=connectionls.cursor() #sequencePlot([xaxis,yaxis],cls,("T=0.5",filt),orderBy='edMax', # suppressShow=True) #sequencePlot([xaxis,yaxis],chs,("T=0.5",filt),orderBy='edMax', # suppressShow=True) #sequencePlot([xaxis,yaxis],cgam2,("T=0.0",filt),orderBy='edMax',grid=True, # title="Blue - LS220 Green - HS Red - Gamma=2") chs.execute("SELECT DISTINCT T FROM MODELS") Ts=chs.fetchall() print "TS: ", Ts for T in Ts: sequencePlot([xaxis,yaxis],chs, (filt,"T>%s" % (T[0] -.01), "T<%s" % (T[0] +.01)), orderBy='edMax', suppressShow=True,grid=True,title="HS EOS", c=matplotlib.pyplot.cm.jet(T[0]/40.0)) sequencePlot([xaxis,yaxis],chs, (filt,), orderBy='edMax', colorBy="T", suppressShow=False,grid=True,title="HS EOS", s=0.01)
cls = connectionls.cursor() #sequencePlot([xaxis,yaxis],cls,("T=0.5",filt),orderBy='edMax', # suppressShow=True) #sequencePlot([xaxis,yaxis],chs,("T=0.5",filt),orderBy='edMax', # suppressShow=True) #sequencePlot([xaxis,yaxis],cgam2,("T=0.0",filt),orderBy='edMax',grid=True, # title="Blue - LS220 Green - HS Red - Gamma=2") chs.execute("SELECT DISTINCT T FROM MODELS") Ts = chs.fetchall() print "TS: ", Ts for T in Ts: sequencePlot([xaxis, yaxis], chs, (filt, "T>%s" % (T[0] - .01), "T<%s" % (T[0] + .01)), orderBy='edMax', suppressShow=True, grid=True, title="HS EOS", c=matplotlib.pyplot.cm.jet(T[0] / 40.0)) sequencePlot([xaxis, yaxis], chs, (filt, ), orderBy='edMax', colorBy="T", suppressShow=False, grid=True, title="HS EOS", s=0.01)
resultTov = [] for T in cTov.execute("SELECT DISTINCT T FROM models"): T = T[0] print T answer = edOfMaxMass(equalsFiltersFromDict({'T': T}), connectionTov) resultTov.append((-.10, T, answer[plotVar])) databaseFile = '/home/jeff/work/rotNSruns/allRuns3-20-13.db' connection = sqlite3.connect(databaseFile) c = connection.cursor() sequencePlot(["edMax", "baryMass"], c, equalsFiltersFromDict({ 'a': 0.8, 'RedMax': 0.0 }), "rpoe", grid=True) c.execute("SELECT DISTINCT a,T FROM models") grid = c.fetchall() c.execute("SELECT DISTINCT T FROM MODELS") Ts = c.fetchall() c.execute("SELECT DISTINCT a FROM MODELS") ahs = c.fetchall() resultS = [] resultT = []