mass_dependent_ia_model
================================================

adds mass dependence to IA, with a linear equation

+-------------+-------------------------------------------------------------------+
| File        | intrinsic_alignments/mass_dependent_ia/mass_dependent_ia_model.py |
+-------------+-------------------------------------------------------------------+
| Attribution | Marika Asgari                                                     |
+-------------+-------------------------------------------------------------------+
|             | Benjamin Joachimi                                                 |
+-------------+-------------------------------------------------------------------+
|             | Konrad Kuijken                                                    |
+-------------+-------------------------------------------------------------------+
| URL         |                                                                   |
+-------------+-------------------------------------------------------------------+

Takes input IA Cls, Cl^ij_II and Cl^ij_GI. Then creates 
new Cl^ij_II = [f_r^(i) * (<M_*>^(i) / M_{*, pivot})^beta] * [f_r^(j) * (<M_*>^(j) / M_{*, pivot})^beta]  * Cl^ij_II
new Cl^ij_GI = [f_r^(j) * (<M_*>^(j) / M_{*, pivot})^beta] * Cl^ij_GI
new Cl^ij_IG = [f_r^(i) * (<M_*>^(i) / M_{*, pivot})^beta] * Cl^ij_IG


Assumptions
-----------

 - Input model only produces E-modes.
 - Currently only works for cosmic shear IA



Setup Parameters
----------------

.. list-table::
   :header-rows: 1

   * - Name
     - Type
     - Default
     - Description

   * - suffix
     - str
     - 
     - suffix used to distinguish between different data vectors
   * - new_suffix
     - str
     - 
     - If provided will save the resulting IA Cls using this suffix


Input values
----------------

.. list-table::
   :header-rows: 1

   * - Section
     - Name
     - Type
     - Default
     - Description

   * - intrinsic_alignment_parameters_{suffix}
     - beta
     - real
     - 
     - Mass dependence power for early-type galaxies
   * - 
     - log10_M_pivot
     - real
     - 
     - log10 of pivot halo mass (M_sun/h), can be derive from mean halo mass of early-type galaxies
   * - 
     - log10_M_mean_{i}
     - real
     - 
     - log10 of mean halo mass (M_sun/h) of early-type galaxies per redshif bin from catalogues
   * - 
     - f_r_{i}
     - real
     - 
     - Fraction of early-type galaxies per redshift bin, can be calculated using T_B from catalogues
   * - shear_cl_ii
     - nbin_a
     - int
     - 
     - Number of tomographic bins
   * - 
     - bin_{i}_{j}
     - real 1d
     - 
     - Intrinsic-intrinsic angular spectra for pairs of i,j values C_ell
   * - shear_cl_gi
     - bin_{i}_{j}
     - real 1d
     - 
     - Shear-intrinsic angular spectra for pairs of i,j values C_ell
   * - 
     - bin_{j}_{i}
     - real 1d
     - 
     - Intrinsic-shear angular spectra for pairs of i,j values C_ell


Output values
----------------


.. list-table:: Output values
   :header-rows: 1

   * - Section
     - Name
     - Type
     - Description

   * - shear_cl_ii
     - nbin_a
     - int
     - Number of tomographic bins
   * - 
     - bin_{i}_{j}
     - real 1d
     - Intrinsic-intrinsic angular spectra for pairs of i,j values C_ell
   * - shear_cl_gi
     - bin_{i}_{j}
     - real 1d
     - Shear-intrinsic angular spectra for pairs of i,j values C_ell
   * - 
     - bin_{j}_{i}
     - real 1d
     - Intrinsic-shear angular spectra for pairs of i,j values C_ell