from __future__ import division
from cctbx import miller
from cctbx import crystal

def run():
  crystal_symmetry = crystal.symmetry((10,20,30, 90,90,90), "P212121")
  ms = miller.build_set(
    crystal_symmetry = crystal_symmetry,
    anomalous_flag   = True,
    d_min            = 2,
    d_max            = 10)
  crystal_symmetry.show_summary()  
  print "Resolution range:", ms.d_max_min()
  print "Number of reflections:", ms.indices().size()

if __name__ == '__main__':
  run()