mosesdecoder/mert/SentenceLevelScorer.cpp

//
//  SentenceLevelScorer.cpp
//  mert_lib
//
//  Created by Hieu Hoang on 22/06/2012.
//  Copyright 2012 __MyCompanyName__. All rights reserved.
//

#include "SentenceLevelScorer.h"

#include <iostream>
#include <boost/spirit/home/support/detail/lexer/runtime_error.hpp>

using namespace std;

namespace MosesTuning
{

SentenceLevelScorer::SentenceLevelScorer(const string& name, const string& config)
    : Scorer(name, config),
      m_regularisationStrategy(REG_NONE),
      m_regularisationWindow(0) {
  Init();
}

SentenceLevelScorer::~SentenceLevelScorer() {}

void SentenceLevelScorer::Init() {
    // Configure regularisation.
    static string KEY_TYPE = "regtype";
    static string KEY_WINDOW = "regwin";
    static string KEY_CASE = "case";
    static string TYPE_NONE = "none";
    static string TYPE_AVERAGE = "average";
    static string TYPE_MINIMUM = "min";
    static string TRUE = "true";
    static string FALSE = "false";

    const string type = getConfig(KEY_TYPE, TYPE_NONE);
    if (type == TYPE_NONE) {
      m_regularisationStrategy = REG_NONE;
    } else if (type == TYPE_AVERAGE) {
      m_regularisationStrategy = REG_AVERAGE;
    } else if (type == TYPE_MINIMUM) {
      m_regularisationStrategy = REG_MINIMUM;
    } else {
      throw boost::lexer::runtime_error("Unknown scorer regularisation strategy: " + type);
    }
    cerr << "Using scorer regularisation strategy: " << type << endl;

    const string window = getConfig(KEY_WINDOW, "0");
    m_regularisationWindow = atoi(window.c_str());
    cerr << "Using scorer regularisation window: " << m_regularisationWindow << endl;

    const string preservecase = getConfig(KEY_CASE, TRUE);
    if (preservecase == TRUE) {
      m_enable_preserve_case = true;
    } else if (preservecase == FALSE) {
      m_enable_preserve_case = false;
    }
    cerr << "Using case preservation: " << m_enable_preserve_case << endl;
}

void  SentenceLevelScorer::score(const candidates_t& candidates, const diffs_t& diffs,
                                 statscores_t& scores)
{
  //cout << "*******SentenceLevelScorer::score" << endl;
  if (!m_score_data) {
    throw runtime_error("Score data not loaded");
  }
  //calculate the score for the candidates
  if (m_score_data->size() == 0) {
    throw runtime_error("Score data is empty");
  }
  if (candidates.size() == 0) {
    throw runtime_error("No candidates supplied");
  }
  const int numCounts = m_score_data->get(0,candidates[0]).size();
  vector<float> totals(numCounts);
  for (size_t i = 0; i < candidates.size(); ++i) {
    //cout << " i " << i << " candi " << candidates[i] ;
    ScoreStats stats = m_score_data->get(i,candidates[i]);
    if (stats.size() != totals.size()) {
      stringstream msg;
      msg << "Statistics for (" << "," << candidates[i] << ") have incorrect "
      << "number of fields. Found: " << stats.size() << " Expected: "
      << totals.size();
      throw runtime_error(msg.str());
    }
    //Add up scores for all sentences, would normally be just one score
    for (size_t k = 0; k < totals.size(); ++k) {
      totals[k] += stats.get(k);
      //cout << " stats " << stats.get(k) ;
    }
    //cout << endl;
  }
  //take average
  for (size_t k = 0; k < totals.size(); ++k) {
    //cout << "totals = " << totals[k] << endl;
    //cout << "cand = " << candidates.size() << endl;
    totals[k] /= candidates.size();
    //cout << "finaltotals = " << totals[k] << endl;
  }

  scores.push_back(calculateScore(totals));

  candidates_t last_candidates(candidates);
  //apply each of the diffs, and get new scores
  for (size_t i = 0; i < diffs.size(); ++i) {
    for (size_t j = 0; j < diffs[i].size(); ++j) {
      const size_t sid = diffs[i][j].first;
      const size_t nid = diffs[i][j].second;
      //cout << "sid = " << sid << endl;
      //cout << "nid = " << nid << endl;
      const size_t last_nid = last_candidates[sid];
      for (size_t k  = 0; k < totals.size(); ++k) {
        const float diff = m_score_data->get(sid,nid).get(k)
                           - m_score_data->get(sid,last_nid).get(k);
        //cout << "diff = " << diff << endl;
        totals[k] += diff/candidates.size();
        //cout << "totals = " << totals[k] << endl;
      }
      last_candidates[sid] = nid;
    }
    scores.push_back(calculateScore(totals));
  }

  //regularisation. This can either be none, or the min or average as described in
  //Cer, Jurafsky and Manning at WMT08
  if (m_regularisationStrategy == REG_NONE || m_regularisationWindow <= 0) {
    //no regularisation
    return;
  }

  //window size specifies the +/- in each direction
  statscores_t raw_scores(scores);//copy scores
  for (size_t i = 0; i < scores.size(); ++i) {
    size_t start = 0;
    if (i >= m_regularisationWindow) {
      start = i - m_regularisationWindow;
    }
    const size_t end = min(scores.size(), i + m_regularisationWindow+1);
    if (m_regularisationStrategy == REG_AVERAGE) {
      scores[i] = score_average(raw_scores, start, end);
    } else {
      scores[i] = score_min(raw_scores, start, end);
    }
  }
}

}