QuarticSymmBreaking.h

#ifndef gridripper_phys_scalar_QuarticSymmBreaking_h
#define gridripper_phys_scalar_QuarticSymmBreaking_h

#include "ScalarPotential.h"
#include <gridripper/lang/IllegalArgumentException.h>

namespace gridripper {

class Parameters;

namespace phys { namespace scalar {

using namespace std;

class QuarticSymmBreaking: public ScalarPotential
{
private:
    GReal_t lambda;
    GReal_t sigmaSq;
    GReal_t C;

public:
    QuarticSymmBreaking(const string& str, const Parameters& p)
        throw(IllegalArgumentException&);

    GReal_t eval_V(GReal_t psi, GReal_t T) {
        double psisq = psi*psi;
        GReal_t w = psisq - sigmaSq;
        return 0.5*(lambda*w*w + C*T*T*psisq);
    }

    GReal_t eval_dVdpsi(GReal_t psi, GReal_t T) {
        double w = psi*psi - sigmaSq;
        return 2*lambda*w*psi + C*T*T*psi;
    }

    GReal_t get_lambda() const {
        return lambda;
    }

    GReal_t get_C() const {
        return C;
    }
};

} } } // namespace gridripper::phys::scalar

#endif /* gridripper_phys_scalar_QuarticSymmBreaking_h */