Versions Compared

Key

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

...

Code Block
titleList of events from Meng
collapsetrue
Liang, Mengning   
Tue 2/28/2017 4:09 PM
To:   Dubrovin, Mikhail;  

148 – double hit in equator – 3 peaks
162 – no peaks
525 – double hit in equator - 6 peaks
1456 – should find multiple arc peaks
1521 – double hit equator – 4 total peaks
2551 – no peaks
I could not find parameters to find the right peaks in these events – can see by eye where the peaks are. 
Thanks,
Meng

 

2017-03-02 Peak selection parameters

EQU

Code Block
titlePeak reconstruction/selection parameters
peaks_equ = alg_equ.peak_finder_v4r2(nda, thr_low=40, thr_high=60, rank=4, r0=6, dr=2)
alg_equ.set_peak_selection_pars(npix_min=3, npix_max=1000, amax_thr=0, atot_thr=500, son_min=4)

Image Removed

Note

For further histograms discard peaks with pk.r<105

 

Image RemovedImage Removed

Note

For further histograms discard peaks with pk.atot<600

Image RemovedImage RemovedImage Removed

Image RemovedImage Removed

Image RemovedImage RemovedImage Removed

Image Removed

Image RemovedImage RemovedImage Removed

Image RemovedImage RemovedImage Removed

Tentative peak selection criteria

Code Block
        if pk.r<100     : return False
        if pk.r>454     : return False
        if pk.atot<600  : return False
        if pk.rms>25    : return False
        if pk.bkgd>15   : return False
        if pk.bkgd<-15  : return False

pk.son ?
sigma_hov ?

ARC

Code Block
alg_arc.set_peak_selection_pars(npix_min=3, npix_max=1000, amax_thr=0, atot_thr=500, son_min=4)
peaks_arc = alg_arc.peak_finder_v4r2(nda, thr_low=40, thr_high=60, rank=8, r0=10, dr=2)

Image Removed

Note

Further discard peaks with pk.atot<1000

Image RemovedImage RemovedImage Removed

Image RemovedImage Removed

Image RemovedImage RemovedImage Removed

Image RemovedImage RemovedImage Removed

Image RemovedImage Removed

Image RemovedImage Removed

Image RemovedImage Removed

 

Tentative peak selection criteria

 

Code Block
titleTentative peak selection criteria
        if pk.atot<1000 : return False
        if pk.r<435     : return False
        if pk.r>443     : return False
        if pk.rms>25    : return False
        if pk.bkgd>15   : return False
        if pk.bkgd<-15  : return False
pk.son ?
sigma_hov ?