示例#1
0
        public void TestMaxEntDeconvoluter()
        {
            const string rawFileFolder = @"\\proto-11\MSXML_Cache\PBF_Gen_1_214\2015_4";
            const string fname         = "WHIM2_LoHi_T2DD_HCD_GF07_02";
            var          rawFile       = string.Format(@"{0}\{1}.pbf", rawFileFolder, fname);
            var          ms1ft         = string.Format(@"\\protoapps\UserData\Jungkap\CompRef\lowRes\{0}.ms1ft", fname);

            var run         = PbfLcMsRun.GetLcMsRun(rawFile, 1.4826, 0);
            var ms1ScanNums = run.GetMs1ScanVector();

            var featureFinder = new LcMsPeakMatrixLowResolution(run);

            foreach (var scan in ms1ScanNums)
            {
                var fts = featureFinder.DetectMs1Features(scan);
                //Console.WriteLine("{0}\t{1}",scan, fts.Count);
            }

            var features = featureFinder.GetLcMsFeatures();

            var writer = new StreamWriter(ms1ft);
            var id     = 1;

            writer.WriteLine("FeatureID\tMinScan\tMaxScan\tMinCharge\tMaxCharge\tMonoMass\tAbundance\tRepScan\tMaxElutionTime\tElutionLength\tLikelihoodRatio");
            foreach (var feature in features.OrderBy(f => f.Mass))
            {
                writer.WriteLine("{0}\t{1}\t{2}\t{3}\t{4}\t{5}\t{6}\t{7}\t{8}\t{9}\t{10}", id, feature.MinScanNum, feature.MaxScanNum, feature.MinCharge,
                                 feature.MaxCharge, feature.Mass, feature.Abundance, feature.RepresentativeScanNum, feature.MinElutionTime, feature.MaxElutionTime, 0);
                id++;
            }
            writer.Close();
        }
示例#2
0
        public void TestMaxEntDeconvoluter()
        {

            const string rawFileFolder = @"\\proto-11\MSXML_Cache\PBF_Gen_1_214\2015_4";
            const string fname = "WHIM2_LoHi_T2DD_HCD_GF07_02";
            string rawFile = string.Format(@"{0}\{1}.pbf", rawFileFolder, fname);
            string ms1ft = string.Format(@"\\protoapps\UserData\Jungkap\CompRef\lowRes\{0}.ms1ft", fname);

            var run = PbfLcMsRun.GetLcMsRun(rawFile, 1.4826, 0);
            var ms1ScanNums = run.GetMs1ScanVector();
            
            var featureFinder = new LcMsPeakMatrixLowResolution(run);

            foreach (var scan in ms1ScanNums)
            {
                var fts = featureFinder.DetectMs1Features(scan);
                //Console.WriteLine("{0}\t{1}",scan, fts.Count);
            }

            var features = featureFinder.GetLcMsFeatures();
            
            var writer = new StreamWriter(ms1ft);
            var id = 1;
            writer.WriteLine("FeatureID\tMinScan\tMaxScan\tMinCharge\tMaxCharge\tMonoMass\tAbundance\tRepScan\tMaxElutionTime\tElutionLength\tLikelihoodRatio");
            foreach (var feature in features.OrderBy(f => f.Mass))
            {
                writer.WriteLine("{0}\t{1}\t{2}\t{3}\t{4}\t{5}\t{6}\t{7}\t{8}\t{9}\t{10}", id, feature.MinScanNum, feature.MaxScanNum, feature.MinCharge,
                    feature.MaxCharge, feature.Mass,  feature.Abundance, feature.RepresentativeScanNum, feature.MinElutionTime, feature.MaxElutionTime, 0);
                id++;
            }            
            writer.Close();
        }