Sensitivity study with a D and B mesons modular simulation code of heavy flavor RAA and azimuthal anisotropies based on beam energy, initial conditions, hadronization, and suppression mechanisms

Heavy flavor probes provide important information about the in-medium properties of the quark-gluon plasma produced in heavy-ion collisions. In this work, we investigate the effects of 2D+1 event-by-event fluctuating hydrodynamic backgrounds on the nuclear suppression factor and momentum anisotropies of heavy flavor mesons and nonphotonic electrons. Using the state-of-the-art D and B mesons modular simulation code (called "DAB-MOD"), we perform a systematic comparison of different transport equations in the same background, including a few energy-loss models - with and without energy-loss fluctuations - and a relativistic Langevin model with different drag parametrizations. We present the resulting D and B mesons RAA, v2, v3, and v4 as well as multiparticle cumulants, in AuAu collisions at sNN=200GeV and PbPb collisions at sNN=2.76TeV and sNN=5.02TeV, and compare them to the available experimental data. The v2{4}/v2{2} ratio, which is known to be a powerful probe of the initial conditions and flow fluctuations in the soft sector, is also studied in the context of heavy flavor. We also investigate the correlations between the transverse anisotropies of heavy mesons and all charged particles to better understand how heavy quarks couple to the hydrodynamically expanding quark-gluon plasma. We study the influence that different initial conditions and the implementation of heavy-light quark coalescence has on our results.