def find_roi(mca, left, right, energy=False):
"""
This procedure finds the index number of the ROI with a specified
left and right channel number.
Parameters:
-----------
* mca: mca object
* left: Left channel number (or energy) of this ROI
* right: Right channel number (or energy) of this ROI
* energy: Set this flag to True to indicate that Left and Right are
in units of energy rather than channel number.
Output:
-------
* Returns the index of the specified ROI,
-1 if the ROI was not found.
Example:
--------
>>index = find_roi(mca, 100, 200)
"""
if energy:
kws = dict(offset=mca.offset, slope=mca.slope,
quad=mca.quad, clip=len(mca.data))
left = energy_to_channel(left, **kws)
right = energy_to_channel(right, **kws)
index = 0
for roi in self.rois:
if left == roi.left and right == roi.right:
return index
index = index + 1
return -1
def med_copy_rois(med, source_mca=0):
"""
This procedure copies the ROIs defined for one Mca in the Med to all of
the other Mcas.
Parameters:
-----------
* source_mca: The index number of the Mca from which the ROIs are to
be copied. This number ranges from 0 to self.n_detectors-1.
The default is the first Mca (index=0).
Notes:
------
The ROIs are copied by their position
in energy rather than in channels. This is very useful when
copying ROIs when the calibration parameters for each Mca in
the Med are not identical.
"""
units = "channel"
if energy:
units = "keV"
kws = dict(offset = med.mca[source_mca].offset,
slope = med.mca[source_mca].slope,
quad = med.mca[source_mca].quad,
clip = med.mca[source_mca].clip)
left, right, label, bgr_width = [], [], [], []
for roi in med.mca[source_mca].rois:
left.append(channel_to_energy(roi.left, **kws))
right.append(channel_to_energy(roi.right, **kws))
label.append(roi.label)
bgr_width.append(roi.bgr)
for j in range(med.n_detectors):
off = med.mca[j].offset
slope = med.mca[j].slope
quad = med.mca[j].quad
clip = len(med.mca[j].data)
med.mca[j].rois = []
for k in range(len(left)):
l = energy_to_channel(left[k],offset=off,slope=slope,quad=quad,clip=clip)
r = energy_to_channel(right[k],offset=off,slope=slope,quad=quad,clip=clip)
roi = ROI(left=int(l),right=int(r), label=label[k], bgr_width=bgr_width[k])
