Double check the use of weights in polarization maps

[lucabaldini]

Dear Luca,

Thank you very much for adding the weighting application and the plot.

I’ve updated the ixpeobssim and tried the xpaddmodfweight app. It works well and generates new events files with a W_MODF column!

For a quick test, I used the following lines to generate Stokes cubes and polarization maps, as I did before. But I cannot see any difference between the obtained images and those without weighting. I am still checking and will try to figure out the problems. I understand that you are busy in recent days. In case you can quickly find a problem with my script, please kindly let me know. 

def make_pmap_cube():
        """Create the Stokes cubes.
        """
        file_list = get_sim_files()
        suffix='%s_%skev_%sarcs_pmap_cube' %(E_BINNING[0], E_BINNING[1], pixsize)
        pipeline.xpbin(*file_list, algorithm='PMAPCUBE', ebinalg='LIST', 
                       ebinning=E_BINNING,npix=npix, overwrite=True,
                       weights=True, weightcol='W_MODF',
                       suffix=pipeline.suffix(suffix))

‌

def plot_pmap_cube():
        """Open the PMAP files and run the standard plot.
        """
        file_list = get_pmap_cube_files()
        pmap_cube = xBinnedPolarizationMapCube.from_file_list(file_list)
        NUM_SIGMA=3.0
        pmap_cube.plot_polarization_angle(num_sigma=NUM_SIGMA)

Thanks,

Ping

Original url: https://bitbucket.org/ixpesw/ixpeobssim/issues/540
Created on 2022-01-27T09:10:13.927897+00:00