IncoherentSum.h

Go to the documentation of this file.

#ifndef AMPGEN_INCOHERENTSUM_H
#define AMPGEN_INCOHERENTSUM_H
 
#include <complex>
#include <iomanip>
#include <string>
#include <vector>
 
#include "AmpGen/EventList.h"
#include "AmpGen/CoherentSum.h"
#include "AmpGen/Integrator.h"
#include "AmpGen/MsgService.h"
#include "AmpGen/Event.h"
#include "AmpGen/Types.h"
 
namespace AmpGen {
  class EventType;
  class FitFraction;
  class LinearErrorPropagator;
  class MinuitParameterSet;

  class IncoherentSum : public CoherentSum {
  public:
    IncoherentSum(const EventType &eventType, const AmpGen::MinuitParameterSet &mps, const std::string &prefix = "Inco");

    real_t prob(const Event &evt) const;
    real_t operator()(const Event &evt) const { return prob(evt); }
    real_v operator()(const real_v *, const unsigned) const;
#if ENABLE_AVX
    double operator()(const double *, const unsigned) const;
#endif

    real_t prob_unnormalised(const Event &evt) const;

    real_t norm() const;
    complex_t norm(const size_t &i, const size_t &j) { return i == j ? m_normalisations.get(i, i) : 0; }
    complex_t norm(const size_t &i) { return m_normalisations.get(i, i); }
    real_t norm(const Bilinears &norms) const;
    std::vector<FitFraction> fitFractions(const LinearErrorPropagator &linProp);
 
    real_t prob_unnormalisedNoCache(const Event &evt) const;
    void debug(const Event &evt, const std::string &nameMustContain = "");
 
    std::function<real_t(const Event &)> evaluator(const EventList_type * = nullptr) const;
    KeyedFunctors<double(Event)> componentEvaluator(const EventList_type * = nullptr) const;
 
    // void prepare();
  };
} // namespace AmpGen
 
#endif