FFISGW2Base.hh

Go to the documentation of this file.

///
/// @file  FFISGW2Base.hh
/// @brief Hammer base class for ISGW2 form factors
///

//**** This file is a part of the HAMMER library
//**** Copyright (C) 2016 - 2020 The HAMMER Collaboration
//**** HAMMER is licensed under version 3 of the GPL; see COPYING for details
//**** Please note the MCnet academic guidelines; see GUIDELINES for details

// -*- C++ -*-
#ifndef HAMMER_FF_ISGW2_BASE
#define HAMMER_FF_ISGW2_BASE

#include "Hammer/FormFactorBase.hh"

namespace Hammer {

    /// @brief Base class for ISGW2 form factors
    /// implementation matched to EvtGen
    ///
    /// Provides the rate contents as tensor...
    ///
    /// @ingroup FormFactors
    class FFISGW2Base : public FormFactorBase {

    public:
        FFISGW2Base();

        FFISGW2Base(const FFISGW2Base& other) = default;
        FFISGW2Base& operator=(const FFISGW2Base& other) = delete;
        FFISGW2Base(FFISGW2Base&& other) = delete;
        FFISGW2Base& operator=(FFISGW2Base&& other) = delete;
        virtual ~FFISGW2Base() = default;

        /// @brief method to evaluate the object on a specific particle set
        /// @param[in] parent  the parent Particle
        /// @param[in] daughters the daughters (and grand-daughters, if necessary) Particle list
        /// @param[in] references the parent Particle siblings (necessary e.g. for helicity amplitude phase conventions)
        virtual void eval(const Particle& parent, const ParticleList& daughters,
                          const ParticleList& references);

    protected:
        /// @brief
        /// @param[in] z
        /// @return
        double GetGammaji(double z) const;

        /// @brief
        /// @param[in] mq1
        /// @param[in] mq2
        /// @return
        double Getas(double mq1, double mq2) const;

        /// @brief
        /// @param[in] m
        /// @return
        double Getas(double m) const;

        /// @brief
        /// @param[in] point
        /// @param[in] masses
        /// @return
        virtual void evalAtPSPoint(const std::vector<double>& point, const std::vector<double>& masses = {}) = 0;

        /// @brief
        virtual void defineSettings() = 0;
        
    protected:
        void addRefs() const; 
        
    protected:
        double msb;
        double msd;
        double bb2;
        double mbb;
        double nf;
        double cf;
        double msq;
        double bx2;
        double mbx;
        double nfp;
    };

} // namespace Hammer

#endif