Bug #5036
closedlosing particles when inverting the particle stack
0%
Description
I have found problems while using myami/trunk the number of particles differ between stacks when changing from invert to no invert. The number of particles decrease by approximately 2%. I have tried this multiple times on different data sets. When switching back to myami/3.2 this problem goes away. When I run the very first stack using myami/3.2 instead of myami/trunk I gain 3,000 particles.
Here are example of two command lines
/opt/applications/myami/trunk/bin/makestack2.py \
--single=start.hed --fromstackid=7 --no-invert --normalize-method=edgenorm \
--boxsize=384 --bin=2 --filetype=imagic --forceInsert --description=stack12 \
--runname=stack12 --rundir=/gpfs/group/em/appion/kubalek/17may30d/stacks/stack12 \
--preset=ed-a-DW --commit --projectid=503 --session=17may30d \
--no-rejects --no-wait --continue --expid=19529 --jobtype=makestack2
[Tue Jun 13 15:39:30 2017]
finished run of makestack2
The number of particles in the above stack =61,018
running the above stack with only the density inverted
[ kubalek@emb0633: Wed Jun 21 14:13:13 2017 ]
/opt/applications/myami/trunk/bin/makestack2.py \
--single=start.hed --fromstackid=13 --invert --normalize-method=edgenorm \
--boxsize=384 --bin=2 --filetype=imagic --forceInsert --description='test redo of stack12' \
--runname=stack21 --rundir=/gpfs/group/em/appion/kubalek/17may30d/stacks/stack21 \
--preset=ed-a-DW --commit --projectid=503 --session=17may30d \
--no-rejects --no-wait --continue --expid=19529 --jobtype=makestack2
the number in this stack now 58,746
Updated by Anchi Cheng over 7 years ago
- Subject changed from myami/trunk to losing particles when inverting the particle stack
- Status changed from New to Assigned
- Assignee changed from Sargis Dallakyan to Neil Voss
- Affected Version changed from Appion/Leginon 3.2 to Appion/Leginon 3.3
Updated by Neil Voss over 7 years ago
Hi Liz, I do not have access to your data, but I notice in stack12 you have "--fromstackid=7", but in stack21 you have "--fromstackid=13". Does stack12 have id #13? If you want the same data, I would keep the fromstackid the same.
Updated by Neil Voss over 7 years ago
- Assignee changed from Neil Voss to Elizabeth Wilson-Kubalek
Updated by Anchi Cheng about 6 years ago
- Status changed from Assigned to Won't Fix or Won't Do
user error