ShipFieldPar.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
 
#include "ShipFieldPar.h"
 
#include <iostream>
 
#include "FairParamList.h"
#include "ShipConstField.h"
using std::cerr;
using std::cout;
using std::endl;
 
// ------   Constructor   --------------------------------------------------
ShipFieldPar::ShipFieldPar(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.)
 
{}
// -------------------------------------------------------------------------
 
ShipFieldPar::ShipFieldPar()
    : 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   ---------------------------------------------------
ShipFieldPar::~ShipFieldPar() {}
// -------------------------------------------------------------------------
 
// ------   Put parameters   -----------------------------------------------
void ShipFieldPar::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 ShipFieldPar::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 ShipFieldPar::SetParameters(FairField* field) {
  if (!field) {
    cerr << "-W- ShipFieldPar::SetParameters: Empty field pointer!" << endl;
    return;
  }
 
  fType = field->GetType();
 
  if (fType == 0) {  // constant field
    ShipConstField* fieldConst = dynamic_cast<ShipConstField*>(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- ShipFieldPar::SetParameters: Unknown field type " << fType
         << "!" << endl;
    fBx = fBy = fBz = 0.;
    fXmin = fXmax = fYmin = fYmax = fZmin = fZmax = 0.;
    fMapName = "";
    fPosX = fPosY = fPosZ = fScale = 0.;
  }
 
  return;
}
// -------------------------------------------------------------------------