Merge branch 'pr/aukeheerdink/75'

Author: douweschulte

rustyms/src/aminoacid/aminoacid.rs

@@ -373,6 +373,35 @@ impl AminoAcid {
         Self::Valine,
     ];
 
+    /// All amino acids (including I/L/J/B/Z but excluding X)
+    pub const ALL_AMINO_ACIDS: &'static [Self] = &[
+        Self::Alanine,
+        Self::AmbiguousAsparagine,
+        Self::AmbiguousGlutamine,
+        Self::AmbiguousLeucine,
+        Self::Arginine,
+        Self::Asparagine,
+        Self::AsparticAcid,
+        Self::Cysteine,
+        Self::GlutamicAcid,
+        Self::Glutamine,
+        Self::Glycine,
+        Self::Histidine,
+        Self::Isoleucine,
+        Self::Leucine,
+        Self::Lysine,
+        Self::Methionine,
+        Self::Phenylalanine,
+        Self::Proline,
+        Self::Pyrrolysine,
+        Self::Selenocysteine,
+        Self::Serine,
+        Self::Threonine,
+        Self::Tryptophan,
+        Self::Tyrosine,
+        Self::Valine,
+    ];
+
     // TODO: generalise over used storage type, so using molecularformula, monoisotopic mass, or average mass, also make sure that AAs can return these numbers in a const fashion
     #[expect(clippy::too_many_lines, clippy::too_many_arguments)]
     pub(crate) fn fragments(

rustyms/src/peptidoform/peptidoform.rs

@@ -1508,7 +1508,12 @@ impl<Complexity: AtMax<Linear>> Peptidoform<Complexity> {
     }
 
     /// Digest this sequence with the given protease and the given maximal number of missed cleavages.
-    pub fn digest(&self, protease: &Protease, max_missed_cleavages: usize) -> Vec<Self> {
+    pub fn digest(
+        &self,
+        protease: &Protease,
+        max_missed_cleavages: usize,
+        size_range: impl RangeBounds<usize>,
+    ) -> Vec<Self> {
         let mut sites = vec![0];
         sites.extend_from_slice(&protease.match_locations(&self.sequence));
         sites.push(self.len());
@@ -1517,7 +1522,9 @@ impl<Complexity: AtMax<Linear>> Peptidoform<Complexity> {
 
         for (index, start) in sites.iter().enumerate() {
             for end in sites.iter().skip(index + 1).take(max_missed_cleavages + 1) {
-                result.push(self.sub_peptide((*start)..*end));
+                if size_range.contains(&(end - start)) {
+                    result.push(self.sub_peptide((*start)..*end));
+                }
             }
         }
         result