int main (void)
{
cout << "Hello world for data processing 2" << endl;
DataProcessing toDo;
toDo.Test();
cout << "Press enter to finish..." << endl;
getchar();
return 0;
}
int main(int argc, char *argv[])
{
//QTimer timerWriter;
QCoreApplication a(argc, argv);
DataProcessing processing;
QFile file(":/testFiles/origin.txt");
file.open(QIODevice::ReadOnly);
QString block = file.readAll();
QString blockOrigin = processing.getBlockOrigin(QDateTime::fromString("2015-01-20 11:30:39.7","yyyy-MM-dd hh:mm:ss.z"), QDateTime::fromString("2015-01-20 11:30:40.6","yyyy-MM-dd hh:mm:ss.z"));
QString blockPick = processing.getBlockPick(QDateTime::fromString("2015-01-20 11:30:54.5","yyyy-MM-dd hh:mm:ss.z"), QDateTime::fromString("2015-01-20 11:30:56.4","yyyy-MM-dd hh:mm:ss.z"));
std::cout << blockOrigin.toStdString() << std::endl;
std::set<DATEBLOCK> blockDateTime = processing.getDateTimeBlocks(block);
QList<ANIMATIONBLOCK> blockSecuence = processing.getSecuence(blockDateTime);
//SimulationPlanner planner("prueba",QDir(QDir::home ()),blockSecuence,blockSecuence);
a.exec ();
}
void MainWindow::on_pushButton_3_clicked()
{
if(myData.vector.size() != 0 && ui->label_3->text() != "")
{
myData.vector.remove(myData.vector.indexOf(ui->label_3->text()));
myData.Write(mFileName);
ui->label_3->setText("");
}
}
MainWindow::MainWindow(QWidget *parent) :
QMainWindow(parent),
ui(new Ui::MainWindow)
{
ui->setupUi(this);
if(myData.vector.size() == 0){
myData.labels.append(ui->label);
myData.labels.append(ui->label_2);
myData.labels.append(ui->label_3);
myData.labels.append(ui->label_4);
myData.Read(mFileName);
}
}
ExitCodes main_(int , const char**)
{
String in = getStringOption_("in");
String out_meta = getStringOption_("out");
String in_cached = in + ".cached";
String out_cached = out_meta + ".cached";
bool convert_back = getFlag_("convert_back");
if (!convert_back)
{
MapType exp;
CachedmzML cacher;
MzMLFile f;
cacher.setLogType(log_type_);
f.setLogType(log_type_);
f.load(in,exp);
cacher.writeMemdump(exp, out_cached);
DataProcessing dp;
std::set<DataProcessing::ProcessingAction> actions;
actions.insert(DataProcessing::FORMAT_CONVERSION);
dp.setProcessingActions(actions);
dp.setMetaValue("cached_data", "true");
for (Size i=0; i<exp.size(); ++i)
{
exp[i].getDataProcessing().push_back(dp);
}
std::vector<MSChromatogram<ChromatogramPeak> > chromatograms = exp.getChromatograms();
for (Size i=0; i<chromatograms.size(); ++i)
{
chromatograms[i].getDataProcessing().push_back(dp);
}
exp.setChromatograms(chromatograms);
cacher.writeMetadata(exp, out_meta);
}
else
{
MzMLFile f;
MapType meta_exp;
CachedmzML cacher;
MapType exp_reading;
cacher.setLogType(log_type_);
f.setLogType(log_type_);
f.load(in,meta_exp);
cacher.readMemdump(exp_reading, in_cached);
std::cout << " read back, got " << exp_reading.size() << " spectra " << exp_reading.getChromatograms().size() << " chromats " << std::endl;
{
for (Size i=0; i<meta_exp.size(); ++i)
{
for (Size j = 0; j < meta_exp[i].getDataProcessing().size(); j++)
{
DataProcessing& dp = meta_exp[i].getDataProcessing()[j];
if (dp.metaValueExists("cached_data"))
{
dp.removeMetaValue("cached_data");
}
}
}
std::vector<MSChromatogram<ChromatogramPeak> > chromatograms = meta_exp.getChromatograms();
for (Size i=0; i<chromatograms.size(); ++i)
{
for (Size j = 0; j < chromatograms[i].getDataProcessing().size(); j++)
{
DataProcessing& dp = chromatograms[i].getDataProcessing()[j];
if (dp.metaValueExists("cached_data"))
{
dp.removeMetaValue("cached_data");
}
}
}
}
if (meta_exp.size() != exp_reading.size())
{
std::cerr << " Both experiments need to have the same size!";
}
for (Size i=0; i<exp_reading.size(); ++i)
{
for (Size j = 0; j < exp_reading[i].size(); j++)
{
meta_exp[i].push_back(exp_reading[i][j]);
}
}
std::vector<MSChromatogram<ChromatogramPeak> > chromatograms = exp_reading.getChromatograms();
std::vector<MSChromatogram<ChromatogramPeak> > old_chromatograms = meta_exp.getChromatograms();
for (Size i=0; i<chromatograms.size(); ++i)
{
for (Size j = 0; j < chromatograms[i].size(); j++)
{
old_chromatograms[i].push_back(chromatograms[i][j]);
}
}
meta_exp.setChromatograms(old_chromatograms);
f.store(out_meta,meta_exp);
}
return EXECUTION_OK;
}
void ChromatogramExtractor::return_chromatogram(std::vector< OpenSwath::ChromatogramPtr > & chromatograms,
std::vector< ChromatogramExtractor::ExtractionCoordinates > & coordinates,
OpenMS::TargetedExperiment & transition_exp_used, SpectrumSettings settings,
std::vector<OpenMS::MSChromatogram<> > & output_chromatograms, bool ms1) const
{
typedef std::map<String, const ReactionMonitoringTransition* > TransitionMapType;
TransitionMapType trans_map;
for (Size i = 0; i < transition_exp_used.getTransitions().size(); i++)
{
trans_map[transition_exp_used.getTransitions()[i].getNativeID()] = &transition_exp_used.getTransitions()[i];
}
for (Size i = 0; i < chromatograms.size(); i++)
{
const OpenSwath::ChromatogramPtr & chromptr = chromatograms[i];
const ChromatogramExtractor::ExtractionCoordinates & coord = coordinates[i];
TargetedExperiment::Peptide pep;
OpenMS::ReactionMonitoringTransition transition;
OpenMS::MSChromatogram<> chrom;
// copy data
OpenSwathDataAccessHelper::convertToOpenMSChromatogram(chrom, chromptr);
chrom.setNativeID(coord.id);
// Create precursor and set
// 1) the target m/z
// 2) the isolation window (upper/lower)
// 3) the peptide sequence
Precursor prec;
if (ms1)
{
pep = transition_exp_used.getPeptideByRef(coord.id);
prec.setMZ(coord.mz);
chrom.setChromatogramType(ChromatogramSettings::BASEPEAK_CHROMATOGRAM);
}
else
{
transition = (*trans_map[coord.id]);
pep = transition_exp_used.getPeptideByRef(transition.getPeptideRef());
prec.setMZ(transition.getPrecursorMZ());
if (settings.getPrecursors().size() > 0)
{
prec.setIsolationWindowLowerOffset(settings.getPrecursors()[0].getIsolationWindowLowerOffset());
prec.setIsolationWindowUpperOffset(settings.getPrecursors()[0].getIsolationWindowUpperOffset());
}
// Create product and set its m/z
Product prod;
prod.setMZ(transition.getProductMZ());
chrom.setProduct(prod);
chrom.setChromatogramType(ChromatogramSettings::SELECTED_REACTION_MONITORING_CHROMATOGRAM);
}
prec.setMetaValue("peptide_sequence", pep.sequence);
chrom.setPrecursor(prec);
// Set the rest of the meta-data
chrom.setInstrumentSettings(settings.getInstrumentSettings());
chrom.setAcquisitionInfo(settings.getAcquisitionInfo());
chrom.setSourceFile(settings.getSourceFile());
for (Size i = 0; i < settings.getDataProcessing().size(); ++i)
{
DataProcessing dp = settings.getDataProcessing()[i];
dp.setMetaValue("performed_on_spectra", "true");
chrom.getDataProcessing().push_back(dp);
}
output_chromatograms.push_back(chrom);
}
}
DateTime time;
time.set("2000-10-09 08:07:40");
DataProcessing* ptr = 0;
DataProcessing* nullPointer = 0;
START_SECTION(DataProcessing())
ptr = new DataProcessing();
TEST_NOT_EQUAL(ptr, nullPointer)
END_SECTION
START_SECTION(~DataProcessing())
delete ptr;
END_SECTION
START_SECTION(const DateTime& getCompletionTime() const)
DataProcessing tmp;
TEST_EQUAL(tmp.getCompletionTime().get(),"0000-00-00 00:00:00");
END_SECTION
START_SECTION(void setCompletionTime(const DateTime& completion_time))
DataProcessing tmp;
tmp.setCompletionTime(time);
TEST_EQUAL(tmp.getCompletionTime()==time,true);
END_SECTION
START_SECTION(Software& getSoftware())
DataProcessing tmp;
TEST_EQUAL(tmp.getSoftware()==Software(),true)
END_SECTION
START_SECTION(const Software& getSoftware() const)
