Hi,
I'm having difficulty running phenix.real_space_correlation. I get the following error when comparing a pdb to a ccp4 map:
phenix.real_space_correlation ../17may31a_cryoSPARC_C1.pdb ../17may31a_cryoSPARC_C1.ccp4
-------------------------------------------------------------------------------
Input PDB file name: ../17may31a_cryoSPARC_C1.pdb
*************************************************
Number of scatterers: 95640
At special positions: 0
Unit cell: (592.359, 589.659, 590.506, 90, 90, 90)
Space group: P 1 (No. 1)
-------------------------------------------------------------------------------
Input map file name: ../17may31a_cryoSPARC_C1.ccp4
**************************************************
Map type: CCP4-format
Traceback (most recent call last):
File "/opt/X-tal_suite/phenix-1.12-2829/build/../modules/cctbx_project/mmtbx/command_line/real_space_correlation.py", line 10, in <module>
command_name = "phenix.real_space_correlation")
File "/opt/X-tal_suite/phenix-1.12-2829/modules/cctbx_project/mmtbx/real_space_correlation.py", line 335, in cmd_run
(fobs_handle.file_name, map_handle.file_name))
NameError: global name 'fobs_handle' is not defined
I get the following error when comparing a pdb to the reflection files.
phenix.real_space_correlation ../17may31a_cryoSPARC_C1.pdb ../17may31a_cryoSPARC_C1.mtz
-------------------------------------------------------------------------------
Input PDB file name: ../17may31a_cryoSPARC_C1.pdb
*************************************************
Number of scatterers: 95640
At special positions: 0
Unit cell: (592.359, 589.659, 590.506, 90, 90, 90)
Space group: P 1 (No. 1)
-------------------------------------------------------------------------------
Input reflection file name: ../17may31a_cryoSPARC_C1.mtz
********************************************************
Miller array info: ../17may31a_cryoSPARC_C1.mtz:F-obs,SIGF-obs
Observation type: xray.amplitude
Type of data: double, size=762857
Type of sigmas: double, size=762857
Number of Miller indices: 762857
Anomalous flag: False
Unit cell: (278.528, 278.528, 278.528, 90, 90, 90)
Space group: P 1 (No. 1)
Systematic absences: 0
Centric reflections: 0
Resolution range: 278.528 3.90017
Completeness in resolution range: 1
Completeness with d_max=infinity: 1
Wavelength: 1.0000
Traceback (most recent call last):
File "/opt/X-tal_suite/phenix-1.12-2829/build/../modules/cctbx_project/mmtbx/command_line/real_space_correlation.py", line 10, in <module>
command_name = "phenix.real_space_correlation")
File "/opt/X-tal_suite/phenix-1.12-2829/modules/cctbx_project/mmtbx/real_space_correlation.py", line 354, in cmd_run
r_free_flags = r_free_flags)
File "/opt/X-tal_suite/phenix-1.12-2829/modules/cctbx_project/mmtbx/utils/__init__.py", line 1572, in fmodel_simple
ro = outliers_rejection,om=optimize_mask)
File "/opt/X-tal_suite/phenix-1.12-2829/modules/cctbx_project/mmtbx/utils/__init__.py", line 1559, in get_fmodel
twin_law = tl)
File "/opt/X-tal_suite/phenix-1.12-2829/modules/cctbx_project/mmtbx/utils/__init__.py", line 1507, in fmodel_manager
n_resolution_bins_output = n_resolution_bins_output)
File "/opt/X-tal_suite/phenix-1.12-2829/modules/cctbx_project/mmtbx/f_model/f_model.py", line 439, in __init__
if(f_calc is None): f_calc = self.compute_f_calc(miller_array=f_obs)
File "/opt/X-tal_suite/phenix-1.12-2829/modules/cctbx_project/mmtbx/f_model/f_model.py", line 526, in compute_f_calc
exp_table_one_over_step_size = p.exp_table_one_over_step_size)
File "/opt/X-tal_suite/phenix-1.12-2829/modules/cctbx_project/cctbx/miller/__init__.py", line 1515, in structure_factors_from_scatterers
algorithm=algorithm)
File "/opt/X-tal_suite/phenix-1.12-2829/modules/cctbx_project/cctbx/xray/structure_factors/from_scatterers.py", line 53, in __call__
algorithm=algorithm) # passing algorithm allows f to decide on CPU/GPU implementation
File "/opt/X-tal_suite/phenix-1.12-2829/modules/cctbx_project/cctbx/xray/structure_factors/from_scatterers_fft.py", line 18, in __init__
manager, xray_structure, miller_set, algorithm="fft")
File "/opt/X-tal_suite/phenix-1.12-2829/modules/cctbx_project/cctbx/xray/structure_factors/manager.py", line 233, in __init__
assert xray_structure.unit_cell().is_similar_to(miller_set.unit_cell())
AssertionError
Any suggestions would be appreciated.
Best wishes,
Reza