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)
Exemple #2
0
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')
Exemple #3
0

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')
Exemple #4
0
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)