DecayingParticle.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_DecayingParticle_h
#define yap_DecayingParticle_h

#include "fwd/AmplitudeComponent.h"
#include "fwd/BlattWeisskopf.h"
#include "fwd/DataAccessor.h"
#include "fwd/DecayChannel.h"
#include "fwd/DecayTree.h"
#include "fwd/FreeAmplitude.h"
#include "fwd/Model.h"
#include "fwd/ParticleCombination.h"
#include "fwd/PhaseSpaceFactor.h"
#include "fwd/QuantumNumbers.h"
#include "fwd/SpinAmplitude.h"

#include "Particle.h"

#include <memory>

namespace yap {

class DecayingParticle : public Particle
{
protected:

    DecayingParticle(std::string name, const QuantumNumbers& q, double radialSize, std::shared_ptr<PhaseSpaceFactorFactory> phsp_factory);

public:

    static std::shared_ptr<DecayingParticle> create(std::string name, const QuantumNumbers& q, double radialSize, std::shared_ptr<PhaseSpaceFactorFactory> phsp_factory = DefaultPHSPFactory)
    { return std::shared_ptr<DecayingParticle>(new DecayingParticle(name, q, radialSize, phsp_factory)); }

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

    virtual bool consistent() const override;

    virtual void checkDecayChannel(const DecayChannel& c) const;

    virtual std::shared_ptr<DecayChannel> addChannel(std::shared_ptr<DecayChannel> c);

    std::shared_ptr<DecayChannel> addChannel(const ParticleVector& daughters);

    template <typename ... Types>
    std::shared_ptr<DecayChannel> addChannel(std::shared_ptr<Particle> A, std::shared_ptr<Particle> B, Types ... other_daughters)
    { ParticleVector V{A, B, other_daughters...}; return addChannel(V); }


    const DecayChannelVector& channels() const
    { return Channels_;}

    std::shared_ptr<RealParameter> radialSize()
    { return RadialSize_; }

    const BlattWeisskopfMap& blattWeisskopfs() const
    { return BlattWeisskopfs_; }


    void printDecayChain() const
    { printDecayChainLevel(0); }

    const Model* model() const override;

    friend Model;

    friend DecayChannel;

protected:

    virtual void addParticleCombination(const std::shared_ptr<ParticleCombination>& c) override;

    virtual void pruneParticleCombinations() override;

    virtual void registerWithModel() override;

    void fixSolitaryFreeAmplitudes();

    void printDecayChainLevel(int level) const;

    virtual void modifyDecayTree(DecayTree& dt);

    virtual std::shared_ptr<PhaseSpaceFactor> phaseSpaceFactor(const DecayChannel& dc, const SpinAmplitude& sa);

    std::shared_ptr<PhaseSpaceFactorFactory> phaseSpaceFactorFactory()
    { return PhaseSpaceFactorFactory_; }

private:

    DecayChannelVector Channels_;

    BlattWeisskopfMap BlattWeisskopfs_;

    std::shared_ptr<RealParameter> RadialSize_;

    DecayTreeVectorMap DecayTrees_;

    std::shared_ptr<PhaseSpaceFactorFactory> PhaseSpaceFactorFactory_;

};
    
std::string to_string(const DecayTreeVectorMap& m_dtv_map);

DecayTreeSet decay_trees(const DecayingParticle& dp);

template <typename Last, typename ... Predicates>
DecayTreeSet decay_trees(const DecayingParticle& dp, Last p, Predicates ... P)
{ return filter(decay_trees(dp), p, P...); }

template <typename Last, typename ... Predicates>
DecayTreeSet::value_type decay_tree(const DecayingParticle& dp, Last p, Predicates ... P)
{ return lone_elt(filter(decay_trees(dp), p, P...)); }

FreeAmplitudeSet free_amplitudes(const DecayingParticle& dp);

template <typename Last, typename ... UnaryPredicates>
FreeAmplitudeSet free_amplitudes(const DecayingParticle& dp, Last p, UnaryPredicates ... P)
{ return filter(free_amplitudes(dp), p, P...); }

template <typename Last, typename ... UnaryPredicates>
FreeAmplitudeSet::value_type free_amplitude(const DecayingParticle& dp, Last p, UnaryPredicates ... P)
{ return lone_elt(free_amplitudes(dp, p, P...)); }

ParticleSet particles(DecayingParticle& dp);

template <typename Last, typename ... UnaryPredicates>
DecayChannelSet decay_channels(const DecayingParticle& dp, Last p, UnaryPredicates ... P)
{ return filter(DecayChannelSet(dp.channels().begin(), dp.channels().end()), p, P...); }

template <typename Last, typename ... UnaryPredicates>
DecayChannelSet::value_type decay_channel(const DecayingParticle& dp, Last p, UnaryPredicates ... P)
{ return lone_elt(decay_channels(dp, p, P...)); }

}

#endif