def calculateRowIntensity(row):
intensity = []
for column in range(detG.xdim):
xB, yB, zB = samG.generate_random_points_within_sample(
NG, slit_xG, cdf_xG, slit_yG, cdf_yG)
xD, yD, zD = detG.generate_random_points_within_pixel(NG, row, column)
beam_path = beam_path_within_sample(samG, xB, yB, zB, xD, yD, zD)
intensity.append(np.mean(np.exp(-beam_path / samG.mu)))
return row, intensity
def calculateIntensity(input):
i, j = input
intensity = []
for detPixel in range(i, i + j):
row, column = detPixel / det.xdim, detPixel % det.xdim
xB, yB, zB = sam.generate_random_points_within_sample(N, slit_x, cdf_x, slit_y, cdf_y)
xD, yD, zD = det.generate_random_points_within_pixel(N, row, column)
beam_path = beam_path_within_sample(sam, xB, yB, zB, xD, yD, zD)
intensity.append(np.mean(np.exp(-beam_path / sam.mu)))
return intensity