filep.close()
xc = numpy.mean(xlist)
yc = numpy.mean(ylist)
zc = numpy.mean(zlist)
for i in range(len(xlist)):
xlist[i] = xlist[i] - xc
ylist[i] = ylist[i] - yc
zlist[i] = zlist[i] - zc
nanoparticles_init = []
botx, topx, boty, topy, botz, topz = \
nanoparticle.file_to_nanoparticle_list(nanopfile, nanoparticles_init)
print >> sys.stderr, "Read done"
nanoparticles = []
for selectedid in range(len(nanoparticles_init)):
nanop = nanoparticles_init[selectedid]
pcx, pcy, pcz = nanop.get_center()
# seleziona solo particelle interne
distx = min(math.fabs(pcx-botx), math.fabs(pcx-topx))
disty = min(math.fabs(pcy-boty), math.fabs(pcy-topy))
distz = min(math.fabs(pcz-botz), math.fabs(pcz-topz))
#print min(distx, disty, distz), 2.0 * common.R
import vtk
# no mi interessano le intersezioni
nanoparticle.POINTINSIDEDIM = 0
filename = "nanoparticle_final_config.txt"
if (len(sys.argv)) == 2:
filename = sys.argv[1]
nanoparticles = []
print ('Start reading '+filename+' ...')
botx, topx, boty, topy, botz, topz = \
nanoparticle.file_to_nanoparticle_list(filename, nanoparticles)
print ('Done')
if (botx >= topx) or (boty >= topy) or \
(boty >= topy):
print ('Error Invalid BOX\n')
exit()
# nanoparticle
minbox_x = 100000.0
maxbox_x =-100000.0
minbox_y = 100000.0
maxbox_y =-100000.0
minbox_z = 100000.0
maxbox_z =-100000.0
plt.savefig("curve_" + id + ".png", bbox_inches='tight')
# too big in case of many traps
#fullsetoftraps = []
if args.filmfile != "":
# non mi inter)essano le intersezioni
nanoparticle.POINTINSIDEDIM = 0
nanoparticles = []
print >> sys.stderr, "Reading nanoparticles file"
botx, topx, boty, topy, botz, topz = \
nanoparticle.file_to_nanoparticle_list(args.filmfile, \
nanoparticles)
if (botx >= topx) or (boty >= topy) or (boty >= topy):
print "Error Invalid BOX"
exit()
sources = []
counter = 0
for nanop in nanoparticles:
p1, p2, tetha = nanop.get_rotation_info()
cx, cy, cz = nanop.get_center()
counter = counter + 1
