ParticleFactory.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_ParticleFactory_h
#define yap_ParticleFactory_h

#include "fwd/DecayingParticle.h"
#include "fwd/FinalStateParticle.h"
#include "fwd/MassShape.h"
#include "fwd/ParticleFactory.h"
#include "fwd/PhaseSpaceFactor.h"
#include "fwd/Resonance.h"

#include "QuantumNumbers.h"

#include <memory>
#include <string>
#include <vector>

namespace yap {

struct ParticleTableEntry : public QuantumNumbers {
    ParticleTableEntry(int pdg = 0, std::string name = "", QuantumNumbers q = QuantumNumbers(), double mass = -1, std::vector<double> parameters = {});

    bool consistent() const override;

    int PDG;

    std::string Name;

    double Mass;

    std::vector<double> MassShapeParameters;
};

class ParticleFactory
{
public:

    using value_type = ParticleTableEntry;

    using iterator = ParticleTableMap::iterator;

    std::shared_ptr<FinalStateParticle> fsp(int PDG) const;

    std::shared_ptr<DecayingParticle> decayingParticle(int PDG, double radialSize, std::shared_ptr<PhaseSpaceFactorFactory> phsp_factory = DefaultPHSPFactory) const;

    std::shared_ptr<Resonance> resonance(int PDG, double radialSize, std::shared_ptr<MassShape> massShape, std::shared_ptr<PhaseSpaceFactorFactory> phsp_factory = DefaultPHSPFactory) const;

    ParticleFactory& operator+=(const ParticleFactory& rhs);


    const ParticleTableEntry& operator[](int PDG) const;

    const ParticleTableEntry& operator[](std::string name) const
    { return (*this)[pdgCode(name)]; }

    const QuantumNumbers& quantumNumbers(int PDG) const
    { return static_cast<const QuantumNumbers&>((*this)[PDG]); }

    const QuantumNumbers& quantumNumbers(std::string name) const
    { return static_cast<const QuantumNumbers&>((*this)[name]); }

    std::pair<ParticleTableMap::iterator, bool> insert(const ParticleTableEntry& entry);

    ParticleTableMap::iterator insert(ParticleTableMap::iterator hint, const ParticleTableEntry& entry);

    friend std::insert_iterator<ParticleFactory> inserter(ParticleFactory& F)
    { return std::insert_iterator<ParticleFactory>(F, F.ParticleTable_.end()); }

    // find PDG number by particle name
    // \return PDG code number
    // \param name Particle name as listed in particle table
    int pdgCode(std::string name) const;


private:

    ParticleTableMap ParticleTable_;
};

}

#endif