Моделирование процессов в физике частиц (7 сем) / Model_M

Свободный пробег в кусочно однородной среде

xi − точки пересечениятраектории с границами ;μi между xi−1 и xi

τ= −ln(γ)

τ= τ − μ0 x0 ; τ = τ − μ1 (x1 − x0); ... τ = τ − μk (xk − xk −1)

τ< 0 ξ = xk + μτk

Свободный пробег в кусочно однородной среде

double ksi, tau, x[ ], mu[ ]; int k, n = ;

k = 0; tau = -ln(gamma()) - mu[0]*x[0]; while(k < n && tau >= 0){

k = k + 1; tau = tau – mu[k]*(x[k] – x[k-1]);

}

if(tau > 0) { } // out !!!

ksi = x[k] + tau / mu[k]; // результат

Пример

Установка из 2-х частей длинами вдоль траектории l1 , l2 3 процесса взаимодействия μ(11) , μ(12), μ(13); μ(21), μ(22), μ(23)

Моделировать точку взаимодействия xint и номер процесса.

-------------------------------------------------------------------------------

μ1 = μ(11) + μ(12) + μ(13)

μ2 = μ(21) + μ(22) + μ(23)

τ= −ln (γ1)

Макроскопическое сечение

μ макроскопическое сечение[обратная длина]

wk относительная массовая доля