DecayTreeVectorIntegral.h

Go to the documentation of this file.

/*  YAP - Yet another PWA toolkit
    Copyright 2015, Technische Universitaet Muenchen,
    Authors: Daniel Greenwald, Johannes Rauch

    This program is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with this program.  If not, see <http://www.gnu.org/licenses/>.
 */


#ifndef yap_DecayTreeVectorIntegral_h
#define yap_DecayTreeVectorIntegral_h

#include "IntegralElement.h"

#include "fwd/DecayTreeVectorIntegral.h"

#include "fwd/DecayTree.h"
#include "fwd/Model.h"

#include <array>
#include <complex>
#include <map>

namespace yap {

class DecayTreeVectorIntegral
{
public:

    explicit DecayTreeVectorIntegral(const DecayTreeVector& dtv);

    const RealIntegralElement integral(unsigned i) const;

    const RealIntegralElement integral(unsigned i, unsigned j) const;

    const DecayTreeVector& decayTrees() const
    { return DecayTrees_; }

    const Model* model() const;

    const RealIntegralElementVector& diagonals() const
    { return Diagonals_; }

    const ComplexIntegralElementMatrix& offDiagonals() const
    { return OffDiagonals_; }

    const ComplexIntegralElement component(unsigned i, unsigned j) const;

    // addition operator
    DecayTreeVectorIntegral& operator+=(const DecayTreeVectorIntegral& rhs);

    // multiplication operator
    DecayTreeVectorIntegral& operator*=(double rhs);

    friend class Integrator;

protected:

    DecayTreeVectorIntegral& reset();

private:

    DecayTreeVector DecayTrees_;

    RealIntegralElementVector Diagonals_;

    ComplexIntegralElementMatrix OffDiagonals_;

};

const RealIntegralElement integral(const DecayTreeVectorIntegral& dtvi);

const ComplexIntegralElementMatrix cached_integrals(const DecayTreeVectorIntegral& dtvi);

const ComplexIntegralElementMatrix integrals(const DecayTreeVectorIntegral& dtvi);

const RealIntegralElementVector diagonal_integrals(const DecayTreeVectorIntegral& dtvi);

const RealIntegralElementVector fit_fractions(const DecayTreeVectorIntegral& dtvi);


}

#endif