def METHOD_NAME(self, fs=30, functional_channel=0, nbaseline_frames=30, k=200, nchannels=2):
    # MOTION CORRECTION
    wfield_cli._motion(str(self.data_path), nchannels=nchannels, plot_ext='.png')
    # COMPUTE AVERAGE FOR BASELINE
    wfield_cli._baseline(str(self.data_path), nbaseline_frames, nchannels=nchannels)
    # DATA REDUCTION
    wfield_cli._decompose(str(self.data_path), k=k, nchannels=nchannels)
    # HAEMODYNAMIC CORRECTION
    # check if it is 2 channel
    dat = wfield_cli.load_stack(str(self.data_path), nchannels=nchannels)
    if dat.shape[1] == 2:
        del dat
        wfield_cli._hemocorrect(str(self.data_path), fs=fs, functional_channel=functional_channel, plot_ext='.png')