AtConstField.cxx

Go to the documentation of this file.

/********************************************************************************
 *    Copyright (C) 2014 GSI Helmholtzzentrum fuer Schwerionenforschung GmbH    *
 *                                                                              *
 *              This software is distributed under the terms of the             *
 *         GNU Lesser General Public Licence version 3 (LGPL) version 3,        *
 *                  copied verbatim in the file "LICENSE"                       *
 ********************************************************************************/
 
// -------------------------------------------------------------------------
// -----                    AtConstField source file                  -----
// -----                Created 06/01/04  by M. Al/Turany              -----
// -----                Redesign 13/02/06  by V. Friese                -----
// -------------------------------------------------------------------------
#include "AtConstField.h"
 
#include "AtFieldPar.h"
 
#include <FairField.h>
 
#include <TString.h>
 
#include <iomanip>
#include <iostream>
 
using std::cerr;
using std::cout;
using std::endl;
using std::setw;
 
// -----   Default constructor   -------------------------------------------
AtConstField::AtConstField()
   : FairField(), fXmin(0.), fXmax(0.), fYmin(0.), fYmax(0.), fZmin(0.), fZmax(0.), fBx(0.), fBy(0.), fBz(0.)
{
   fType = 0;
}
// -------------------------------------------------------------------------
 
// -----   Standard constructor   ------------------------------------------
AtConstField::AtConstField(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   -------------------------------
AtConstField::AtConstField(AtFieldPar *fieldPar)
   : FairField(), fXmin(0.), fXmax(0.), fYmin(0.), fYmax(0.), fZmin(0.), fZmax(0.), fBx(0.), fBy(0.), fBz(0.)
{
   if (!fieldPar) {
      cerr << "-W- AtConstField::AtConstField: 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   ----------------------------------------------------
AtConstField::~AtConstField() = default;
// -------------------------------------------------------------------------
 
// -----   Set field region   ----------------------------------------------
void AtConstField::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 AtConstField::SetField(Double_t bX, Double_t bY, Double_t bZ)
{
   fBx = bX;
   fBy = bY;
   fBz = bZ;
}
// -------------------------------------------------------------------------
 
// -----   Get x component of field   --------------------------------------
Double_t AtConstField::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 AtConstField::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 AtConstField::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 AtConstField::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;
}
// -------------------------------------------------------------------------
 
ClassImp(AtConstField)