The irradiation of sub-regions of the parotid has been linked to xerostomia development in patients with head and neck cancer (HNC). In this study, we compared the xerostomia classification performance of radiomics features calculated on clinically relevant and de novo sub-regions of the parotid glands of HNC patients. All patients (N = 117) were treated with TomoTherapy in 30-35 fractions of 2-2.167 Gy per fraction with daily mega-voltage-CT (MVCT) acquisition for image-guidance purposes. Radiomics features (N = 123) were extracted from daily MVCTs for the whole parotid gland and nine sub-regions. The changes in feature values after each complete week of treatment were considered as predictors of xerostomia (CTCAEv4.03, grade ≥ 2) at 6 and 12 months. Combinations of predictors were generated following the removal of statistically redundant information and stepwise selection. The classification performance of the logistic regression models was evaluated on train and test sets of patients using the Area Under the Curve (AUC) associated with the different sub-regions at each week of treatment and benchmarked with the performance of models solely using dose and toxicity at baseline. In this study, radiomics-based models predicted xerostomia better than standard clinical predictors. Models combining dose to the parotid and xerostomia scores at baseline yielded an AUC test of 0.63 and 0.61 for xerostomia prediction at 6 and 12 months after radiotherapy while models based on radiomics features extracted from the whole parotid yielded a maximum AUC test of 0.67 and 0.75, respectively. Overall, across sub-regions, maximum AUC test was 0.76 and 0.80 for xerostomia prediction at 6 and 12 months. Within the first two weeks of treatment, the cranial part of the parotid systematically yielded the highest AUC test . Our results indicate that variations of radiomics features calculated on sub-regions of the parotid glands can lead to earlier and improved prediction of xerostomia in HNC patients.