Method getLogScore_

src/openms/source/ANALYSIS/DECHARGING/ILPDCWrapper.cpp:479–516

  } // !compute_slice

  double ILPDCWrapper::getLogScore_(const PairsType::value_type& pair, const FeatureMap& fm) const
  {
    double score;
    String e;
    if (getenv("M") != nullptr)
      e = String(getenv("M"));
    if (e.empty())
    {
      //std::cout << "1";
      score = pair.getCompomer().getLogP();
      /*double charge_enhance = 0;

      if (pairs[i].getCharge(0) == fm[pairs[i].getElementIndex(0)].getCharge())
          charge_enhance += log(0.9); else charge_enhance += log(0.1);

      if (pairs[i].getCharge(1) == fm[pairs[i].getElementIndex(1)].getCharge())
          charge_enhance += log(0.9); else charge_enhance += log(0.1);

      score += charge_enhance;
      */

    }
    else
    {
      //std::cout << "2";
      double rt_diff =  fabs(fm[pair.getElementIndex(0)].getRT() - fm[pair.getElementIndex(1)].getRT());
      // enhance correct charge
      double charge_enhance = ((pair.getCharge(0) == fm[pair.getElementIndex(0)].getCharge())
                              &&
                               (pair.getCharge(1) == fm[pair.getElementIndex(1)].getCharge()))
                              ? 100 : 1;
      score = charge_enhance * (1 / (pair.getMassDiff() + 1) + 1 / (rt_diff + 1));
    }

    //std::cout << "logscore: " << score << "\n";

    return score;
  }

}