Welcome @veeus18!
The main variation your listed steps have from the current recommendations are the response function estimation and FOD calculation algorithm. The dhollander algorithm is generally recommended for response function estimation. "dwi2fod csd invokes the original CSD algorithm, which involves a “soft” non-negativity constraint (negative values are penalised) rather than the “hard” non-negativity constraint (negative values are forbidden) of the MSMT CSD algorithm. Despite the naming, the “multi-shell” algorithm can still be used with what’s considered “single-shell” data: with two unique b-values the data can be decomposed into two tissues, typically WM and CSF.
As far as connectome construction using a volume-based atlas, take a look at my BIDS App and see if your target parcellation is already supported there. I’ve not yet been able to finalise and release a decent volume of changes, so it’s a little out of date, but new parcellations are not too difficult to add. Otherwise you can check out the relevant wiki post which explains the underlying logic & steps.
Cheers
Rob