PWIZ_API_DECL void initializeTiny(MSData& msd) { msd.id = "urn:lsid:psidev.info:mzML.instanceDocuments.tiny.pwiz"; // cvList msd.cvs = defaultCVList(); // fileDescription FileContent& fc = msd.fileDescription.fileContent; fc.set(MS_MSn_spectrum); fc.set(MS_centroid_spectrum); SourceFilePtr sfp(new SourceFile); sfp->id = "tiny1.yep"; sfp->name = "tiny1.yep"; sfp->location = "file://*****:*****@higglesworth.edu"); // paramGroupList ParamGroupPtr pg1(new ParamGroup); pg1->id = "CommonMS1SpectrumParams"; pg1->set(MS_MS1_spectrum); pg1->set(MS_positive_scan); msd.paramGroupPtrs.push_back(pg1); ParamGroupPtr pg2(new ParamGroup); pg2->id = "CommonMS2SpectrumParams"; pg2->set(MS_MSn_spectrum); pg2->set(MS_negative_scan); msd.paramGroupPtrs.push_back(pg2); // sampleList SamplePtr samplePtr(new Sample); samplePtr->id = "20090101 - Sample 1"; samplePtr->name = "Sample 1"; msd.samplePtrs.push_back(samplePtr); // instrumentConfigurationList InstrumentConfigurationPtr instrumentConfigurationPtr(new InstrumentConfiguration("LCQ Deca")); instrumentConfigurationPtr->set(MS_LCQ_Deca); instrumentConfigurationPtr->set(MS_instrument_serial_number,"23433"); instrumentConfigurationPtr->componentList.push_back(Component(MS_nanoelectrospray, 1)); instrumentConfigurationPtr->componentList.push_back(Component(MS_quadrupole_ion_trap, 2)); instrumentConfigurationPtr->componentList.push_back(Component(MS_electron_multiplier, 3)); SoftwarePtr softwareCompassXtract(new Software); softwareCompassXtract->id = "CompassXtract"; softwareCompassXtract->set(MS_CompassXtract); softwareCompassXtract->version = "2.0.5"; instrumentConfigurationPtr->softwarePtr = softwareCompassXtract; msd.instrumentConfigurationPtrs.push_back(instrumentConfigurationPtr); // softwareList SoftwarePtr softwareBioworks(new Software); softwareBioworks->id = "Bioworks"; softwareBioworks->set(MS_Bioworks); softwareBioworks->version = "3.3.1 sp1"; SoftwarePtr softwarepwiz(new Software); softwarepwiz->id = "pwiz"; softwarepwiz->set(MS_pwiz); softwarepwiz->version = "1.0"; msd.softwarePtrs.push_back(softwareBioworks); msd.softwarePtrs.push_back(softwarepwiz); msd.softwarePtrs.push_back(softwareCompassXtract); // dataProcessingList DataProcessingPtr dpCompassXtract(new DataProcessing); dpCompassXtract->id = "CompassXtract processing"; ProcessingMethod procCXT; procCXT.order = 1; procCXT.softwarePtr = softwareCompassXtract; procCXT.set(MS_deisotoping); procCXT.set(MS_charge_deconvolution); procCXT.set(MS_peak_picking); dpCompassXtract->processingMethods.push_back(procCXT); DataProcessingPtr dppwiz(new DataProcessing("pwiz_processing")); ProcessingMethod procpwiz; procpwiz.order = 2; procpwiz.softwarePtr = softwarepwiz; procpwiz.set(MS_Conversion_to_mzML); dppwiz->processingMethods.push_back(procpwiz); msd.dataProcessingPtrs.push_back(dpCompassXtract); //msd.dataProcessingPtrs.push_back(dppwiz); ScanSettingsPtr as1(new ScanSettings("tiny scan settings")); as1->sourceFilePtrs.push_back(sfp_parameters); Target t1; t1.set(MS_selected_ion_m_z, 1000, MS_m_z); Target t2; t2.set(MS_selected_ion_m_z, 1200, MS_m_z); as1->targets.push_back(t1); as1->targets.push_back(t2); msd.scanSettingsPtrs.push_back(as1); // run msd.run.id = "Experiment 1"; msd.run.defaultInstrumentConfigurationPtr = instrumentConfigurationPtr; msd.run.samplePtr = samplePtr; msd.run.startTimeStamp = "2007-06-27T15:23:45.00035"; msd.run.defaultSourceFilePtr = sfp; shared_ptr<SpectrumListSimple> spectrumList(new SpectrumListSimple); msd.run.spectrumListPtr = spectrumList; spectrumList->dp = dppwiz; spectrumList->spectra.push_back(SpectrumPtr(new Spectrum)); spectrumList->spectra.push_back(SpectrumPtr(new Spectrum)); spectrumList->spectra.push_back(SpectrumPtr(new Spectrum)); spectrumList->spectra.push_back(SpectrumPtr(new Spectrum)); spectrumList->spectra.push_back(SpectrumPtr(new Spectrum)); Spectrum& s19 = *spectrumList->spectra[0]; s19.id = "scan=19"; s19.index = 0; s19.set(MS_ms_level, 1); s19.set(MS_centroid_spectrum); s19.set(MS_lowest_observed_m_z, 400.39, MS_m_z); s19.set(MS_highest_observed_m_z, 1795.56, MS_m_z); s19.set(MS_base_peak_m_z, 445.347, MS_m_z); s19.set(MS_base_peak_intensity, 120053, MS_number_of_detector_counts); s19.set(MS_total_ion_current, 1.66755e+007); s19.paramGroupPtrs.push_back(pg1); s19.scanList.scans.push_back(Scan()); s19.scanList.set(MS_no_combination); Scan& s19scan = s19.scanList.scans.back(); s19scan.instrumentConfigurationPtr = instrumentConfigurationPtr; s19scan.set(MS_scan_start_time, 5.890500, UO_minute); s19scan.set(MS_filter_string, "+ c NSI Full ms [ 400.00-1800.00]"); s19scan.set(MS_preset_scan_configuration, 3); s19scan.scanWindows.resize(1); ScanWindow& window = s19.scanList.scans.back().scanWindows.front(); window.set(MS_scan_window_lower_limit, 400.000000, MS_m_z); window.set(MS_scan_window_upper_limit, 1800.000000, MS_m_z); BinaryDataArrayPtr s19_mz(new BinaryDataArray); s19_mz->dataProcessingPtr = dpCompassXtract; s19_mz->set(MS_m_z_array, "", MS_m_z); s19_mz->data.resize(15); for (int i=0; i<15; i++) s19_mz->data[i] = i; BinaryDataArrayPtr s19_intensity(new BinaryDataArray); s19_intensity->dataProcessingPtr = dpCompassXtract; s19_intensity->set(MS_intensity_array, "", MS_number_of_detector_counts); s19_intensity->data.resize(15); for (int i=0; i<15; i++) s19_intensity->data[i] = 15-i; s19.binaryDataArrayPtrs.push_back(s19_mz); s19.binaryDataArrayPtrs.push_back(s19_intensity); s19.defaultArrayLength = s19_mz->data.size(); Spectrum& s20 = *spectrumList->spectra[1]; s20.id = "scan=20"; s20.index = 1; s20.paramGroupPtrs.push_back(pg2); s20.set(MS_ms_level, 2); s20.set(MS_profile_spectrum); s20.set(MS_lowest_observed_m_z, 320.39, MS_m_z); s20.set(MS_highest_observed_m_z, 1003.56, MS_m_z); s20.set(MS_base_peak_m_z, 456.347, MS_m_z); s20.set(MS_base_peak_intensity, 23433, MS_number_of_detector_counts); s20.set(MS_total_ion_current, 1.66755e+007); s20.precursors.resize(1); Precursor& precursor = s20.precursors.front(); precursor.spectrumID= s19.id; precursor.isolationWindow.set(MS_isolation_window_target_m_z, 445.3, MS_m_z); precursor.isolationWindow.set(MS_isolation_window_lower_offset, .5, MS_m_z); precursor.isolationWindow.set(MS_isolation_window_upper_offset, .5, MS_m_z); precursor.selectedIons.resize(1); precursor.selectedIons[0].set(MS_selected_ion_m_z, 445.34, MS_m_z); precursor.selectedIons[0].set(MS_peak_intensity, 120053, MS_number_of_detector_counts); precursor.selectedIons[0].set(MS_charge_state, 2); precursor.activation.set(MS_collision_induced_dissociation); precursor.activation.set(MS_collision_energy, 35.00, UO_electronvolt); s20.scanList.scans.push_back(Scan()); s20.scanList.set(MS_no_combination); Scan& s20scan = s20.scanList.scans.back(); s20scan.instrumentConfigurationPtr = instrumentConfigurationPtr; s20scan.set(MS_scan_start_time, 5.990500, UO_minute); s20scan.set(MS_filter_string, "+ c d Full ms2 [email protected] [ 110.00-905.00]"); s20scan.set(MS_preset_scan_configuration, 4); s20scan.scanWindows.resize(1); ScanWindow& window2 = s20scan.scanWindows.front(); window2.set(MS_scan_window_lower_limit, 110.000000, MS_m_z); window2.set(MS_scan_window_upper_limit, 905.000000, MS_m_z); BinaryDataArrayPtr s20_mz(new BinaryDataArray); s20_mz->dataProcessingPtr = dpCompassXtract; s20_mz->set(MS_m_z_array, "", MS_m_z); s20_mz->data.resize(10); for (int i=0; i<10; i++) s20_mz->data[i] = i*2; BinaryDataArrayPtr s20_intensity(new BinaryDataArray); s20_intensity->dataProcessingPtr = dpCompassXtract; s20_intensity->set(MS_intensity_array, "", MS_number_of_detector_counts); s20_intensity->data.resize(10); for (int i=0; i<10; i++) s20_intensity->data[i] = (10-i)*2; s20.binaryDataArrayPtrs.push_back(s20_mz); s20.binaryDataArrayPtrs.push_back(s20_intensity); s20.defaultArrayLength = s20_mz->data.size(); // spectrum with no data Spectrum& s21 = *spectrumList->spectra[2]; s21.id = "scan=21"; s21.index = 2; s21.paramGroupPtrs.push_back(pg1); s21.set(MS_ms_level, 1); s21.set(MS_centroid_spectrum); s21.userParams.push_back(UserParam("example", "spectrum with no data")); s21.setMZIntensityArrays(vector<double>(), vector<double>(), MS_number_of_detector_counts); s21.scanList.scans.push_back(Scan()); s21.scanList.scans.back().instrumentConfigurationPtr = instrumentConfigurationPtr; s21.scanList.set(MS_no_combination); // Cover ETD, ETD+SA, and ECD precursor activation mode usage Spectrum& s22 = *spectrumList->spectra[3]; s22.id = "scan=22"; s22.index = 3; s22.paramGroupPtrs.push_back(pg2); s22.set(MS_ms_level, 2); s22.set(MS_profile_spectrum); s22.set(MS_lowest_observed_m_z, 320.39, MS_m_z); s22.set(MS_highest_observed_m_z, 1003.56, MS_m_z); s22.set(MS_base_peak_m_z, 456.347, MS_m_z); s22.set(MS_base_peak_intensity, 23433, MS_number_of_detector_counts); s22.set(MS_total_ion_current, 1.66755e+007); s22.precursors.resize(1); Precursor& precursor22 = s22.precursors.front(); precursor22.spectrumID= s19.id; precursor22.isolationWindow.set(MS_isolation_window_target_m_z, 545.3, MS_m_z); precursor22.isolationWindow.set(MS_isolation_window_lower_offset, .5, MS_m_z); precursor22.isolationWindow.set(MS_isolation_window_upper_offset, .5, MS_m_z); precursor22.selectedIons.resize(1); precursor22.selectedIons[0].set(MS_selected_ion_m_z, 545.34, MS_m_z); precursor22.selectedIons[0].set(MS_peak_intensity, 120053, MS_number_of_detector_counts); precursor22.selectedIons[0].set(MS_charge_state, 2); precursor22.activation.set(MS_ETD); precursor22.activation.set(MS_CID); precursor22.activation.set(MS_collision_energy, 60.00, UO_electronvolt); s22.scanList.scans.push_back(Scan()); s22.scanList.set(MS_no_combination); Scan& s22scan = s22.scanList.scans.back(); s22scan.instrumentConfigurationPtr = instrumentConfigurationPtr; s22scan.set(MS_scan_start_time, 6.5, UO_minute); s22scan.set(MS_filter_string, "+ c d Full ms2 [email protected] [ 110.00-905.00]"); s22scan.set(MS_preset_scan_configuration, 4); s22scan.scanWindows.resize(1); window2 = s22scan.scanWindows.front(); window2.set(MS_scan_window_lower_limit, 110.000000, MS_m_z); window2.set(MS_scan_window_upper_limit, 905.000000, MS_m_z); BinaryDataArrayPtr s22_mz(new BinaryDataArray); s22_mz->dataProcessingPtr = dpCompassXtract; s22_mz->set(MS_m_z_array, "", MS_m_z); s22_mz->data.resize(10); for (int i=0; i<10; i++) s22_mz->data[i] = i*2; BinaryDataArrayPtr s22_intensity(new BinaryDataArray); s22_intensity->dataProcessingPtr = dpCompassXtract; s22_intensity->set(MS_intensity_array, "", MS_number_of_detector_counts); s22_intensity->data.resize(10); for (int i=0; i<10; i++) s22_intensity->data[i] = (10-i)*2; s22.binaryDataArrayPtrs.push_back(s22_mz); s22.binaryDataArrayPtrs.push_back(s22_intensity); s22.defaultArrayLength = s22_mz->data.size(); // spectrum with MALDI spot information Spectrum& s23 = *spectrumList->spectra[4]; s23.id = "sample=1 period=1 cycle=23 experiment=1"; s23.index = 4; s23.spotID = "A1,42x42,4242x4242"; s23.sourceFilePtr = sfp2; s23.set(MS_ms_level, 1); s23.set(MS_centroid_spectrum); s23.set(MS_lowest_observed_m_z, 142.39, MS_m_z); s23.set(MS_highest_observed_m_z, 942.56, MS_m_z); s23.set(MS_base_peak_m_z, 422.42, MS_m_z); s23.set(MS_base_peak_intensity, 42, MS_number_of_detector_counts); s23.set(MS_total_ion_current, 4200); s23.userParams.push_back(UserParam("alternate source file", "to test a different nativeID format")); s23.paramGroupPtrs.push_back(pg1); s23.scanList.scans.push_back(Scan()); s23.scanList.set(MS_no_combination); Scan& s23scan = s23.scanList.scans.back(); s23scan.instrumentConfigurationPtr = instrumentConfigurationPtr; s23scan.set(MS_scan_start_time, 42.0500, UO_second); s23scan.set(MS_filter_string, "+ c MALDI Full ms [100.00-1000.00]"); s23scan.scanWindows.resize(1); ScanWindow& window3 = s23scan.scanWindows.front(); window3.set(MS_scan_window_lower_limit, 100.000000, MS_m_z); window3.set(MS_scan_window_upper_limit, 1000.000000, MS_m_z); s23.binaryDataArrayPtrs.push_back(s19_mz); s23.binaryDataArrayPtrs.push_back(s19_intensity); s23.defaultArrayLength = s19_mz->data.size(); // chromatograms shared_ptr<ChromatogramListSimple> chromatogramList(new ChromatogramListSimple); msd.run.chromatogramListPtr = chromatogramList; chromatogramList->dp = dppwiz; chromatogramList->chromatograms.push_back(ChromatogramPtr(new Chromatogram)); chromatogramList->chromatograms.push_back(ChromatogramPtr(new Chromatogram)); Chromatogram& tic = *chromatogramList->chromatograms[0]; tic.id = "tic"; tic.index = 0; tic.defaultArrayLength = 15; tic.dataProcessingPtr = dpCompassXtract; tic.set(MS_total_ion_current_chromatogram); BinaryDataArrayPtr tic_time(new BinaryDataArray); tic_time->dataProcessingPtr = dppwiz; tic_time->set(MS_time_array, "", UO_second); tic_time->data.resize(15); for (int i=0; i<15; i++) tic_time->data[i] = i; BinaryDataArrayPtr tic_intensity(new BinaryDataArray); tic_intensity->dataProcessingPtr = dppwiz; tic_intensity->set(MS_intensity_array, "", MS_number_of_detector_counts); tic_intensity->data.resize(15); for (int i=0; i<15; i++) tic_intensity->data[i] = 15-i; tic.binaryDataArrayPtrs.push_back(tic_time); tic.binaryDataArrayPtrs.push_back(tic_intensity); Chromatogram& sic = *chromatogramList->chromatograms[1]; sic.id = "sic"; sic.index = 1; sic.defaultArrayLength = 10; sic.dataProcessingPtr = dppwiz; sic.set(MS_selected_ion_current_chromatogram); sic.precursor.isolationWindow.set(MS_isolation_window_target_m_z, 456.7, MS_m_z); sic.precursor.activation.set(MS_CID); sic.product.isolationWindow.set(MS_isolation_window_target_m_z, 678.9, MS_m_z); BinaryDataArrayPtr sic_time(new BinaryDataArray); sic_time->dataProcessingPtr = dppwiz; sic_time->set(MS_time_array, "", UO_second); sic_time->data.resize(10); for (int i=0; i<10; i++) sic_time->data[i] = i; BinaryDataArrayPtr sic_intensity(new BinaryDataArray); sic_intensity->dataProcessingPtr = dppwiz; sic_intensity->set(MS_intensity_array, "", MS_number_of_detector_counts); sic_intensity->data.resize(10); for (int i=0; i<10; i++) sic_intensity->data[i] = 10-i; sic.binaryDataArrayPtrs.push_back(sic_time); sic.binaryDataArrayPtrs.push_back(sic_intensity); } // initializeTiny()
PWIZ_API_DECL void addMIAPEExampleMetadata(TraData& td) { ContactPtr contactPtr(new Contact("JQP")); contactPtr->set(MS_contact_name, "John Q. Public"); contactPtr->set(MS_contact_affiliation, "Department of Redundancy Department"); contactPtr->set(MS_contact_address, "1600 Pennsylvania Ave."); td.contactPtrs.push_back(contactPtr); Publication publication; publication.id = "Al_et_al"; publication.set(MS_PubMed_identifier, 123456); td.publications.push_back(publication); InstrumentPtr lcqInstrumentPtr(new Instrument("LCQ")); lcqInstrumentPtr->set(MS_LCQ_Deca); lcqInstrumentPtr->set(MS_instrument_serial_number,"23433"); td.instrumentPtrs.push_back(lcqInstrumentPtr); InstrumentPtr qtrapInstrumentPtr(new Instrument("QTRAP")); qtrapInstrumentPtr->set(MS_4000_QTRAP); td.instrumentPtrs.push_back(qtrapInstrumentPtr); SoftwarePtr softwarePtr(new Software("Xcalibur")); softwarePtr->set(MS_Xcalibur); softwarePtr->version = "2.0.5"; SoftwarePtr softwareMaRiMba(new Software("MaRiMba")); softwareMaRiMba->set(MS_MaRiMba); softwareMaRiMba->version = "0.5"; SoftwarePtr softwarepwiz(new Software("pwiz")); softwarepwiz->set(MS_pwiz); softwarepwiz->version = "1.0"; td.softwarePtrs.push_back(softwareMaRiMba); td.softwarePtrs.push_back(softwarepwiz); td.softwarePtrs.push_back(softwarePtr); ProteinPtr proteinPtr(new Protein("Q123")); proteinPtr->set(MS_protein_accession, "Q123"); proteinPtr->sequence = "ABCD"; proteinPtr->set(MS_protein_name, "A short protein."); proteinPtr->set(MS_molecular_mass, 12345, UO_dalton); td.proteinPtrs.push_back(proteinPtr); PeptidePtr peptide1Ptr(new Peptide("Pep1")); peptide1Ptr->sequence = "AB"; peptide1Ptr->set(MS_theoretical_mass, 1234, UO_dalton); peptide1Ptr->proteinPtrs.push_back(proteinPtr); peptide1Ptr->modifications.push_back(Modification()); peptide1Ptr->modifications.back().set(UNIMOD_Methylmalonylation); peptide1Ptr->modifications.back().location = 1; peptide1Ptr->modifications.back().monoisotopicMassDelta = 123; peptide1Ptr->retentionTimes.push_back(RetentionTime()); peptide1Ptr->retentionTimes.back().set(MS_predicted_retention_time, 42, UO_minute); peptide1Ptr->retentionTimes.back().softwarePtr = softwareMaRiMba; peptide1Ptr->evidence.set(MS_confident_peptide_sequence_number, 6); td.peptidePtrs.push_back(peptide1Ptr); CompoundPtr compound1Ptr(new Compound("Cmp1")); compound1Ptr->set(MS_theoretical_mass, 1234, UO_dalton); compound1Ptr->retentionTimes.push_back(RetentionTime()); compound1Ptr->retentionTimes.back().set(MS_predicted_retention_time, 42, UO_minute); compound1Ptr->retentionTimes.back().softwarePtr = softwareMaRiMba; td.compoundPtrs.push_back(compound1Ptr); Interpretation interpretation1; interpretation1.set(MS_frag__y_ion); interpretation1.set(MS_product_ion_series_ordinal, 8); interpretation1.set(MS_product_ion_m_z_delta, 0.03); Interpretation interpretation2; interpretation2.set(MS_frag__b_ion___H2O); interpretation2.set(MS_product_ion_series_ordinal, 9); interpretation2.set(MS_product_ion_m_z_delta, -0.43); Validation qtrapValidation; qtrapValidation.set(MS_transition_optimized_on_specified_instrument); qtrapValidation.set(MS_4000_QTRAP); qtrapValidation.set(MS_peak_intensity, 4072, MS_percent_of_base_peak_times_100); qtrapValidation.set(MS_peak_intensity_rank, 2); qtrapValidation.set(MS_peak_targeting_suitability_rank, 1); Configuration qtrapConfiguration; qtrapConfiguration.instrumentPtr = qtrapInstrumentPtr; qtrapConfiguration.contactPtr = contactPtr; qtrapConfiguration.set(MS_dwell_time, 0.12, UO_second); qtrapConfiguration.set(MS_collision_gas, "argon"); qtrapConfiguration.set(MS_collision_gas_pressure, 12, UO_pascal); qtrapConfiguration.set(MS_collision_energy, 26, UO_electronvolt); qtrapConfiguration.set(MS_cone_voltage, 1200, UO_volt); qtrapConfiguration.set(MS_interchannel_delay, 0.1, UO_second); qtrapConfiguration.set(MS_tube_lens_voltage, 23, UO_volt); qtrapConfiguration.validations.push_back(qtrapValidation); Transition& tra0 = td.transitions[0]; tra0.peptidePtr = peptide1Ptr; tra0.precursor.set(MS_charge_state, 2); tra0.product.set(MS_charge_state, 1); tra0.prediction.softwarePtr = softwareMaRiMba; tra0.prediction.set(MS_transition_purported_from_an_MS_MS_spectrum_on_a_different__specified_instrument); tra0.prediction.set(MS_linear_ion_trap); tra0.prediction.set(MS_peak_intensity, 10000, UO_electronvolt); tra0.prediction.set(MS_peak_intensity_rank, 1); tra0.prediction.set(MS_peak_targeting_suitability_rank, 1); tra0.retentionTime.set(MS_local_retention_time, 40.02, UO_minute); tra0.retentionTime.set(MS_retention_time_window_lower_offset, 3, UO_minute); tra0.retentionTime.set(MS_retention_time_window_upper_offset, 3, UO_minute); tra0.configurationList.push_back(qtrapConfiguration); Transition& tra1 = td.transitions[1]; tra1.compoundPtr = compound1Ptr; tra1.precursor.set(MS_charge_state, 2); tra1.product.set(MS_charge_state, 1); Configuration lcqConfiguration; lcqConfiguration.instrumentPtr = lcqInstrumentPtr; lcqConfiguration.contactPtr = contactPtr; lcqConfiguration.set(MS_collision_energy, 42, UO_electronvolt); Target& tar0 = td.targets.targetIncludeList[0]; tar0.compoundPtr = compound1Ptr; tar0.precursor.set(MS_charge_state, 2); tar0.retentionTime.set(MS_local_retention_time, 12.3, UO_minute); tar0.retentionTime.set(MS_retention_time_window_lower_offset, 1, UO_minute); tar0.retentionTime.set(MS_retention_time_window_upper_offset, 2, UO_minute); tar0.configurationList.push_back(lcqConfiguration); td.targets.targetExcludeList[0].peptidePtr = peptide1Ptr; td.targets.targetExcludeList[0].precursor.set(MS_charge_state, 2); } // addMIAPEExampleMetadata()