Flatte.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_Flatte_h
#define yap_Flatte_h

#include "fwd/Flatte.h"

#include "fwd/DataPartition.h"
#include "fwd/DataPoint.h"
#include "fwd/FinalStateParticle.h"
#include "fwd/Parameter.h"
#include "fwd/ParticleCombination.h"
#include "fwd/StatusManager.h"

#include "MassShapeWithNominalMass.h"

#include <array>
#include <complex>
#include <memory>

namespace yap {

class Flatte : public MassShapeWithNominalMass
{
public:

    Flatte(double m = -1) : MassShapeWithNominalMass(m) {}

    void add(FlatteChannel fc);

    const std::vector<FlatteChannel>& channels() const
    { return FlatteChannels_; }

    virtual void checkDecayChannel(const DecayChannel& c) const override;

    virtual bool consistent() const override;

protected:

    virtual void calculateT(DataPartition& D, const std::shared_ptr<ParticleCombination>& pc, unsigned si) const override;

private:

    std::vector<FlatteChannel> FlatteChannels_;

};

struct FlatteChannel {
    std::shared_ptr<RealParameter> Coupling;
    
    std::array<std::shared_ptr<FinalStateParticle>, 2> Particles;
    
    FlatteChannel(std::shared_ptr<RealParameter> coupling, FinalStateParticle& A, FinalStateParticle& B);
    
    FlatteChannel(double coupling, FinalStateParticle& A, FinalStateParticle& B) :
    FlatteChannel(std::make_shared<RealParameter>(coupling), A, B) {}
};

}

#endif