Simulation in ALICE

ALICE, the experiment dedicated to the study of heavy ion collisions at the LHC, uses an object-oriented framework for simulation, reconstruction and analysis (AliRoot) based on ROOT. Here, we describe the general ALICE simulation strategy and those components of the framework related to simulation. Two main requirements have driven the development of the simulation components. First, the possibility to run different transport codes with the same user code for geometry and detector response has led to the development of the Virtual Monte Carlo concept. Second, simulation has to provide tools to efficiently study events ranging from low-multiplicity pp collisions to Pb-Pb collisions with up to 80000 primary particles per event. This has led to the development of a variety of collaborating generator classes and specific classes for event merging.Comment: Talk from the 2003 Computing in High Energy and Nuclear Physics (CHEP03), La Jolla, Ca, USA, March 2003, 6 pages, LaTeX, 5 eps figures. PSN TUMT00

ALICE experience with GEANT4

Since its release in 1999, the LHC experiments have been evaluating GEANT4 in view of adopting it as a replacement for the obsolescent GEANT3 transport MonteCarlo. The ALICE collaboration has decided to perform a detailed physics validation of elementary hadronic processes against experimental data already used in international benchmarks. In one test, proton interactions on different nuclear targets have been simulated, and the distribution of outgoing particles has been compared to data. In a second test, penetration of quasi-monoenergetic low energy neutrons through a thick shielding has been simulated and again compared to experimental data. In parallel, an effort has been put on the integration of GEANT4 in the AliRoot framework. An overview of the present status of ALICE GEANT4 simulation and the remaining problems will be presented. This document will describe in detail the results of these tests, together with the improvements that the GEANT4 team has made to the program as a result of the feedback received from the ALICE collaboration. We will also describe the remaining problems that have been communicated to GEANT4 but not yet addressed.Comment: 8 pages, 12 figures, for the CHEP03 conference proceeding

Rayleigh scattering and atomic dynamics in dissipative optical lattices

We investigate Rayleigh scattering in dissipative optical lattices. In particular, following recent proposals [S. Guibal et al., Phys. Rev. Lett. 78, 4709 (1997); C. Jurczak et al., Phys. Rev. Lett. 77, 1727 (1996)], we study whether the Rayleigh resonance originates from the diffraction on a density grating and is therefore a probe of transport of atoms in optical lattices. It turns out that this is not the case: the Rayleigh line is instead a measure of the cooling rate, while spatial diffusion contributes to the scattering spectrum with a much broader resonance

Normal fault earthquakes or graviquakes

Earthquakes are dissipation of energy throughout elastic waves. Canonically is the elastic energy accumulated during the interseismic period. However, in crustal extensional settings, gravity is the main energy source for hangingwall fault collapsing. Gravitational potential is about 100 times larger than the observed magnitude, far more than enough to explain the earthquake. Therefore, normal faults have a different mechanism of energy accumulation and dissipation (graviquakes) with respect to other tectonic settings (strike-slip and contractional), where elastic energy allows motion even against gravity. The bigger the involved volume, the larger is their magnitude. The steeper the normal fault, the larger is the vertical displacement and the larger is the seismic energy released. Normal faults activate preferentially at about 60° but they can be shallower in low friction rocks. In low static friction rocks, the fault may partly creep dissipating gravitational energy without releasing great amount of seismic energy. The maximum volume involved by graviquakes is smaller than the other tectonic settings, being the activated fault at most about three times the hypocentre depth, explaining their higher b-value and the lower magnitude of the largest recorded events. Having different phenomenology, graviquakes show peculiar precursor

Do regions matter in ALICE?: Social relationships and data exchanges in the Grid

This study aims at investigating the impact of regional affiliations of centres on the organisation of collaborations within the Distributed Computing ALICE infrastructure, based on social networks methods. A self-administered questionnaire was sent to all centre managers about support, email interactions and wished collaborations in the infrastructure. Several additional measures, stemming from technical observations were collected, such as bandwidth, data transfers and Internet Round Trip Time (RTT) were also included. Information for 50 centres were considered (about 70% response rate). Empirical analysis shows that despite the centralisation on CERN, the network is highly organised by regions. The results are discussed in the light of policy and efficiency issue

Shear-free perfect fluids with a solenoidal electric curvature

We prove that the vorticity or the expansion vanishes for any shear-free perfect fluid solution of the Einstein field equations where the pressure satisfies a barotropic equation of state and the spatial divergence of the electric part of the Weyl tensor is zero.Comment: 9 page

The ALICE Off-Line Strategy: A Successful Migration to OO

The ALICE Experiment has chosen to start developing its software directly in OO, using the services of the ROOT system, which is ALICE's official candidate for the common LHC framework. This had lead to the definition of a complete environment (AliRoot) where the software developed by the different experimental groups is being integrated. Different test-benches for I/O and Simulation have been set up based on real production code. This allows early assessment of technology, both software and hardware in a realistic production environment. Different codes, such as GEANT3, GEANT4 and FLUKA, or the reconstruction algorithms by the physicists developing the detectors, have been easily integrated in the framework, that has shown to be both evolutive and modular.The ALICE Collaboration has adopted this setup and we are now successfully migrating the users into it. This talk describes the AliRoot environment and its future evolution

Cold atom realizations of Brownian motors

Brownian motors are devices which "rectify" Brownian motion, i.e. they can generate a current of particles out of unbiased fluctuations. Brownian motors are important for the understanding of molecular motors, and are also promising for the realization of new nanolelectronic devices. Among the different systems that can be used to study Brownian motors, cold atoms in optical lattices are quite an unusual one: there is no thermal bath and both the potential and the fluctuations are determined by laser fields. In this article recent experimental implementations of Brownian motors using cold atoms in optical lattices are reviewed

Stochastic resonance in periodic potentials: realization in a dissipative optical lattice

We have observed the phenomenon of stochastic resonance on the Brillouin propagation modes of a dissipative optical lattice. Such a mode has been excited by applying a moving potential modulation with phase velocity equal to the velocity of the mode. Its amplitude has been characterized by the center-of-mass (CM) velocity of the atomic cloud. At Brillouin resonance, we studied the CM-velocity as a function of the optical pumping rate at a given depth of the potential wells. We have observed a resonant dependence of the CM velocity on the optical pumping rate, corresponding to the noise strength. This corresponds to the experimental observation of stochastic resonance in a periodic potential in the low-damping regime

Shearfree perfect fluids with solenoidal magnetic curvature and a gamma-law equation of state

We show that shearfree perfect fluids obeying an equation of state p=(gamma -1) mu are non-rotating or non-expanding under the assumption that the spatial divergence of the magnetic part of the Weyl tensor is zero.Comment: 11 page