Пример #1
0
        Index_1=-1.504042874,
        Index_2=1.891184873,
        E_break=e_break,
        beta=0.1,
        e0=200)

    smooth_soft = SmoothBrokenPowerLaw(
        Index_1=+0.2,
        Index_2=1.891184873,
        E_break=e_break,
        beta=0.1,
        e0=200)
    smooth_soft.set_prefactor(smooth_hard(e_break),e_break)

    plaw=PowerLaw(index=1.891184873, e0=200)
    plaw.set_prefactor(smooth_hard(e_break), e_break)

    if spectrum == 'PowerLaw':
        W44.model = plaw
    elif spectrum == 'SmoothBrokenPowerLawHard':
        W44.model = smooth_hard
    elif spectrum == 'SmoothBrokenPowerLawSoft':
        W44.model = smooth_soft


elif source == 'IC443':
    """
    SmoothBrokenPowerlaw Spectrum is from:
        /u/gl/funk/data3/ExtendedSources/NewAnalysis/gtlike/IC443/SmoothBrokenPowerlaw/IC443_fitted_BINNED_freeMore_BSP_lowEnergy.minos.xml

        <spectrum type="SmoothBrokenPowerLaw">