+
+// x-sections for pp @ 7 TeV:
+// -charmonia: 4pi integral of fit function for inclusive J/psi dsigma/dy LHC data
+// gives 60 mub; so sigma_prompt = 54 mub, while Ref = R.Vogt_arXiv:1003.3497 (Table 2)
+// gives 35 mub. Below we use sigma_direct from the Ref scaled by the factor 54/35.
+// -bottomonia: 4pi integral of fit function for inclusive Upsilon1S dsigma/dy LHC data
+// gives 0.56 mub, sigmas for 2S & 3S obtained using LHCb data for ratios 2S/1S & 3S/1S
+// -ccbar & bbbar: NLO pQCD computations - http://www-alice.gsi.de/ana/MNR/results.html
+ fCMSEnergyTeVArray[0] = 7.00;
+ fSigmaReactionArray[0] = 0.070;
+ fSigmaJPsiArray[0] = 33.6e-6;
+ fSigmaChic1Array[0] = 32.6e-6;
+ fSigmaChic2Array[0] = 53.8e-6;
+ fSigmaPsiPArray[0] = 7.6e-6;
+ fSigmaUpsilonArray[0] = 0.56e-6;
+ fSigmaUpsilonPArray[0] = 0.18e-6;
+ fSigmaUpsilonPPArray[0] = 0.08e-6;
+ fSigmaCCbarArray[0] = 6.91e-3;
+ fSigmaBBbarArray[0] = 0.232e-3;
+
+//x-sections for pp @ 10 TeV: charmonia and bottomonia from 14 TeV numbers
+// scaled down according to ccbar and bbbar cross-sections
+ fCMSEnergyTeVArray[1] = 10.00;
+ fSigmaReactionArray[1] = 0.070;
+ fSigmaJPsiArray[1] = 26.06e-6;
+ fSigmaChic1Array[1] = 25.18e-6;
+ fSigmaChic2Array[1] = 41.58e-6;
+ fSigmaPsiPArray[1] = 5.88e-6;
+ fSigmaUpsilonArray[1] = 0.658e-6;
+ fSigmaUpsilonPArray[1] = 0.218e-6;
+ fSigmaUpsilonPPArray[1] = 0.122e-6;
+ fSigmaCCbarArray[1] = 8.9e-3;
+ fSigmaBBbarArray[1] = 0.33e-3;
+
+//x-sections for pp @ 14 TeV: charmonia from hep-ph/0311048 Tab.9, page 19,
+// bottomonium from hep-ph/0311048 Tab.9, page 19 taken into account that
+// feed-down from chib is included
+ fCMSEnergyTeVArray[2] = 14.00;
+ fSigmaReactionArray[2] = 0.070;
+ fSigmaJPsiArray[2] = 32.9e-6;
+ fSigmaChic1Array[2] = 31.8e-6;
+ fSigmaChic2Array[2] = 52.5e-6;
+ fSigmaPsiPArray[2] = 7.43e-6;
+ fSigmaUpsilonArray[2] = 0.989e-6;
+ fSigmaUpsilonPArray[2] = 0.502e-6;
+ fSigmaUpsilonPPArray[2] = 0.228e-6;
+ fSigmaCCbarArray[2] = 11.2e-3;
+ fSigmaBBbarArray[2] = 0.445e-3;
+
+// x-sections for Min. Bias p-Pb & Pb-p @ 2.76 TeV: charmonia and bottomonia
+// from 7 TeV numbers scaled according to pQCD ccbar and bbbar x-sections
+// and with Glauber scaling
+ fCMSEnergyTeVArray[3] = 2.00; // for 2.76 TeV
+ fSigmaReactionArray[3] = 2.10;
+ fSigmaJPsiArray[3] = 3.54e-3; // 208*0.507*33.6e-6
+ fSigmaChic1Array[3] = 3.44e-3;
+ fSigmaChic2Array[3] = 5.66e-3;
+ fSigmaPsiPArray[3] = 0.80e-3;
+ fSigmaUpsilonArray[3] = 0.045e-3; // 208*0.384*0.56e-6
+ fSigmaUpsilonPArray[3] = 0.014e-3;
+ fSigmaUpsilonPPArray[3] = 0.006e-3;
+ fSigmaCCbarArray[3] = 0.73; // 208*3.50e-3
+ fSigmaBBbarArray[3] = 0.019; // 208*0.089e-3
+
+ fCMSEnergyTeVArray[4] = -fCMSEnergyTeVArray[3];
+ fSigmaReactionArray[4] = fSigmaReactionArray[3];
+ fSigmaJPsiArray[4] = fSigmaJPsiArray[3];
+ fSigmaChic1Array[4] = fSigmaChic1Array[3];
+ fSigmaChic2Array[4] = fSigmaChic2Array[3];
+ fSigmaPsiPArray[4] = fSigmaPsiPArray[3];
+ fSigmaUpsilonArray[4] = fSigmaUpsilonArray[3];
+ fSigmaUpsilonPArray[4] = fSigmaUpsilonPArray[3];
+ fSigmaUpsilonPPArray[4] = fSigmaUpsilonPPArray[3];
+ fSigmaCCbarArray[4] = fSigmaCCbarArray[3];
+ fSigmaBBbarArray[4] = fSigmaBBbarArray[3];
+
+// x-sections for Min. Bias p-Pb & Pb-p @ 4.4 TeV: charmonia and bottomonia
+// from 7 TeV numbers scaled according to pQCD ccbar and bbbar x-sections
+// and with Glauber scaling
+ fCMSEnergyTeVArray[5] = 4.00; // for 4.4 TeV
+ fSigmaReactionArray[5] = 2.10;
+ fSigmaJPsiArray[5] = 5.00e-3; // 208*0.715*33.6e-6
+ fSigmaChic1Array[5] = 4.86e-3;
+ fSigmaChic2Array[5] = 7.99e-3;
+ fSigmaPsiPArray[5] = 1.12e-3;
+ fSigmaUpsilonArray[5] = 0.074e-3; // 208*0.629*0.56e-6
+ fSigmaUpsilonPArray[5] = 0.023e-3;
+ fSigmaUpsilonPPArray[5] = 0.010e-3;
+ fSigmaCCbarArray[5] = 1.03; // 208*4.94e-3
+ fSigmaBBbarArray[5] = 0.030; // 208*0.146e-3
+
+ fCMSEnergyTeVArray[6] = -fCMSEnergyTeVArray[5];
+ fSigmaReactionArray[6] = fSigmaReactionArray[5];
+ fSigmaJPsiArray[6] = fSigmaJPsiArray[5];
+ fSigmaChic1Array[6] = fSigmaChic1Array[5];
+ fSigmaChic2Array[6] = fSigmaChic2Array[5];
+ fSigmaPsiPArray[6] = fSigmaPsiPArray[5];
+ fSigmaUpsilonArray[6] = fSigmaUpsilonArray[5];
+ fSigmaUpsilonPArray[6] = fSigmaUpsilonPArray[5];
+ fSigmaUpsilonPPArray[6] = fSigmaUpsilonPPArray[5];
+ fSigmaCCbarArray[6] = fSigmaCCbarArray[5];
+ fSigmaBBbarArray[6] = fSigmaBBbarArray[5];
+
+// x-sections for Min. Bias p-Pb & Pb-p @ 5.03 TeV: charmonia and bottomonia
+// from 7 TeV numbers scaled according to pQCD ccbar and bbbar x-sections
+// and with Glauber scaling
+ fCMSEnergyTeVArray[7] = 5.00; // for 5.03 TeV
+ fSigmaReactionArray[7] = 2.10;
+ fSigmaJPsiArray[7] = 5.50e-3; // 208*0.787*33.6e-6
+ fSigmaChic1Array[7] = 5.35e-3;
+ fSigmaChic2Array[7] = 8.79e-3;
+ fSigmaPsiPArray[7] = 1.23e-3;
+ fSigmaUpsilonArray[7] = 0.083e-3; // 208*0.716*0.56e-6
+ fSigmaUpsilonPArray[7] = 0.026e-3;
+ fSigmaUpsilonPPArray[7] = 0.011e-3;
+ fSigmaCCbarArray[7] = 1.13; // 208*5.44e-3
+ fSigmaBBbarArray[7] = 0.035; // 208*0.166e-3
+
+ fCMSEnergyTeVArray[8] = -fCMSEnergyTeVArray[7];
+ fSigmaReactionArray[8] = fSigmaReactionArray[7];
+ fSigmaJPsiArray[8] = fSigmaJPsiArray[7];
+ fSigmaChic1Array[8] = fSigmaChic1Array[7];
+ fSigmaChic2Array[8] = fSigmaChic2Array[7];
+ fSigmaPsiPArray[8] = fSigmaPsiPArray[7];
+ fSigmaUpsilonArray[8] = fSigmaUpsilonArray[7];
+ fSigmaUpsilonPArray[8] = fSigmaUpsilonPArray[7];
+ fSigmaUpsilonPPArray[8] = fSigmaUpsilonPPArray[7];
+ fSigmaCCbarArray[8] = fSigmaCCbarArray[7];
+ fSigmaBBbarArray[8] = fSigmaBBbarArray[7];
+
+// x-sections for Min. Bias p-Pb & Pb-p @ 8.8 TeV: charmonia and bottomonia
+// from 7 TeV numbers scaled according to pQCD ccbar and bbbar x-sections
+// and with Glauber scaling
+ fCMSEnergyTeVArray[9] = 9.00; // for 8.8 TeV
+ fSigmaReactionArray[9] = 2.10;
+ fSigmaJPsiArray[9] = 8.19e-3; // 208*1.172*33.6e-6
+ fSigmaChic1Array[9] = 7.95e-3;
+ fSigmaChic2Array[9] = 13.1e-3;
+ fSigmaPsiPArray[9] = 1.85e-3;
+ fSigmaUpsilonArray[9] = 0.146e-3; // 208*1.25*0.56e-6
+ fSigmaUpsilonPArray[9] = 0.047e-3;
+ fSigmaUpsilonPPArray[9] = 0.021e-3;
+ fSigmaCCbarArray[9] = 1.68; // 208*8.1e-3
+ fSigmaBBbarArray[9] = 0.061; // 208*0.29e-3
+
+ fCMSEnergyTeVArray[10] = -fCMSEnergyTeVArray[9];
+ fSigmaReactionArray[10] = fSigmaReactionArray[9];
+ fSigmaJPsiArray[10] = fSigmaJPsiArray[9];
+ fSigmaChic1Array[10] = fSigmaChic1Array[9];
+ fSigmaChic2Array[10] = fSigmaChic2Array[9];
+ fSigmaPsiPArray[10] = fSigmaPsiPArray[9];
+ fSigmaUpsilonArray[10] = fSigmaUpsilonArray[9];
+ fSigmaUpsilonPArray[10] = fSigmaUpsilonPArray[9];
+ fSigmaUpsilonPPArray[10] = fSigmaUpsilonPPArray[9];
+ fSigmaCCbarArray[10] = fSigmaCCbarArray[9];
+ fSigmaBBbarArray[10] = fSigmaBBbarArray[9];
+
+// x-sections for Min. Bias Pb-Pb @ 2.76 TeV: charmonia and bottomonia
+// from 7 TeV numbers scaled according to pQCD ccbar and bbbar x-sections
+// and with Glauber scaling
+ fCMSEnergyTeVArray[11] = 3.00; // for 2.76 TeV
+ fSigmaReactionArray[11] = 7.65;
+ fSigmaJPsiArray[11] = 0.737; // 208*208*0.507*33.6e-6
+ fSigmaChic1Array[11] = 0.715;
+ fSigmaChic2Array[11] = 1.179;
+ fSigmaPsiPArray[11] = 0.166;
+ fSigmaUpsilonArray[11] = 0.0093; // 208*208*0.384*0.56e-6
+ fSigmaUpsilonPArray[11] = 0.0030;
+ fSigmaUpsilonPPArray[11] = 0.0013;
+ fSigmaCCbarArray[11] = 151.; // 208*208*3.50e-3
+ fSigmaBBbarArray[11] = 3.8; // 208*208*0.089e-3
+