Purpose: To investigate the effect of patient specific vessel cooling on head and neck hyperthermia treatment planning (HTP). Methods and materials: Twelve patients undergoing radiotherapy were scanned using computed tomography (CT), magnetic resonance imaging (MRI) and contrast enhanced MR angiography (CEMRA). 3D patient models were constructed using the CT and MRI data. The arterial vessel tree was constructed from the MRA images using the ‘graph-cut’ method, combining information from Frangi vesselness filtering and region growing, and the results were validated against manually placed markers in/outside the vessels. Patient specific HTP was performed and the change in thermal distribution prediction caused by arterial cooling was evaluated by adding discrete vasculature (DIVA) modeling to the Pennes bioheat equation (PBHE). Results: Inclusion of arterial cooling showed a relevant impact, i.e., DIVA modeling predicts a decreased treatment quality by on average 0.19 °C (T90), 0.32 °C (T50) and 0.35 °C (T20) that is robust against variations in the inflow blood rate (|ΔT| < 0.01 °C). In three cases, where the major vessels transverse target volume, notable drops (|ΔT| > 0.5 °C) were observed. Conclusion: Addition of patient-specific DIVA into the thermal modeling can significantly change predicted treatment quality. In cases where clinically detectable vessels pass the heated region, we advise to perform DIVA modeling.