Page History
...
New peakfinder more than doubles the number of signal events in the distance distribution.
Equ region
Peak selection in equatorial region
...
| Code Block |
|---|
def peakIsSelectedEqu() :
"""Apply peak selection criteria to each peak from file
"""
if sp.sonc<9 : return False
if sp.atot<1800 : return False
if sp.r<100 : return False
if sp.r>450 : return False
return True |
Event selection
...
| Code Block |
|---|
def eventIsSelected() :
"""Apply selection criteria to entire event, based on list of peaks
"""
sp.event_is_selected = False
if sp.count_arc_pks_sel > 2 : return False
if sp.count_equ_pks_sel > 5 : return False
#if sp.count_equ_pks_sel < 0 : return False
# Require all peaks in the narrow range of dphi[deg]
for pk in sp.lst_equ_evt_peaks :
if math.fabs(pk.dphicmb) > 6 : return False
sp.event_is_selected = True
return True |
Plots
...
for raw peak data
Delta phi
delta phi for peaks around 0°, 180°, and combined distribution
Angles phi and beta
Evaluated for 2-peak events
Fitted for 2-peak events - these histograms are almost identical to evaluated
Fitted for >2-peak events
Fitted for ≥2-peak events
References
- class ImgAlgos.PyAlgos - collection of n-d array processing algorithms
- class Detector.PyDetector - direct access to detector data and calibration n-d arrays
Overview
Content Tools