Chi2Estimator.h

Go to the documentation of this file.

#ifndef AMPGEN_CHI2ESTIMATOR_H
#define AMPGEN_CHI2ESTIMATOR_H
 
#include <functional>
#include <string>
#include <vector>
 
#include "AmpGen/BinDT.h"
 
#if ENABLE_AVX
#include "AmpGen/EventListSIMD.h"
#else 
#include "AmpGen/EventList.h"
#endif
 
 
namespace AmpGen
{
  class EventType; 
  class Event;
 
  class Chi2Estimator
  {
    #if ENABLE_AVX
    typedef EventListSIMD EventList_type;
    #else
    typedef EventList     EventList_type;
    #endif
    public:
    template <typename... argument_types> 
    Chi2Estimator( const EventList_type& dataEvents, const EventList_type& mcEvents,
                   const std::function<double( const Event& )>& fcn, 
                   const argument_types&... args ) : m_binning(dataEvents.begin(), dataEvents.end(), ArgumentPack(args...) ) 
    {
      doChi2(dataEvents, mcEvents, fcn);
    }
 
    double chi2() const;
    double nBins() const; 
    void writeBinningToFile( const std::string& filename ); 
    void doChi2( const EventList_type& dataEvents, const EventList_type& mcEvents,
                 const std::function<double( const Event& )>& fcn );
  private: 
    double  m_chi2;
    size_t  m_nBins;
    BinDT   m_binning;
  };
} // namespace AmpGen
 
#endif /* end of include guard: AMPGEN_CHI2ESTIMATOR_H */