def oneImage(self,framenumber):
self.reporters[framenumber] = []
from dxtbx.format.Registry import Registry
filename = self.phil_params.distl.image[framenumber]
reader = Registry.find(filename)
img = reader(filename)
detector = img.get_detector()
beam = img.get_beam()
S0 = beam.get_s0()
data = img.get_raw_data()
scan = img.get_scan()
print scan
if scan is None:
print "No scan"
RR = (0,1)
else:
print scan.get_oscillation()
RR = scan.get_oscillation_range()
from spotfinder.dxtbx_toolbox import Distl
sfall = Distl(params = self.phil_params, detector = detector, beam = beam, data = data)
resolutions = flex.double()
spotlist = []
from dials.model.data import ReflectionList,Reflection
reflections = ReflectionList()
for ip,panel in enumerate(detector):
for spot in sfall.finderlist[ip].spots:
resolutions.append( panel.get_resolution_at_pixel(S0, (spot.ctr_mass_x(), spot.ctr_mass_y())) )
spotlist.append(spot)
refl = Reflection()
refl.panel_number = ip
refl.centroid_position = (spot.ctr_mass_x(), spot.ctr_mass_y(),0.0)
refl.centroid_variance = (0.5,0.5,0.0)
reflections.append(refl)
selection = (resolutions>0.0)
if self.phil_params.distl.res.outer is not None:
selection = (selection and (resolutions>self.phil_params.distl.res.outer))
if self.phil_params.distl.res.inner is not None:
selection = (selection and (resolutions<self.phil_params.distl.res.inner))
reflections = reflections.select(selection)
return dict(detector=detector, beam=beam, reflections=reflections, scan = scan,
gonio = img.get_goniometer())
def oneImage(self, framenumber):
self.reporters[framenumber] = []
import dxtbx.format.Registry
filename = self.phil_params.distl.image[framenumber]
reader = dxtbx.format.Registry.get_format_class_for_file(filename)
img = reader(filename)
detector = img.get_detector()
beam = img.get_beam()
S0 = beam.get_s0()
data = img.get_raw_data()
scan = img.get_scan()
print(scan)
if scan is None:
print("No scan")
RR = (0, 1)
else:
print(scan.get_oscillation())
RR = scan.get_oscillation_range()
from spotfinder.dxtbx_toolbox import Distl
sfall = Distl(params=self.phil_params,
detector=detector,
beam=beam,
data=data)
resolutions = flex.double()
spotlist = []
from dials.model.data import ReflectionList, Reflection
reflections = ReflectionList()
for ip, panel in enumerate(detector):
for spot in sfall.finderlist[ip].spots:
resolutions.append(
panel.get_resolution_at_pixel(
S0, (spot.ctr_mass_x(), spot.ctr_mass_y())))
spotlist.append(spot)
refl = Reflection()
refl.panel_number = ip
refl.centroid_position = (spot.ctr_mass_x(), spot.ctr_mass_y(),
0.0)
refl.centroid_variance = (0.5, 0.5, 0.0)
reflections.append(refl)
selection = (resolutions > 0.0)
if self.phil_params.distl.res.outer is not None:
selection = (selection
and (resolutions > self.phil_params.distl.res.outer))
if self.phil_params.distl.res.inner is not None:
selection = (selection
and (resolutions < self.phil_params.distl.res.inner))
reflections = reflections.select(selection)
return dict(detector=detector,
beam=beam,
reflections=reflections,
scan=scan,
gonio=img.get_goniometer())
