BendBase.h

Go to the documentation of this file.

#ifndef CLASSIC_BENDBASE_H
#define CLASSIC_BENDBASE_H
 
#include "AbsBeamline/Component.h"
#include "Fields/Definitions.h"
#include "Physics/Units.h"
 
#include <string>
#include <vector>
 
class BendBase: public Component {
public:
    BendBase();
    BendBase(const BendBase&);
    BendBase(const std::string&);
 
    virtual bool bends() const;
 
    double getChordLength() const;
    virtual void setBendAngle(double angle);
    double getBendAngle() const;
    virtual void setEntranceAngle(double entranceAngle);
    double getEntranceAngle() const;
    virtual double getExitAngle() const = 0;
    void setFullGap(double);
    double getFullGap() const;
 
    virtual void setDesignEnergy(const double& energy, bool changeable = true);
    double getDesignEnergy() const;
    std::vector<Vector_t> getDesignPath() const;
 
    void setFieldAmplitude(double k0, double k0s);
    double getFieldAmplitude() const;
 
    void setFieldMapFN(std::string fileName);
    std::string getFieldMapFN() const;
 
protected:
    double calcDesignRadius(double fieldAmplitude) const;
    double calcFieldAmplitude(double radius) const;
    double calcBendAngle(double chordLength, double radius) const;
    double calcDesignRadius(double chordLength, double angle) const;
    double calcGamma() const;
    double calcBetaGamma() const;
 
    double chordLength_m;
    double angle_m;         
    double entranceAngle_m; 
    Fieldmap fieldmap_m;      
    const bool fast_m = false; 
 
    double gap_m; 
 
    double designEnergy_m; 
    bool designEnergyChangeable_m;
    std::vector<Vector_t> refTrajMap_m;
 
    double fieldAmplitudeX_m; 
    double fieldAmplitudeY_m; 
 
    double fieldAmplitude_m;  
 
    std::string fileName_m;
};
 
inline
bool BendBase::bends() const {
    return true;
}
 
inline
double BendBase::getChordLength() const {
    return chordLength_m;
}
 
inline
void BendBase::setBendAngle(double angle) {
    angle_m = angle;
}
 
inline
double BendBase::getBendAngle() const {
    return angle_m;
}
 
inline
void BendBase::setEntranceAngle(double angle)
{
    entranceAngle_m = angle;
}
 
inline
double BendBase::getEntranceAngle() const {
    return entranceAngle_m;
}
 
inline
void BendBase::setFullGap(double gap) {
    gap_m = std::abs(gap);
}
 
inline
double BendBase::getFullGap() const {
    return gap_m;
}
 
inline
void BendBase::setDesignEnergy(const double& energy, bool changeable) {
    if (designEnergyChangeable_m) {
        designEnergy_m = std::abs(energy) * Units::MeV2eV;
        designEnergyChangeable_m = changeable;
    }
}
 
inline
double BendBase::getDesignEnergy() const {
    return designEnergy_m;
}
 
inline
double BendBase::getFieldAmplitude() const
{
    return fieldAmplitude_m;
}
 
inline
void BendBase::setFieldMapFN(std::string fileName) {
    fileName_m = fileName;
}
 
inline
std::string BendBase::getFieldMapFN() const {
    return fileName_m;
}
 
#endif