Tools/CIC/cic/Event.cc

Go to the documentation of this file.

// \author Darko Veberic
// \date 2014

#include <cic/Event.h>
#include <utl/Math.h>

using namespace std;


namespace cic {

  vector<Event>
  Bootstrap(const vector<Event>& events)
  {
    const unsigned int n = events.size();
    vector<Event> ret;
    ret.reserve(n);
    for (unsigned int i = 0; i < n; ++i)
      ret.push_back(events[utl::gRandom3->Integer(n)]);
    return ret;
  }


  vector<Event>
  ResampleWithErrors(const vector<Event>& events)
  {
    vector<Event> ret;
    ret.reserve(events.size());
    for (const Event& e : events) {
      const double th = utl::gRandom3->Gaus(e.GetTheta(), e.fThetaErr);
      if (0 <= th && th <= M_PI) {
        const double lnS1000 = utl::gRandom3->Gaus(e.fLnS1000, e.fLnS1000Err);
        ret.push_back(Event(utl::Sin2(th), e.fThetaErr, lnS1000, e.fLnS1000Err));
      }
    }
    return ret;
  }

}