ISiTGR Readme

This version of ISiTGR is for the December 2013 version of CosmoMC

Introduction

ISiTGR is an integrated set of modified modules for the software package CosmoMC for use in testing whether observational data is consistent with general relativity on cosmological scales. This latest version of the code has been updated to allow for the consideration of non-flat universes. It incorporates modifications to the codes: CAMB, CosmoMC, the ISW-galaxy cross correlation likelihood code of Ho et al, and two of our independently developed weak lensing likelihood codes: one for the CFHTLenS weak lensing tomography of Heymans et al; and one for the refined COSMOS 3D weak lensing tomography of Schrabback et al.

To use ISiTGR you must first download and install CosmoMC (12/2013 version). You can find download instructions, system requirements, and setup instructions for CosmoMC in the CosmoMC ReadMe. Then, simply copy the contents of the unzipped ISiTGR folder to the CosmoMC folder, edit the Makefile for your compilers, compile (extra compiling step here), and you are ready to go.

Basic usage of ISiTGR is the same as CosmoMC. See below for extra options in the params.ini file and added organizational features in ISiTGR.

Differences between ISiTGR and CosmoMC

Other than the changes to the code to incorporate the parameters used to test general relativy as described in our papers: arXiv:1109.4583 and arXiv:1205.2422, the following additional changes have been made

  • ini files

    In both ISiTGR and ISiTGR_BIN the following additions have been made to the test.ini and associated files.
    • To run the ISW-galaxy cross correlation likelihood code the line use_ISW = T has been added to batch1/likelihood_batch1.ini. Set this to F if you do not want to use this likelihood.
    • To use the CFHTLenS weak lensing tomography data set simply make sure the line DEFAULT(batch1/CFHTLens.ini) is uncommented in test.ini. This file enables use of the full CFHTLenS data sample as described in Heymans et al. If you would like to use one of the other galaxy samples instead you can change that line to:
      • "Early Type," Red sample - DEFAULT(batch1/CFHTLens_red.ini)
      • "Late Type," Blue sample - DEFAULT(batch1/CFHTLens_blu.ini)
      • "Optomized Early Type" sample - DEFAULT(batch1/CFHTLens_rfbb.ini)
      These files each contain various settings to enable the different data sets as well as the default parameter ranges for the CFHTLenS intrinsic alignment (IA) nuissance parameter, ACFHTLenS, which allows this parameter to vary from -10 to 10. If you would like to fix this parameter to 0 and thus ignore the (IA) signal you can add the following line at the end of test.ini:

      param[CFHTLensA] = 0 0 0 0 0

    • To use the COSMOS 3D weak lensing tomography data set simply make sure the line DEFAULT(batch1/WeakLen.ini) is uncommented in test.ini. This file includes settings to vary the weak lensing nuissance parameter fz which accounts for a 10% uncertainty in the redshift distribution of galaxies in the the sixth redshift bin of the refined COSMOS data.

    Additional parameters and options for testing general relativity: ISiTGR

    The following options have been included in batch1/params_TGR_defaults.ini

    • The option: Use_R_Function = lets you decide which parameters for testing GR you want to evolve using the funcitonal form shown in our paper. Setting this option to T evolves the parameters Q and R, while setting this option to F evolves Q and Σ (previously called D).
    • For varying the parameters used to test GR we add the following parameter lines, please refer to our paper for a detailed explanation of each of these parameters.

      param[Q0] = 1 0 10 0.05 0.05 ## 1 1 1 0 0

      param[Qinf] 1 1 1 0 0 ## 1 0 10 0.05 0.05

      param[Sigma0] = 1 0 10 0.05 0.05 ## 1 1 1 0 0

      param[Sigmainf] = 1 1 1 0 0 ## 1 0 10 0.05 0.05

      param[kc] = 0.01 0.01 0.01 0 0

      param[s] = 0 0 3 -1 1 ## 0 0 0 0 0 ##

    • The option: Scale_Dependent = lets you decide whether the evolution of the parmaters is scale dependent. Setting this option to T enforces scale dependence, and the parameters ending in inf above should be varied or changed to whatever value you want the parameters to take on small scales.

    Additional parameters and options for testing general relativity: ISiTGR_BIN

    The following options have been included in batch1/params_TGR_defaults.ini

    • For varying the parameters used to test GR we add the following parameter lines, please refer to our paper for a detailed explanation of each of these parameters.

      param[Q1] = 1 0 10 0.05 0.05

      param[Q2] = 1 0 10 0.05 0.05

      param[Q3] = 1 0 10 0.05 0.05

      param[Q4] = 1 0 10 0.05 0.05

      param[Sigma1] = 1 0 10 0.05 0.05

      param[Sigma2] = 1 0 10 0.05 0.05

      param[Sigma3] = 1 0 10 0.05 0.05

      param[Sigma4] = 1 0 10 0.05 0.05

      param[kc] = 0.01 0.01 0.01 0 0

    • The option: z_grid_spacing = defines the size of your bins in redshift bins.
    • The option: Do_Exponential_Binning = lets you choose how you want to transition between scale bins. Setting this option to T sets transitions between scale bins to behave as an exponential function with a decay constant kc (This is the hybrid method described in our paper). Setting this option to F sets transitions between scale bins to behave as a hyberbolic tangent functions (near step functions) with the bins divided at k = kc.
  • ISW-galaxy cross correlation likelihood code

    To implement the ISW-galaxy cross correlations likelihood code of Ho et al, the following files have been added to source/: lrg_2dCl.f90 and iswdata.f90 . The files jl.f90 and lrg_pk.f90 have been removed in this latest version of the code.
    The folder bdndz_code which contains c-codes called by the routine Iswlnlike in iswdata.f90 has been added to the root directory. Before running either version of ISiTGR for the first time, these c-codes must be compiled. This can be done by changing the line TCC = to the appropriate compiler call, and running make within the bdndz_code directory.
    When Iswlnlike calls the c-program dndz.x it passes input files and output file names to the program. By default, these files are saved in the IO directory within the bdndz_code directory, however you can enable writing of these files to a different directory by setting the environmental variable TMPDIR.
  • Weak Lensing Likelihood code

    The various weak lensing likelihood codes and the associated subroutines are contained in CFHTLens.f90, AngularCorr.f90, WeakLen.f90, and AngularCorrSchr.f90 which are located in the source directory.
    The first two of these files are for the CFHTLenS weak lensing tomography of Heymans et al, while the latter two are for the data of Schrabback et al.
    All of these files are based upon the files contained in the original COSMOS weak lensing likelihood code of Lesgourges et al.
  • Other Stuff

    We have also added the files fid.ini and batch1/*_fid.ini for testing the code. These files contain settings which fix the Planck 2013 cosmology and runs the vaious likelihoods in order to test the code.

Referencing ISiTGR

We would ask that when using ISiTGR or a modified version of it, you cite: our papers (arXiv:1109.4583 and arXiv:1205.2422), this website; the original CAMB paper; the original CosmoMC paper; and the papers on original ISW-galaxy cross correlation likelihood by Ho et al and Hirata et al. When using the CFHTLenS weak lensing tomography likelihood, please cite our paper, Heymans et al, and follow the CFHTLenS publication guidlines given here. If the refined HST COSMOS 3D weak lensing tomography likelihood in our code is used, the papers by Schrabback et al., Scoville et al., and Ilbert et al should be cited.  Additionally please cite the use of any other datasets already included in the original version of CosmoMC.

Version History

1.2 (Released 01/25/15, this version): Updates for Dec. 2013 version of CosmoMC, new CFHTLenS likelihood module.

1.1.1 (Released 05/29/13): Includes bug fixes (Special thanks to Ana Caramete and Lucia Popa).

1.1: Updated code to allow for non-flat universes when varying the MG parameters.

1.02.2 (Released 05/29/13): Includes bug fixes.

1.02: Updated to January 2012 version of CosmoMC.

1.01.2 (Released 05/29/13): Includes bug fixes.

1.01: Updated to October 2011 version of CosmoMC. BAO module now comes included in CosmoMC (Thanks to Antony Lewis).

1.0.1 (Released 05/29/13): Includes bug fixes.

1.0: Initial release, for August 2011 version of CosmoMC.

Contact

If you have any questions or comments about this readme, or you want further details on the changes made to CosmoMC or CAMB in ISiTGR, please feel free to email Jason Dossett at jason.dossett@brera.inaf.it