 p, chipq (-q, ..., q); bcdev=2500 m
 1, 10.026513098524001+0.0j, 0.0+0.0j, -10.026513098524001+0.0j
 2, 0.0+0.0j, 0.0+0.0j, -28.359261614488254+0.0j, 0.0+0.0j, 0.0+0.0j
