  Galactic PeVatrons are sources that can accelerate cosmic rays to PeV
energies. The high-energy cosmic rays are expected to interact with the
surrounding ambient material or radiation, resulting in the production of gamma
rays and neutrinos. To optimize for the detection of such associated production
of gamma rays and neutrinos for a given source morphology and spectrum, a
multi-messenger analysis that combines gamma rays and neutrinos is required. In
this study, we use the Multi-Mission Maximum Likelihood framework (3ML) with
IceCube Maximum Likelihood Analysis software (i3mla) and HAWC Accelerated
Likelihood (HAL) to search for a correlation between 22 known gamma-ray sources
from the third HAWC gamma-ray catalog and 14 years of IceCube track-like data.
No significant neutrino emission from the direction of the HAWC sources was
found. We report the best-fit gamma-ray model and 90% CL neutrino flux limit
from the 22 sources. From the neutrino flux limit, we conclude that the
gamma-ray emission from five of the sources can not be produced purely from
hadronic interactions. We report the limit for the fraction of gamma rays
produced by hadronic interactions for these five sources.
