31 #include "gsl/gsl_sf_pow_int.h"
33 using namespace endfieldmodel;
36 const std::string &
name):
39 planarArcGeometry_m(1.0, 1.0),
46 maxOrderX_m(right.maxOrderX_m),
47 recursion_m(right.recursion_m),
48 planarArcGeometry_m(right.planarArcGeometry_m),
49 angle_m(right.angle_m) {
62 double alpha =
atan(R[2] / (R[0] + radius ));
63 if (alpha != 0.0 &&
angle_m != 0.0) {
64 R[0] = R[2] /
sin(alpha) - radius;
65 R[2] = radius *
alpha;
76 B[0] = Bx *
cos(theta) - Bs *
sin(theta);
77 B[2] = Bx *
sin(theta) + Bs *
cos(theta);
83 while (maxOrder >= N) {
113 for (std::size_t j = 0;
117 for (std::size_t i = 0;
123 func += (
recursion_m.at(n).evaluatePolynomial(x / rho, i, j)
125 / gsl_sf_pow_int(rho, 2 * n - i - 2 * j);
128 func *= gsl_sf_pow_int(-1.0, n);
~MultipoleTCurvedConstRadius()
Interface for basic beam line object.
std::vector< polynomial::RecursionRelation > recursion_m
void truncate(std::size_t highestXorder)
Tps< T > sin(const Tps< T > &x)
Sine.
virtual void transformCoords(Vector_t &R) override
virtual double getFn(const std::size_t &n, const double &x, const double &s) override
std::size_t getTransMaxOrder() const
MultipoleTCurvedConstRadius(const std::string &name)
constexpr double alpha
The fine structure constant, no dimension.
virtual double getScaleFactor(const double &x, const double &s) override
PETE_TUTree< FnArcTan, typename T::PETE_Expr_t > atan(const PETE_Expr< T > &l)
PartBunchBase< double, 3 > * RefPartBunch_m
virtual ElementBase * clone() const override
double getTransDeriv(const std::size_t &n, const double &x)
void resizeX(const std::size_t &xDerivatives)
Tps< T > cos(const Tps< T > &x)
Cosine.
double getFringeDeriv(const std::size_t &n, const double &s)
virtual double getRadius(const double &s) override
virtual void setMaxOrder(const std::size_t &maxOrder)
virtual void setMaxOrder(const std::size_t &maxOrder) override
virtual void transformBField(Vector_t &B, const Vector_t &R) override