Astroparticle experiments such as IceCube or MAGIC require a deconvolution of their measured data with respect to the response function of the detector to provide the distributions of interest, e.g. energy spectra. In this paper, appropriate uncertainty limits that also allow to draw conclusions on the geometric shape of the underlying distribution are determined using bootstrap methods, which are frequently applied in statistical applications. Bootstrap is a collective term for resampling methods that can be employed to approximate unknown probability distributions or features thereof. A clear advantage of bootstrap methods is their wide range of applicability. For instance, they yield reliable results, even if the usual normality assumption is violated. The use, meaning and construction of uncertainty limits to any user-specific confidence level in the form of confidence intervals and levels are discussed. The precise algorithms for the implementation of these methods, applicable for any deconvolution algorithm, are given. The proposed methods are applied to Monte Carlo simulations to show their feasibility and their precision in comparison to the statistical uncertainties calculated with the deconvolution software TRUEE.

N-body simulations are widely used to simulate the dynamical evolution of a variety of systems, among them star clusters. Much of our understanding of their evolution rests on the results of such direct N-body simulations. They provide insight in the structural evolution of star clusters, as well as into the occurrence of stellar exotica. Although the major pure N-body codes STARLAB/KIRA and NBODY4 are widely used for a range of applications, there is no thorough comparison study yet. Here we thoroughly compare basic quantities as derived from simulations performed either with STARLAB/KIRA or NBODY4. We construct a large number of star cluster models for various stellar mass function settings (but without stellar/binary evolution, primordial binaries, external tidal fields etc), evolve them in parallel with STARLAB/KIRA and NBODY4, analyse them in a consistent way and compare the averaged results quantitatively. For this quantitative comparison we develop a bootstrap algorithm for functional dependencies. We find an overall excellent agreement between the codes, both for the clusters' structural and energy parameters as well as for the properties of the dynamically created binaries. However, we identify small differences, like in the energy conservation before core collapse and the energies of escaping stars, which deserve further studies. Our results reassure the comparability and the possibility to combine results from these two major N-body codes, at least for the purely dynamical models (i.e. without stellar/binary evolution) we performed. (abridged)

A complete map of the 3D distribution of molecular (CO) gas was constructed using a realistic dynamical model of the gas flow in the barred potential of the Milky Way. The map shows two prominent spiral arms starting at the bar ends connecting smoothly to the 4-armed spiral pattern observed in the atomic hydrogen gas in the outer Galaxy. Unlike previous attempts, our new map uncovers the gas distribution in the bar region of the Galaxy and the far side of the disk. For the first time, we can follow spiral arms in gas as they pass behind the galactic centre.

We present a new model of the three-dimensional distribution of molecular gas in the Milky Way Galaxy, based on CO line data. Our analysis is based on a gas-flow simulation of the inner Galaxy using smoothed-particle hydrodynamics (SPH) using a realistic barred gravitional potential derived from the observed COBE/DIRBE near-IR light distribution. The gas model prescribes the gas orbits much better than a simple circular rotation model and is highly constrained by observations, but it cannot predict local details. In this study, we provide a 3D map of the observed molecular gas distribution using the velocity field from the SPH model. A comparison with studies of the Galactic Center region suggests that the main structures are reproduced but somewhat stretched along the line-of-sight, probably on account of limited resolution of the underlying SPH simulation. The gas model will be publicly available and may prove useful in a number of applications, among them the analysis of diffuse gamma-ray emission as measured with GLAST.

The luminosity functions (LFs) of star cluster systems (i.e. the number of clusters per luminosity interval) are vital diagnostics to probe the conditions of star cluster formation. Early studies have revealed a clear dichotomy between old globular clusters and young clusters, with the former characterised by Gaussian-shaped LFs, and the latter following a power law. Recently, this view was challenged by studies of galaxy merger remnants and post-starburst galaxies. In this paper we re-evaluate the young ($\lta$ few hundreds of Myrs, with the majority $\lta$ few tens of Myrs) star cluster system in the ongoing spiral-spiral major merger system NGC 4038/39, the "Antennae" galaxies. The Antennae galaxies represent a very active and complex star-forming environment, which hampers cluster selection and photometry as well as the determination of observational completeness fractions. A main issue of concern is the large number of bright young stars contained in most earlier studies, which we carefully exclude from our cluster sample by accurately determining the source sizes. The resulting LFs are fitted both with Gaussian and with power-law distributions, taking into account both the observational completeness fractions and photometric errors, and compared using a likelihood ratio test. The likelihood ratio results are rigidly evaluated using Monte Carlo simulations. We perform a number of additional tests, e.g. with subsets of the total sample, all confirming our main result: that a Gaussian distribution fits the observed LFs of clusters in this preferentially very young cluster system significantly better than a power-law distribution, at a (statistical) error probability of less than 0.5 per cent.

We build a stellar-dynamical model of the Milky Way barred bulge and disk, using a newly implemented adaptive particle method. The underlying mass model has been previously shown to match the Galactic near-infrared surface brightness as well as gas-kinematic observations. Here we show that the new stellar-dynamical model also matches the observed stellar kinematics in several bulge fields, and that its distribution of microlensing event timescales reproduces the observed timescale distribution of the {\it MACHO} experiment with a reasonable stellar mass function. The model is therefore an excellent basis for further studies of the Milky Way. We also predict the observational consequences of this mass function for parallax shifted events.

We discuss the systematic uncertainties inherent to analyses of observed (broad-band) Spectral Energy Distributions (SEDs) of star clusters with evolutionary synthesis models. We investigate the effects caused by restricting oneself to a limited number of available passbands, choices of various passband combinations, finite observational errors, non-continuous model input parameter values, and restrictions in parameter space allowed during analysis. Starting from a complete set of UBVRIJH passbands (respectively their Hubble Space Telescope/WFPC2 equivalents) we investigate to which extent clusters with different combinations of age, metallicity, internal extinction and mass can or cannot be disentangled in the various evolutionary stages throughout their lifetimes and what are the most useful passbands required to resolve the ambiguities. We find the U and B bands to be of the highest significance, while the V band and near-infrared data provide additional constraints. A code is presented that makes use of luminosities of a star cluster system in all of the possibly available passbands, and tries to find ranges of allowed age-metallicity-extinction-mass combinations for individual members of star cluster systems. Numerous tests and examples are presented. We show the importance of good photometric accuracies and of determining the cluster parameters independently without any prior assumptions.

We analyse multi-wavelength Hubble Space Telescope (HST) observations of a large number of star clusters in the nearby (post-) starburst dwarf galaxy NGC 1569. Their spectral energy distributions (SEDs) cover at least the wavelength range from U to I in equivalent HST filters, in most cases supplemented by near-infrared data. Using our most up-to-date evolutionary synthesis models of the Goettingen GALEV code we determine ages, metallicities, extinction values and masses for each individual cluster robustly and independently. We confirm the youth of most of these objects. The majority were formed in a very intense starburst starting around 25 Myr ago. While there are two prominent ``super star clusters'' present in this galaxy, with masses of (5-16) * 10^5 Msun, almost all remaining clusters are significantly less massive than an average Milky Way globular cluster, and are generally consistent with open cluster-type objects. We determine the cluster mass function from individual cluster masses derived by scaling the model SEDs of known mass to the observed cluster SEDs for each individual cluster. We find signs of a change in the cluster mass function as the burst proceeds, which we attribute to the special conditions of star cluster formation in this starburst dwarf galaxy environment.

We present new gas flow models for the Milky Way inside the solar circle. Using a mass model derived from the COBE/DIRBE maps and clump giant star counts, and using a parametric model for the spiral arm pattern in the disk, we calculate the gas flow and compare with 12CO observations. We find that models with 4 spiral arms fit the observations better then 2-armed models. We also find that models with separate pattern speeds for the bar and spiral arms can explain the gas flow in the bar corotation region better than single-pattern speed models.

In astrophysical (inverse) regression problems it is an important task to decide whether a given parametric model describes the observational data sufficiently well or whether a non-parametric modelling becomes necessary. However, in contrast to common practice this cannot be decided by solely comparing the quality of fit due to possible over-fitting by the non-parametric method. Therefore, in this paper we present a resampling algorithm which allows to decide whether deviations between a parametric and a non-parametric model are systematic or due to noise. The algorithm is based on a statistical comparison of the corresponding residuals, under the assumption of the parametric model as well as under violation of this assumption. This yields a graphical tool for a robust decision making of parametric versus non-parametric modelling. Moreover, our approach can be used for the selection of the most proper model among several competitors (model selection). The methods are illustrated by the problem of recovering the luminosity density in the Milky Way [MW] from near-infrared [NIR] surface brightness data of the DIRBE experiment on board of the COBE satellite. Among the parametric models investigated one with 4-armed spiral structure performs best. In this model the Sagittarius-Carina arm and its counter-arm are significantly weaker than the other pair of arms. Furthermore, we find statistical evidence for an improvement over a range of parametric models with different spiral structure morphologies by a non-parametric model of Bissantz & Gerhard (2002).

We present new gas flow models for the Milky Way inside the solar circle. To this end we use SPH simulations in gravitational potentials determined from the NIR luminosity distribution (including spiral arms) which are based on the COBE/DIRBE maps. Gas flows in models which include massive spiral arms clearly match the observed 12CO lvplot better than if the potential does not include spiral structure. Besides single pattern speed models we investigate models with separate pattern speeds for the bar and spiral arms. The most important difference is that in the latter case the gas spiral arms go through the bar corotation region, keeping the gas aligned with the arms there. In the (l,v) plot this results in characteristic regions which appear to be nearly void of gas.In single pattern speed models these regions are filled with gas because the spiral arms dissolve in the bar corotation region. Comparing with the 12CO data we find evidence for separate pattern speeds in the Milky Way.From a series of models the preferred range for the bar pattern speed is Om_p=60\pm5 /Gyr, corresponding to corotation at 3.4\pm0.3kpc. The spiral pattern speed is less well constrained, but our preferred value is Om_sp\approx 20 /Gyr. A further series of gas models is computed for different bar angles, using separately determined luminosity models and gravitational potentials in each case. We find acceptable gas models for 20<=\phibar<=25. The model with (\phibar=20, Om_p=60 /Gyr, Om_sp=20 /Gyr) gives an excellent fit to the spiral arm ridges in the observed (l,v) plot.

A common problem in physics is to fit regression data by a parametric class of functions, and to decide whether a certain functional form allows for a good fit of the data. Common goodness of fit methods are based on the calculation of the distribution of certain statistical quantities under the assumption that the model under consideration {\it holds true}. This proceeding bears methodological flaws, e.g. a good ``fit'' - albeit the model is wrong - might be due to over-fitting, or to the fact that the chosen statistical criterion is not powerful enough against the present particular deviation between model and true regression function. This causes particular difficulties when models with different numbers of parameters are to be compared. Therefore the number of parameters is often penalised additionally. We provide a methodology which circumvents these problems to some extent. It is based on the consideration of the error distribution of the goodness of fit criterion under a broad range of possible models - and not only under the assumption that a given model holds true. We present a graphical method to decide for the most evident model from a range of parametric models of the data. The method allows to quantify statistical evidence {\it for} the model (up to some distance between model and true regression function) and not only {\it absence of evidence} against, as common goodness of fit methods do. Finally we apply our method to the problem of recovering the luminosity density of the Milky Way from a de-reddened {\it COBE/DIRBE} L-band map. We present statistical evidence for flaring of the stellar disc inside the solar circle.

A new model for the luminosity distribution in the inner Milky Way is found, using a non-parametric penalized maximum-likelihood algorithm to deproject a dereddened COBE/DIRBE L-band map of the inner Galaxy. The model is also constrained by the apparent magnitude (line-of-sight) distributions of clump giant stars in certain bulge fields. An important new feature is the inclusion of a spiral arm model in the disk. Spiral arms make the model appear broader on the sky, thus our bar is more elongated than in previous eight-fold symmetric models. They also lead to a smoother disk model interior to the Sun. The bar length is approx. 3.5kpc and its axis ratios are 1:(0.3-0.4):0.3, independent of whether the spiral arm model is 4-armed or 2-armed. The larger elongation in the plane makes it possible to reproduce the observed clump giant distributions as well. For the bar angle we estimate 15^o <= phi <= 30^o, with the best models obtained for 20^o <= phi <= 25^o. We use our reference model to predict a microlensing optical depth map towards the bulge, normalising its mass by the observed terminal velocity curve. For clump giant sources at (l,b)=(3.9,-3.8) we find tau_{-6}:=tau/10^{-6}=1.27, within 1.8*sigma of the new MACHO measurement given by Popowski et al. The value for all sources at (l,b)=(2.68,-3.35) is tau_{-6}=1.1, still >3*sigma away from the published MACHO DIA value. The dispersion of these tau_{-6} values within our models is 10%. Because the distribution of sources is well-fit by the NIR model, increasing the predicted optical depths by >20% will be difficult. Thus the high value of the measured clump giant optical depth argues for a near-maximal disk in the Milky Way. (abridged)

The assumption that a parametric class of functions fits the data structure sufficiently well is common in fitting curves and surfaces to regression data. One then derives a parameter estimate resulting from a least squares fit, say, and in a second step various kinds of chi^2 goodness of fit measures, to assess whether the deviation between data and estimated surface is due to random noise and not to systematic departures from the model. In this paper we show that commonly-used chi^2-measures are invalid in regression models, particularly when inhomogeneous noise is present. Instead we present a bootstrap algorithm which is applicable in problems described by noisy versions of Fredholm integral equations. of the first kind. We apply the suggested method to the problem of recovering the luminosity density in the Milky Way from data of the DIRBE experiment on board the COBE satellite.

We generalize to elliptical models the argument of Kuijken (1997), which connects the microlensing optical depth towards the Galactic bulge to the Galactic rotation curve. When applied to the latest value from the MACHO collaboration for the optical depth for microlensing of bulge sources, the argument implies that the Galactic bar cannot plausibly reconcile the measured values of the optical depth, the rotation curve and the local mass density. Either there is a problem with the interpretation of the microlensing data, or our line of sight to the Galactic centre is highly atypical in that it passes through a massive structure that wraps only a small distance around the Galactic centre.

We examine the left-right asymmetry in the cleaned COBE/DIRBE near-infrared data of the inner Galaxy and show (i) that the Galactic bar is probably not seen very nearly end-on, and (ii) that even if it is, it is not highly elongated. The assumption of constant mass-to-light ratio is used to derive simulated terminal-velocity plots for the ISM from our model luminosity distributions. By comparing these plots with observed terminal velocities we determine the mass-to-light ratio of the near-IR bulge and disk. Assuming that all this mass contributes to gravitational microlensing we compute optical depths $\tau$ for microlensing in Galactic-centre fields. For three models with bar major axis between $10\deg-25\deg$ from the Sun-Galactic Center line, the resulting optical depths in Baade's window lie in the range $0.83\times10^{-6} \lta \tau \lta 0.89\times10^{-6}$ for main-sequence stars and $1.2\times10^{-6} \lta \tau \lta 1.3\times10^{-6}$ for red-clump giants. We discuss a number of uncertainties including possible variations of the near-infrared mass-to-light ratio. We conclude that, although the values predicted from analyzing the COBE and gas velocity data are inconsistent at the $2-2.5\sigma$ level with recent observational determinations of $\tau$, we believe they should be taken seriously.