{'fraction_non_alphanumeric': 0.03326238677424145, 'fraction_numerical': 0.026151231665264233, 'mean_word_length': 4.594941382811779, 'pattern_counts': {'":': 0, '<': 15, '<?xml version=': 0, '>': 7, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
The ACS Survey of Galactic Globular Clusters. VII. Relative Ages

The ACS Survey of Galactic Globular Clusters is a Hubble Space Telescope (HST) Treasury program designed to provide a new large, deep and homogeneous photometric database. Based on observations from this program, we have measured precise relative ages for a sample of 64 Galactic globular clusters by comparing the relative position of the clusters' main sequence turn offs, using main-sequence fitting to cross-compare clusters within the sample. This method provides relative ages to a formal precision of 2-7%. We demonstrate that the calculated relative ages are independent of the choice of theoretical model. We find that the Galactic globular cluster sample can be divided into two groups -- a population of old clusters with an age dispersion of ~5% and no age-metallicity relation, and a group of younger clusters with an age-metallicity relation similar to that of the globular clusters associated with the Sagittarius dwarf galaxy. These results are consistent with the Milky Way halo having formed in two phases or processes. The first one would be compatible with a rapid (<0.8 Gyr) assembling process of the halo, in which the clusters in the old group were formed. Also these clusters could have been formed before reionization in dwarf galaxies that would later merge to build the Milky Way halo as predicted by Lambda-CDM cosmology. However, the galactocentric metallicity gradient shown by these clusters seems difficult to reconcile with the latter. As for the younger clusters, it is very tempting to argue that their origin is related to their formation within Milky Way satellite galaxies that were later accreted, but the origin of the age-metallicity relation remains unclear.


Introduction
One of the keys to understanding the structure and evolution of the Milky Way is the ability to divide its stars and clusters into separate Galactic Populations as first hinted at by Oort (1926) and formally introduced by Baade (1944). In the early 1950s, Sandage published the first deep color−magnitude diagrams (CMDs) of the northern Galactic globular clusters (GGCs) M3 and M92 (Sandage 1953;Arp, Baum, & Sandage 1953) and discovered the faint sequence of stars similar to the main sequence of Population I. It was almost immediately realized that GGCs are ancient compared to the stars in the Solar Neighborhood.
Over the last half century, it has become clear that GGCs are astronomical fossils providing important information on the formation and evolution of the Milky Way. One characteristic of GGCs that makes them so valuable is that they are the oldest Galactic objects for which reliable ages can be measured.
GGC ages are of great interest to a variety of cosmological and comogonical issues since their absolute ages place a lower limit on the age of the Universe. Unfortunately, absolute age measurements are still compromised by a number of uncertainties, particularly GGC distances and foreground reddenings as well as metallicities, bolometric corrections and some poorly known aspects of the stellar evolution. Relative ages are almost (although not totally) free of these effects and can still reveal fundamental information about Galaxy formation mechanisms and time scales. For this reason they have become of paramount importance in the context of the study of GGC.
The study of GGC ages has evolved apace with improvements in instrumentation, in our theoretical understanding of stellar evolution and with the advance of computational facilities. Reliable age determination requires deep photometry that reaches at least the main−sequence turn off (MSTO), together with a good understanding of the stellar evolution theoretical modeling. During the first half of the 20 th century, the available instrumentation did not allow observations deeper than the horizontal branch of most GGCs. Over the past three decades, however, low quantum efficiency non−linear photographic imaging has given way to digital imaging using high quantum efficiency linear CCDs. This has resulted in ever deeper and more precise cluster CMDs, that reached the MSTO and allowed the application of age dating techniques. First relative age studies were carried out by Gratton (1985), Peterson (1987), Sarajedini & King (1989) and Sarajedini & Demarque (1990), as representative examples. However, uncertainties in both the input physics (Chaboyer et al. 1998) and the color−T eff transformation of the theoretical models into the observational plane (Buonanno et al. 1998) constitute a significant source of errors in the final relative age determination, particularly when ages of clusters with different metallicities are compared. In the last decade, many efforts have been made in GGC relative ages, improving our understanding with each technological or theoretical advance (e. g. Chaboyer, Demarque, & Sarajedini 1996;Richer et al. 1996;Stetson et al. 1996;Salaris & Weiss 1998;Buonanno et al. 1998;Rosenberg et al. 1999;VandenBerg 2000;De Angeli et al. 2005).
Relative ages can be estimated to high precision by measuring the position of the MSTO relative to CMD features that have little or no age dependence (Stetson et al. 1996;Sarajedini et al. 1997). However, a key aspect of any relative age experiment is having a large homogenous database of globular cluster photometry. Such databases have only recently appeared in the literature but have significantly advanced our understanding of Milky Way formation and evolution. The largest dedicated efforts are those of Rosenberg et al. (2000a,b), a ground-based survey of 35 nearby GGCs (Rosenberg et al. 1999), andDe Angeli et al. (2005), an analysis of 55 clusters from Piotto et al. (2002)'s HST snapshot catalog. These and other recent relative age studies (Buonanno et al. 1998;Salaris & Weiss 1998;Rosenberg et al. 1999;VandenBerg 2000;De Angeli et al. 2005) consistently point towards a scenario in which metal−poor clusters are largely coeval but metal−rich clusters show some age dispersion, ∼1.0 Gyr (rms), and a possible age−metallicity relation. No correlation between age and Galactocentric distance has been conclusively demonstrated.
In this paper, we present an analysis of a uniform set of HST/ACS data for 64 GGCs 2 . This dataset, previously introduced by Sarajedini et al. (2007), has produced well-defined CMDs for all of the target clusters from the tip of the red giant branch (RGB) to at least ∼6.5 magnitudes below the MSTO . The determination of GGC ages was a principle goal of the HST/ACS survey and exposures times were defined to maximize the S/N ratio at the MSTO. This makes the survey an excellent resource for the evaluation of relative ages.
Due to the difficulty in determining the position of the MSTO 3 to the high precision required for age studies, many different methods have been developed in previous relative ages works (see Meissner, & Weiss 2006, for a complete discussion). The two most used relative age indicators are those based on the color difference between the MSTO and the RGB at a given magnitude level (horizontal method), and the magnitude difference between the MSTO and the zero−age HB 4 (vertical method). VandenBerg, Bolte, & Stetson (1990) were pioneers in measuring GGC relative ages making use of the horizontal method. They essentially derived relative ages comparing the color difference between the MSTO and the location of a well defined point in the lower RGB, located 2.5 magnitudes brighter than a point in the upper main−sequence (MS) that is 0.05 magnitudes redder than the MSTO. However, horizontal method results are sensitive to the assumed value of the mixing length parameter and to the color−T eff transformation of the theoretical models into the observational plane. As a result, relative ages determined using the horizontal method are model dependent. Potentially, more reliable age indicators are those related to the brightness of the MSTO, in particular the vertical method. First surveys of GGC relative ages using the vertical method were carried out by Gratton (1985) and Peterson (1987). They measured relative ages from a sample of 26 and 41 GGCs, respectively, collected from the available literature. Later on, Sarajedini & King (1989) used the brightness difference between the HB and the MSTO to estimate the ages of 31 GGCs. Nevertheless, accurate determination of the zero−age HB is delicate, in particular for the clusters at the extreme boundaries of the metallicity distribution, which generally have a red (or blue) HB and no RR Lyrae stars. Moreover, the dependence of the zero−age HB luminosity on the cluster metallicity is still controversial.
In this paper, we use an alternative approach. Our data reaches at least ∼6.5 magnitudes below the MSTO, providing a very well defined MS. This allows, for the first time, the use of relative MS fitting for a large number of clusters. As it will be shown later, the relative MS fitting procedure used in this study does essentially the following. The MS of two clusters, with similar metallicities, are superimposed in the CMD by shifting one of the cluster's mean ridge line in both color and magnitude (explained in detail in section 4). This way, any differences between the two cluster's distance and reddening are compensated, and an intrinsic MSTO magnitude difference is obtained. This relative MSTO brightness is then used to derive relative ages. This allows us to use what is, in essence, an improved version of the traditional vertical method that substitutes the well-defined MS for the contentious zero−age HB.
We also take the additional step of comparing our observational measures to updated theoretical models. We have primarily used the most recent state−of−art models of Dotter et al. (2007) (hereafter D07), transformed to the observational plane with an empirical Color−T eff transformation. However, we also compare our results to the models of Pietrinferni et al. (2004) (hereafter P04), Bertelli et al. (1994) (hereafter B94) and Girardi et al. (2000) (hereafter G00). We find that the various theoretical models agree closely on the MSTO absolute magnitude and color, resulting in negligible model dependence.
This paper is organized as follows. Section 2 describes the observations and data reduction pipeline, and Section 3 presents the globular cluster database. The mean ridge line fitting and MS−fitting procedures used in this study are described in Section 4. The analysis of the uncertainties is presented in Section 5. Finally, relative age results and a discussion are presented in Sections 6 and 7, respectively. The paper ends with a summary of conclusions (Section 8).

Observations and data reduction
Observations were performed with the ACS/WFC instrument on board the HST . A total of 65 GGCs were observed in the F 606W (∼ V ) and F 814W (∼ I) filters. Due to its high differential reddening, Pal 2 is not considered here.
Details on the observations, data reduction and photometry are given in Sarajedini et al. (2007) and Anderson et al. (2008). For most globular clusters, the observation and data reduction pipeline produces ∼12 magnitudes of precise photometry, extending from nearly the tip of the RGB to several magnitudes below the MSTO. In general, the photometry reaches a depth corresponding to an approximate stellar mass of ∼ 0.2 m ⊙ . This uniform and deep photometry offers a database of unprecedented quality, opening a window for new data analysis approaches. In particular, cluster MSs are so well defined that it is possible, for the first time, to apply relative MS−fitting techniques to a large number of clusters within the same database.

The database
The position of the MSTO depends not only upon age but also upon metallicity. This effect requires the use of a large and homogeneous metallicity database in order to determine accurate relative ages. To account for this, we have adopted the metallicities listed in the Rutledge et al. (1997) catalog 5 over both the Carretta & Gratton (1997) (hereafter CG) and the Zinn & West (1984) (hereafter ZW) metallicity scales. This catalog is maybe not as precise as a true compilation of metallicities based on high dispersion spectra, but it has the the advantage of being homogeneous. For clusters not listed in Rutledge et al. (1997), the original metallicity listed in the ZW catalog has been transformed to the CG metallicity scale using the following equation (Carretta et al. 2001 Lyngå 7, E3, and Pal 1 metallicities are not listed in the previous catalogs. Lyngå7 and E3 metallicities in the CG scale have been calculated from the ZW metallicities of Tavarez & Friel (1995) and Harris (1996), respectively. For Pal 1, we used the measure of Rosenberg et al. (1998). Table 1 presents the list of our target GGCs, grouped by metallicity. Column 3 gives the cluster identity while columns 4 and 5 present the adopted metallicities in the ZW and CG scales, respectively. Column 6 lists the Galactocentric distance (r GC , in kpc), taken from Harris (1996).

Mean ridge lines and main sequence fitting method
In this section, we describe the general properties of the cluster mean ridge line fitting procedure. We then use the derived ridge lines to perform relative MS−fitting between clusters.

Mean ridge lines and main sequence turn-off apparent magnitude
In order to determine accurate relative ages for each cluster, a precise determination of the MS and RGB mean ridge line is required. With this aim, we have created a new software package for the ACS program. It is worth mentioning that mean ridge lines are determined in the observational plane (F 606W -F 814W , F 814W ) instead of (F 606W -F 814W , F 606W ), as the SGB appears more vertical in the CMD and the following procedure produces more accurate results. For the rest of the analysis, the (F 606W -F 814W , F 606W ) plane is used.
For each globular cluster, the mean ridge line is determined in two steps ( Figure 1). First, color histograms are computed for stars grouped in 0.4 magnitude (F 814W ) wide bins. The value of 0.4 magnitudes was determined experimentally as a compromise between two effects -the low resolution produced by larger bins and the noisy ridge line produced by smaller bins. Color histograms are constructed along the RGB and MS with a step of 0.04 magnitudes (a tenth of the bin size) in F814W. That is, a moving bin of fixed width (0.4 magnitudes), with successive steps of 0.04 magnitudes, is used. A preliminary mean ridge line is then derived with (F 606W − F 814W ) colors set at each histogram's maximum, and the F 814W magnitude of the bin center.
Once the preliminary mean ridge line has been determined, the fitting process is iterated three more times. However, we now use rotated histograms. For each F 814W value, the tangential angle of the mean ridge line is determined. This angle is used to rotate the reference system on which the histogram is computed. The X−axis of the new reference system is perpendicular to the preliminary mean ridge line at the considered F 814W magnitude. The histogram is then computed perpendicular to the preliminary mean ridge line and considering stars in the new reference system's Y−axis interval [-0.04, 0.04]. The histogram's maximum is determined, defining a point of the new mean ridge line in the rotated reference system, (max ,0). Note that this point shows the correction to be applied to the preliminary mean ridge line. The coordinates of this point are de −rotated and its location in the original reference system, the CMD, is calculated. The process is repeated for all the F 814W values to derive the revised mean ridge line, which is then smoothed to derive the final mean ridge line. Figure 2 shows some examples of mean ridge line fitting results. Four selected clusters covering the metallicity range are shown.
The empirical MSTO is the bluest point on the mean ridge line. It is worth mentioning that the final mean ridge line consists of a discrete number of points, it is not a continuous line. For this reason, in order to determine the MSTO position, a continuous spline has been fitted to the main−sequence−to−subgiant−branch portion of the final mean ridge line. Its bluest point is then adopted as the empirical MSTO. As an example, the position of the MSTO is marked with an open circle in Figure 2. In order to increase the precision of the MSTO position, the mean ridge lines have been determined ten times for each cluster, with the starting point of the histograms offset in steps of 0.004 magnitudes. The MSTO color and apparent magnitude have been determined as the mean of the ten different measurements.
This method provides a consistently defined mean ridge line, minimizing the effects of the binary population and differential reddening (see section 5). Additionally, the rotated coordinate system ensures that the best fit is derived around the steeply curved MSTO and sub giant branch region of the CMD. For high metallicity clusters, in particular, the sub giant branch is nearly horizontal and the more commonly used mean ridge line determination method (drawing an approximate ridge line, rendering the sequence to a nearly vertical line and then taking histograms in small ranges of magnitude) lacks the resolution around the sub giant branch needed for precise age measurement.

Main sequence fitting method
The CMD location of the faint MS of a GGC is independent of its age, but highly dependent upon its metallicity. For this reason, we have divided the HST/ACS sample into six metallicity groups (see Table 1) using the following divisions: [Fe/H] = -0.3, -0.8, -1.1, -1.3, -1.5, -1,8 and -2.3. For each group, a reference cluster has been chosen: NGC 6304, NGC 6723, NGC 6981, NGC 6681, NGC 6101 and NGC4590, respectively. Reference clusters have been selected to have low differential reddening, low field star contamination, and a well defined MS. Using the derived mean ridge lines, we performed MS− and RGB−fitting between each cluster in each metallicity group and the corresponding reference cluster. Figure 3 shows examples of the fitting, in which the mean ridge lines of clusters with −1.3 ≤ [F e/H] CG < −1.1 (solid lines) have been shifted in both magnitude and color to fit the reference cluster, NGC 6981 (dashed line). The fitting was performed in a least squares fashion and taking into account two CMD regions that have little dependence upon cluster age. The two intervals are [(M TO F606W − 2.5) < M F606W < (M TO F606W − 1.5)] and [(M TO F606W + 1.5) < M F606W < (M TO F606W + 3.0)], M TO F606W being the magnitude of the reference cluster's MSTO. These two intervals have been defined making use of D07 theoretical isochrones. In particular, D07 isochrones with similar metallicity and different ages were superimposed on the same CMD. Based on visual inspection, these two particular regions were found to have little dependence upon cluster age, and were adopted as the optimum intervals for the MS fitting procedure. These regions are shaded in Figure 3. This fitting procedure produces unequivocal results, with no degeneracy between distance modulus and reddening.
The MS−fitting provides relative MSTO magnitudes for clusters inside each metallicity group. In order to calculate relative MSTOs between two consecutive metallicity groups, we applied the MS−fitting technique between the corresponding reference clusters. This was repeated for every two adjacent metallicity groups, producing a unified sample of relative MSTO magnitudes for all the analyzed clusters.
In order to compare relative MSTO measures to model predictions, we had to calculate the absolute magnitudes of the turnoffs. As the relative MSTO magnitudes are known, the MSTO absolute magnitude is only needed for one cluster. We chose the well-measured NGC 6752 cluster, adopting a subdwarf-based distance modulus of (m − M) V = 13.24 ± 0.08 from Gratton et al. (2003). We calculate an MSTO magnitude of M T O F 606W = 4.02 ± 0.08 for NGC 6752.
We performed similar calculations using the NGC 0104 and NGC 6397 distance moduli listed in Gratton et al. (2003) of (m − M) V = 13.50 ± 0.08 and 12.58 ± 0.08 , respectively. We calculate MSTO absolute magnitudes of M T O F 606W = 4.03±0.08 and M T O F 606W = 3.77±0.08 for NGC 0104 and NGC 6397, respectively. As the relative MSTOs are known, NGC 6752's MSTO absolute magnitude can be also derived independently from NGC 0104 and NGC 6397 (3.72 ± 0.08 and 3.86 ± 0.08, respectively). The mean of the three NGC 6752 MSTO absolute magnitude determinations is M T O F 606W = 3.87 ± 0.15 and we adopt this as the final MSTO zeropoint.
We note that any error in the MSTO magnitude zeropoint would affect absolute age determination, but not relative age measurements. For this reason, the zeropoint uncertainty is not taken into consideration during relative age determination.
Using the adopted zeropoint, we calculate MSTO absolute magnitudes for our entire sample, which are listed in column 4 of Table 4. The listed uncertainties are determined in next section.

Main sequence turn off uncertainties
Several effects must be taken into account when deriving the uncertainty associated with the measured MSTO magnitudes (M TO F606W ). Among these, we address the cluster differential reddening, the binary fraction and total number of cluster stars. The total number of stars affects the MSTO magnitude measurements in two ways. The first one is related to the total number of stars close to the MSTO that can be used to measure its magnitude: the smaller this number, the more uncertain the MSTO magnitude will be. On the other hand, a large number of cluster stars implies more severe crowding, and therefore less accurate magnitude measurement. To these uncertainties, we need to add the uncertainty from the MS−fitting and the mean ridge line determination procedure itself. The final uncertainty associated to the MSTO magnitude, σ 2 (M TO F606W ), is: where σ 2 RED is the variance due to the differential reddening, σ 2 BIN due to the binary star population, σ 2 NUM due to the number of cluster stars, σ 2 CROW due to the crowding, σ 2 MRL is the MSTO's magitude variance from the mean ridge line determination, and σ 2 MSF from the MS−fitting method. It is worth noting that field star contamination does not affect the MSTO determination. Our fields cover the central region in each cluster and the ratio of cluster stars to field stars is very high. Moreover, as noted above, the mean ridge line fitting method fits the dominant population, ignoring both field stars and blue stragglers.
As it has been described in Section 4.1, each MSTO's magnitude has been determined as the mean of the MSTO's magnitude obtained from ten different mean ridge line fittings. Here we assume that the variance associated with the mean ridge line determination, σ 2 MRL , is the variance of the previous ten measurements. Typical values of σ MRL are in the range ∼0.01-0.04 mag.
With the aim of performing an external cross-check for the σ 2 MRL , the MSTO's magnitude variance from to the mean ridge line determination, we have considered four clusters having low reddening and a large number of stars (NGC 0104, NGC 5024, NGC 5272 and NGC7089). These clusters' CMDs have been randomly divided into 5 partial CMDs. A mean ridge line has been fitted to each partial CMD and its MSTO's magnitude has been measured. Finally, σ/ √ 5 − 1 for these five MSTO's measurements has been computed for each cluster, where σ is the standard deviation of the five measurements.
Results are shown in Table 2. Column 1 shows the cluster's name, column 2 lists σ MRL and column 3 lists σ cross−check = σ/ √ 5 − 1 derived in this external cross-check. It can be seen that the uncertainties obtained here are equivalent to the derived σ MRL . This test indicates that the adopted σ MRL is a good estimation.
For σ 2 MSF , we adopted the uncertainty calculated by the least squares main sequence fitting procedure. Typical values of σ MSF are in the interval ∼0.01-0.06 mag.
To estimate the uncertainties induced by differential reddening, binary star population, and total number of cluster stars, as well as by differences between the measured MSTO and the actual GGC's MSTO, we generated over 250 synthetic CMDs using the IAC−star (Aparicio & Gallart 2004) synthetic CMD program. Stellar populations are calculated on a star-by-star basis. L, T eff , and g of each star are computed by direct bilogarithmic interpolation in both the metallicity and age grids of a library of stellar evolution tracks, metallicity and age following continuous distributions. All the stellar evolution phases are covered accurately, and mass loss is computed during the RGB and the asymptotic giant branch phases (see Aparicio & Gallart 2004, for details). A set of different star formation rates, initial mass functions and metallicity laws are allowed, and binary star content can be computed. The result is a synthetic CMD with a smooth realistic stellar distribution.
Synthetic CMDs were generated with an age of 13 Gyr, [M/H] = −1.2, differential reddening E(F 606W − F 814W ) values between 0.0 and 0.2, binary fractions (with 0.5 < m 2 /m 1 < 1.0) between 0.0 and 0.7, and a number of stars brighter than M F606W = 6 (approximately two mag. below the MSTO) between 1,000 and 50,000. The P04 stellar evolution library was used to generate the stellar libraries. For each combination of differential reddening, binary fraction and number of stars, 10 different CMDs were generated. Observational errors were simulated in all the synthetic CMDs by applying a gaussian dispersion to the synthetic CMD, with σ taken from the dispersion of the observed CMD as a function of the magnitude.
Using the same mean ridge line fitting method used for the observational data, MSTO's magnitudes were measured for the synthetic CMD. The results are shown in Figure 4. Grey points show the MSTO measurements from the 10 different CMDs generated for each experiment, black points show the mean, and the error bars correspond to the standard deviation of the 10 measurements. Solid line represents the input MSTO magnitude, that is, the magnitude of the MSTO obtained directly from the input isochrone. It can be seen that a systematic effect on the measured MSTO's magnitudes exists. Output MSTOs (measured) are systematically ∼0.02 mag. fainter than input ones. This would have an effect on the measured absolute ages, but as this work is focused on relative ages, the variance of the difference between the input and output MSTO's magnitudes, σ 2 IO , will be also considered. The left panel of Figure 4 shows the M TO F606W measurements for the synthetic CMDs with no differential reddening, no binary stars and varying the number of stars brighter than M F606W = 6 from 1,000 to 50,000. With the aim of determining σ 2 NUM , a gaussian has been generated for each black point, being the mean and the sigma of the gaussian equal to the value of the corresponding point and its sigma, respectively. Finally, all gaussians have been added together and the variance of this sum of gaussians has been adopted as σ 2 NUM . Since most of the target clusters have more than 1,000 stars brighter than F 606W = 6, we adopt σ NUM = 1.6 × 10 −2 mag. (dotted lines). Estimated this way, σ 2 NUM also includes the σ 2 IO contribution.
The central panel of Figure 4 shows the effect of differential reddening upon the estimated M TO F606W value. For a synthetic CMD, the reddening has been simulated by shifting every star along the reddening vector a quantity between zero and a maximum differential reddening value, E(F 606W − F 814W ) max . We generated a set of nine synthetic CMDs with 50,000 stars brighter than F 606W = 6 and different values of E(F 606W − F 814W ) max from 0.0 to 0.2 magnitudes. It is worth noting that the black points' variance also includes the contribution from σ 2 IO . In order to avoid duplicity, the contribution of σ 2 IO must be eliminated. As these synthetic diagrams have 50,000 stars brighter than F 606W = 6, their σ 2 IO is equal to the variance of the point in the left panel corresponding to a number of stars brighter than F 606W = 6 equal to 50,000, that is σ IO = 9.7 × 10 −3 mag. In order to compute the effect of differential reddening upon the estimated M TO F606W value, σ 2 IO has been quadratically subtracted to each point's variance in the central panel of Figure 4. The resulting variance accounts for the differential reddening effect only. Then a gaussian has been generated for each black point, being the mean and the sigma of the gaussian equal to the value of the corresponding point and the obtained quadratically subtracted sigma, respectively. Finally, all gaussians have been added together and the variance of this sum of gaussians has been adopted as σ 2 RED . As we expect most of the target clusters to have differential reddening lower than 0.15 magnitudes, we adopt σ RED = 2.0 × 10 −2 mag.
The right panel of Figure 4 shows the M TO F606W measurements of the synthetic CMDs with 50,000 stars brighter than M F606W = 6, no differential reddening, and varying fractions of binary stars. It can be seen that the fraction of binary stars affects M TO F606W , making the measured MSTO fainter for clusters with a higher binary fractions. Following the same procedure as before, σ 2 IO has been quadratically subtracted to each point's variance in the right panel of Figure 4. Then a gaussian has been generated for each black point, being the mean and the sigma of the gaussian equal to the value of the corresponding point and the sigma obtained by subtracting σ 2 IO to the rms of the point, respectively. Finally, all gaussians have been added together and the variance of this sum of gaussians has been adopted as σ 2 BIN . In a parallel work, we are measuring the fraction of binary stars in each cluster in our database (Milone et al. 2008b). For most clusters, it is smaller than 20%. Accordingly, we consider the fraction of binary stars interval [0.0, 0.2] which results in σ BIN = 5.9 × 10 −3 mag.
To estimate the effect of cluster crowding on the measured MSTO magnitude, we used the artificial−star tests described in Anderson et al. (2008). In principle, crowding would produce some systematic shift between output and input magnitudes of artificial stars. We are interested in evaluating if these shifts significantly change as a function of crowding. Since different clusters have different crowding levels, determining the artificial star magnitude shifts in several of them will provide the information we need. We will use the NGC 0104, NGC 1851, NGC 2808, NGC 5139, NGC 5286, NGC 6171 and E3 clusters. This sample contains both one of the most and one of the least crowded clusters, namely NGC 2808 and E3, respectively.
Each cluster has associated a table containing input and output positions and magnitudes of the artificial stars. For each cluster, we selected artificial stars with input F 814W (I) magnitudes within ±0.25 mag of the MSTO. For most clusters, ∼ 4, 000 artificial stars fulfill this condition. For each of them, we calculated the quantity δI = I out − I in , where I out and I in are the recovered and input magnitudes, respectively. To have a representative estimate of δI and also of σ(δI) for each cluster, we have divided the artificial stars in 10 subsamples and calculated med(δI) for each of them. The average of medians, < med(δI) >, and the mean square root, σ[med(δI)], of the 10 med(δI) values have been computed for each cluster. Standard errors σ n = σ[med(δI)]/ (n − 1), with n = 10, are obtained. These can be considered as an estimate of the internal errors of med(δI) of each cluster. Average values < med(δI) > run from 3.0 ×10 −4 mag for E3, to −4.5 ×10 −3 mag for NGC2808, while σ n is always of the order of 1.0 × 10 −4 mag, indicating that the med(δI) are significant. In summary, we can say that ∼ 5 × 10 −3 mag is the maximum shift introduced by crowding in our MSTO measurements. Finally, we can estimate the uncertainty produced by crowding that we were looking for, σ CROW , as the root mean square of the < med(δI) > values of each cluster. This results in σ CROW = 1.5 × 10 −3 mag. Table 3 lists a summary of the different contributions to M TO F606W uncertainties. The MSTO's magnitude variances produced by the cluster number of stars, differential reddening, binary stars content and crowding, evaluated in this section using simulations, have been quadratically added and the resulting sigma (0.026 mag.) is also shown in Table 3. Finally, M TO F606W uncertainty has been evaluated for each cluster adding quadratically 0.026 mag. to the obtained σ 2 MRL and σ 2 MSF (equation 2), and the results are listed in Table 4 (column 4). It is worth mentioning that the total uncertainties tend to increase both at low and high metallicities. This is due to the error propagation in the σ 2 MSF contribution, as we are using the intermediate metallicity cluster NGC 6752 for the MSTO absolute magnitude zeropoint.

Results
In this section, we present the relative ages derived for the sample of 64 GGCs. We compare the derived ages to the theoretical isochrones of D07 using both the ZW and CG abundance scales. In the second part of this section, we test the MS−fitting procedure and check the self−consistency of the relative age determination. The relative age results are then compared with those derived using the stellar evolution libraries of P04, B94, and G00. Finally, we present a comparison of our results with previous work.

Relative ages
Relative ages were calculated using a stellar evolution library (D07) Table 4, respectively. We finally estimate the age of each cluster based on the interpolated curves.
While this procedure provides absolute ages for all of our clusters, it must be noted that absolute ages depend on the theoretical model as well as the adopted MSTO magnitude zeropoint. A much more detailed analysis is required to derive reliable absolute ages, and this will be done in a forthcoming paper. This paper concentrates only on relative ages, which are much less dependent on the MSTO zero point.
In order to derive relative ages, we divided the target clusters into three groups: low−metallicity Using the D07 models, we compute the mean age of the low−metallicity group as 12.80 ± 0.17 Gyr. The rms scatter of the 13 low−metallicity GGCs is equal to 0.6 Gyr. The absolute age of each cluster was then divided by this mean age, obtaining what we will call normalized age from here on. Columns 5 and 6 of Table 4 list the normalized ages from the D07 models using the ZW and CG metallicity scales, respectively. The adopted ± 0.05 dex uncertainty for the α-element enhancement translates into an uncertainty of ∼ 0.012 on the final relative age. With the aim of taking into consideration this effect on the relative ages, this quantity has been added quadratically to obtain the final relative age uncertainties listed in Table 4.

Testing self−consistency
In this section, the self−consistency of the MS−fitting procedure and relative age determination is tested. In particular, we are interested in quantifying the validity of matching up clusters with different metallicities on the grounds of theoretical stellar evolution models.
We created a set of 13 Gyr synthetic isochrones using the D07 models with metallicities similar to those of the selected reference clusters NGC 6304, NGC 6723, NGC 6981, NGC 6681, NGC 6101 and NGC4590. That is, isochrones with [M/H] = -0.5, -0.8, -1.0, -1.1, -1.5 and -1.8 have been considered. We applied the same MS−fitting procedure to the isochrones as for the observational data (section 4.2). The [M/H] = -1.1 isochrone, located in the center of the metallicity interval, was chosen for the MSTO zeropoint determination. Figure 6 summarizes the test results. Upper panel shows the MS−fitting for the 13 Gyr synthetic isochrones with the same metallicities as the adopted reference clusters. Lower panel shows the obtained relative ages results. The MSTO of the [M/H]=-1.1 isochrone, that has been adopted as MSTO zeropoint, has been encircled. Error bars represent the uncertainties derived from the MS−fitting, σ MSF . It can be seen that the main sequence fitting procedure is self consistent, and that the derived ages are similar to the input ages with a typical uncertainty of less than 2% over the entire metallicity range. We thus conclude that our metallicity grouping should not induce a significant bias in the relative age estimation. Figure 7 shows the theoretical MSTO magnitudes derived using P04, B94, and G00 stellar evolution libraries. D07 model predictions have been also plotted for comparison. Lines represent M TO F606W for different ages as a function of global metallicity, in 1 Gyr (solid lines) and 0.5 Gyr (dashed lines) age steps. Points represent our measured MSTO magnitudes as a function of metallicity in the CG scale. In the case of G00 models, metallicities lower than Z=0.0004 have been extrapolated.

Comparison to other stellar evolution libraries
Using these theoretical grids, normalized ages have been derived following the procedure described in section 6.1. The corresponding average age of the low metallicity group, as derived from each isochrone set, has been used for normalization.
Columns 7, 8 and 9 of Table 4 list the normalized ages derived from the P04, B94 and G00 stellar evolution libraries, respectively, using the CG metallicity scale. Figure 8 shows the difference between the normalized ages from the D07 models, and those from the P04, B94 and G00 models, adopting the CG metallicity scale. The error bars are the quadrature of the relative age's uncertainties form the two corresponding models. Interestingly, the relatives ages are the same, within the uncertainties, independent of the adopted theoretical library.  (2005) ages are consistent with ours within the error bars and no appreciable age trend is observed. On the other hand, some trend is seen when comparing our results to those by Rosenberg et al. (1999): Rosenberg et al. (1999) find marginally older ages for cluster with −1.4 <[M/H]< −0.9 and younger ages for high metallicity clusters. This discrepancy is in part due to the different methods applied to measure the observed cluster parameters (see discussion in De Angeli et al. 2005, and in particular their Figure 7), and in part to the different evolutionary models adopted for the age determination (De Angeli et al. 2005, uses the P04 models, which we have shown to give consistent results with our adopted D07 models) .

The special case of NGC 0288 and NGC 0362
In this study MSTO's magnitude is measured for each GGCs in our database, and then obtained magnitudes are transformed into ages using a set of theoretical stellar evolution models. During this transformation, canonical chemical abundances are assumed for all GGCs. We are aware of the fact that if a particular cluster has different chemical abundances, for example, an anomalous CNO content, then its age determination could be affected by this effect. It is worth mentioning that age uncertainties could be underestimated in this study because of the uncertainties coming from the chemical inhomogeneities which, to date, are impossible to quantify.
In particular, our results indicate that NGC 0288 and NGC 0362 have the same age within ± 0.9 Gyr, while previous works (Sarajedini & Demarque 1990;Bellazzini at al. 2001) have found NGC 0288 to be 1-2 Gyr older than NGC 0362. The relative ages from these studies come from a comparison of all the evolutionary sequences (MS, SGB, HB, etc.), and not on a measurement of the absolute magnitude of the MSTO, as in the present paper. Therefore, we are planning a forthcoming paper in which we will, among other things, consider the case of these two clusters in more detail. In this context, it is worth mentioning that the present study does not take into consideration the distribution of stars along the HB, or finer features such as the RGB tip luminosity, the RGB bump, AGB bump and multiple populations.
A similar argument could be used for the case of NGC 5272 and NGC 6205, for which we obtain identical age, while previous works (VandenBerg, Bolte, & Stetson 1990;Chaboyer, Demarque, & Sarajedini 1996) have found NGC 6205 to be 1-2 Gyr older than NGC 5272. . The mean and dispersion of the low−metallicity clusters' age are marked with solid and dashed lines, respectively. For each of the three metallicity groups, mean age and rms is indicated. Figure 11 shows the same data as Figure 10, but for the ZW metallicity scale.

Discussion
Overall, we notice an increase in age dispersion (which is model independent) with the metallicity. However, we can look at the results shown in the upper panels of Figures 10  and 11 in another way. They show two branches in the age−metallicity relation of GGCs. On one hand, a 'young' branch showing a clear age−metallicity relation, age decreasing at higher metallicities. On the other hand, a largely coeval 'old' branch. The 'old' branch shows an apparent relative age dispersion of ∼0.05 and no age−metallicity relation. It is worth mentioning that while the age−metallicity relation may be model dependent (though it is somehow reassuring that the four most recent theoretical libraries provide consistent results on this respect) the age dispersion−metallicity relation is not model dependent. From now on, we will divide the clusters in two groups that we call 'young' and 'old'. GGCs with total metallicity higher than -1.4 and normalized age younger than 0.95, that is one sigma smaller than the mean of the low−metallicity clusters' age, are considered members of the 'young' group. The remaining GGCs are considered members of the 'old' group. The only exception is NGC 6441, a very high metallicity cluster which is much closer to the old branch than the young one, so it is included in the old group. Column 10 in Table 4 lists whether the target cluster belongs to the old or young group.
An age trend with galactocentric distance is seen in the lower panels of Figures 10  and 11. The fraction of young group's clusters increases significantly as the galactocentric distance increases. As a result, age dispersion increases at increasing galactocentric distance. Other authors have commented on the increased dispersion in ages at larger galactocentric distances. (e. g. Richer et al. 1996, which is one of the first to state the empirical problem clearly). Also in these panels, a clear metallicity trend with galactocentric distance can be observed for the old group of clusters, with an increasing metallicity at lower galactocentric distances. This trend is explicitly shown in Figure 12, which plots GGC's [M/H] in the CG metallicity scale as a function of the galactocentric distance. Open and filled circles represent clusters in the young and old groups, respectively. In the following, the age dispersion of the old clusters group, the age−metallicity relation of the young clusters group and the age trend with galactocentric distance are discussed in detail.

The age dispersion of the old clusters group
Figures 10 and 11 indicate that the mean normalized age of the clusters in the old group is the same for all metallicities. The relative age dispersion of this group is ∼0.05. If an absolute age of 12.8 Gyr (from the D07 models) is assumed as the mean age of these clusters, age dispersion would be ∼0.6 Gyr. We note that if we assume that all sources of uncertainty have been taken into account in the error bars, then at least part of this dispersion must be real. The average normalized age uncertainty among the clusters in the old group is 0.04, which is approximately 0.5 Gyr. Subtracting this quadratically from the observed dispersion, we find an intrinsic dispersion of 0.03 in relative age for the clusters in the old group. Using 12.8 Gyr for the age of this group, the intrinsic dispersion is ∼0.4 Gyr.
These results are consistent with a scenario in which the old group of clusters formed within a fast assembling process of the halo, lasting ∼0.8 Gyr and which should also account for an increasing metallicity towards the center of the Galaxy. This implies that the chemical enrichment of the protogalactic cloud was faster than the assembling timescale. Regarding the nature of such assembling process, it is interesting to note the following. Wilkinson & Evans (1999) estimate the mass and scalelength of the Milky Way dark matter halo to be ∼ 1.9×10 12 M ⊙ and 170 kpc, respectively. The free-fall time of a homogeneous sphere of those mass and radius is ∼ 0.84 Gyr, similar to the age range we find here for the old clusters. In other words, the age dispersion of old globular clusters is not in contradiction with the formation from the colapse of a single protosystem, resembling the model proposed by Eggen, Lynden-Bell, & Sandage (1962). On another side, the standard ΛCDM scenario would predict a general star formation, including star clusters, in protogalaxy building blocks of mass ∼ 10 8 M ⊙ before reionization (Moore et al. 2006). These blocks merge together afterwards to form larger galaxies, and would, in principle, be compatible with the old, coeval cluster population that we observe. More difficult seems to account for the metallicity gradient (Fig. 12) under this scenario. A more detailed analysis is necessary, but it is important to note that any successful galaxy formation scenario must account for the existence of a large nubmer of old, coeval globular clusters in the Milky Way.

The Age−metallicity relation of the young clusters group
Clusters in the young group, however, show a clear age−metallicity relationship, with younger clusters being more metal rich than older ones. The nature of this "young branch" can perhaps be better understood by looking at the GGCs associated with known or putative accreted dwarf galaxies: Sagittarius, Canis Major and Monoceros. Ibata, Gilmore, & Irwin (1995) discovered the accreted Sagittarius dwarf galaxy and its GC system (Terzan 7, Terzan 8, Arp 2, and M54). Dinescu et al. (2000) and Bellazzini et al. (2003) have since then argued, respectively, that Pal 12 and NGC 4147 are part of the extended tidal stream.
The Monoceros Ring was discovered by Newberg et al. (2002), and it has been pro-posed to be the tidal stream of the accreted Monoceros galaxy. Crane et al. (2003) and Frinchaboy et al. (2004) identified five possible GC candidates of this accreted galaxy: NGC 2298, NGC 2808, NGC 5286, Pal 1, and BH 176. In addition, Martin et al. (2004) identify four GGCs -NGC 1851, NGC 1904, NGC 2298 and NGC 2808 -possibly associated with a stellar overdensity in Canis Major, which they suspect to be a dwarf galaxy. However, the nature of these two structures is controversial. The existence of the Canis Major dwarf galaxy has been seriously questioned by Momany et al. (2004Momany et al. ( , 2006. In Momany et al. (2006), the Canis Major stellar over−density has been completely accounted for as the effect of the Milky Way disk warp (an ubiquitous property of all massive galaxies). In addition, Momany et al. (2006) have shown that also the Monoceros Ring seems to reflect the signature of another Galactic disk property; the flaring (the increase in scale-height as a function of Galactocentric distance) of the outer disk. Figures 13A and 13B show the age−metallicity relation, using the CG metallicity scale 6 . Sagittarius, Monoceros and Canis Major's clusters have been marked. It is interesting that, for metallicities higher than [M/H] = -1.4, all of Sagittarius, Monoceros and Canis Major's accreted clusters fall in the young group, with the exception of NGC 5286, and they follow the same age−metallicity relation as the rest of the young group's clusters. This result suggests a different origin for the old and young groups of clusters. It is worth mentioning that there is a number of clusters, other than the Sagittarius, Monoceros and Canis Major's clusters among the young sample.
Another interesting question is related to the recent discovery that GGCs are not simple single stellar populations. The high mass and multiple stellar populations of ωCen has led to speculation that it is the remnant nucleus of an accreted Milky Way satellite galaxy (Lee et al. 1999;Villanova et al. 2007). But ωCen is not a unique case. NGC 2808 (Piotto et al. 2007), NGC 1851 (Milone et al. 2008a), NGC 6388 (Piotto 2008) and NGC 6441 (Caloi & D'Antona 2007) are four massive clusters hosting multiple stellar populations. Freeman (1993) suggested that GGCs form as the nuclei of dwarf galaxies in the early universe, and are accreted as their host galaxies merge onto larger structures. In this context, these peculiar GGCs would be good candidates to be the remnants of accreted satellites. Multiple population clusters have been also marked in Figure 13A.
We note that multiple stellar population clusters have broadened MSTOs and at least a fraction of the stars in these clusters have rather anomalous chemical compositions. For this reason, we should be hesitant in adopting their relative ages obtained in this study. A much more detailed analysis is needed for these peculiar clusters. Villanova et al. (2007) and Cassisi et al. (2008) present particularly illuminating discussions of the difficulties in estimating ages for the multiple populations in ωCen and NGC 1851, respectively. Figure 13B shows the same Figure 13A, but here GGCs with r GC > 10 kpc have are represented with filled circles, while those with r GC < 10 kpc are plotted with open circles. Interestingly, all intermediate− and high−metallicity clusters with r GC > 10 kpc fall on the young group of clusters.
In order to characterize the young group's age−metallicity relation, a least squares fit has been performed, and results are also shown in Figure 13B. The age−metallicity relation observed in the young group of clusters can be described by the equation: This relation is plotted with a solid line in Figure 13B. The dispersion of the clusters' relative ages with respect to the previous equation is also 0.05, and it is also represented in Figure 13B (dashed lines). Following a similar argument as for the old group, if an absolute age of 12.8 Gyr is assumed as the mean of the clusters in the old group, the young group's age dispersion with respect to the previous age−metallicity relation would be ∼0.6 Gyr. The average normalized age uncertainty among the clusters in the young group is 0.04, which is approximately 0.5 Gyr. Subtracting this quadratically from the observed dispersion, we find an intrinsic dispersion with respect to its age−metallicity relation of 0.03 in relative age, or ∼0.4 Gyr.
In summary, our results are consistent with the clusters of the young group having been formed in a different process spanning an interval of time as long as 0.45 in relative age, or ∼6 Gyr, and resulting in a group of GCs with a clear age−metallicity relation. It is very tempting to argue that the origin of this second group of clusters is related to their formation within Milky Way satellite galaxies that were later accreted. However, the reason why the age−metallicity relation of the GGCs associated with the Sagittarius dwarf galaxy, as well as Monoceros and Canis Major, is similar to that of the rest of 'accreted' clusters should be investigated. Perhaps all of the young group clusters share the same origin, which would require an association between Sagittarius, Monoceros and Canis Major, or alternatively, all dwarf galaxy systems may share a common (or very similar) star formation and metal enrichment history, and hence a common age−metallicity relation, which seems unlikely.

Age trend with galactocentric distance
With the aim of analyzing in more detail the observed age trend with galactocentric distance, GGCs in the old and young groups are considered separately. Figure 13C shows normalized ages versus galactocentric distance (same as Figure 10, lower panel), but now clusters in the old group are represented with open circles and clusters in the young group with filled circles. Besides, Sagittarius, Monoceros and Canis Major's clusters are shown as encircled points.
It can be seen that if the old group is considered, an age trend with galactocentric distance is not observed. The age's variance remains constant with the galactocentric distance for this group.
On the other hand if the young group is taken into account, it can be seen that the age's variance increases with galactocentric distance. In order to determine if this variance increase depends either on the galactocentric distance or the metallicity, a three dimensional principal−component (PC) analysis has been done considering the clusters in the young group only. The obtained eigenvalues are 2.03, 0.83 and 0.14. The first PC, of eigenvalue 2.03, accounts for 68% of the total variance. According to the commonly used average or eigenvalue-one criterion, only this PC should be retained. Figure 14 shows the relations between the original data and PCs projected on the two first PCs plane. On the one hand it can be seen that the first PC is strongly correlated with [M/H] and age, but weakly with r GC , while the second PC carries on most of the r GC variance. On another hand it is clear that young cluster age dispersion can be explained by an age−metallicity relation and that age is not related to r GC .
In fact, the effect of metallicity on age of clusters in the young group can be eliminated by subtracting the fitted age−metallicity relation (equation 4) to the cluster's ages. By doing this, age−metallicity 'corrected' ages are obtained. Figure 13D shows the same as Figure 13C, but here age−metallicity 'corrected' clusters are represented with solid circles. It can be clearly seen that most of the age dispersion observed in Figure 13C disappears, and that age is not related to r GC In summary, it can be concluded that a significant part of the age dispersion present in Figure 13C for r GC ≤ 15 kpc and most of that for r GC > 15 can be explained by a strong age−metallicity relation for the clusters of the young group.
It is worth mentioning that our database is biased in distance. No outer−halo clusters, with galactocentric distances larger than ∼20 kpc, have been considered for this study. So the analysis of the galactic halo's age structure is limited by the galactocentric distance range of the database. Given that outer−halo clusters (Pal 3, Pal 4, Pal 14, Eridanus, AM-1) have younger ages (Sarajedini 1997;Stetson et al. 1999;Dotter et al. 2008), it is possible that our age range might be a lower limit on the actual range.

Conclusions
Normalized ages have been derived for a sample of 64 GGCs using the stellar evolution models of D07 and both ZW and CG metallicity scales. We have also performed an analysis using the stellar models of P04, B94 and G00. The result is the most extensive and precise database of normalized ages so far produced . Our results are: • We find that we are able to measure relative ages to a formal precision of 2-7%. The relative ages are independent of the choice of theoretical model. Four independent sets of stellar evolution libraries (D07, P04, B94, G00) produce essentially identical results.
• We find that the GGCs fall into two well-defined groups. The first one represents a population of old clusters that have the same age, a dispersion of 5% in relative age and no age−metallicity relation. If we assume an absolute age of 12.8 Gyr, the absolute age dispersion of the old clusters is ∼0.6 Gyr. Accounting for the measurement uncertainties, we obtain an intrinsic age dispersion of ∼0.4 Gyr. The second group of clusters shows a clear age−metallicity relation, with young clusters being more metal rich than older ones. Also in this case, the intrinsic age dispersion with respect to its age−metallicity relation is ∼0.4 Gyr. Besides, there is no age dispersion−metallicity relation if we look at the two samples (old and young) separately.
• These results are consistent with a scenario in which the formation of the Milky Way CCG system took place in two phases. The first one produced clusters with ages within a range of ∼0.8 Gyr. This age range is compatible with the timescale for the collapse of a protogalaxy of the same mass and scalelength as that of the Milky Way dark matter halo. In other words, the age dispersion of the old group of globular clusters is not in contradiction with the formation from the colapse of a single protosystem, resembling the model proposed by Eggen, Lynden-Bell, & Sandage (1962). The standard ΛCDM scenario would also predict a significant star formation in protogalactic building blocks before reionization, which, after merging, would produce an old, coeval population of globular clusters in a large galaxy like the Milky Way. However, to account for the galactocentric metallicity gradient of this group seems difficult in the context of this mechanism. The second phase spanned a time interval as long as ∼6 Gyr and resulted in a group of GGCs with a clear age−metallicity relation. It is very tempting to argue that the origin of this second group of clusters is related to their formation within Milky Way satellite galaxies that were later accreted. However, this would not account for the fact that the clusters of this group share the same age−metallicity relation.
• We find that the increasing dispersion in age with galactocentric distance is explained by the age−metallicity relation of GGCs.
• These results are independent of the assumed metallicity scale.
• It is worth mentioning that our database is biased by distance. No outer halo clusters, with galactocentric distances larger than ∼20 kpc, have been considered. The analysis of the galactic halo's age structure is therefore limited by the galactocentric distance range of the database.    measurements over synthetic CMDs with 50,000 stars brighter than M F606W = 6, no differential reddening, and varying the fraction of binary stars. Grey points show the MSTO measurements for the 10 different CMDs generated in each case, and black points show the average and standard deviation of each of the 10 measurement sets. Solid line represents the MSTO's magnitude measured directly on the isochrone. The Central panel shows the same as left panel, but now the fraction of binary stars is 0.0 and the maximum differential reddening is varying between 0.0 and 0.2. The left panel shows the same, but now the fraction of binary stars is 0.0, the differential reddening is 0.0 and the number of stars brighter then M F606W = 6 is varying between 1,000 to 50,000. Dashed lines represent the adopted ±σ in each case. See text for details.   for different ages as a function of global metallicity, in 1 Gyr (solid lines) and 0.5 Gyr (dashed lines) age steps. The lower isochrone corresponds to an age of 15 Gyr and the upper one to 6 Gyr. In the case of G00 models, metallicities lower than Z=0.0004 have been extrapolated. Points represent observational data. Fig. 8.-Comparison of normalized ages derived from different sets of models in the CG metallicity scale. The difference between normalized ages derived using P04, B94 and G00 models, respectively, and the D07 models are shown. It can be seen that relative ages derived using the main sequence fitting method are independent from the adopted model.  , and as a function of galactocentric distance (r GC kpc, lower pannel). These results have been derived using the D07 stellar evolution library. Open circles, filled triangles and filled circles represent GCs within the low−, intermediate− and high−metallicity groups, respectively. For each of the three metallicity groups, mean age and rms is indicated. See text for details.   and normalized ages versus galactocentric distance (C and D). In A, clusters associated with Sagittarius, Monoceros and Canis Major have been marked, together with multiple stellar population clusters. In B, GGCs with r GC < 10 kpc are represented with open circles, while those with r GC > 10 kpc are plotted with filled circles. A least squares fit to the young group's age−metallicity relation is also plotted. In C, Sagittarius, Monoceros and Canis Major' clusters are also marked. Finally, D shows the same as C, but here ages of clusters in the young group have been age−metallicity 'corrected' using the age−metallicity relation shown in B. See text for details.