phi = (time - epoch + .25 * period) % period
dur = clip.trapFit.duration_hrs / 24
plt.figure()
plt.plot(phi[isbad], isbad[isbad] * -1 * clip.trapFit.depth_frac, 'ro')
plt.plot(phi, bestmodel, 'b.')
plt.xlim([clip.thermal.phimin - dur, clip.thermal.phimax + dur])
plt.ylim([-1.1 * clip.trapFit.depth_frac, 0.05 * clip.trapFit.depth_frac])
plt.title(epic)
#plt.figure()
#plt.plot(time,isbad,'rx')*-1*clip.trapFit.depth_frac
#plt.plot(time,bestmodel,'b-.')
#plt.title(epic)
#%%
import dave.pipeline.gather as g
cliploc = '/soc/nfs/so-nfs/dave/c6-v2/'
sfile = '/soc/nfs/so-nfs/dave/c6-v2/c6Candidates.txt'
clipar = np.loadtxt(sfile, dtype='str', usecols=[0], skiprows=0)
cliplist = list()
for v in clipar:
cliplist.append("%s/%s/c%s-06.clip" % (cliploc, v, v))
#%%
import dave.susanplay.thermalFlags as tf
fun = tf.getFractionThermal
(epics, values) = g.gatherFunctionB(cliplist[0:20], fun)
import dave.pipeline.exporter as ex
cliploc = '/soc/nfs/so-nfs/dave/c6-v2/'
listfile = "%s/%s" % (cliploc, 'cliplist')
clipar = np.loadtxt(listfile, dtype='str')
cliplist = list()
for v in clipar:
cliplist.append("%s/%s" % (cliploc, v))
#%%
func = ex.writeTableLine
(epics, table) = gr.gatherFunctionB(cliplist, func)
outtable = '/soc/nfs/so-nfs/dave/c6-v2/c6results.txt'
fid = open(outtable, 'w')
for e in table:
fid.write("%s\n" % e)
fid.close()
#%%
func = ex.writeCandidates
(epics, table) = gr.gatherFunctionB(cliplist, func)
outtable = '/soc/nfs/so-nfs/dave/c6-v2/c6Candidates.txt'
fid = open(outtable, 'w')
cliploc='/soc/nfs/so-nfs/dave/c6-v2/'
listfile="%s/%s" % (cliploc,'cliplist')
clipar=np.loadtxt(listfile,dtype='str')
cliplist=list()
for v in clipar:
cliplist.append("%s/%s" % (cliploc,v))
#%%
func=ex.writeTableLine
(epics,table)=gr.gatherFunctionB(cliplist, func)
outtable='/soc/nfs/so-nfs/dave/c6-v2/c6results.txt'
fid=open(outtable,'w')
for e in table:
fid.write("%s\n" % e)
fid.close()
#%%
func=ex.writeCandidates
(epics,table)=gr.gatherFunctionB(cliplist,func)
outtable='/soc/nfs/so-nfs/dave/c6-v2/c6Candidates.txt'
fid=open(outtable,'w')
phi = (time-epoch + .25*period) % period
dur=clip.trapFit.duration_hrs/24;
plt.figure()
plt.plot(phi[isbad],isbad[isbad]*-1*clip.trapFit.depth_frac,'ro')
plt.plot(phi,bestmodel,'b.')
plt.xlim([clip.thermal.phimin-dur,clip.thermal.phimax+dur])
plt.ylim([-1.1*clip.trapFit.depth_frac,0.05*clip.trapFit.depth_frac])
plt.title(epic)
#plt.figure()
#plt.plot(time,isbad,'rx')*-1*clip.trapFit.depth_frac
#plt.plot(time,bestmodel,'b-.')
#plt.title(epic)
#%%
import dave.pipeline.gather as g
cliploc='/soc/nfs/so-nfs/dave/c6-v2/'
sfile='/soc/nfs/so-nfs/dave/c6-v2/c6Candidates.txt'
clipar=np.loadtxt(sfile,dtype='str',usecols=[0],skiprows=0)
cliplist=list()
for v in clipar:
cliplist.append("%s/%s/c%s-06.clip" % (cliploc,v,v))
#%%
import dave.susanplay.thermalFlags as tf
fun=tf.getFractionThermal
(epics,values)=g.gatherFunctionB(cliplist[0:20],fun)