PICOcode.AnalysisModules package
Subpackages
- PICOcode.AnalysisModules.dataset_modules package
- Submodules
- PICOcode.AnalysisModules.dataset_modules.dataset_processor_base module
- PICOcode.AnalysisModules.dataset_modules.handscanning_data module
- PICOcode.AnalysisModules.dataset_modules.level2_acoustic_analysis module
- PICOcode.AnalysisModules.dataset_modules.level2_dytran_analysis module
- Module contents
- PICOcode.AnalysisModules.event_modules package
- Subpackages
- Submodules
- PICOcode.AnalysisModules.event_modules.dytran_analysis module
DytranAnalysisBaseDytranAnalysisBase.HeaviQuad()DytranAnalysisBase.QuartFit()DytranAnalysisBase.blockDownsample()DytranAnalysisBase.blockSmooth()DytranAnalysisBase.find_t0()DytranAnalysisBase.fit_dytran()DytranAnalysisBase.fit_dytran_heaviQuad()DytranAnalysisBase.fit_dytran_new()DytranAnalysisBase.get_dytran_list()DytranAnalysisBase.run()DytranAnalysisBase.sum_dytran()
- PICOcode.AnalysisModules.event_modules.fastdaq_analysis module
- PICOcode.AnalysisModules.event_modules.get_history module
- PICOcode.AnalysisModules.event_modules.processor_base module
- Module contents
- PICOcode.AnalysisModules.run_modules package
Submodules
PICOcode.AnalysisModules.pico40l_cuts module
- apply_common_cuts(run_type_fractions=None, cut_dict={'AP:alpha:gt': 18.17414536944306, 'AP:bulk:gt': 0, 'AP:neutron:lt': 18.17414536944306, 'R2:bulk:ge': 0, 'R2:bulk:le': 19321, 'R2:wall:ge': 19321, 'Z:bulk:ge': -400, 'Z:bulk:le': 0, 'Z:wall:ge': -400, 'Z:wall:le': 0, 'abs(PT4-pset):bulk:le': 1, 'abs(piezo_t0_v1_global-dytran_t0):bulk:lt': 0.02, 'dytranC:single:lt': 1.2, 'dytran_aprime:multiples:ge': [0.0727, 0.078], 'dytran_aprime:singles:ge': [0.02, 0.02], 'dytran_aprime:singles:le': [0.0727, 0.078], 'dytran_combo_sum_param:bulk:gt': 0, 'dytran_cprime:bulk:ge': [-0.0013, 0.0018], 'dytran_cprime:singles:le': [-0.00077, 0.0027], 'dytran_cprime:wall:le': [-0.0013, 0.0018], 'piezo_t0_v1_global:bulk:gt': -999, 'pset:bulk:eq': 30, 'run_type:alpha:eq': [99], 'run_type:neutron:eq': [91, 95], 'te:bulk:ge': 5, 'trigger_camera:bulk:ge': 4, 'trigger_mask:bulk:eq': 12, 'trigger_plc:pressure:eq': 16, 'trigger_route:bulk:eq': 'trigger_camera'}, recon_path='master', recon=None, include_multis=False, random=False, other_cuts=None, verbose=False, debug=False)
Cut on a recon file, and return a numpy cut for the events which pass the cuts.
By default, this applies cuts to return bulk events.
- Parameters:
- run_type_fractions: dict(int: int) or None, default=None
None, or a dictionary with keys that are run types, and values that are either: - between 0 and 1, inclusive: fraction of all events passing cuts
which are returned in the cut array.
Greater than 1: return up to the value’s number of events which pass cuts.
If random is False, then the first N values are returned. If run_type_fractions is None, make no cuts on run types.
- cut_dict: dict(str: (int, float, list)),
default=cuts (defined in file)
Cuts to apply. Must have the same keys as the cuts dictionary above.
- recon_path: str
Path to the recon file.
- recon: None or PICOcode.ReconFile
Recon file on which to apply cuts. If None, read the file at recon_path and apply cuts to it. Useful if a recon file has been modified prior to applying cuts.
- include_multis: bool, default=False
Whether to include multi-bubble events in the returned cut array.
- random: bool, default=False
Whether to randomize the results. If True and the run_type_fraction of each run_type is less than 1 (or less than the number of events passing cuts), then the returned array has a random sample of events passing cuts.
- other_cuts: array of bool or None, default=None
Other cuts to apply to the recon data. Useful when using random=True, otherwise identical to using get_bulk_events(…) & other_cuts. If None, then does not apply additional cuts.
- verbose: bool, default=False
Print the number of events passing cuts before returning.
- debug: bool, default=False
Print the number of events passing after each cut. Note the numbers printed are with all previous cuts applied as well, so each step will have the same or fewer events.
- Returns:
- np.ndarray
An array which contains True on events passing cuts (up to the fraction or number of events specified in run_type_fractions).
Notes
The returned array may be used in conjunction with other cuts, e.g. >>> cut = get_bulk_events() & (recon.AP < cuts[“AP_neutron_lt”])
If include_multis is True, the returned array only includes the first entry in the recon file for that event, since the cut on ibub==1 is still made. This is intentional to avoid double-counting of multi-bubble events.
- get_bulk_events(run_type_fractions=None, cut_dict={'AP:alpha:gt': 18.17414536944306, 'AP:bulk:gt': 0, 'AP:neutron:lt': 18.17414536944306, 'R2:bulk:ge': 0, 'R2:bulk:le': 19321, 'R2:wall:ge': 19321, 'Z:bulk:ge': -400, 'Z:bulk:le': 0, 'Z:wall:ge': -400, 'Z:wall:le': 0, 'abs(PT4-pset):bulk:le': 1, 'abs(piezo_t0_v1_global-dytran_t0):bulk:lt': 0.02, 'dytranC:single:lt': 1.2, 'dytran_aprime:multiples:ge': [0.0727, 0.078], 'dytran_aprime:singles:ge': [0.02, 0.02], 'dytran_aprime:singles:le': [0.0727, 0.078], 'dytran_combo_sum_param:bulk:gt': 0, 'dytran_cprime:bulk:ge': [-0.0013, 0.0018], 'dytran_cprime:singles:le': [-0.00077, 0.0027], 'dytran_cprime:wall:le': [-0.0013, 0.0018], 'piezo_t0_v1_global:bulk:gt': -999, 'pset:bulk:eq': 30, 'run_type:alpha:eq': [99], 'run_type:neutron:eq': [91, 95], 'te:bulk:ge': 5, 'trigger_camera:bulk:ge': 4, 'trigger_mask:bulk:eq': 12, 'trigger_plc:pressure:eq': 16, 'trigger_route:bulk:eq': 'trigger_camera'}, recon_path='master', recon=None, include_multis=False, random=False, other_cuts=None, verbose=False, debug=False)
- get_multiples_events(run_type_fractions=None, cut_dict={'AP:alpha:gt': 18.17414536944306, 'AP:bulk:gt': 0, 'AP:neutron:lt': 18.17414536944306, 'R2:bulk:ge': 0, 'R2:bulk:le': 19321, 'R2:wall:ge': 19321, 'Z:bulk:ge': -400, 'Z:bulk:le': 0, 'Z:wall:ge': -400, 'Z:wall:le': 0, 'abs(PT4-pset):bulk:le': 1, 'abs(piezo_t0_v1_global-dytran_t0):bulk:lt': 0.02, 'dytranC:single:lt': 1.2, 'dytran_aprime:multiples:ge': [0.0727, 0.078], 'dytran_aprime:singles:ge': [0.02, 0.02], 'dytran_aprime:singles:le': [0.0727, 0.078], 'dytran_combo_sum_param:bulk:gt': 0, 'dytran_cprime:bulk:ge': [-0.0013, 0.0018], 'dytran_cprime:singles:le': [-0.00077, 0.0027], 'dytran_cprime:wall:le': [-0.0013, 0.0018], 'piezo_t0_v1_global:bulk:gt': -999, 'pset:bulk:eq': 30, 'run_type:alpha:eq': [99], 'run_type:neutron:eq': [91, 95], 'te:bulk:ge': 5, 'trigger_camera:bulk:ge': 4, 'trigger_mask:bulk:eq': 12, 'trigger_plc:pressure:eq': 16, 'trigger_route:bulk:eq': 'trigger_camera'}, recon_path='master', recon=None, random=False, other_cuts=None, verbose=False, debug=False)
- get_singles_events(run_type_fractions=None, cut_dict={'AP:alpha:gt': 18.17414536944306, 'AP:bulk:gt': 0, 'AP:neutron:lt': 18.17414536944306, 'R2:bulk:ge': 0, 'R2:bulk:le': 19321, 'R2:wall:ge': 19321, 'Z:bulk:ge': -400, 'Z:bulk:le': 0, 'Z:wall:ge': -400, 'Z:wall:le': 0, 'abs(PT4-pset):bulk:le': 1, 'abs(piezo_t0_v1_global-dytran_t0):bulk:lt': 0.02, 'dytranC:single:lt': 1.2, 'dytran_aprime:multiples:ge': [0.0727, 0.078], 'dytran_aprime:singles:ge': [0.02, 0.02], 'dytran_aprime:singles:le': [0.0727, 0.078], 'dytran_combo_sum_param:bulk:gt': 0, 'dytran_cprime:bulk:ge': [-0.0013, 0.0018], 'dytran_cprime:singles:le': [-0.00077, 0.0027], 'dytran_cprime:wall:le': [-0.0013, 0.0018], 'piezo_t0_v1_global:bulk:gt': -999, 'pset:bulk:eq': 30, 'run_type:alpha:eq': [99], 'run_type:neutron:eq': [91, 95], 'te:bulk:ge': 5, 'trigger_camera:bulk:ge': 4, 'trigger_mask:bulk:eq': 12, 'trigger_plc:pressure:eq': 16, 'trigger_route:bulk:eq': 'trigger_camera'}, recon_path='master', recon=None, random=False, other_cuts=None, verbose=False, debug=False)
- get_wall_events(run_type_fractions=None, cut_dict={'AP:alpha:gt': 18.17414536944306, 'AP:bulk:gt': 0, 'AP:neutron:lt': 18.17414536944306, 'R2:bulk:ge': 0, 'R2:bulk:le': 19321, 'R2:wall:ge': 19321, 'Z:bulk:ge': -400, 'Z:bulk:le': 0, 'Z:wall:ge': -400, 'Z:wall:le': 0, 'abs(PT4-pset):bulk:le': 1, 'abs(piezo_t0_v1_global-dytran_t0):bulk:lt': 0.02, 'dytranC:single:lt': 1.2, 'dytran_aprime:multiples:ge': [0.0727, 0.078], 'dytran_aprime:singles:ge': [0.02, 0.02], 'dytran_aprime:singles:le': [0.0727, 0.078], 'dytran_combo_sum_param:bulk:gt': 0, 'dytran_cprime:bulk:ge': [-0.0013, 0.0018], 'dytran_cprime:singles:le': [-0.00077, 0.0027], 'dytran_cprime:wall:le': [-0.0013, 0.0018], 'piezo_t0_v1_global:bulk:gt': -999, 'pset:bulk:eq': 30, 'run_type:alpha:eq': [99], 'run_type:neutron:eq': [91, 95], 'te:bulk:ge': 5, 'trigger_camera:bulk:ge': 4, 'trigger_mask:bulk:eq': 12, 'trigger_plc:pressure:eq': 16, 'trigger_route:bulk:eq': 'trigger_camera'}, recon_path='master', recon=None, include_multis=False, random=False, other_cuts=None, verbose=False, debug=False)