AtFieldPar.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"                       *
 ********************************************************************************/
 
// -------------------------------------------------------------------------
// -----                    AtFieldPar  file                           -----
// -----                Created 26/03/14  by M. Al-Turany              -----
// -------------------------------------------------------------------------
 
#include "AtFieldPar.h"
 
#include "AtConstField.h"
 
#include <FairField.h>
#include <FairParGenericSet.h>
#include <FairParamList.h>
 
#include <iostream>
 
using std::cerr;
using std::cout;
using std::endl;
 
const int kMaxLen = 2048;
 
// ------   Constructor   --------------------------------------------------
AtFieldPar::AtFieldPar(const char *name, const char *title, const char *context)
   : FairParGenericSet(name, title, context), fType(-1), fXmin(0.), fXmax(0.), fYmin(0.), fYmax(0.), fZmin(0.),
     fZmax(0.), fBx(0.), fBy(0.), fBz(0.), fMapName(""), fPosX(0.), fPosY(0.), fPosZ(0.), fScale(0.), fPeak(0.),
     fMiddle(0.)
 
{
}
// -------------------------------------------------------------------------
 
AtFieldPar::AtFieldPar()
   : FairParGenericSet(), fType(-1), fXmin(0.), fXmax(0.), fYmin(0.), fYmax(0.), fZmin(0.), fZmax(0.), fBx(0.), fBy(0.),
     fBz(0.), fMapName(""), fPosX(0.), fPosY(0.), fPosZ(0.), fScale(0.), fPeak(0.), fMiddle(0.)
 
{
}
// -------------------------------------------------------------------------
 
// ------   Destructor   ---------------------------------------------------
AtFieldPar::~AtFieldPar() = default;
// -------------------------------------------------------------------------
 
// ------   Put parameters   -----------------------------------------------
void AtFieldPar::putParams(FairParamList *list)
{
 
   if (!list)
      return;
 
   list->add("Field Type", fType);
 
   if (fType == 0) { // constant field
      list->add("Field min x", fXmin);
      list->add("Field max x", fXmax);
      list->add("Field min y", fYmin);
      list->add("Field max y", fYmax);
      list->add("Field min z", fZmin);
      list->add("Field max z", fZmax);
      list->add("Field Bx", fBx);
      list->add("Field By", fBy);
      list->add("Field Bz", fBz);
   } else if (fType >= 1 && fType <= kMaxFieldMapType) { //
      list->add("Field Peak", fPeak);
      list->add("Field Middle", fMiddle);
   } else if (fType >= 2 && fType <= kMaxFieldMapType) { // field map
      list->add("Field map name", fMapName);
      list->add("Field x position", fPosX);
      list->add("Field y position", fPosY);
      list->add("Field z position", fPosZ);
      list->add("Field scaling factor", fScale);
   }
}
// -------------------------------------------------------------------------
 
// --------   Get parameters   ---------------------------------------------
Bool_t AtFieldPar::getParams(FairParamList *list)
{
 
   if (!list)
      return kFALSE;
 
   if (!list->fill("Field Type", &fType))
      return kFALSE;
 
   if (fType == 0) { // constant field
      if (!list->fill("Field min x", &fXmin))
         return kFALSE;
      if (!list->fill("Field max x", &fXmax))
         return kFALSE;
      if (!list->fill("Field min y", &fYmin))
         return kFALSE;
      if (!list->fill("Field max y", &fYmax))
         return kFALSE;
      if (!list->fill("Field min z", &fZmin))
         return kFALSE;
      if (!list->fill("Field max z", &fZmax))
         return kFALSE;
      if (!list->fill("Field Bx", &fBx))
         return kFALSE;
      if (!list->fill("Field By", &fBy))
         return kFALSE;
      if (!list->fill("Field Bz", &fBz))
         return kFALSE;
 
   } else if (fType >= 1 && fType <= kMaxFieldMapType) {
 
      if (!list->fill("Field Peak", &fPeak))
         return kFALSE;
      if (!list->fill("Field Middle", &fMiddle))
         return kFALSE;
 
   } else if (fType >= 2 && fType <= kMaxFieldMapType) { // field map
      Text_t mapName[80];
      if (!list->fill("Field map name", mapName, 80))
         return kFALSE;
      fMapName = mapName;
      if (!list->fill("Field x position", &fPosX))
         return kFALSE;
      if (!list->fill("Field y position", &fPosY))
         return kFALSE;
      if (!list->fill("Field z position", &fPosZ))
         return kFALSE;
      if (!list->fill("Field scaling factor", &fScale))
         return kFALSE;
   }
 
   return kTRUE;
}
// -------------------------------------------------------------------------
 
void AtFieldPar::SetParameters(FairField *field)
{
 
   if (!field) {
      cerr << "-W- AtFieldPar::SetParameters: Empty field pointer!" << endl;
      return;
   }
 
   fType = field->GetType();
 
   if (fType == 0) { // constant field
      auto *fieldConst = dynamic_cast<AtConstField *>(field);
      fBx = fieldConst->GetBx();
      fBy = fieldConst->GetBy();
      fBz = fieldConst->GetBz();
      fXmin = fieldConst->GetXmin();
      fXmax = fieldConst->GetXmax();
      fYmin = fieldConst->GetYmin();
      fYmax = fieldConst->GetYmax();
      fZmin = fieldConst->GetZmin();
      fZmax = fieldConst->GetZmax();
      fMapName = "";
      fPosX = fPosY = fPosZ = fScale = 0.;
   } else {
      cerr << "-W- AtFieldPar::SetParameters: Unknown field type " << fType << "!" << endl;
      fBx = fBy = fBz = 0.;
      fXmin = fXmax = fYmin = fYmax = fZmin = fZmax = 0.;
      fMapName = "";
      fPosX = fPosY = fPosZ = fScale = 0.;
   }
 
   return;
}
// -------------------------------------------------------------------------
 
ClassImp(AtFieldPar)