solar_modulation module

Module with class to simulate effect of solar modulation on antiproton flux; contains class for the force field approximation, and can be extended with additional classes with same attributes for different solar modulation models.

class pbarlike.solar_modulation.ForceFieldApprox(data, propagation_config, drn, verbose=False)

Bases: object

Contains methods to help calculate the solar modulated antiproton flux measured at the top of the atmosphere (TOA), which depends on solar modulation parameters.

Note

  • Always add new classes with same attributes for compatibility with rest of the code.

  • Solar modulation parameters are profiled over when using user-defined propagation parameters; and marginalized over when using the multinest sample.

Parameters:

  • propagation_config (object) – Object of class Propagation

  • drn (object) – Object of class DRNet

TOA_sim(phi_DM_LIS, errors='correlated')

Calculates primary and secondary antiproton flux after solar modulation according to the force field approximation.

nuisance_estimation(chi2)

Finds solar modulation parameters that minimize the given chi-squared function.

solar_mod(phi_LIS, V, Z=-1.0, A=1.0, m=0.9382720881604903)

Modulates the given CR flux depending on the given solar modulation potential, charge and mass of the incoming nuclei; defaults are set for antiproton nuclei.

Parameters:

  • phi_LIS (array) – Antiproton flux at the local interstellar (LIS) region.

  • V (float) – Solar modulation potential

  • Z (float) – Atomic number of CR nuclei

  • A (float) – Mass number of CR nuclei

  • m (float) – Mass of CR nuclei

Returns:

Solar modulated CR flux

Return type:

array

solar_modulation_margdparams(phi_LIS)

Calculates solar modulated flux for solar modulation parameters contained in the multinest sample.

solmod_c_profiledparams(phi_LIS)
solmod_uc_profiledparams(phi_LIS)