from herschel.hifi.pipeline.wbs import *
hdf=AccessDataFrameTask()(obsid=obsid,apid=1030)
spectra=WbsSpectrumDataset(hdf) 
DoWbsScanCountTask()(spectra) #step 1
calPixel=MkWbsBadPixelsTask()(spectra) #step 2
MkWbsBadPixelsTask()(spectra,calPixel) #step 3 
DoWbsDarkTask()(spectra) #step 4 
DoWbsNonlinTask()(spectra) #step 5 
zeroCal=MkWbsZeroTask()(spectra) #step 6 
DoWbsZeroTask()(spectra,zeroCal) #step 7
freqCal=MkWbsFreqTask(spectra) #step 8
DoWbsFreqTask()(spectra,freqCal) #step 9
att=MkWbsFluxAttenTask()(spectra) #step 10
DoWbsSubbandsTask()(spectra)