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