OpalVacuum.cpp

Go to the documentation of this file.

//
// Class OpalVacuum
//   Defines the VACUUM element and its attributes.
//
// Copyright (c) 2018 - 2021, Pedro Calvo, CIEMAT, Spain
// All rights reserved
//
// Implemented as part of the PhD thesis
// "Optimizing the radioisotope production of the novel AMIT
// superconducting weak focusing cyclotron"
//
// 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/OpalVacuum.h"
#include "Attributes/Attributes.h"
#include "BeamlineCore/VacuumRep.h"
#include "Structure/ParticleMatterInteraction.h"
 
OpalVacuum::OpalVacuum()
    : OpalElement(
              SIZE, "VACUUM",
              "The \"VACUUM\" element defines the vacuum conditions "
              "for beam stripping interactions."),
      parmatint_m(nullptr) {
    itsAttr[GAS] = Attributes::makePredefinedString(
            "GAS", "The composition of residual gas", {"AIR", "H2"});
 
    itsAttr[PRESSURE] = Attributes::makeReal("PRESSURE", " Pressure in the accelerator, [mbar]");
 
    itsAttr[PMAPFN] = Attributes::makeString("PMAPFN", "Filename for the Pressure fieldmap");
 
    itsAttr[PSCALE] = Attributes::makeReal("PSCALE", "Scale factor for the P-field", 1.0);
 
    itsAttr[TEMPERATURE] =
            Attributes::makeReal("TEMPERATURE", " Temperature of the accelerator, [K]");
 
    itsAttr[STOP] = Attributes::makeBool(
            "STOP", "Option whether stop tracking after beam stripping. Default: true", true);
 
    registerOwnership();
 
    setElement(new VacuumRep("VACUUM"));
}
 
OpalVacuum::OpalVacuum(const std::string& name, OpalVacuum* parent)
    : OpalElement(name, parent), parmatint_m(nullptr) {
    setElement(new VacuumRep(name));
}
 
OpalVacuum::~OpalVacuum() { delete parmatint_m; }
 
OpalVacuum* OpalVacuum::clone(const std::string& name) { return new OpalVacuum(name, this); }
 
void OpalVacuum::update() {
    OpalElement::update();
 
    VacuumRep* vac = dynamic_cast<VacuumRep*>(getElement());
 
    double length      = Attributes::getReal(itsAttr[LENGTH]);
    std::string gas    = Attributes::getString(itsAttr[GAS]);
    double pressure    = Attributes::getReal(itsAttr[PRESSURE]);
    std::string pmap   = Attributes::getString(itsAttr[PMAPFN]);
    double pscale      = Attributes::getReal(itsAttr[PSCALE]);
    double temperature = Attributes::getReal(itsAttr[TEMPERATURE]);
    bool stop          = Attributes::getBool(itsAttr[STOP]);
 
    vac->setElementLength(length);
    vac->setResidualGas(gas);
    vac->setPressure(pressure);
    vac->setPressureMapFN(pmap);
    vac->setPScale(pscale);
    vac->setTemperature(temperature);
    vac->setStop(stop);
 
    // Transmit "unknown" attributes.
    OpalElement::updateUnknown(vac);
}