Researcher profile

M. J. Betancourt

M. J. Betancourt contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
33works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

33 published item(s)

preprint2015arXiv

Optimizing The Integrator Step Size for Hamiltonian Monte Carlo

Hamiltonian Monte Carlo can provide powerful inference in complex statistical problems, but ultimately its performance is sensitive to various tuning parameters. In this paper we use the underlying geometry of Hamiltonian Monte Carlo to construct a universal optimization criteria for tuning the step size of the symplectic integrator crucial to any implementation of the algorithm as well as diagnostics to monitor for any signs of invalidity. An immediate outcome of this result is that the suggested target average acceptance probability of 0.651 can be relaxed to $0.6 \lesssim a \lesssim 0.9$ with larger values more robust in practice.

preprint2015arXiv

The Fundamental Incompatibility of Hamiltonian Monte Carlo and Data Subsampling

Leveraging the coherent exploration of Hamiltonian flow, Hamiltonian Monte Carlo produces computationally efficient Monte Carlo estimators, even with respect to complex and high-dimensional target distributions. When confronted with data-intensive applications, however, the algorithm may be too expensive to implement, leaving us to consider the utility of approximations such as data subsampling. In this paper I demonstrate how data subsampling fundamentally compromises the efficient exploration of Hamiltonian flow and hence the scalable performance of Hamiltonian Monte Carlo itself.

preprint2014arXiv

Measurement of Charge Multiplicity Asymmetry Correlations in High Energy Nucleus-Nucleus Collisions at 200 GeV

A study is reported of the same- and opposite-sign charge-dependent azimuthal correlations with respect to the event plane in Au+Au collisions at 200 GeV. The charge multiplicity asymmetries between the up/down and left/right hemispheres relative to the event plane are utilized. The contributions from statistical fluctuations and detector effects were subtracted from the (co-)variance of the observed charge multiplicity asymmetries. In the mid- to most-central collisions, the same- (opposite-) sign pairs are preferentially emitted in back-to-back (aligned on the same-side) directions. The charge separation across the event plane, measured by the difference, $Δ$, between the like- and unlike-sign up/down $-$ left/right correlations, is largest near the event plane. The difference is found to be proportional to the event-by-event final-state particle ellipticity (via the observed second-order harmonic $v^{\rm obs}_{2}$), where $Δ=(1.3\pm1.4({\rm stat})^{+4.0}_{-1.0}({\rm syst}))\times10^{-5}+(3.2\pm0.2({\rm stat})^{+0.4}_{-0.3}({\rm syst}))\times10^{-3}v^{\rm obs}_{2}$ for 20-40% Au+Au collisions. The implications for the proposed chiral magnetic effect are discussed.

preprint2014arXiv

The Geometric Foundations of Hamiltonian Monte Carlo

Although Hamiltonian Monte Carlo has proven an empirical success, the lack of a rigorous theoretical understanding of the algorithm has in many ways impeded both principled developments of the method and use of the algorithm in practice. In this paper we develop the formal foundations of the algorithm through the construction of measures on smooth manifolds, and demonstrate how the theory naturally identifies efficient implementations and motivates promising generalizations.

preprint2013arXiv

A General Metric for Riemannian Manifold Hamiltonian Monte Carlo

Markov Chain Monte Carlo (MCMC) is an invaluable means of inference with complicated models, and Hamiltonian Monte Carlo, in particular Riemannian Manifold Hamiltonian Monte Carlo (RMHMC), has demonstrated impressive success in many challenging problems. Current RMHMC implementations, however, rely on a Riemannian metric that limits their application to analytically-convenient models. In this paper I propose a new metric for RMHMC without these limitations and verify its success on a distribution that emulates many hierarchical and latent models.

preprint2013arXiv

Cruising The Simplex: Hamiltonian Monte Carlo and the Dirichlet Distribution

Due to its constrained support, the Dirichlet distribution is uniquely suited to many applications. The constraints that make it powerful, however, can also hinder practical implementations, particularly those utilizing Markov Chain Monte Carlo (MCMC) techniques such as Hamiltonian Monte Carlo. I demonstrate a series of transformations that reshape the canonical Dirichlet distribution into a form much more amenable to MCMC algorithms.

preprint2013arXiv

Fluctuations of charge separation perpendicular to the event plane and local parity violation in sqrt(sNN)=200 GeV Au+Au collisions at RHIC

Previous experimental results based on data (15 million events) collected by the STAR detector at RHIC suggest event-by-event charge separation fluctuations perpendicular to the event plane in non-central heavy-ion collisions. Here we present the correlator previously used split into its two component parts to reveal correlations parallel and perpendicular to the event plane. The results are from a high statistics 200 GeV Au+Au collisions data set (57 million events) collected by the STAR experiment. We explicitly count units of charge separation from which we find clear evidence for more charge separation fluctuations perpendicular than parallel to the event plane. We also employ a modified correlator to study the possible P-even background in same and opposite charge correlations, and find that the P-even background may largely be explained by momentum conservation and collective motion.

preprint2013arXiv

Freeze-out Dynamics via Charged Kaon Femtoscopy in sqrt(sNN)=200 GeV Central Au+Au Collisions

We present measurements of three-dimensional correlation functions of like-sign low transverse momentum kaon pairs from sqrt(sNN)=200 GeV Au+Au collisions. A Cartesian surface-spherical harmonic decomposition technique was used to extract the kaon source function. The latter was found to have a three-dimensional Gaussian shape and can be adequately reproduced by Therminator event generator simulations with resonance contributions taken into account. Compared to the pion one, the kaon source function is generally narrower and does not have the long tail along the pair transverse momentum direction. The kaon Gaussian radii display a monotonic decrease with increasing transverse mass m_T over the interval of 0.55<=m_T<=1.15 GeV/c^2. While the kaon radii are adequately described by the m_T-scaling in the outward and sideward directions, in the longitudinal direction the lowest m_T value exceeds the expectations from a pure hydrodynamical model prediction.

preprint2013arXiv

Generalizing the No-U-Turn Sampler to Riemannian Manifolds

Hamiltonian Monte Carlo provides efficient Markov transitions at the expense of introducing two free parameters: a step size and total integration time. Because the step size controls discretization error it can be readily tuned to achieve certain accuracy criteria, but the total integration time is left unconstrained. Recently Hoffman and Gelman proposed a criterion for tuning the integration time in certain systems with their No U-Turn Sampler, or NUTS. In this paper I investigate the dynamical basis for the success of NUTS and generalize it to Riemannian Manifold Hamiltonian Monte Carlo.

preprint2013arXiv

Hamiltonian Monte Carlo for Hierarchical Models

Hierarchical modeling provides a framework for modeling the complex interactions typical of problems in applied statistics. By capturing these relationships, however, hierarchical models also introduce distinctive pathologies that quickly limit the efficiency of most common methods of in- ference. In this paper we explore the use of Hamiltonian Monte Carlo for hierarchical models and demonstrate how the algorithm can overcome those pathologies in practical applications.

preprint2013arXiv

Single Spin Asymmetry $A_N$ in Polarized Proton-Proton Elastic Scattering at $\sqrt{s}=200$ GeV

We report a high precision measurement of the transverse single spin asymmetry $A_N$ at the center of mass energy $\sqrt{s}=200$ GeV in elastic proton-proton scattering by the STAR experiment at RHIC. The $A_N$ was measured in the four-momentum transfer squared $t$ range $0.003 \leqslant |t| \leqslant 0.035$ $\GeVcSq$, the region of a significant interference between the electromagnetic and hadronic scattering amplitudes. The measured values of $A_N$ and its $t$-dependence are consistent with a vanishing hadronic spin-flip amplitude, thus providing strong constraints on the ratio of the single spin-flip to the non-flip amplitudes. Since the hadronic amplitude is dominated by the Pomeron amplitude at this $\sqrt{s}$, we conclude that this measurement addresses the question about the presence of a hadronic spin flip due to the Pomeron exchange in polarized proton-proton elastic scattering.

preprint2013arXiv

System size dependence of transverse momentum correlations at $\sqrt{s_{NN}}=$ 62.4 and 200 GeV at the BNL Relativistic Heavy Ion Collider

We present a study of the average transverse momentum ($p_t$) fluctuations and $p_t$ correlations for charged particles produced in Cu+Cu collisions at midrapidity for $\sqrt{s_{NN}} =$ 62.4 and 200 GeV. These results are compared with those published for Au+Au collisions at the same energies, to explore the system size dependence. In addition to the collision energy and system size dependence, the $p_t$ correlation results have been studied as functions of the collision centralities, the ranges in $p_t$, the pseudorapidity $η$, and the azimuthal angle $ϕ$. The square root of the measured $p_t$ correlations when scaled by mean $p_t$ is found to be independent of both colliding beam energy and system size studied. Transport-based model calculations are found to have a better quantitative agreement with the measurements compared to models which incorporate only jetlike correlations.

preprint2013arXiv

Third Harmonic Flow of Charged Particles in Au+Au Collisions at sqrtsNN = 200 GeV

We report measurements of the third harmonic coefficient of the azimuthal anisotropy, v_3, known as triangular flow. The analysis is for charged particles in Au+Au collisions at sqrtsNN = 200 GeV, based on data from the STAR experiment at the Relativistic Heavy Ion Collider. Two-particle correlations as a function of their pseudorapidity separation are fit with narrow and wide Gaussians. Measurements of triangular flow are extracted from the wide Gaussian, from two-particle cumulants with a pseudorapidity gap, and also from event plane analysis methods with a large pseudorapidity gap between the particles and the event plane. These results are reported as a function of transverse momentum and centrality. A large dependence on the pseudorapidity gap is found. Results are compared with other experiments and model calculations.

preprint2012arXiv

Anomalous centrality evolution of two-particle angular correlations from Au-Au collisions at $\sqrt{s_{\rm NN}}$ = 62 and 200 GeV

We present two-dimensional (2D) two-particle angular correlations on relative pseudorapidity $η$ and azimuth $ϕ$ for charged particles from Au-Au collisions at $\sqrt{s_{\rm NN}} = 62$ and 200 GeV with transverse momentum $p_t \geq 0.15$ GeV/$c$, $|η| \leq 1$ and $2π$ azimuth. Observed correlations include a {same-side} (relative azimuth $< π/2$) 2D peak, a closely-related away-side azimuth dipole, and an azimuth quadrupole conventionally associated with elliptic flow. The same-side 2D peak and away-side dipole are explained by semihard parton scattering and fragmentation (minijets) in proton-proton and peripheral nucleus-nucleus collisions. Those structures follow N-N binary-collision scaling in Au-Au collisions until mid-centrality where a transition to a qualitatively different centrality trend occurs within a small centrality interval. Above the transition point the number of same-side and away-side correlated pairs increases rapidly {relative to} binary-collision scaling, the $η$ width of the same-side 2D peak also increases rapidly ($η$ elongation) and the $ϕ$ width actually decreases significantly. Those centrality trends are more remarkable when contrasted with expectations of jet quenching in a dense medium. Observed centrality trends are compared to {\sc hijing} predictions and to the expected trends for semihard parton scattering and fragmentation in a thermalized opaque medium. We are unable to reconcile a semihard parton scattering and fragmentation origin for the observed correlation structure and centrality trends with heavy ion collision scenarios which invoke rapid parton thermalization. On the other hand, if the collision system is effectively opaque to few-GeV partons the observations reported here would be inconsistent with a minijet picture.

preprint2012arXiv

Di-electron spectrum at mid-rapidity in $p+p$ collisions at $\sqrt{s} = 200$ GeV

We report on mid-rapidity mass spectrum of di-electrons and cross sections of pseudoscalar and vector mesons via $e^{+}e^{-}$ decays, from $\sqrt{s} = 200$ GeV $p+p$ collisions, measured by the large acceptance experiment STAR at RHIC. The ratio of the di-electron continuum to the combinatorial background is larger than 10% over the entire mass range. Simulations of di-electrons from light-meson decays and heavy-flavor decays (charmonium and open charm correlation) are found to describe the data. The extracted $ω\rightarrow e^{+}e^{-}$ invariant yields are consistent with previous measurements. The mid-rapidity yields ($dN/dy$) of $ϕ$ and $J/ψ$ are extracted through their di-electron decay channels and are consistent with the previous measurements of $ϕ\rightarrow K^{+}K^{-}$ and $J/ψ\rightarrow e^{+}e^{-}$. Our results suggest a new upper limit of the branching ratio of the $η\rightarrow e^{+}e^{-}$ of $1.7\times10^{-5}$ at 90% confidence level.

preprint2012arXiv

Longitudinal and transverse spin asymmetries for inclusive jet production at mid-rapidity in polarized p+p collisions at sqrt{s}=200 GeV

We report STAR measurements of the longitudinal double-spin asymmetry A_LL, the transverse single-spin asymmetry A_N, and the transverse double-spin asymmetries A_Sigma and A_TT for inclusive jet production at mid-rapidity in polarized p+p collisions at a center-of-mass energy of sqrt{s} = 200 GeV. The data represent integrated luminosities of 7.6 /pb with longitudinal polarization and 1.8 /pb with transverse polarization, with 50-55% beam polarization, and were recorded in 2005 and 2006. No evidence is found for the existence of statistically significant jet A_N, A_Sigma, or A_TT at mid-rapidity. Recent model calculations indicate the A_N results may provide new limits on the gluon Sivers distribution in the proton. The asymmetry A_LL significantly improves the knowledge of gluon polarization in the nucleon.

preprint2012arXiv

Measurements of $D^{0}$ and $D^{*}$ Production in $p$ + $p$ Collisions at $\sqrt{s}$ = 200 GeV

We report measurements of charmed-hadron ($D^{0}$, $D^{*}$) production cross sections at mid-rapidity in $p$ + $p$ collisions at a center-of-mass energy of 200 GeV by the STAR experiment. Charmed hadrons were reconstructed via the hadronic decays $D^{0}\rightarrow K^{-}π^{+}$, $D^{*+}\rightarrow D^{0}π^{+}\rightarrow K^{-}π^{+}π^{+}$ and their charge conjugates, covering the $p_T$ range of 0.6$-$2.0 GeV/$c$ and 2.0$-$6.0 GeV/$c$ for $D^{0}$ and $D^{*+}$, respectively. From this analysis, the charm-pair production cross section at mid-rapidity is $dσ/dy|_{y=0}^{c\bar{c}}$ = 170 $\pm$ 45 (stat.) $^{+38}_{-59}$ (sys.) $μ$b. The extracted charm-pair cross section is compared to perturbative QCD calculations. The transverse momentum differential cross section is found to be consistent with the upper bound of a Fixed-Order Next-to-Leading Logarithm calculation.

preprint2012arXiv

Strangeness Enhancement in Cu+Cu and Au+Au Collisions at \sqrt{s_{NN}} = 200 GeV

We report new STAR measurements of mid-rapidity yields for the $Λ$, $\barΛ$, $K^{0}_{S}$, $Ξ^{-}$, $\barΞ^{+}$, $Ω^{-}$, $\barΩ^{+}$ particles in Cu+Cu collisions at \sNN{200}, and mid-rapidity yields for the $Λ$, $\barΛ$, $K^{0}_{S}$ particles in Au+Au at \sNN{200}. We show that at a given number of participating nucleons, the production of strange hadrons is higher in Cu+Cu collisions than in Au+Au collisions at the same center-of-mass energy. We find that aspects of the enhancement factors for all particles can be described by a parameterization based on the fraction of participants that undergo multiple collisions.

preprint2012arXiv

Transverse Single-Spin Asymmetry and Cross-Section for pi0 and eta Mesons at Large Feynman-x in Polarized p+p Collisions at sqrt(s)=200 GeV

Measurements of the differential cross-section and the transverse single-spin asymmetry, A_N, vs. x_F for pi0 and eta mesons are reported for 0.4 < x_F < 0.75 at an average pseudorapidity of 3.68. A data sample of approximately 6.3 pb^{-1} was analyzed, which was recorded during p+p collisions at sqrt{s} = 200 GeV by the STAR experiment at RHIC. The average transverse beam polarization was 56%. The cross-section for pi0 is consistent with a perturbative QCD prediction, and the eta/pi0 cross-section ratio agrees with previous mid-rapidity measurements. For 0.55 < x_F < 0.75, A_N for eta (0.210 +- 0.056) is 2.2 standard deviations larger than A_N for pi0 (0.081 +- 0.016).

preprint2011arXiv

Directed and elliptic flow of charged particles in Cu+Cu collisions at $\sqrt{\bm {s_{NN}}} =$ 22.4 GeV

This paper reports results for directed flow $v_{1}$ and elliptic flow $v_{2}$ of charged particles in Cu+Cu collisions at $\sqrt{s_{NN}}=$ 22.4 GeV at the Relativistic Heavy Ion Collider. The measurements are for the 0-60% most central collisions, using charged particles observed in the STAR detector. Our measurements extend to 22.4 GeV Cu+Cu collisions the prior observation that $v_1$ is independent of the system size at 62.4 and 200 GeV, and also extend the scaling of $v_1$ with $η/y_{\rm beam}$ to this system. The measured $v_2(p_T)$ in Cu+Cu collisions is similar for $\sqrt{s_{NN}} = 22.4-200$ GeV. We also report a comparison with results from transport model (UrQMD and AMPT) calculations. The model results do not agree quantitatively with the measured $v_1(η), v_2(p_T)$ and $v_2(η)$.

preprint2011arXiv

Energy and system-size dependence of two- and four-particle $v_2$ measurements in heavy-ion collisions at RHIC and their implications on flow fluctuations and nonflow

We present STAR measurements of azimuthal anisotropy by means of the two- and four-particle cumulants $v_2$ ($v_2\{2\}$ and $v_2\{4\}$) for Au+Au and Cu+Cu collisions at center of mass energies $\sqrt{s_{_{\mathrm{NN}}}} = 62.4$ and 200 GeV. The difference between $v_2\{2\}^2$ and $v_2\{4\}^2$ is related to $v_{2}$ fluctuations ($σ_{v_2}$) and nonflow $(δ_{2})$. We present an upper limit to $σ_{v_2}/v_{2}$. Following the assumption that eccentricity fluctuations $σ_ε$ dominate $v_2$ fluctuations $\frac{σ_{v_2}}{v_2} \approx \frac{σ_ε}ε$ we deduce the nonflow implied for several models of eccentricity fluctuations that would be required for consistency with $v_2\{2\}$ and $v_2\{4\}$. We also present results on the ratio of $v_2$ to eccentricity.

preprint2011arXiv

Measurement of the $W \to e ν$ and $Z/γ^* \to e^+e^-$ Production Cross Sections at Mid-rapidity in Proton-Proton Collisions at $\sqrt{s}$ = 500 GeV

We report measurements of the charge-separated $W^{+(-)} \to e^{+(-)} + ν_e(\barν_e)$ and $Z/γ^* \to e^+e^-$ production cross sections at mid-rapidity in proton-proton collisions at $\sqrt{s}$ = 500 GeV. These results are based on 13.2 pb$^{-1}$ of data recorded in 2009 by the STAR detector at RHIC. Production cross sections for W bosons that decay via the $e ν$ channel were measured to be $σ(pp \to W^+ X) \cdot BR(W^+ \to e^+ ν_e)$ = 117.3 \pm 5.9(stat) \pm 6.2(syst) \pm 15.2(lumi) pb, and $σ(pp \to W^- X) \cdot BR(W^- \to e^- \barν_e)$ = 43.3 \pm 4.6(stat) \pm 3.4(syst) \pm 5.6(lumi) pb. For $Z/γ^*$ production, $σ(pp \to Z/γ^* X) \cdot BR(Z/γ^* \to e^+ e^-)$ = 7.7 \pm 2.1(stat) $^{+0.5}_{-0.9}$(syst) \pm 1.0(lumi) pb for di-lepton invariant masses $m_{e^+e^-}$ between 70 and 110 GeV/$c^2$. First measurements of the W cross section ratio, $σ(pp \to W^+ X) / σ(pp \to W^- X)$, at $\sqrt{s}$ = 500 GeV are also reported. Theoretical predictions, calculated using recent parton distribution functions, are found to agree with the measured cross sections.

preprint2011arXiv

Pion femtoscopy in p+p collisions at sqrt(s)=200 GeV

The STAR Collaboration at RHIC has measured two-pion correlation functions from p+p collisions at sqrt(s)=200 GeV. Spatial scales are extracted via a femtoscopic analysis of the correlations, though this analysis is complicated by the presence of strong non-femtoscopic effects. Our results are put into the context of the world dataset of femtoscopy in hadron-hadron collisions. We present the first direct comparison of femtoscopy in p+p and heavy ion collisions, under identical analysis and detector conditions.

preprint2011arXiv

Studies of di-jet survival and surface emission bias in Au+Au collisions via angular correlations with respect to back-to-back leading hadrons

We report first results from an analysis based on a new multi-hadron correlation technique, exploring jet-medium interactions and di-jet surface emission bias at RHIC. Pairs of back-to-back high transverse momentum hadrons are used for triggers to study associated hadron distributions. In contrast with two- and three-particle correlations with a single trigger with similar kinematic selections, the associated hadron distribution of both trigger sides reveals no modification in either relative pseudo-rapidity or relative azimuthal angle from d+Au to central Au+Au collisions. We determine associated hadron yields and spectra as well as production rates for such correlated back-to-back triggers to gain additional insights on medium properties.

preprint2011arXiv

System size and energy dependence of near-side di-hadron correlations

Two-particle azimuthal ($Δϕ$) and pseudorapidity ($Δη$) correlations using a trigger particle with large transverse momentum ($p_T$) in $d$+Au, Cu+Cu and Au+Au collisions at $\sqrt{s_{NN}}$ =\xspace 62.4 GeV and 200~GeV from the STAR experiment at RHIC are presented. The \ns correlation is separated into a jet-like component, narrow in both $Δϕ$ and $Δη$, and the ridge, narrow in $Δϕ$ but broad in $Δη$. Both components are studied as a function of collision centrality, and the jet-like correlation is studied as a function of the trigger and associated $p_T$. The behavior of the jet-like component is remarkably consistent for different collision systems, suggesting it is produced by fragmentation. The width of the jet-like correlation is found to increase with the system size. The ridge, previously observed in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV, is also found in Cu+Cu collisions and in collisions at $\sqrt{s_{NN}}$ =\xspace 62.4 GeV, but is found to be substantially smaller at $\sqrt{s_{NN}}$ =\xspace 62.4 GeV than at $\sqrt{s_{NN}}$ = 200 GeV for the same average number of participants ($ \langle N_{\mathrm{part}}\rangle$). Measurements of the ridge are compared to models.

preprint2010arXiv

A Bayesian Approach To Histogram Comparison

Determining if two histograms are consistent, whether they have been drawn from the same underlying distribution or not, is a common problem in physics. Existing approaches are not only limited in power but also inapplicable to histograms filled with importance weights, a common feature of Monte Carlo simulations. From a Bayesian perspective, the comparison between a single underlying distribution and two underlying distributions is readily solved within the context of model comparison. I introduce an implementation of Bayesian model comparison to the problem, including the extension to importance sampling.

preprint2010arXiv

An Experimental Exploration of the QCD Phase Diagram: The Search for the Critical Point and the Onset of De-confinement

The QCD phase diagram lies at the heart of what the RHIC Physics Program is all about. While RHIC has been operating very successfully at or close to its maximum energy for almost a decade, it has become clear that this collider can also be operated at lower energies down to 5 GeV without extensive upgrades. An exploration of the full region of beam energies available at the RHIC facility is imperative. The STAR detector, due to its large uniform acceptance and excellent particle identification capabilities, is uniquely positioned to carry out this program in depth and detail. The first exploratory beam energy scan (BES) run at RHIC took place in 2010 (Run 10), since several STAR upgrades, most importantly a full barrel Time of Flight detector, are now completed which add new capabilities important for the interesting physics at BES energies. In this document we discuss current proposed measurements, with estimations of the accuracy of the measurements given an assumed event count at each beam energy.

preprint2010arXiv

Azimuthal di-hadron correlations in d+Au and Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV from STAR

Yields, correlation shapes, and mean transverse momenta \pt{} of charged particles associated with intermediate to high-\pt{} trigger particles ($2.5 < \pt < 10$ \GeVc) in d+Au and Au+Au collisions at $\snn=200$ GeV are presented. For associated particles at higher $\pt \gtrsim 2.5$ \GeVc, narrow correlation peaks are seen in d+Au and Au+Au, indicating that the main production mechanism is jet fragmentation. At lower associated particle $\pt < 2$ \GeVc, a large enhancement of the near- ($\dphi \sim 0$) and away-side ($\dphi \sim π$) associated yields is found, together with a strong broadening of the away-side azimuthal distributions in Au+Au collisions compared to d+Au measurements, suggesting that other particle production mechanisms play a role. This is further supported by the observed significant softening of the away-side associated particle yield distribution at $\dphi \sim π$ in central Au+Au collisions.

preprint2010arXiv

Higher Moments of Net-proton Multiplicity Distributions at RHIC

We report the first measurements of the kurtosis (κ), skewness (S) and variance (σ^2) of net-proton multiplicity (N_p - N_pbar) distributions at midrapidity for Au+Au collisions at \sqrt(s_NN) = 19.6, 62.4, and 200 GeV corresponding to baryon chemical potentials (μ_B) between 200 - 20 MeV. Our measurements of the products κσ^2 and S σ, which can be related to theoretical calculations sensitive to baryon number susceptibilities and long range correlations, are constant as functions of collision centrality. We compare these products with results from lattice QCD and various models without a critical point and study the \sqrt(s_NN) dependence of κσ^2. From the measurements at the three beam energies, we find no evidence for a critical point in the QCD phase diagram for μ_B below 200 MeV.

preprint2010arXiv

K*0 production in Cu+Cu and Au+Au collisions at \sqrt{s_NN} = 62.4 GeV and 200 GeV

We report on K*0 production at mid-rapidity in Au+Au and Cu+Cu collisions at \sqrt{s_{NN}} = 62.4 and 200 GeV collected by the Solenoid Tracker at RHIC (STAR) detector. The K*0 is reconstructed via the hadronic decays K*0 \to K+ pi- and \bar{K*0} \to K-pi+. Transverse momentum, pT, spectra are measured over a range of pT extending from 0.2 GeV/c to 5 GeV/c. The center of mass energy and system size dependence of the rapidity density, dN/dy, and the average transverse momentum, <pT>, are presented. The measured N(K*0)/N(K) and N(ϕ)/N(K*0) ratios favor the dominance of re-scattering of decay daughters of K*0 over the hadronic regeneration for the K*0 production. In the intermediate pT region (2.0 < pT < 4.0 GeV/c), the elliptic flow parameter, v2, and the nuclear modification factor, RCP, agree with the expectations from the quark coalescence model of particle production.

preprint2010arXiv

Nested Sampling with Constrained Hamiltonian Monte Carlo

Nested sampling is a powerful approach to Bayesian inference ultimately limited by the computationally demanding task of sampling from a heavily constrained probability distribution. An effective algorithm in its own right, Hamiltonian Monte Carlo is readily adapted to efficiently sample from any smooth, constrained distribution. Utilizing this constrained Hamiltonian Monte Carlo, I introduce a general implementation of the nested sampling algorithm.

preprint2010arXiv

Upsilon cross section in p+p collisions at sqrt(s) = 200 GeV

We report on a measurement of the Upsilon(1S+2S+3S) -> e+e- cross section at midrapidity in p+p collisions at sqrt(s)=200 GeV. We find the cross section to be 114 +/- 38 (stat.) +23,-24 (syst.) pb. Perturbative QCD calculations at next-to-leading order in the Color Evaporation Model are in agreement with our measurement, while calculations in the Color Singlet Model underestimate it by 2 sigma. Our result is consistent with the trend seen in world data as a function of the center-of-mass energy of the collision and extends the availability of Upsilon data to RHIC energies. The dielectron continuum in the invariant mass range near the Upsilon is also studied to obtain a combined cross section of Drell-Yan plus (b b-bar) -> e+e-.

preprint2009arXiv

Pion Interferometry in Au+Au and Cu+Cu Collisions at $\sqrt{s_{\rm{NN}}}$ = 62.4 and 200 GeV

We present a systematic analysis of two-pion interferometry in Au+Au collisions at $\sqrt{s_{\rm{NN}}}$ = 62.4 GeV and Cu+Cu collisions at $\sqrt{s_{\rm{NN}}}$ = 62.4 and 200 GeV using the STAR detector at RHIC. The multiplicity and transverse momentum dependences of the extracted correlation lengths (radii) are studied. The scaling with charged particle multiplicity of the apparent system volume at final interaction is studied for the RHIC energy domain. The multiplicity scaling of the measured correlation radii is found to be independent of colliding system and collision energy.