MatrixElement.h

Go to the documentation of this file.

#ifndef AMPGEN_TRANSITIONMATRIXELEMENT_H
#define AMPGEN_TRANSITIONMATRIXELEMENT_H
 
#include "AmpGen/AmplitudeRules.h"
#include "AmpGen/CompiledExpression.h"
#include "AmpGen/Particle.h"
 
namespace AmpGen {
  class MatrixElement : public CompiledExpression<void(complex_v *, const size_t *, const real_t *, const real_v *)> {
  public:
    using amp_type = CompiledExpression<void(complex_v *, const size_t *, const real_t *, const real_v *)>;
    MatrixElement() = default;
    MatrixElement(const Particle &dt, const TotalCoupling &coupling, const amp_type &amp);
    MatrixElement(const Particle &dt, const TotalCoupling &coupling, const MinuitParameterSet &mps, const std::map<std::string, unsigned> &evtFormat,
                  const bool &debugThis = false);
 
    const std::vector<complex_v> operator()(const Event &event) const;
    template <class... arg_types> auto operator()(arg_types... args) const { return amp_type::operator()(args...); }
    void debug(const Event &event) const;
    const std::string decayDescriptor() const { return decayTree.decayDescriptor(); }
    DebugSymbols db;
    Particle decayTree;
    TotalCoupling coupling;
    complex_t coefficient;
    bool workToDo = {false};
    unsigned size = {1}; // size of the returned object
  };
  std::vector<size_t> processIndex(const std::vector<MatrixElement> &tm, const std::string &label);
  size_t findIndex(const std::vector<MatrixElement> &tm, const std::string &decayDescriptor);
  std::vector<size_t> findIndices(const std::vector<MatrixElement> &tm, const std::string &decayDescriptor);
}
#endif