Azimuthal jet flavor tomography with CUJET2.0 of nuclear collisions at RHIC and LHC

A perturbative QCD based jet tomographic Monte Carlo model, CUJET2.0, is presented to predict jet quenching observables in relativistic heavy ion collisions at RHIC/BNL and LHC/CERN energies. This model generalizes the DGLV theory of flavor dependent radiative energy loss by including multi-scale running strong coupling effects. It generalizes CUJET1.0 by computing jet path integrations though more realistic 2+1D transverse and longitudinally expanding viscous hydrodynamical fields contrained by fits to low $p_T$ flow data. The CUJET2.0 output depends on three control parameters, $(\alpha_{max},f_E,f_M)$, corresponding to an assumed upper bound on the vacuum running coupling in the infrared and two chromo-electric and magnetic QGP screening mass scales $(f_E \mu(T), f_M \mu(T))$ where $\mu(T)$ is the 1-loop Debye mass. We compare numerical results as a function of $\alpha_{max}$ for pure and deformed HTL dynamically enhanced scattering cases corresponding to $(f_E=1,2, f_M=0)$ to data of the nuclear modification factor, $R^f_{AA}(p_T,\phi; \sqrt{s}, b)$ for jet fragment flavors $f=\pi,D, B, e$ at $\sqrt{s}=0.2-2.76$ ATeV c.m. energies per nucleon pair and with impact parameter $b=2.4, 7.5$ fm. A $\chi^2$ analysis is presented and shows that $R^\pi_{AA}$ data from RHIC and LHC are consistent with CUJET2.0 at the $\chi^2/d.o.f< 2$ level for $\alpha_{max}=0.23-0.30$. The corresponding $\hat{q}(E_{jet}, T)/T^3$ effective jet transport coefficient field of this model is computed to facilitate comparison to other jet tomographic models in the literature. The predicted elliptic asymmetry, $v_2(p_T;\sqrt{s},b)$ is, however, found to significantly underestimated relative to RHIC and LHC data. We find the $\chi^2_{v_2}$ analysis shows that $v_2$ is very sensitive to allowing even as little as 10\% variations of the path averaged $\alpha_{max}$ along in and out of reaction plane paths.