calculate_gammas_secondary function

Gamma factors for fitting secondary models