def findInjTransits(clip):
"""
Now that the light curve has a transit in it, run the rest of the pipeline.
dpp.lppMetricTask
dpp.modshiftTask dpp.measureDiffImgCentroidsTask dpp.dispositionTask
dpp.saveClip
"""
clip=dpp.detrendDataTask(clip)
clip=dpp.fblsTask(clip)
clip=dpp.trapezoidFitTask(clip)
clip=dpp.dispositionTask(clip)
clip=dpp.lppMetricTask(clip)
clip=dpp.saveClip(clip)
return clip
acti=i+span[0];
cfg['campaign'] = int(rundata[acti,1])
print epicid,cfg['campaign']
clip=mS.runOne(epicid, cfg)
#plt.plot(clip.bls.bls_search_periods,clip.bls.convolved_bls,'r')
try:
f=interp1d(clip.bls.bls_search_periods,clip.bls.convolved_bls)
ibls=f(periods)
#bls[n,:]=ibls;
clip.bls['ibls']=ibls;
clip.bls['iperiods']=periods
clip.bls['kepmag']=rundata[acti,2]
clip=pipe.saveClip(clip)
n=n+1;
except AttributeError:
pass
#np.save('/home/sthomp/DAVE/playK2/ch42-bls.npy', bls)
#%%
#Code to gather up the bls spectra into a 2d array belongs here.
import dave.pipeline.clipboard as cb
listdir='/home/sthomp/daveOutput/bls/list.txt'
clipnames=np.loadtxt(listdir,dtype='string')
bls=np.zeros((len(clipnames),25000))
kepmag=np.zeros((len(clipnames),1))
for i,v in enumerate(clipnames):
c=cb.loadClipboard(v)
bls[i,:]=c.bls.ibls;
