Particle.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_Particle_h
#define yap_Particle_h

#include "fwd/Particle.h"

#include "fwd/DataPoint.h"
#include "fwd/FinalStateParticle.h"
#include "fwd/Model.h"
#include "fwd/Parameter.h"
#include "fwd/Particle.h"
#include "fwd/Spin.h"
#include "fwd/StatusManager.h"

#include "ParticleCombination.h"
#include "QuantumNumbers.h"

#include <iostream>
#include <memory>
#include <string>

namespace yap {

class Particle :
    public std::enable_shared_from_this<Particle>
{
protected:

    Particle(std::string name, const QuantumNumbers& q)
        : std::enable_shared_from_this<Particle>(),
        QuantumNumbers_(q), Name_(name) {}

public:

    virtual bool consistent() const;


    const QuantumNumbers& quantumNumbers() const
    { return QuantumNumbers_; }

    const std::string& name() const
    { return Name_; }

    std::string& name()
    { return Name_; }

    virtual const Model* model() const = 0;

    const ParticleCombinationSet& particleCombinations() const
    { return ParticleCombinations_; }


    friend class DecayChannel;

protected:

    virtual void addParticleCombination(const std::shared_ptr<ParticleCombination>& pc) = 0;

    virtual void pruneParticleCombinations();

    virtual void registerWithModel() = 0;

private:

    QuantumNumbers QuantumNumbers_;

    std::string Name_;

    ParticleCombinationSet ParticleCombinations_;

};

const SpinVector spins(const ParticleVector& v);

const bool decays_to_full_final_state(const Particle& p);

inline std::string to_string(const Particle& p)
{ return p.name() + "(" + to_string(p.quantumNumbers()) + ")"; }

inline std::string to_string(const ParticleVector& p)
{
    return std::accumulate(p.begin(), p.end(), std::string(""),
                           [](std::string& s, const ParticleVector::value_type& p)
                           {return s += " " + p->name();}).erase(0, 1);
}

inline std::ostream& operator<<(std::ostream& os, const Particle& p)
{ os << to_string(p); return os; }

}

#endif