EvtRareLbToLllFF.hh

Go to the documentation of this file.

/***********************************************************************
* Copyright 1998-2020 CERN for the benefit of the EvtGen authors       *
*                                                                      *
* This file is part of EvtGen.                                         *
*                                                                      *
* EvtGen 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.                                  *
*                                                                      *
* EvtGen 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 EvtGen.  If not, see <https://www.gnu.org/licenses/>.     *
***********************************************************************/

#ifndef EVTRARELBTOLLLFF_HH
#define EVTRARELBTOLLLFF_HH 1

// Include files

#include "EvtGenBase/EvtIdSet.hh"
#include "EvtGenBase/EvtParticle.hh"

#include "EvtGenModels/EvtRareLbToLllFFBase.hh"

#include <array>
#include <map>
#include <memory>
#include <string>

class EvtRareLbToLllFF : public EvtRareLbToLllFFBase {
  public:
    class FormFactorDependence final {
      public:
        FormFactorDependence();

        FormFactorDependence( const double al, const double ap );

        FormFactorDependence( const double a0, const double a2, const double a4,
                              const double al, const double ap );

        FormFactorDependence( const FormFactorDependence& other );

        FormFactorDependence* clone() const;

        void param( const double al, const double ap );

        void param( const double a0, const double a2, const double a4,
                    const double al, const double ap );

        double a0_;
        double a2_;
        double a4_;
        double al_;
        double ap_;
    };

    class FormFactorSet final {
      public:
        FormFactorSet();

        FormFactorSet( const FormFactorSet& other );

        EvtRareLbToLllFF::FormFactorDependence F1;
        EvtRareLbToLllFF::FormFactorDependence F2;
        EvtRareLbToLllFF::FormFactorDependence F3;
        EvtRareLbToLllFF::FormFactorDependence F4;

        EvtRareLbToLllFF::FormFactorDependence G1;
        EvtRareLbToLllFF::FormFactorDependence G2;
        EvtRareLbToLllFF::FormFactorDependence G3;
        EvtRareLbToLllFF::FormFactorDependence G4;

        EvtRareLbToLllFF::FormFactorDependence H1;
        EvtRareLbToLllFF::FormFactorDependence H2;
        EvtRareLbToLllFF::FormFactorDependence H3;
        EvtRareLbToLllFF::FormFactorDependence H4;
        EvtRareLbToLllFF::FormFactorDependence H5;
        EvtRareLbToLllFF::FormFactorDependence H6;
    };

    void init() override;

    void getFF( EvtParticle* parent, EvtParticle* lambda,
                EvtRareLbToLllFFBase::FormFactors& FF ) override;

  private:
    double func( const double p, EvtRareLbToLllFF::FormFactorDependence& dep );

    std::array<std::unique_ptr<EvtRareLbToLllFF::FormFactorSet>, 2> FF_;
    std::map<int, EvtRareLbToLllFF::FormFactorSet*> FFMap_;

    void DiracFF( EvtParticle* parent, EvtParticle* lambda,
                  EvtRareLbToLllFF::FormFactorSet& FFset,
                  EvtRareLbToLllFF::FormFactors& FF );

    void RaritaSchwingerFF( EvtParticle* parent, EvtParticle* lambda,
                            EvtRareLbToLllFF::FormFactorSet& FFset,
                            EvtRareLbToLllFF::FormFactors& FF );
};

#endif    // EVTRARELBTOLLLFF_HH