spectrum.cpp

Go to the documentation of this file.

/**
 * \file pappsomspp/core/processing/cbor/mzcbor/spectrum.h
 * \date 24/11/2025
 * \author Olivier Langella
 * \brief PSI spectrum object for mzML/mzCBOR
 */
 
/*******************************************************************************
 * Copyright (c) 2025 Olivier Langella <Olivier.Langella@universite-paris-saclay.fr>.
 *
 * This file is part of PAPPSOms-tools.
 *
 *     PAPPSOms-tools is free software: you can redistribute it and/or modify
 *     it under the terms of the GNU General Public License as published by
 *     the Free Software Foundation, either version 3 of the License, or
 *     (at your option) any later version.
 *
 *     PAPPSOms-tools is distributed in the hope that it will be useful,
 *     but WITHOUT ANY WARRANTY; without even the implied warranty of
 *     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *     GNU General Public License for more details.
 *
 *     You should have received a copy of the GNU General Public License
 *     along with PAPPSOms-tools.  If not, see <http://www.gnu.org/licenses/>.
 *
 ******************************************************************************/
 
 
#include "spectrum.h"
#include "pappsomspp/core/exception/exceptionnotfound.h"
 
 
void
pappso::cbor::mzcbor::Spectrum::fromCbor(CborStreamReader &reader, bool want_binary_data)
{
  QString txt_value;
  cborBinaryDataArray.clear();
  precursorList.clear();
  defaultArrayLength = 0;
  index              = 0;
  cvParamMap.clear();
  scanListCvParamMap.clear();
  scanCvParamMap.clear();
 
  if(reader.enterContainer())
    {
      while(reader.hasNext() && (!reader.isInvalid()))
        {
          if(reader.isString())
            {
              reader.decodeString(txt_value);
              qDebug() << txt_value;
              if(txt_value == "id")
                {
                  reader.decodeString(txt_value);
                  qDebug() << txt_value;
                  id = txt_value;
                }
              else if(txt_value == "index")
                {
                  index = reader.toUnsignedInteger();
                  reader.next();
                }
              else if(txt_value == "defaultArrayLength")
                {
                  defaultArrayLength = reader.toUnsignedInteger();
                  reader.next();
                }
              else if(txt_value == "cvParam")
                {
                  cvParamMap = CvParam::getCvParamsMapFromCbor(reader);
                }
              else if(txt_value == "scanList")
                {
                  // mpa_cborReader->next();
                  reader.enterContainer();
 
                  while(reader.hasNext() && (!reader.isInvalid()))
                    {
                      if(reader.isString())
                        {
                          if(reader.decodeString(txt_value))
                            {
                              qDebug() << txt_value;
                              if(txt_value == "cvParam")
                                {
                                  scanListCvParamMap = CvParam::getCvParamsMapFromCbor(reader);
                                }
                              else if(txt_value == "scan")
                                {
                                  // scan array
                                  reader.enterContainer();
                                  while(reader.hasNext())
                                    {
                                      auto new_cv = readScanCvParams(reader);
                                      scanCvParamMap.insert(new_cv.begin(), new_cv.end());
                                    }
                                  reader.leaveContainer();
                                }
                            }
                          else
                            {
                              reader.next();
                            }
                        }
                      else
                        {
                          reader.next();
                        }
                    }
                  reader.leaveContainer();
                }
              else if(txt_value == "precursorList")
                {
 
                  reader.enterContainer();
 
                  while(reader.hasNext() && (!reader.isInvalid()))
                    {
                      if(reader.isString())
                        {
                          if(reader.decodeString(txt_value))
                            {
                              qDebug() << txt_value;
 
                              if(txt_value == "cvParam")
                                {
                                  precursorListCvParamMap = CvParam::getCvParamsMapFromCbor(reader);
                                }
                              else if(txt_value == "precursor")
                                {
 
                                  reader
                                    .enterContainer(); //                                     array
 
                                  while(reader.hasNext() && (!reader.isInvalid()))
                                    {
                                      // scan array
                                      Precursor precursor;
                                      precursor.fromCbor(reader);
                                      precursorList.emplace_back(precursor);
                                    }
                                  reader
                                    .leaveContainer(); //                                     array
 
                                  // reader.next();
                                }
                              else
                                {
                                  reader.next();
                                }
                            }
                          else
                            {
                              reader.next();
                            }
                        }
                      else
                        {
                          reader.next();
                        }
                    }
                  reader.leaveContainer();
                }
              else if(txt_value == "binaryDataArray")
                {
                  if(want_binary_data)
                    {
 
                      reader.enterContainer(); // start array
 
                      while(reader.hasNext())
                        {
 
                          BinaryDataArray data_array;
                          data_array.fromCbor(reader);
                          binaryDataArrayList.emplace_back(data_array);
                        }
 
                      reader.leaveContainer(); // end array
                    }
                  else
                    {
                      reader.next();
                    }
                }
              else
                {
                  reader.next();
                }
            }
          // reader.next();
        }
 
      reader.leaveContainer();
    }
}
 
 
std::map<QString, pappso::cbor::mzcbor::CvParam>
pappso::cbor::mzcbor::Spectrum::readScanCvParams(CborStreamReader &reader)
{
  std::map<QString, CvParam> cv_params;
  reader.enterContainer();
  QString txt_value;
 
  while(reader.hasNext() && (!reader.isInvalid()))
    {
      if(reader.isString())
        {
          if(reader.decodeString(txt_value))
            {
              qDebug() << txt_value;
              if(txt_value == "cvParam")
                {
                  cv_params = CvParam::getCvParamsMapFromCbor(reader);
                }
              else
                {
                  reader.next();
                }
            }
        }
      else
        {
          reader.next();
        }
    }
 
  reader.leaveContainer();
  return cv_params;
}
 
 
uint
pappso::cbor::mzcbor::Spectrum::getMsLevel() const
{
  //        <cvParam cvRef="MS" accession="MS:1000511" value="1" name="ms level" />
 
  auto it = cvParamMap.find("MS:1000511");
  if(it == cvParamMap.end())
    {
      throw pappso::ExceptionNotFound(QObject::tr("MS level not found in cvParam map"));
    }
  qDebug() << it->first << " " << it->second.valueInt;
  return it->second.valueInt;
}
 
double
pappso::cbor::mzcbor::Spectrum::getRtInSeconds() const
{
  //            <cvParam cvRef="MS" accession="MS:1000016" value="0.0048758833" name="scan start
  //            time" unitAccession="UO:0000031" unitName="minute" unitCvRef="UO" />
 
  auto it = scanCvParamMap.find("MS:1000016");
  if(it == scanCvParamMap.end())
    {
      throw pappso::ExceptionNotFound(QObject::tr("retention time not found in cvParam map"));
    }
  double rt_seconds = it->second.valueDouble;
  if(it->second.unitAccession == "UO:0000031")
    {
      // // minutes
      rt_seconds = rt_seconds * 60;
    }
  return rt_seconds;
}
 
void
pappso::cbor::mzcbor::Spectrum::decodeTrace(pappso::Trace &trace)
{
  std::vector<pappso_double> xVector;
  std::vector<pappso_double> yVector;
  for(auto &binary_array : binaryDataArrayList)
    {
      if(binary_array.isIntensity())
        {
          binary_array.decodeVector(defaultArrayLength, yVector);
        }
      else if(binary_array.isMz())
        {
          binary_array.decodeVector(defaultArrayLength, xVector);
        }
      else
        {
 
          throw pappso::ExceptionNotFound(QObject::tr(
            "error in pappso::cbor::mzcbor::Spectrum::decodeTrace: binary vector type not found"));
        }
    }
 
  trace.initialize(xVector, yVector);
}