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;