This program deals with hadronic afterburner output files and perform various analysis on hadronic observables
Single particle spectra for charged hadrons and identified particles are collected. The anisotropic flow coefficients from two particle correlation are also analyzed. As well as event-plane correlation, vn_distribution, rn ratio for flow factorizatin breaking.
The HBT correlation between identified particle pair is originated from the quantum statistic of the final wave function. The Monte-Carlo samples of the emitted particles do not have this correlation built in. So we need to add it by hand.
The HBT correlation function can be approximatly computed as,
C(q, K) - 1 = \frac{\int d^4 x \int d^4 y s(x, p1) s(y, p2) cos(q \cdot (x - y)}{\int d^4 x \int d^4 y s(x, p1) s(y, p2)},
where K = 0.5(p1 + p2) and q = p1 - p2.
In the program, we introduce oversampling parameter to increase statistics by first grouping multiply events together to a single big event and then form the pairs. The number of the pairs will increase by oversampling factor.
Pairs from mixed event are form by first rotating the mixture event by a random angle which eliminate the unwanted correlations.
The output correlation function is in 3-d q space (q_out, q_side, q_long).
Particle yield distrbution from a given pT and rapidity cut can be collected from UrQMD outputs