Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

  • psana-cxif5315-r0169-cspad-ds2-NDArrDropletFinder.py.txt (download) - python script, which used *.cfg file to produce all ndarrays in the event store, gets these ndarrays, does preliminary processing, saves list of peaks in the output text file, draw events (if requested).
    This script can be run interactively or in batch by commands:

    Code Block
    // interactive command:
    psana -c psana-cxif5315-r0169-cspad-ds2-NDArrDropletFinder.cfg exp=cxif5315:run=169
    
    
    // batch job submission command:
    bsub -q psfehq -o log-r169.log python psana-cxif5315-r0169-cspad-ds2-NDArrDropletFinder.py

    At least 2 parameters need to be adjusted in this script for each mode:

    Code Block
        do_plot = True             - should be = False in the batch mode
        #events_max  = 10000000    - for batch mode loop over all events
        events_max  = 1000         - for loop over 1000 events in interactive mode

    This job produces text file with all found peaks with name like work/peaks-2015-05-09-23:52:30-cxif5315-r0169-all.txt, which can be used for further analysis

  • proc-peaks-from-file.py (download) - python script for processing of the text file with list of peaks.
    Input file name and output file path prefix should be set at the end of this script:

    Code Block
        procPeaksFromFiles (('work/peaks-2015-05-09-23:52:30-cxif5315-r0169-all.txt',), '2015-05-14-figs-cxif5315/plot-cxif5315-r0169-all')

    After that this script can be executed by the command:

    Code Block
    python proc-peaks-from-file.py

    which loops over peaks in file, does analysis, and plots requested histograms at the end.

E-mails

2015-05-04 from Kartik Ayyer:

Note
Regarding my code, I have made a git repo public here:
The code is in C. It is a complete mess and is not optimized as I have been tinkering with it. The main file is peakalign.c, and the functions you asked for are calc_phi(), calc_beta() and calc_merge().
My workflow has been Cheetah -> cxi file output -> Align -> Signal and background RZ merge -> Subtract background from signal. I have been posting my results on DESY confluence, which is unfortunately private, but Meng has access to it.

 

2015-05-04 from Meng:
Note
fraser1.m - matlab code for doing the reciprocal space mapping.

 

References

...