ShipConstField.cxx

Go to the documentation of this file.

// SPDX-License-Identifier: LGPL-3.0-or-later
// SPDX-FileCopyrightText: Copyright CERN for the benefit of the SHiP
// Collaboration
 
// -------------------------------------------------------------------------
// -----                    ShipConstField source file                  -----
// -----                Created 06/01/04  by M. Al/Turany              -----
// -----                Redesign 13/02/06  by V. Friese                -----
// -------------------------------------------------------------------------
#include "ShipConstField.h"
 
#include <iomanip>
#include <iostream>
 
#include "ShipFieldPar.h"
 
using std::cerr;
using std::cout;
using std::endl;
using std::setw;
 
// -----   Default constructor   -------------------------------------------
ShipConstField::ShipConstField()
    : FairField(),
      fXmin(0.),
      fXmax(0.),
      fYmin(0.),
      fYmax(0.),
      fZmin(0.),
      fZmax(0.),
      fBx(0.),
      fBy(0.),
      fBz(0.) {
  fType = 0;
}
// -------------------------------------------------------------------------
 
// -----   Standard constructor   ------------------------------------------
ShipConstField::ShipConstField(const char* name, Double_t xMin, Double_t xMax,
                               Double_t yMin, Double_t yMax, Double_t zMin,
                               Double_t zMax, Double_t bX, Double_t bY,
                               Double_t bZ)
    : FairField(name),
      fXmin(xMin),
      fXmax(xMax),
      fYmin(yMin),
      fYmax(yMax),
      fZmin(zMin),
      fZmax(zMax),
      fBx(bX),
      fBy(bY),
      fBz(bZ) {
  fType = 0;
}
// -------------------------------------------------------------------------
 
// --------   Constructor from CbmFieldPar   -------------------------------
ShipConstField::ShipConstField(ShipFieldPar* fieldPar)
    : FairField(),
      fXmin(0.),
      fXmax(0.),
      fYmin(0.),
      fYmax(0.),
      fZmin(0.),
      fZmax(0.),
      fBx(0.),
      fBy(0.),
      fBz(0.) {
  if (!fieldPar) {
    cerr << "-W- ShipConstField::ShipConstField: empty parameter container!"
         << endl;
    fType = 0;
  } else {
    fXmin = fieldPar->GetXmin();
    fXmax = fieldPar->GetXmax();
    fYmin = fieldPar->GetYmin();
    fYmax = fieldPar->GetYmax();
    fZmin = fieldPar->GetZmin();
    fZmax = fieldPar->GetZmax();
    fBx = fieldPar->GetBx();
    fBy = fieldPar->GetBy();
    fBz = fieldPar->GetBz();
    fType = fieldPar->GetType();
  }
}
// -------------------------------------------------------------------------
 
// -----   Destructor   ----------------------------------------------------
ShipConstField::~ShipConstField() {}
// -------------------------------------------------------------------------
 
// -----   Set field region   ----------------------------------------------
void ShipConstField::SetFieldRegion(Double_t xMin, Double_t xMax, Double_t yMin,
                                    Double_t yMax, Double_t zMin,
                                    Double_t zMax) {
  fXmin = xMin;
  fXmax = xMax;
  fYmin = yMin;
  fYmax = yMax;
  fZmin = zMin;
  fZmax = zMax;
}
// -------------------------------------------------------------------------
 
// -----   Set field values   ----------------------------------------------
void ShipConstField::SetField(Double_t bX, Double_t bY, Double_t bZ) {
  fBx = bX;
  fBy = bY;
  fBz = bZ;
}
// -------------------------------------------------------------------------
 
// -----   Get x component of field   --------------------------------------
Double_t ShipConstField::GetBx(Double_t x, Double_t y, Double_t z) {
  if (x < fXmin || x > fXmax || y < fYmin || y > fYmax || z < fZmin ||
      z > fZmax)
    return 0.;
  return fBx;
}
// -------------------------------------------------------------------------
 
// -----   Get y component of field   --------------------------------------
Double_t ShipConstField::GetBy(Double_t x, Double_t y, Double_t z) {
  if (x < fXmin || x > fXmax || y < fYmin || y > fYmax || z < fZmin ||
      z > fZmax)
    return 0.;
  return fBy;
}
// -------------------------------------------------------------------------
 
// -----   Get z component of field   --------------------------------------
Double_t ShipConstField::GetBz(Double_t x, Double_t y, Double_t z) {
  if (x < fXmin || x > fXmax || y < fYmin || y > fYmax || z < fZmin ||
      z > fZmax)
    return 0.;
  return fBz;
}
// -------------------------------------------------------------------------
 
// -----   Screen output   -------------------------------------------------
void ShipConstField::Print() {
  cout << "======================================================" << endl;
  cout << "----  " << fTitle << " : " << fName << endl;
  cout << "----" << endl;
  cout << "----  Field type    : constant" << endl;
  cout << "----" << endl;
  cout << "----  Field regions : " << endl;
  cout << "----        x = " << setw(4) << fXmin << " to " << setw(4) << fXmax
       << " cm" << endl;
  cout << "----        y = " << setw(4) << fYmin << " to " << setw(4) << fYmax
       << " cm" << endl;
  cout << "----        z = " << setw(4) << fZmin << " to " << setw(4) << fZmax
       << " cm" << endl;
  cout.precision(4);
  cout << "----  B = ( " << fBx << ", " << fBy << ", " << fBz << " ) kG"
       << endl;
  cout << "======================================================" << endl;
}
// -------------------------------------------------------------------------