The Offline Software Framework of the Pierre Auger Observatory

The Pierre Auger Observatory is designed to unveil the nature and the origins of the highest energy cosmic rays. The large and geographically dispersed collaboration of physicists and the wide-ranging collection of simulation and reconstruction tasks pose some special challenges for the offline analysis software. We have designed and implemented a general purpose framework which allows collaborators to contribute algorithms and sequencing instructions to build up the variety of applications they require. The framework includes machinery to manage these user codes, to organize the abundance of user-contributed configuration files, to facilitate multi-format file handling, and to provide access to event and time-dependent detector information which can reside in various data sources. A number of utilities are also provided, including a novel geometry package which allows manipulation of abstract geometrical objects independent of coordinate system choice. The framework is implemented in C++, and takes advantage of object oriented design and common open source tools, while keeping the user side simple enough for C++ novices to learn in a reasonable time. The distribution system incorporates unit and acceptance testing in order to support rapid development of both the core framework and contributed user code.Comment: 4 pages, 2 figures, presented at IEEE NSS/MIC, 23-29 October 2005, Puerto Ric

The Multifunctional Sorting Protein PACS-2 Regulates SIRT1-Mediated Deacetylation of p53 to Modulate p21-Dependent Cell-Cycle Arrest

SIRT1 regulates the DNA damage response by deacetylating p53, thereby repressing p53 transcriptional output. Here, we demonstrate that the sortingprotein PACS-2 regulates SIRT1-mediated deacetylation of p53 to modulate the DNA damage response. PACS-2 knockdown cells failed to efficiently undergo p53-induced cell-cycle arrest in response to DNA damage. Accordingly, p53 acetylation was reduced both in PACS-2 knockdown cells and thymocytes from Pacs-2-/- mice, thereby blunting induction of the cyclin-dependent kinase inhibitor p21 (CDKN1A). The SIRT1 inhibitor EX-527 or SIRT1 knockdown restored p53 acetylation and p21 induction as well as p21-dependent cell-cycle arrest in PACS-2 knockdown cells. Traffickingstudies revealed that cytoplasmic PACS-2 shuttled to the nucleus, where it interacted with SIRT1 andrepressed SIRT1-mediated p53 deacetylation. Correspondingly, invitro assays demonstrated that PACS-2 directly inhibited SIRT1-catalyzed p53 deacetylation. Together, these findings identify PACS-2 as an invivo mediator of the SIRT1-p53-p21 axis that modulates the DNA damage response

Resolving galaxies in time and space: II: Uncertainties in the spectral synthesis of datacubes

In a companion paper we have presented many products derived from the application of the spectral synthesis code STARLIGHT to datacubes from the CALIFA survey, including 2D maps of stellar population properties and 1D averages in the temporal and spatial dimensions. Here we evaluate the uncertainties in these products. Uncertainties due to noise and spectral shape calibration errors and to the synthesis method are investigated by means of a suite of simulations based on 1638 CALIFA spectra for NGC 2916, with perturbations amplitudes gauged in terms of the expected errors. A separate study was conducted to assess uncertainties related to the choice of evolutionary synthesis models. We compare results obtained with the Bruzual & Charlot models, a preliminary update of them, and a combination of spectra derived from the Granada and MILES models. About 100k CALIFA spectra are used in this comparison. Noise and shape-related errors at the level expected for CALIFA propagate to 0.10-0.15 dex uncertainties in stellar masses, mean ages and metallicities. Uncertainties in A_V increase from 0.06 mag in the case of random noise to 0.16 mag for shape errors. Higher order products such as SFHs are more uncertain, but still relatively stable. Due to the large number statistics of datacubes, spatial averaging reduces uncertainties while preserving information on the history and structure of stellar populations. Radial profiles of global properties, as well as SFHs averaged over different regions are much more stable than for individual spaxels. Uncertainties related to the choice of base models are larger than those associated with data and method. Differences in mean age, mass and metallicity are ~ 0.15 to 0.25 dex, and 0.1 mag in A_V. Spectral residuals are ~ 1% on average, but with systematic features of up to 4%. The origin of these features is discussed. (Abridged)Comment: A&A, accepte

Observational constraints to boxy/peanut bulge formation time

Boxy/peanut bulges are considered to be part of the same stellar structure as bars and both could be linked through the buckling instability. The Milky Way is our closest example. The goal of this letter is determining if the mass assembly of the different components leaves an imprint in their stellar populations allowing to estimate the time of bar formation and its evolution. To this aim we use integral field spectroscopy to derive the stellar age distributions, SADs, along the bar and disc of NGC 6032. The analysis shows clearly different SADs for the different bar areas. There is an underlying old (>=12 Gyr) stellar population for the whole galaxy. The bulge shows star formation happening at all times. The inner bar structure shows stars of ages older than 6 Gyrs with a deficit of younger populations. The outer bar region presents a SAD similar to that of the disc. To interpret our results, we use a generic numerical simulation of a barred galaxy. Thus, we constrain, for the first time, the epoch of bar formation, the buckling instability period and the posterior growth from disc material. We establish that the bar of NGC 6032 is old, formed around 10 Gyr ago while the buckling phase possibly happened around 8 Gyr ago. All these results point towards bars being long-lasting even in the presence of gas.Comment: Accepted for publication in MNRAS Letter

The Offline Software Framework of the Pierre Auger Observatory

To be published in the ProceedingsInternational audienceThe Pierre Auger Observatory is designed to unveil the nature and the origins of the highest energy cosmic rays. The large and geographically dispersed collaboration of physicists and the wide-ranging collection of simulation and reconstruction tasks pose some special challenges for the offline analysis software. We have designed and implemented a general purpose framework which allows collaborators to contribute algorithms and sequencing instructions to build up the variety of applications they require. The framework includes machinery to manage these user codes, to organize the abundance of user-contributed configuration files, to facilitate multi-format file handling, and to provide access to event and time-dependent detector information which can reside in various data sources. A number of utilities are also provided, including a novel geometry package which allows manipulation of abstract geometrical objects independent of coordinate system choice. The framework is implemented in C++, and takes advantage of object oriented design and common open source tools, while keeping the user side simple enough for C++ novices to learn in a reasonable time. The distribution system incorporates unit and acceptance testing in order to support rapid development of both the core framework and contributed user code

Ionized gas kinematics of galaxies in the CALIFA survey I: Velocity fields, kinematic parameters of the dominant component, and presence of kinematically distinct gaseous systems

This work provides an overall characterization of the kinematic behavior of the ionized gas of the galaxies included in the Calar Alto Legacy Integral field Area (CALIFA), offering kinematic clues to potential users of this survey for including kinematical criteria for specific studies. From the first 200 galaxies observed by CALIFA, we present the 2D kinematic view of the 177 galaxies satisfying a gas detection threshold. After removing the stellar contribution, we used the cross-correlation technique to obtain the radial velocity of the dominant gaseous component. The main kinematic parameters were directly derived from the radial velocities with no assumptions on the internal motions. Evidence of the presence of several gaseous components with different kinematics were detected by using [OIII] profiles. Most objects in the sample show regular velocity fields, although the ionized-gas kinematics are rarely consistent with simple coplanar circular motions. 35% of the objects present evidence of a displacement between the photometric and kinematic centers larger than the original spaxel radii. Only 17% of the objects in the sample exhibit kinematic lopsidedness when comparing receding and approaching sides of the velocity fields, but most of them are interacting galaxies exhibiting nuclear activity. Early-type galaxies in the sample present clear photometric-kinematic misaligments. There is evidence of asymmetries in the emission line profiles suggesting the presence of kinematically distinct gaseous components at different distances from the nucleus. This work constitutes the first determination of the ionized gas kinematics of the galaxies observed in the CALIFA survey. The derived velocity fields, the reported kinematic peculiarities and the identification of the presence of several gaseous components might be used as additional criteria for selecting galaxies for specific studies.Comment: 38 pages, 16 figures, 4 tables. Paper accepted for publication in A&

Discovering a Light Higgs Boson with Light

We evaluate the prospects for detecting a non-standard light Higgs boson with a significant branching ratio to two photons, in Run II of the Fermilab Tevatron. We derive the reach for several channels: $2\gamma$ inclusive, $2\gamma+1$ jet and $2\gamma+2$ jets. We present the expected Run II limits on the branching ratio of $h\to\gamma\gamma$ as a function of the Higgs mass, for the case of ``bosonic'', as well as ``topcolor'' Higgs bosons.Comment: 11 pages, LaTeX, 7 figures, 4 tables, uses aipproc2.sty, contributed to the Physics at Run II Workshop, analysis redone with optimized cuts and improved background estimate, references adde

The Pierre Auger Observatory offline software

The Pierre Auger Observatory aims to discover the nature and origins of the highest energy cosmic rays. The large number of physicists involved in the project and the diversity of simulation and reconstruction tasks pose a challenge for the offline analysis software, not unlike the challenges confronting software for very large high energy physics experiments. Previously we have reported on the design and implementation of a general purpose but relatively lightweight framework which allows collaborators to contribute algorithms and sequencing instructions to build up the variety of applications they require. In this report, we update the status of this work and describe some of the successes and difficulties encountered over the last few years of use. We explain the machinery used to manage user contributions, to organize the abundance of configuration files, to facilitate multi-format file handling, and to provide access to event and time-dependent detector information residing in various data sources. We also describe the testing procedures used to help maintain stability of the code in the face of a large number of contributions. Foundation classes will also be discussed, including a novel geometry package which allows manipulation of abstract geometrical objects independent of coordinate system choice

Resolving galaxies in time and space: I: Applying STARLIGHT to CALIFA data cubes

Fossil record methods based on spectral synthesis techniques have matured over the past decade, and their application to integrated galaxy spectra fostered substantial advances on the understanding of galaxies and their evolution. Yet, because of the lack of spatial resolution, these studies are limited to a global view, providing no information about the internal physics of galaxies. Motivated by the CALIFA survey, which is gathering Integral Field Spectroscopy over the full optical extent of 600 galaxies, we have developed an end-to-end pipeline which: (i) partitions the observed data cube into Voronoi zones in order to, when necessary and taking due account of correlated errors, increase the S/N, (ii) extracts spectra, including propagated errors and bad-pixel flags, (iii) feeds the spectra into the STARLIGHT spectral synthesis code, (iv) packs the results for all galaxy zones into a single file, (v) performs a series of post-processing operations, including zone-to-pixel image reconstruction and unpacking the spectral and stellar population properties into multi-dimensional time, metallicity, and spatial coordinates. This paper provides an illustrated description of this whole pipeline and its products. Using data for the nearby spiral NGC 2916 as a show case, we go through each of the steps involved, presenting ways of visualizing and analyzing this manifold. These include 2D maps of properties such as the v-field, stellar extinction, mean ages and metallicities, mass surface densities, star formation rates on different time scales and normalized in different ways, 1D averages in the temporal and spatial dimensions, projections of the stellar light and mass growth (x,y,t) cubes onto radius-age diagrams, etc. The results illustrate the richness of the combination of IFS data with spectral synthesis, providing a glimpse of what is to come from CALIFA and future surveys. (Abridged)Comment: A&A, accepte

Stellar Population gradients in galaxy discs from the CALIFA survey

While studies of gas-phase metallicity gradients in disc galaxies are common, very little has been done in the acquisition of stellar abundance gradients in the same regions. We present here a comparative study of the stellar metallicity and age distributions in a sample of 62 nearly face-on, spiral galaxies with and without bars, using data from the CALIFA survey. We measure the slopes of the gradients and study their relation with other properties of the galaxies. We find that the mean stellar age and metallicity gradients in the disc are shallow and negative. Furthermore, when normalized to the effective radius of the disc, the slope of the stellar population gradients does not correlate with the mass or with the morphological type of the galaxies. Contrary to this, the values of both age and metallicity at $\sim$2.5 scale-lengths correlate with the central velocity dispersion in a similar manner to the central values of the bulges, although bulges show, on average, older ages and higher metallicities than the discs. One of the goals of the present paper is to test the theoretical prediction that non-linear coupling between the bar and the spiral arms is an efficient mechanism for producing radial migrations across significant distances within discs. The process of radial migration should flatten the stellar metallicity gradient with time and, therefore, we would expect flatter stellar metallicity gradients in barred galaxies. However, we do not find any difference in the metallicity or age gradients in galaxies with without bars. We discuss possible scenarios that can lead to this absence of difference.Comment: 24 pages, 17 figures, accepted for publication in A&