ebe_scripts/analysis_multi_particle_correlation_oversample_UrQMDevent.py
def calcualte_vn_2(vn_data_array):
    """
        this function computes vn{2} and its stat. err.
        self correlation is substracted
    """
    vn_data_array = array(vn_data_array)
    nev = len(vn_data_array[:, 0])
    dN = vn_data_array[:, 0]
    dN = dN.reshape(len(dN), 1)
    Qn_array = dN*vn_data_array[:, 1:]
    corr = 1./(dN*(dN - 1.))*(Qn_array*conj(Qn_array) - dN)
    vn_2 = sqrt(real(mean(corr, 0))) + 1e-30
    vn_2_err = std(real(corr), 0)/sqrt(nev)/2./vn_2
    return(vn_2, vn_2_err)

        this function computes vn{2} and its stat. err.
        self correlation is substracted
    
this function computes vn{2} and its stat. err.
self correlation is substracted
this function computes vn { 2 } and its stat . err . self correlation is substracted
def calcualte_vn_2(vn_data_array):
    vn_data_array = array(vn_data_array)
    nev = len(vn_data_array[:, 0])
    dN = vn_data_array[:, 0]
    dN = dN.reshape(len(dN), 1)
    Qn_array = dN*vn_data_array[:, 1:]
    corr = 1./(dN*(dN - 1.))*(Qn_array*conj(Qn_array) - dN)
    vn_2 = sqrt(real(mean(corr, 0))) + 1e-30
    vn_2_err = std(real(corr), 0)/sqrt(nev)/2./vn_2
    return(vn_2, vn_2_err)
def calcualte_vn_2 ( vn_data_array ) : vn_data_array = array ( vn_data_array ) nev = len ( vn_data_array [ : , 0 ] ) dN = vn_data_array [ : , 0 ] dN = dN . reshape ( len ( dN ) , 1 ) Qn_array = dN * vn_data_array [ : , 1 : ] corr = 1. / ( dN * ( dN - 1. ) ) * ( Qn_array * conj ( Qn_array ) - dN ) vn_2 = sqrt ( real ( mean ( corr , 0 ) ) ) + 1e-30 vn_2_err = std ( real ( corr ) , 0 ) / sqrt ( nev ) / 2. / vn_2 return ( vn_2 , vn_2_err )
this function computes vn{2} and its stat.
this function computes vn { 2 } and its stat .
"""
        this function computes vn{2} and its stat. err.
        self correlation is substracted
    """
vn_data_array