Geochemical anomaly separation using the concentration–number (C–N) method at the Haftcheshmeh porphyry system in NW Iran is the aim of this study. We used lithogeochemical data sets to explore Cu, Mo, Au and Re mineralization in gabbroic, dioritic and monzonitic units at the Haftcheshmeh Cu–Mo porphyry system. The obtained results were interpreted using a rather extensive set of information available for each mineralized area, consisting of detailed geological mapping, structural interpretation and alteration data. Threshold values of elemental anomalies for the mineralized zone were computed and compared with the statistical methods based on the data obtained from chemical analyses of samples for the lithological units. Several anomalies at local scale were identified for Cu (40 ppm), Mo (12 ppm), Au (79 ppb), and Re (0.02 ppm), and the results suggest the existence of local Cu anomalies whose magnitude generally is above 500 ppm. The log–log plots show the existence of three stages of Cu and Mo enrichment, and two enrichment stages for Au and Re. The third and most important mineralization event is responsible for presence of Cu at grades above 159 ppm. The identified anomalies in Haftcheshmeh porphyry system, and distribution of the rock types, are mainly gabbrodiorite–monzodiorite, granodiorite and monzodiorite–diorite that have special correlation with Cu–Mo and gabbroic and monzonitic rocks, especially the gabbrodiorite–monzodiorite type, which is of considerable importance. The study shows that these elemental anomalous parts have been concentrated dominantly by potassic and phyllic, argillic and propylitic alterations within the gabbroic, monzonitic and dioritic rocks especially in the gabrodioritic type in certain parts of the area. The results, which were compared with fault distribution patterns, revealed a positive correlation between mineralization in anomalous areas and the faults present in the mineralized system.