def get_array(file_name=None,labels=None): print("Reading from %s" %(file_name)) from iotbx import reflection_file_reader reflection_file = reflection_file_reader.any_reflection_file( file_name=file_name) array_to_use=None if labels: for array in reflection_file.as_miller_arrays(): if ",".join(array.info().labels)==labels: array_to_use=array break else: for array in reflection_file.as_miller_arrays(): if array.is_complex_array() or array.is_xray_amplitude_array() or\ array.is_xray_intensity_array(): array_to_use=array break if not array_to_use: text="" for array in reflection_file.as_miller_arrays(): text+=" %s " %(",".join(array.info().labels)) raise Sorry("Cannot identify array to use...possibilities: %s" %(text)) print("Using the array %s" %(",".join(array_to_use.info().labels))) return array_to_use
def get_array(file_name=None,labels=None): print "Reading from %s" %(file_name) from iotbx import reflection_file_reader reflection_file = reflection_file_reader.any_reflection_file( file_name=file_name) array_to_use=None if labels: for array in reflection_file.as_miller_arrays(): if ",".join(array.info().labels)==labels: array_to_use=array break else: for array in reflection_file.as_miller_arrays(): if array.is_complex_array() or array.is_xray_amplitude_array() or\ array.is_xray_intensity_array(): array_to_use=array break if not array_to_use: text="" for array in reflection_file.as_miller_arrays(): text+=" %s " %(",".join(array.info().labels)) raise Sorry("Cannot identify array to use...possibilities: %s" %(text)) print "Using the array %s" %(",".join(array_to_use.info().labels)) return array_to_use