We apply the recently proposed nonrelativistic EFT framework for double heavy hadrons to the double heavy baryon case. The EFT is build from NRQCD by incorporating the adiabatic expansion between the light quark and the heavy quark pair. At leading order the EFT reduces to the Born-Oppenheimer approximation. The Born-Oppenheimer potentials are obtained from available lattice QCD data. We go beyond the leading order by incorporating the heavy quark spin and angular momentum operators. The corresponding potentials are obtained from their Wilson loop matching expressions by computing these in the short and long distances and interpolating the results for the intermediate region. In the short distance this is done using the multipole expansion. For the long distance we have developed an Effective String Theory with a light fermion constrained to the be in the flux tube string. Several free parameters are obtained by comparing with the spin splittings of double charm baryon lattice data and then used to predict double bottom baryon as well as excited double charm states.