OpalDegrader.cpp

Go to the documentation of this file.

//
// Class OpalDegrader
//   The DEGRADER element.
//
// Copyright (c) 200x - 2023, Paul Scherrer Institut, Villigen PSI, Switzerland
// All rights reserved
//
// This file is part of OPAL.
//
// OPAL 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.
//
// You should have received a copy of the GNU General Public License
// along with OPAL. If not, see <https://www.gnu.org/licenses/>.
//
#include "Elements/OpalDegrader.h"
#include "Attributes/Attributes.h"
#include "BeamlineCore/DegraderRep.h"
#include "Structure/ParticleMatterInteraction.h"


OpalDegrader::OpalDegrader():
    OpalElement(SIZE, "DEGRADER",
                "The \"DEGRADER\" element defines a degrader."),
    parmatint_m(nullptr) {
    itsAttr[XSIZE] = Attributes::makeReal
        ("XSIZE", "Horizontal axis of the transverse elliptical shape [m]", 1E6);

    itsAttr[YSIZE] = Attributes::makeReal
        ("YSIZE", "Vertical axis of the transverse elliptical shape [m]", 1E6);

    registerOwnership();

    setElement(new DegraderRep("DEGRADER"));
}


OpalDegrader::OpalDegrader(const std::string& name, OpalDegrader* parent):
    OpalElement(name, parent),
    parmatint_m(nullptr) {
    setElement(new DegraderRep(name));
}


OpalDegrader::~OpalDegrader() {
    if (parmatint_m)
        delete parmatint_m;
}


OpalDegrader* OpalDegrader::clone(const std::string& name) {
    return new OpalDegrader(name, this);
}


void OpalDegrader::update() {
    OpalElement::update();

    DegraderRep* deg = dynamic_cast<DegraderRep*>(getElement());

    double length = Attributes::getReal(itsAttr[LENGTH]);
    double xsize = Attributes::getReal(itsAttr[XSIZE]);
    double ysize = Attributes::getReal(itsAttr[YSIZE]);

    deg->setElementLength(length);
    deg->setDimensions(xsize, ysize);

    if (itsAttr[PARTICLEMATTERINTERACTION] && parmatint_m == nullptr) {
        const std::string matterDescriptor = Attributes::getString(itsAttr[PARTICLEMATTERINTERACTION]);
        ParticleMatterInteraction* orig = ParticleMatterInteraction::find(matterDescriptor);
        parmatint_m = orig->clone(matterDescriptor);
        parmatint_m->initParticleMatterInteractionHandler(*deg);
        deg->setParticleMatterInteraction(parmatint_m->handler_m);
    }

    // Transmit "unknown" attributes.
    OpalElement::updateUnknown(deg);
}