from qsiprep.workflows.dwi.hmc import init_dwi_model_hmc_wf
wf = init_dwi_model_hmc_wf(modelname='3dSHORE',
                           transform='Affine',
                           num_iters=2)