Commit 92e1b800 authored by Giovanni La Mura's avatar Giovanni La Mura
Browse files

Remove radius weighing from geometric center calculation

parent 7e2b93d4

src/scripts/model_maker.py

+3 −3
Original line number Diff line number Diff line
@@ -679,9 +679,9 @@ def print_model_summary(scatterer, geometry): 
    for i in range(scatterer['nsph']):

        sph_type_index = scatterer['vec_types'][i] - 1

        ros = scatterer['ros'][sph_type_index]

        avgX += (ros * geometry['vec_sph_x'][i])

        avgY += (ros * geometry['vec_sph_y'][i])

        avgZ += (ros * geometry['vec_sph_z'][i])

        avgX += geometry['vec_sph_x'][i]

        avgY += geometry['vec_sph_y'][i]

        avgZ += geometry['vec_sph_z'][i]

        R3tot += math.pow(ros, 3.0)

        if (ros > Rmax):

            Rmax = ros