|
template<typename SpectrumT , typename TransitionT > |
void | pickTransitionGroup (MRMTransitionGroup< SpectrumT, TransitionT > &transition_group) |
| Pick a group of chromatograms belonging to the same peptide. More...
|
|
template<typename SpectrumT , typename TransitionT > |
MRMFeature | createMRMFeature (MRMTransitionGroup< SpectrumT, TransitionT > &transition_group, std::vector< SpectrumT > &picked_chroms, int &chr_idx, int &peak_idx) |
| Create feature from a vector of chromatograms and a specified peak. More...
|
|
template<typename SpectrumT > |
void | remove_overlapping_features (std::vector< SpectrumT > &picked_chroms, double best_left, double best_right) |
| Remove overlapping features. More...
|
|
void | findLargestPeak (std::vector< RichPeakChromatogram > &picked_chroms, int &chr_idx, int &peak_idx) |
| Find largest peak in a vector of chromatograms. More...
|
|
|
| MRMTransitionGroupPicker () |
| Constructor. More...
|
|
| ~MRMTransitionGroupPicker () |
| Destructor. More...
|
|
Public Member Functions inherited from DefaultParamHandler |
| DefaultParamHandler (const String &name) |
| Constructor with name that is displayed in error messages. More...
|
|
| DefaultParamHandler (const DefaultParamHandler &rhs) |
| Copy constructor. More...
|
|
virtual | ~DefaultParamHandler () |
| Destructor. More...
|
|
virtual DefaultParamHandler & | operator= (const DefaultParamHandler &rhs) |
| Assignment operator. More...
|
|
virtual bool | operator== (const DefaultParamHandler &rhs) const |
| Equality operator. More...
|
|
void | setParameters (const Param ¶m) |
| Sets the parameters. More...
|
|
const Param & | getParameters () const |
| Non-mutable access to the parameters. More...
|
|
const Param & | getDefaults () const |
| Non-mutable access to the default parameters. More...
|
|
const String & | getName () const |
| Non-mutable access to the name. More...
|
|
void | setName (const String &name) |
| Mutable access to the name. More...
|
|
const std::vector< String > & | getSubsections () const |
| Non-mutable access to the registered subsections. More...
|
|
|
void | updateMembers_ () |
| Synchronize members with param class. More...
|
|
MRMTransitionGroupPicker & | operator= (const MRMTransitionGroupPicker &rhs) |
| Assignment operator is protected for algorithm. More...
|
|
template<typename SpectrumT , typename TransitionT > |
double | computeQuality_ (MRMTransitionGroup< SpectrumT, TransitionT > &transition_group, std::vector< SpectrumT > &picked_chroms, const int chr_idx, const double best_left, const double best_right, String &outlier) |
| Compute transition group quality (higher score is better) More...
|
|
template<typename SpectrumT > |
void | recalculatePeakBorders_ (std::vector< SpectrumT > &picked_chroms, double &best_left, double &best_right, double max_z) |
| Recalculate the borders of the peak. More...
|
|
double | calculateBgEstimation_ (const RichPeakChromatogram &chromatogram, double best_left, double best_right) |
| Will use the chromatogram to estimate the background noise and then subtract it. More...
|
|
|
template<typename SpectrumT > |
void | prepareMasterContainer_ (const SpectrumT &ref_chromatogram, SpectrumT &master_peak_container, double best_left, double best_right) |
| create an empty master peak container that has the correct mz / RT values set More...
|
|
template<typename SpectrumT > |
SpectrumT | resampleChromatogram_ (const SpectrumT &chromatogram, SpectrumT &master_peak_container, double best_left, double best_right) |
| use the master container from above to resample a chromatogram at those points stored in the master container More...
|
|
Protected Member Functions inherited from DefaultParamHandler |
void | defaultsToParam_ () |
| Updates the parameters after the defaults have been set in the constructor. More...
|
|
The MRMTransitionGroupPicker finds peaks in chromatograms that belong to the same precursors.
It is called through pickTransitionGroup which will accept an MRMTransitionGroup filled with n chromatograms and perform the following steps:
- Step 1: find features (peaks) in individual chromatograms
- Step 2: merge these features to consensus features that span multiple chromatograms
Step 1 is performed by smoothing the individual chromatogram and applying the PeakPickerHiRes.
Step 2 is performed by finding the largest peak overall and use this to create a feature, propagating this through all chromatograms.
double computeQuality_ |
( |
MRMTransitionGroup< SpectrumT, TransitionT > & |
transition_group, |
|
|
std::vector< SpectrumT > & |
picked_chroms, |
|
|
const int |
chr_idx, |
|
|
const double |
best_left, |
|
|
const double |
best_right, |
|
|
String & |
outlier |
|
) |
| |
|
inlineprotected |
Compute transition group quality (higher score is better)
This is only based on the co-elution of the chromatograms and internal consistency without any library information.
For the final score (larger is better), consider these scores:
- missing_peaks (the more peaks are missing, the worse)
- multiple_peaks
- mean of the shapes (1 is very good, 0 is bad)
- mean of the coelutions (0 is good, 1 is ok, above 1 is pretty bad)
These scores are similar to the ones computed by MRMFeatureFinderScoring and a simple sum of these scores is returned.
left_borders / right_borders might not have the same length since we might have peaks missing!!
References MRMTransitionGroup< SpectrumType, TransitionType >::getChromatograms(), MRMTransitionGroup< SpectrumType, TransitionType >::getTransitions(), OpenMS::Constants::k, LOG_DEBUG, mean_and_stddev::mean(), OpenSwath::Scoring::normalizedCrossCorrelation(), and OpenSwath::Scoring::xcorrArrayGetMaxPeak().
void pickTransitionGroup |
( |
MRMTransitionGroup< SpectrumT, TransitionT > & |
transition_group | ) |
|
|
inline |
Pick a group of chromatograms belonging to the same peptide.
Will identify peaks in a set of chromatograms that belong to the same peptide. The chromatograms are given in the MRMTransitionGroup container which also contains the mapping of the chromatograms to their metadata.
The resulting features are added added to the MRMTransitionGroup. Each feature contains the following meta-data:
- PeptideRef
- leftWidth
- rightWidth
- total_xic (fragment trace XIC sum)
- peak_apices_sum
References MRMTransitionGroup< SpectrumType, TransitionType >::addFeature(), MRMTransitionGroup< SpectrumType, TransitionType >::getChromatograms(), Peak2D::getIntensity(), MetaInfoInterface::getMetaValue(), MSSpectrum< PeakT >::isSorted(), OpenMS::Constants::k, PeakPickerMRM::pickChromatogram(), DefaultParamHandler::setParameters(), MSSpectrum< PeakT >::sortByIntensity(), and MSSpectrum< PeakT >::sortByPosition().
Referenced by OpenSwathWorkflow::scoreAllChromatograms().
void recalculatePeakBorders_ |
( |
std::vector< SpectrumT > & |
picked_chroms, |
|
|
double & |
best_left, |
|
|
double & |
best_right, |
|
|
double |
max_z |
|
) |
| |
|
inlineprotected |
Recalculate the borders of the peak.
By collecting all left and right borders of contained peaks, a consensus peak is computed. By looking at the means and standard deviations of all the peak borders it is estimated whether the proposed peak border deviates too much from the consensus one. If the deviation is too high (in this case), then we fall back to the "consensus" (a median here).
References OpenMS::Constants::k, LOG_DEBUG, and OpenMS::Math::mean().
void remove_overlapping_features |
( |
std::vector< SpectrumT > & |
picked_chroms, |
|
|
double |
best_left, |
|
|
double |
best_right |
|
) |
| |
|
inline |
Remove overlapping features.
Remove features that are within the current seed (between best_left and best_right) or overlap with it. An overlapping feature is defined as a feature that has either of its borders within the border of the current peak
Directly adjacent features are allowed, e.g. they can share one border.
References OpenMS::Constants::k.