protected virtual void InitializeProcessingTasks()
{
MSGenerator = MSGeneratorFactory.CreateMSGenerator(Run.MSFileType);
PeakDetector = PeakDetectorFactory.CreatePeakDetector(NewDeconToolsParameters);
Deconvolutor = DeconvolutorFactory.CreateDeconvolutor(NewDeconToolsParameters);
//the new iThrash imports the _peaks.txt file
if (Deconvolutor is DeconTools.Backend.ProcessingTasks.Deconvoluters.InformedThrashDeconvolutor)
{
_deconvolutorRequiresPeaksFile = true;
}
//Will initialize these but whether or not they are used are determined elsewhere
ZeroFiller = new DeconToolsZeroFiller(NewDeconToolsParameters.MiscMSProcessingParameters.ZeroFillingNumZerosToFill);
Smoother = new SavitzkyGolaySmoother(NewDeconToolsParameters.MiscMSProcessingParameters.SavitzkyGolayNumPointsInSmooth,
NewDeconToolsParameters.MiscMSProcessingParameters.SavitzkyGolayOrder);
FitScoreCalculator = new DeconToolsFitScoreCalculator();
ScanResultUpdater = new ScanResultUpdater(NewDeconToolsParameters.ScanBasedWorkflowParameters.ProcessMS2);
ResultValidator = new ResultValidatorTask();
IsosResultExporter = IsosExporterFactory.CreateIsosExporter(Run.ResultCollection.ResultType, ExporterType,
IsosOutputFileName);
ScanResultExporter = ScansExporterFactory.CreateScansExporter(Run.MSFileType, ExporterType,
ScansOutputFileName);
if (!_deconvolutorRequiresPeaksFile)
{
PeakListExporter = PeakListExporterFactory.Create(ExporterType, Run.MSFileType, PeakListExporterTriggerValue,
PeakListOutputFileName);
}
PeakToMSFeatureAssociator = new PeakToMSFeatureAssociator();
}
public void Test1()
{
var peakAssociator = new PeakToMSFeatureAssociator();
Run run = new XCaliburRun2(FileRefs.RawDataMSFiles.OrbitrapStdFile1);
var scanSet = new ScanSet(6005);
run.CurrentScanSet = scanSet;
var msgen = MSGeneratorFactory.CreateMSGenerator(run.MSFileType);
var peakDetector = new DeconToolsPeakDetectorV2(1.3, 2, DeconTools.Backend.Globals.PeakFitType.QUADRATIC, true);
var decon = new HornDeconvolutor();
msgen.Execute(run.ResultCollection);
peakDetector.Execute(run.ResultCollection);
decon.Execute(run.ResultCollection);
peakAssociator.Execute(run.ResultCollection);
var numPeaksAssociatedWithFeatures = 0;
foreach (var msfeature in run.ResultCollection.ResultList)
{
foreach (var peak in msfeature.IsotopicProfile.Peaklist)
{
numPeaksAssociatedWithFeatures++;
}
}
var numPeaksInPeakListWithAssociations = 0;
foreach (MSPeak peak in run.PeakList)
{
Console.WriteLine(peak.MSFeatureID + "\t" + peak.XValue + "\t" + peak.Height + "\t");
if (peak.MSFeatureID != -1)
{
numPeaksInPeakListWithAssociations++;
}
}
//I don't think I can assume this... msfeature1 might overlap and share a peak with msfeature2
//Assert.AreEqual(numPeaksAssociatedWithFeatures, numPeaksInPeakListWithAssociations);
Console.WriteLine("total peaks associated with MSFeatures = " + numPeaksAssociatedWithFeatures);
Console.WriteLine("total peaks in original peaklist that were associated = " + numPeaksInPeakListWithAssociations);
Console.WriteLine("fraction peaks assigned = " + (double)(numPeaksInPeakListWithAssociations / (double)run.PeakList.Count));
}
protected override void InitializeProcessingTasks()
{
MSGenerator = MSGeneratorFactory.CreateMSGenerator(Run.MSFileType);
PeakDetector = PeakDetectorFactory.CreatePeakDetector(NewDeconToolsParameters);
var moreSensitivePeakToBackgroundRatio = NewDeconToolsParameters.PeakDetectorParameters.PeakToBackgroundRatio / 2;
var moreSensitiveSigNoiseThresh = NewDeconToolsParameters.PeakDetectorParameters.SignalToNoiseThreshold;
_moreSensitiveMS1PeakDetector = new DeconToolsPeakDetectorV2(moreSensitivePeakToBackgroundRatio,
moreSensitiveSigNoiseThresh, NewDeconToolsParameters.PeakDetectorParameters.PeakFitType,
NewDeconToolsParameters.PeakDetectorParameters.IsDataThresholded);
_superSensitiveMS1PeakDetector = new DeconToolsPeakDetectorV2(0, 0, NewDeconToolsParameters.PeakDetectorParameters.PeakFitType,
NewDeconToolsParameters.PeakDetectorParameters.IsDataThresholded);
Deconvolutor = DeconvolutorFactory.CreateDeconvolutor(NewDeconToolsParameters);
//Will initialize these but whether or not they are used are determined elsewhere
ZeroFiller = new DeconToolsZeroFiller(NewDeconToolsParameters.MiscMSProcessingParameters.ZeroFillingNumZerosToFill);
Smoother = new SavitzkyGolaySmoother(NewDeconToolsParameters.MiscMSProcessingParameters.SavitzkyGolayNumPointsInSmooth,
NewDeconToolsParameters.MiscMSProcessingParameters.SavitzkyGolayOrder);
FitScoreCalculator = new DeconToolsFitScoreCalculator();
ResultValidator = new ResultValidatorTask();
PeakToMSFeatureAssociator = new PeakToMSFeatureAssociator();
_ms2PeakDetectorForCentroidData = new DeconToolsPeakDetectorV2(0, 0, Globals.PeakFitType.QUADRATIC, true)
{
RawDataType = Globals.RawDataType.Centroided
};
_ms2PeakDetectorForProfileData = new DeconToolsPeakDetectorV2(NewDeconToolsParameters.PeakDetectorParameters.PeakToBackgroundRatio,
NewDeconToolsParameters.PeakDetectorParameters.SignalToNoiseThreshold,
NewDeconToolsParameters.PeakDetectorParameters.PeakFitType,
NewDeconToolsParameters.PeakDetectorParameters.IsDataThresholded);
Check.Ensure(Deconvolutor is ThrashDeconvolutorV2, "Error. Currently the DeconMSn workflow only works with the ThrashV2 deconvolutor. Selected deconvolutor= " + Deconvolutor);
}