# PHOBOS Nuclear Modification Factor in d+Au, tagged p+Au and tagged n+Au
# collisions at 200 GeV/c
#
# 2015-05-19
# 
# 
#                d2Nch(dAu)/deta.dp_t 
# R_dAu = -----------------------------------
#         <Ncoll>(dAu) d2Nch(pbarp)/deta.dp_t
# 
# 
#                d2Nch(pAu)/deta.dp_t 
# R_pAu = -----------------------------------
#         <Ncoll>(pAu) d2Nch(pbarp)/deta.dp_t
# 
# 
#                d2Nch(nAu)/deta.dp_t 
# R_nAu = -----------------------------------
#         <Ncoll>(nAu) d2Nch(pbarp)/deta.dp_t
# 
# 
# R_pnAu = (0.4 * R_pAu) + (0.6 * R_nAu)
# 
# 
# Note that systematic errors scale all points together.
# 
# 
# 
# ----------
# Centrality bin: 70-100% of d+Au cross section
# <Ncoll>(d+Au)          = 2.00250e+00
# <Ncoll>(d+Au) stat err = 3.14761e-03
# <Ncoll>(d+Au) syst err = 3.00000e-01
# 
# 
# Columns: x = p_T (GeV/c), y = R_dAu, statistical error on R_dAu
# 
#   p_T    	   R_dAu   	 stat err
1.50000e-01	6.61171e-01	6.24320e-03
2.50000e-01	6.91226e-01	4.75560e-03
3.50000e-01	7.41091e-01	5.43828e-03
4.50000e-01	7.44109e-01	5.75386e-03
5.50000e-01	7.49973e-01	6.49302e-03
6.50000e-01	7.62249e-01	7.12697e-03
7.50000e-01	7.77900e-01	8.14004e-03
8.50000e-01	7.94962e-01	9.63524e-03
9.50000e-01	8.28468e-01	1.21421e-02
1.05000e+00	8.62651e-01	1.55600e-02
1.15000e+00	8.72615e-01	1.75763e-02
1.25000e+00	8.70301e-01	1.86587e-02
1.35000e+00	9.07216e-01	2.10561e-02
1.45000e+00	9.29579e-01	2.59257e-02
1.55000e+00	9.91418e-01	2.90215e-02
1.65000e+00	9.89654e-01	3.43201e-02
1.75000e+00	1.03631e+00	5.28507e-02
1.85000e+00	9.89804e-01	4.21655e-02
1.95000e+00	9.95756e-01	4.60700e-02
2.10000e+00	1.09347e+00	5.94173e-02
2.30000e+00	1.11524e+00	7.05496e-02
2.50000e+00	1.06021e+00	6.61864e-02
2.70000e+00	1.00134e+00	8.86518e-02
2.90000e+00	1.28159e+00	1.59205e-01
3.12500e+00	1.41420e+00	1.49119e-01
3.37500e+00	9.79826e-01	1.42648e-01
3.62500e+00	1.08757e+00	1.94867e-01
3.87500e+00	1.57946e+00	3.43513e-01
# 
# ----------
# Centrality bin: 70-100% of d+Au cross section
# 
# <Ncoll>(p+Au)          = 1.70586e+00
# <Ncoll>(p+Au) stat err = 1.33009e-03
# <Ncoll>(p+Au) syst err = 3.00000e-01
# 
# <Ncoll>(n+Au)          = 1.70425e+00
# <Ncoll>(n+Au) stat err = 1.32219e-03
# <Ncoll>(n+Au) syst err = 3.00000e-01
# 
# Systematic error on R_pnAu (due to tagging) = 7%
# 
# 
# Columns: x = p_T (GeV/c), y = R_pnAu, statistical error on R_pnAu
# 
#   p_T    	   R_pnAu   	 stat err
1.50000e-01	7.04538e-01	1.02512e-02
2.50000e-01	7.44250e-01	7.73411e-03
3.50000e-01	8.06893e-01	8.79498e-03
4.50000e-01	7.96742e-01	9.30972e-03
5.50000e-01	8.04367e-01	1.03873e-02
6.50000e-01	8.15571e-01	1.21513e-02
7.50000e-01	8.46771e-01	1.33993e-02
8.50000e-01	8.49715e-01	1.60558e-02
9.50000e-01	8.82049e-01	2.10264e-02
1.05000e+00	9.16786e-01	2.61762e-02
1.15000e+00	9.20566e-01	2.83900e-02
1.25000e+00	9.16214e-01	2.89517e-02
1.35000e+00	9.55992e-01	3.29494e-02
1.45000e+00	9.31624e-01	3.79477e-02
1.55000e+00	1.06042e+00	5.15295e-02
1.65000e+00	1.04494e+00	5.79855e-02
1.75000e+00	1.18173e+00	1.07345e-01
1.85000e+00	1.09573e+00	7.75938e-02
1.95000e+00	9.77144e-01	6.55192e-02
2.10000e+00	1.27367e+00	1.36593e-01
2.30000e+00	1.12808e+00	1.04185e-01
2.50000e+00	1.03539e+00	9.39916e-02
2.70000e+00	8.31029e-01	9.95728e-02
2.90000e+00	1.10144e+00	1.51597e-01
3.12500e+00	1.57730e+00	2.80898e-01
3.37500e+00	1.00996e+00	2.43621e-01
3.62500e+00	9.35883e-01	2.46462e-01
3.87500e+00	2.21634e+00	7.19802e-01
# 
# ----------
# Centrality bin: 40-70% of d+Au cross section
# <Ncoll>(d+Au)          = 4.98519e+00
# <Ncoll>(d+Au) stat err = 6.53975e-03
# <Ncoll>(d+Au) syst err = 1.50000e-01
# 
# 
# Columns: x = p_T (GeV/c), y = R_dAu, statistical error on R_dAu
# 
#   p_T    	   R_dAu   	 stat err
1.50000e-01	5.27215e-01	3.49475e-03
2.50000e-01	5.38532e-01	2.82423e-03
3.50000e-01	5.64846e-01	3.13524e-03
4.50000e-01	5.75083e-01	3.37035e-03
5.50000e-01	5.89304e-01	3.71731e-03
6.50000e-01	6.19959e-01	4.33082e-03
7.50000e-01	6.40099e-01	5.08122e-03
8.50000e-01	6.78579e-01	6.13576e-03
9.50000e-01	7.06629e-01	7.34966e-03
1.05000e+00	7.43244e-01	8.79122e-03
1.15000e+00	7.58770e-01	1.04832e-02
1.25000e+00	8.18340e-01	1.26923e-02
1.35000e+00	8.26511e-01	1.48400e-02
1.45000e+00	8.41962e-01	1.73600e-02
1.55000e+00	8.88177e-01	2.07986e-02
1.65000e+00	9.54792e-01	2.48675e-02
1.75000e+00	9.58180e-01	2.82115e-02
1.85000e+00	1.05152e+00	3.37366e-02
1.95000e+00	9.23619e-01	3.57644e-02
2.10000e+00	9.55296e-01	3.16410e-02
2.30000e+00	1.10397e+00	4.35453e-02
2.50000e+00	1.11218e+00	5.58018e-02
2.70000e+00	1.04924e+00	6.92808e-02
2.90000e+00	1.17721e+00	9.00965e-02
3.12500e+00	1.02514e+00	9.38442e-02
3.37500e+00	1.09520e+00	1.29778e-01
3.62500e+00	1.11990e+00	1.63852e-01
3.87500e+00	9.61714e-01	1.94733e-01
# 
# ----------
# Centrality bin: 40-70% of d+Au cross section
# 
# <Ncoll>(p+Au)          = 3.69859e+00
# <Ncoll>(p+Au) stat err = 4.09001e-03
# <Ncoll>(p+Au) syst err = 1.50000e-01
# 
# <Ncoll>(n+Au)          = 3.73205e+00
# <Ncoll>(n+Au) stat err = 4.16696e-03
# <Ncoll>(n+Au) syst err = 1.50000e-01
# 
# Systematic error on R_pnAu (due to tagging) = 7%
# 
# 
# Columns: x = p_T (GeV/c), y = R_pnAu, statistical error on R_pnAu
# 
#   p_T    	   R_pnAu   	 stat err
1.50000e-01	5.23828e-01	8.04660e-03
2.50000e-01	5.51038e-01	6.51419e-03
3.50000e-01	5.80400e-01	7.32137e-03
4.50000e-01	5.87655e-01	7.81475e-03
5.50000e-01	6.08141e-01	8.75810e-03
6.50000e-01	6.37617e-01	1.01884e-02
7.50000e-01	6.62985e-01	1.19904e-02
8.50000e-01	6.64039e-01	1.41449e-02
9.50000e-01	7.23109e-01	1.73967e-02
1.05000e+00	7.41969e-01	2.03060e-02
1.15000e+00	7.61990e-01	2.46071e-02
1.25000e+00	8.46606e-01	3.03313e-02
1.35000e+00	7.85639e-01	3.32692e-02
1.45000e+00	8.48669e-01	4.06807e-02
1.55000e+00	8.35480e-01	4.68164e-02
1.65000e+00	1.06101e+00	6.25846e-02
1.75000e+00	9.23448e-01	6.60803e-02
1.85000e+00	1.12172e+00	8.10731e-02
1.95000e+00	8.79554e-01	8.26653e-02
2.10000e+00	9.58251e-01	7.46697e-02
2.30000e+00	1.05625e+00	9.93084e-02
2.50000e+00	9.87295e-01	1.24341e-01
2.70000e+00	1.38706e+00	1.94250e-01
2.90000e+00	1.28276e+00	2.27346e-01
3.12500e+00	7.51005e-01	1.85689e-01
3.37500e+00	1.23146e+00	3.55175e-01
3.62500e+00	1.38909e+00	4.39212e-01
3.87500e+00	1.61871e+00	6.71792e-01
# 
# ----------
# Centrality bin: 20-40% of d+Au cross section
# <Ncoll>(d+Au)          = 9.42774e+00
# <Ncoll>(d+Au) stat err = 1.01282e-02
# <Ncoll>(d+Au) syst err = 1.00000e-01
# 
# 
# Columns: x = p_T (GeV/c), y = R_dAu, statistical error on R_dAu
# 
#   p_T    	   R_dAu   	 stat err
1.50000e-01	4.23084e-01	2.40418e-03
2.50000e-01	4.23760e-01	1.93056e-03
3.50000e-01	4.49448e-01	2.16229e-03
4.50000e-01	4.59910e-01	2.32469e-03
5.50000e-01	4.75485e-01	2.58340e-03
6.50000e-01	5.03088e-01	3.01296e-03
7.50000e-01	5.33510e-01	3.60079e-03
8.50000e-01	5.63548e-01	4.34044e-03
9.50000e-01	5.95274e-01	5.22323e-03
1.05000e+00	6.44386e-01	6.37773e-03
1.15000e+00	6.60396e-01	7.58339e-03
1.25000e+00	6.96591e-01	9.09946e-03
1.35000e+00	7.37858e-01	1.09151e-02
1.45000e+00	7.67245e-01	1.29260e-02
1.55000e+00	7.99663e-01	1.54295e-02
1.65000e+00	8.17489e-01	1.78660e-02
1.75000e+00	8.88289e-01	2.12902e-02
1.85000e+00	8.91641e-01	2.42159e-02
1.95000e+00	9.34650e-01	2.84515e-02
2.10000e+00	9.52898e-01	2.48273e-02
2.30000e+00	9.48691e-01	3.18443e-02
2.50000e+00	9.67059e-01	4.05313e-02
2.70000e+00	1.03592e+00	5.30547e-02
2.90000e+00	1.04454e+00	6.63145e-02
3.12500e+00	1.08731e+00	7.71874e-02
3.37500e+00	1.13236e+00	1.00827e-01
3.62500e+00	1.28406e+00	1.35861e-01
3.87500e+00	1.07030e+00	1.59823e-01
# 
# ----------
# Centrality bin: 20-40% of d+Au cross section
# 
# <Ncoll>(p+Au)          = 6.72290e+00
# <Ncoll>(p+Au) stat err = 1.27128e-02
# <Ncoll>(p+Au) syst err = 1.00000e-01
# 
# <Ncoll>(n+Au)          = 6.74715e+00
# <Ncoll>(n+Au) stat err = 1.27973e-02
# <Ncoll>(n+Au) syst err = 1.00000e-01
# 
# Systematic error on R_pnAu (due to tagging) = 7%
# 
# 
# Columns: x = p_T (GeV/c), y = R_pnAu, statistical error on R_pnAu
# 
#   p_T    	   R_pnAu   	 stat err
1.50000e-01	4.35708e-01	9.84650e-03
2.50000e-01	4.27388e-01	7.76351e-03
3.50000e-01	4.48461e-01	8.56586e-03
4.50000e-01	4.49989e-01	9.12707e-03
5.50000e-01	4.81660e-01	1.05209e-02
6.50000e-01	5.03005e-01	1.20749e-02
7.50000e-01	5.46479e-01	1.48373e-02
8.50000e-01	5.75905e-01	1.79948e-02
9.50000e-01	5.84708e-01	2.08970e-02
1.05000e+00	6.54091e-01	2.64379e-02
1.15000e+00	5.90334e-01	2.88997e-02
1.25000e+00	6.48389e-01	3.58819e-02
1.35000e+00	6.96392e-01	4.40917e-02
1.45000e+00	7.48098e-01	5.27550e-02
1.55000e+00	8.44914e-01	6.49295e-02
1.65000e+00	7.92490e-01	7.32385e-02
1.75000e+00	8.41514e-01	8.29361e-02
1.85000e+00	8.22797e-01	9.33563e-02
1.95000e+00	9.17855e-01	1.12956e-01
2.10000e+00	9.75397e-01	1.04301e-01
2.30000e+00	1.17020e+00	1.44874e-01
2.50000e+00	1.00848e+00	1.74920e-01
2.70000e+00	1.03561e+00	2.09419e-01
2.90000e+00	4.97812e-01	1.99904e-01
3.12500e+00	7.79151e-01	2.59512e-01
3.37500e+00	1.06352e+00	3.90212e-01
3.62500e+00	1.76677e+00	7.62390e-01
3.87500e+00	1.27953e+00	6.11016e-01
# 
# ----------
# Centrality bin: 0-20% of d+Au cross section
# <Ncoll>(d+Au)          = 1.44910e+01
# <Ncoll>(d+Au) stat err = 1.27472e-02
# <Ncoll>(d+Au) syst err = 7.50000e-02
# 
# 
# Columns: x = p_T (GeV/c), y = R_dAu, statistical error on R_dAu
# 
#   p_T    	   R_dAu   	 stat err
1.50000e-01	3.66567e-01	1.75247e-03
2.50000e-01	3.62558e-01	1.38873e-03
3.50000e-01	3.90904e-01	1.57519e-03
4.50000e-01	3.95534e-01	1.67484e-03
5.50000e-01	4.12852e-01	1.87088e-03
6.50000e-01	4.41239e-01	2.19894e-03
7.50000e-01	4.68256e-01	2.62931e-03
8.50000e-01	5.02508e-01	3.19721e-03
9.50000e-01	5.39061e-01	3.88739e-03
1.05000e+00	5.60632e-01	4.63284e-03
1.15000e+00	6.11901e-01	5.69433e-03
1.25000e+00	6.32775e-01	6.77268e-03
1.35000e+00	6.64530e-01	8.11308e-03
1.45000e+00	7.00180e-01	9.64393e-03
1.55000e+00	7.29794e-01	1.14595e-02
1.65000e+00	7.63877e-01	1.35011e-02
1.75000e+00	7.56713e-01	1.52799e-02
1.85000e+00	7.58568e-01	1.74943e-02
1.95000e+00	8.30900e-01	2.08927e-02
2.10000e+00	8.80785e-01	1.86476e-02
2.30000e+00	9.14845e-01	2.44460e-02
2.50000e+00	9.76592e-01	3.21127e-02
2.70000e+00	8.99174e-01	3.86922e-02
2.90000e+00	9.54979e-01	4.99426e-02
3.12500e+00	1.05013e+00	5.87063e-02
3.37500e+00	1.07094e+00	7.73389e-02
3.62500e+00	1.09168e+00	9.78813e-02
3.87500e+00	1.16408e+00	1.30461e-01
# 
# ----------
# Centrality bin: 0-20% of d+Au cross section
# 
# <Ncoll>(p+Au)          = 8.35508e+00
# <Ncoll>(p+Au) stat err = 3.93495e-02
# <Ncoll>(p+Au) syst err = 7.50000e-02
# 
# <Ncoll>(n+Au)          = 8.35667e+00
# <Ncoll>(n+Au) stat err = 3.92462e-02
# <Ncoll>(n+Au) syst err = 7.50000e-02
# 
# Systematic error on R_pnAu (due to tagging) = 7%
# 
# 
# Columns: x = p_T (GeV/c), y = R_pnAu, statistical error on R_pnAu
# 
#   p_T    	   R_pnAu   	 stat err
1.50000e-01	4.29827e-01	1.56856e-02
2.50000e-01	4.28221e-01	1.23274e-02
3.50000e-01	4.48513e-01	1.35003e-02
4.50000e-01	4.19992e-01	1.38382e-02
5.50000e-01	4.52118e-01	1.53809e-02
6.50000e-01	4.74918e-01	1.86756e-02
7.50000e-01	5.80437e-01	2.39516e-02
8.50000e-01	6.35512e-01	2.97972e-02
9.50000e-01	6.03166e-01	3.38894e-02
1.05000e+00	5.92624e-01	3.94840e-02
1.15000e+00	6.55401e-01	4.73885e-02
1.25000e+00	8.51772e-01	6.53958e-02
1.35000e+00	8.17661e-01	7.75327e-02
1.45000e+00	8.75849e-01	9.09264e-02
1.55000e+00	8.98959e-01	1.07876e-01
1.65000e+00	8.56133e-01	1.17176e-01
1.75000e+00	8.38107e-01	1.28343e-01
1.85000e+00	6.17754e-01	1.41679e-01
1.95000e+00	8.07584e-01	1.68535e-01
2.10000e+00	9.72674e-01	1.66093e-01
2.30000e+00	9.48603e-01	2.11971e-01
2.50000e+00	1.37696e+00	3.17822e-01
2.70000e+00	1.11905e+00	3.31057e-01
2.90000e+00	1.11104e+00	5.03526e-01
3.12500e+00	1.03577e+00	4.73368e-01
3.37500e+00	1.02716e+00	5.99599e-01
3.62500e+00	8.74582e-01	5.06439e-01
3.87500e+00	1.98235e+00	1.65651e+00