diff --git a/README b/README
new file mode 100644
--- /dev/null
+++ b/README
@@ -0,0 +1,25 @@
+*************************************************************************************************
+QCD Uncertainties on Particle Spectra from Dark Matter Annihilation
+Authors: Simone Amoroso, Sascha Caron, Adil Jueid, Roberto Ruiz de Austri, and Peter Skands
+If you use these tables, please cite:
+S. Amoroso et al. arXiv: 1812.07424 [hep-ph], JCAP05(2019)007
+*************************************************************************************************
+
+We provide recipes for particle spectra in dark matter indirect searches in the form of tables. A wide range of dark matter masses from 10 GeV to 100 TeVis covered. We consider 11 primary annihilation channels:
+DM DM -> e+e-, mu+ mu-, tau tau, qq  (q=u,d,s), cc, bb, tt, WW, gg, ZZ, and hh.
+Each file contains 13 columns: the dark matter mass, the Log of energy fraction (kinetic energy of the particle over the dark mass), and normalized spectra for 11 primary channels]. We report tables for the spectra of photons, positrons, electron anti-neutrinos, muon anti-neutrinos and tau anti-neutrinos.The spectra of anti-protons, and anti-atoms will be added soon.
+
+For each particle species, we provide twelve tables which can be found in zip format. The notation of the different tables is given below:
+     1) The table corresponding to the central prediction for the spectra is denoted by 'AtProduction-Hadronization1-$TYPE.dat' with
+        $TYPE=Nuel, Numu, Nuta, Ga, Positrons which refers to the three flavors of neutrinos, photons and positrons respectively.
+     2) There are nine tables corresponding to the different variations of the light quark fragmentation function's parameters. These
+	tables are denoted by 'AtProduction-Hadronization$h-$TYPE.dat' with h=2,..,10.
+     3) The particle spectra corresponding to the variations of the shower evolution scale are denoted by AtProduction-Shower-Var$s-$TYPE.dat
+	with s=1,2  corresponds to 1/2 mu_R and 2 mu_R (mu_R is the shower evolution scale in PYTHIA8). 
+
+Important advices:
+     i)    Uncertainty on the spectra, from hadronization, is obtained from the envelope of all the variations (including the central prediction).
+     ii)   In the variations of the parton shower evolution scale, the parameters of the hadronization function are fixed to their central value.
+     iii)  In principle, showering uncertainties are uncorrelated to hadronization uncertainties. To obtain the full uncertainty, one might combine
+           those uncertainties in quadrature.
+	 
\ No newline at end of file