We present an impedance model of the Fermilab Recycler ring using PyHEADTAIL. The model is constructed by incorporating analytical expressions for the wakefields of beamline components that contribute significantly to impedance. The effects of indirect space charge are included as an inductive impedance. Benchmarking against measured coherent Betatron tune shifts, the impedance model is found to capture 73.4% of observed tune shifts. Our findings serve as a stepping stone for the development of a realistic impedance model crucial for studying impedance-driven instabilities at higher intensity.