Taylor series of a function of a matrix

Taylor series of a function of a matrix

What is the approximation of the following function : $$ \det(({\mathbf
A}+ \epsilon {\mathbf B})^H({\mathbf A}+ \epsilon {\mathbf B}))$$ when
$\epsilon <<1$ ?
(Note that ${\mathbf A}$ and ${\mathbf B}$ are $M \times N $ matrices.)