  Star-forming and starburst galaxies (SFGs and SBGs) are considered to be
powerful emitters of non-thermal gamma-rays and neutrinos. On this regard, the
Fermi-LAT collaboration has found a correlation between the gamma-ray and
infrared luminosities for a sample of local sources. Yet, the physics behind
these non-thermal emissions is still under debate. Using the publicly-available
fermitools, we analyse 15.3 years of gamma-ray between 1-1000 GeV data for 70
sources, 56 of which were not previously detected. We relate this emission to a
physically-motivated model for SBGs in order to constrain Fcal for each source
and then study its correlation with the star formation rate of the sources. We
find at 4sigma level an indication of gamma-ray emission for other two SBGs,
namely M 83 and NGC 1365. By contrast, we find that, with the new description
of background, the significance for the gamma-ray emission of the previously
detected M 33 reduces at $4\sigma$. We also find that the physically-motivated
model for Fcal correctly describes the $\gamma$-ray observations, and it is
crucial to assess the systematic uncertainty on the star formation rate.
Finally, undiscovered sources strongly constraints Fcal at 95\% CL, providing
fundamental information when we interpret the results as common properties of
SFGs and SBGs. Hence, we find that these sources might contribute 12% to the
EGB, while the corresponding diffuse neutrino flux strongly depends on the
spectral index distribution along the source class.
