UnfoldIterative.cxx

Go to the documentation of this file.

#include "CAFAna/Unfold/UnfoldIterative.h"

#include <iostream>

namespace ana
{
  //----------------------------------------------------------------------
  UnfoldIterative::UnfoldIterative(const ReweightableSpectrum& rt,
                                   int nIterations)
    : fRT(rt), fNIterations(nIterations)
  {
    assert(nIterations > 0);
  }

  //----------------------------------------------------------------------
  UnfoldIterative::~UnfoldIterative()
  {
  }

  //----------------------------------------------------------------------
  Spectrum UnfoldIterative::Truth(const Spectrum& reco)
  {
    // Leave the member variable as a record of how the matrix started and work
    // with this copy.
    ReweightableSpectrum rt = fRT;

    rt.ReweightToRecoSpectrum(reco);

    Spectrum truth = rt.WeightingVariable();

    for(int i = 2; i <= fNIterations; ++i){
      rt = fRT; // put it back how it started
      rt.ReweightToTrueSpectrum(truth); // project back the other way
      rt.ReweightToRecoSpectrum(reco); // and then forward again
      truth = rt.WeightingVariable();
    }

    return truth.FakeData(reco.POT());
  }
}