========================================
 AUTHOR:      Antoine Rocher       2019		        
 CONTRIBUTOR: Michel-Andrès Breton 2020		        
 ========================================

This C code implements the CLPT and Gaussian streaming model for redshift-space distortions 
presented in Reid & White 2011 (https://arxiv.org/pdf/1105.4165.pdf) and Wang, Reid 
& White 2014 (https://arxiv.org/pdf/1306.1804.pdf). It includes Alcock-Paczynski
distortions through alpha_perp and alpha_par parameters. The final products are the 
monopole, quadrupole, and hexadecapole moments of the redshift-space correlation function.   

*** Compilation *** 

In the code folder run:
make

To build it as a library: 
make lib

*** CLPT ***

*Need to create a data/ folder to store CLPT predictions*

To get CLPT predictions:
./CLPT powerspectrum_file.dat

This has to be done once per input linear power spectrum.

The outputs in the data/ folder include:

** Xi_L_func.dat 
| Column | Value                                 |
|--------+---------------------------------------|
|      1 | $r$ (in $\mathrm{Mpc\ h}^{-1}$)       |    				
|      2 | $\xi_L	

***Xi_r_CLPT.dat
| Column | Value                                 
|--------+---------------------------------------|
|      1 | $r$ (in $\mathrm{Mpc\ h}^{-1}$)       
|      2 | \xi_m_f0
|      3 | \xi_m_f1 	       
|      4 | \xi_m_f2
|      5 | \xi_m_f1f1 	      
|      6 | \xi_m_f1f2 		
|      7 | \xi_m_f2f2 		

***V_12_CLPT.dat
| Column | Value                                 
|--------+---------------------------------------|
|      1 | $r$ (in $\mathrm{Mpc\ h}^{-1}$)       
|      2 | \v12_f0
|      3 | \v12_f1 	       
|      4 | \v12_f2
|      5 | \v12_f1f1 	      
|      6 | \v12_f1f2 			 

***Sigma_12_CLPT.dat
| Column | Value                                 
|--------+---------------------------------------|
|      1 | $r$ (in $\mathrm{Mpc\ h}^{-1}$)       
|      2 | \s12_par_f0
|      3 | \s12_par_f1 	       
|      4 | \s12_par_f2
|      5 | \s12_par_f1f1 	      
|      6 | \s12_perp_f0
|      7 | \s12_perp_f1 	       
|      8 | \s12_perp_f2
|      9 | \s12_perp_f1f1 	      	 


*** Gaussian Streaming for redshift-space distortions ***

The Gaussian streaming predictions are obtained by running:
./GS parameter_file.txt

The parameter file contains:
smin		  -> minimum s scale
smax		  -> maximum s scale
nbins		  -> number of bins in s between smin and smax							
f (growth rate)   -> linear growth rate of structure parameter				
f1		  -> Lagrangian bias parameter 1
f2		  -> Lagrangian bias parameter 2
sig_shift^2	  -> small-scale velocity dispersion squared 
alpha_perp	  -> Alcock-Paczynski perpendicular distortion parameter
alpha_para	  -> Alcock-Paczynski parallel distortion parameter				
input_dir 	  -> input directory where CLPT predictions are stored

The output is stored in the data/ folder in the file:

*** multipoles_CLPT.dat :
| Column | Value                                 
|--------+---------------------------------------|
|      1 | $r$ (in $\mathrm{Mpc\ h}^{-1}$)       
|      2 | $\xi_0											 
|      3 | $\xi_2 		                               				
|      4 | $\xi_4			                     	       
			                     	       

** Contact
For any question, please send an email to michel-andres.breton@lam.fr