text: 
  battaglia_defaults = {}
  battaglia_defaults['AGN'] = {
      'rho0_A0':4000.,
      'rho0_alpham':0.29,
      'rho0_alphaz':-0.66,
      'alpha_A0':0.88,
      'alpha_alpham':-0.03,
      'alpha_alphaz':0.19,
      'beta_A0':3.83,
      'beta_alpham':0.04,
      'beta_alphaz':-0.025
  }
  battaglia_defaults['SH'] = {
      'rho0_A0':19000.,
      'rho0_alpham':0.09,
      'rho0_alphaz':-0.95,
      'alpha_A0':0.70,
      'alpha_alpham':-0.017,
      'alpha_alphaz':0.27,
      'beta_A0':4.43,
      'beta_alpham':0.005,
      'beta_alphaz':0.037
  }
      
  
  battaglia_defaults['pres'] = {
      'P0_A0':18.1,
      'P0_alpham':0.154,
      'P0_alphaz':-0.758,
      'xc_A0':0.497,
      'xc_alpham':-0.00865,
      'xc_alphaz':0.731,
      'beta_A0':4.35,
      'beta_alpham':0.0393,
      'beta_alphaz':0.415
  }
      
  
  default_params = {
      
      # Mass function
      'st_A': 0.3222,
      'st_a': 0.707,
      'st_p': 0.3,
      'st_deltac': 1.686,
      'sigma2_kmin':1e-4,
      'sigma2_kmax':2000,
      'sigma2_numks':10000,
      'Wkr_taylor_switch':0.01,
  
      # Profiles
      'duffy_A_vir':7.85, # for Mvir
      'duffy_alpha_vir':-0.081,
      'duffy_beta_vir':-0.71,
      'duffy_A_mean':10.14, # for M200rhomeanz
      'duffy_alpha_mean':-0.081,
      'duffy_beta_mean':-1.01,
      'nfw_integral_numxs':40000, # not sufficient
      'nfw_integral_xmax':200,
      'electron_density_profile_integral_numxs':5000,
      'electron_density_profile_integral_xmax':20,
      'electron_pressure_profile_integral_numxs':5000,
      'electron_pressure_profile_integral_xmax':20,
      'battaglia_gas_gamma':-0.2,
      'battaglia_gas_family': 'AGN',
  
      'battaglia_pres_gamma' : -0.3,
      'battaglia_pres_alpha' : 1.,
      'battaglia_pres_family' : 'pres',
      # Power spectra
      'kstar_damping':0.01,
      'default_halofit':'mead',
      
      # Cosmology
      'omch2': 0.1198,
      'ombh2': 0.02225,
      'H0': 67.3,
      'ns': 0.9645,
      'As': 2.2e-9,
      'mnu': 0.0, # NOTE NO NEUTRINOS IN DEFAULT
      'w0': -1.0,
      'tau':0.06,
      'nnu':3.046,
      'wa': 0.,
      'num_massive_neutrinos':3,
      'T_CMB':2.7255e6,
      'parsec': 3.08567758e16,
      'mSun': 1.989e30,
      'thompson_SI': 6.6524e-29,
      'meterToMegaparsec': 3.241e-23,
      'Yp': 0.24,
  
      # HOD
      'hod_A_log10mthresh': 1.0, # This parameter is used to vary log10mthresh
      'hod_sig_log_mstellar': 0.2,
      'hod_alphasat': 1.0,
      'hod_Bsat':9.04,
      'hod_betasat':0.74,
      'hod_Bcut':1.65,
      'hod_betacut':0.59,
      'hod_bisection_search_min_log10mthresh': 7.,
      'hod_bisection_search_max_log10mthresh': 14.,
      'hod_bisection_search_rtol': 1e-4,
      'hod_bisection_search_warn_iter': 20,
  }
  