Hi Muyuan,
Yes I do have the early output, I have truncated parts of the output indicated by "...." for brevity.
I would be happy to send the file to you, do you have a preferred email or otherwise to send the file?
Output:Writing in spt_01
Gather metadata from sets/pent.lst
exclude particles around sets/tail.lst
505 particles, scale by 1.00
excluding 4 out of 48 3d particles. Average distnace 18.4 px
excluding 11 out of 132 3d particles. Average distnace 18.4 px
excluding 7 out of 84 3d particles. Average distnace 18.4 px
excluding 11 out of 132 3d particles. Average distnace 18.4 px
....
....
excluding 15 out of 180 3d particles. Average distnace 18.3 px
excluding 8 out of 96 3d particles. Average distnace 18.3 px
excluding 9 out of 108 3d particles. Average distnace 18.3 px
total 505 out of 6060 particles excluded
Loading 3D particles
Loading 3D particles - 6060 particles from 57 tomograms
Using 3d alignment from the image header.
Using 'orig_class' in image header to assign class
5500/5555
Total 5555 3d particles and 337245 2d particles
Reading previous alignment
Loading particle info from spt_00/particle_info_3d.lst
3d alignment from spt_00/aliptcls3d_06.lst
2d alignment from spt_00/aliptcls2d_06.lst
Scale from old to new alignment 0.250
Sanity check on one tomogram
tomogram Legion_20250619_11_2, tilt id 10
4 particles from old refinement, 4 movement trajectories, 44 particles from new refinement
mean neighbor distance between particles from two refinement: 392.6113
mean subtilt translation 0.6137 from old refinement
mean subtilt translation 0.7361 from new refinement
Compute subtilt translation for new refinement
Legion_20250619_11_2 4 44
Legion_20250619_11_3 11 121
Legion_20250619_11 7 77
Legion_20250619_12_2 11 121
....
....
Legion_20250619_7 4 44
Legion_20250619_8_2 15 165
Legion_20250619_8 8 88
Legion_20250619_9 9 99
Traceback (most recent call last):
File "/home/doulin/miniconda3/envs/eman2/bin/e2spt_gathermeta.py", line 451, in <module>
main()
File "/home/doulin/miniconda3/envs/eman2/bin/e2spt_gathermeta.py", line 356, in main
dxf=[Transform({"type":"xyz","xtilt":x[0],"ytilt":x[1],"ztilt":x[2],"tx":v[0], "ty":v[1]}) for x,v in zip(rts,v1)]
~^^^
IndexError: invalid index to scalar variable.