BmnFieldPar.cxx

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                      BmnFieldPar source file                  -----
// -----                   Created 20/02/06  by V. Friese              -----
// -------------------------------------------------------------------------
#include "BmnFieldPar.h"
 
#include "BmnFieldConst.h"
#include "BmnFieldMap.h"
#include "FairParamList.h"
 
#include <iostream>
using std::cerr;
using std::cout;
using std::endl;
 
const int kMaxLen = 2048;
 
// ------   Constructor   --------------------------------------------------
 
BmnFieldPar::BmnFieldPar(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.)
    , fIsOff(kFALSE)
    , fDistortionFilename("")
    , fParentName("")
    , fTypeOfParent(0)
{}
// -------------------------------------------------------------------------
 
BmnFieldPar::BmnFieldPar()
    : 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.)
    , fIsOff(kFALSE)
    , fDistortionFilename("")
    , fParentName("")
    , fTypeOfParent(0)
{}
// -------------------------------------------------------------------------
 
// ------   Destructor   ---------------------------------------------------
 
BmnFieldPar::~BmnFieldPar() {}
// -------------------------------------------------------------------------
 
// ------   Put parameters   -----------------------------------------------
 
void BmnFieldPar::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);
        list->add("Field is off", fIsOff);
    } else if (fType >= 1 && 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);
        list->add("Field is off", fIsOff);
 
        if (fType == kTypeDistorted) {   // BmnFieldMapDistorted case
            list->add("Field map distortion filename", fDistortionFilename.Data());
            list->add("Field name of parent field", fParentName.Data());
            list->add("Field type of parent field", fTypeOfParent);
        }
    }
}
// -------------------------------------------------------------------------
 
// --------   Get parameters   ---------------------------------------------
 
Bool_t BmnFieldPar::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;
        if (!list->fill("Field is off", &fIsOff))
            return kFALSE;
    } else if (fType >= 1 && 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;
        if (!list->fill("Field is off", &fIsOff))
            return kFALSE;
 
        if (fType == kTypeDistorted) {   // BmnFieldMapDistorted case
            Text_t tmp[kMaxLen];
            fDistortionFilename = "";
            if (!list->fill("Field map distortion filename", tmp, kMaxLen))
                return kFALSE;
            fDistortionFilename = tmp;
            fParentName = "";
            if (!list->fill("Field name of parent field", tmp, kMaxLen))
                return kFALSE;
            fParentName = tmp;
            if (!list->fill("Field type of parent field", &fTypeOfParent))
                return kFALSE;
        }
    }
 
    return kTRUE;
}
// -------------------------------------------------------------------------
 
// ---------   Set parameters from BmnField   ------------------------------
 
void BmnFieldPar::SetParameters(FairField* field)
{
 
    if (!field) {
        cerr << "-W- BmnFieldPar::SetParameters: Empty field pointer!" << endl;
        return;
    }
 
    fType = field->GetType();
 
    if (fType == 0) {   // constant field
        BmnFieldConst* fieldConst = (BmnFieldConst*)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.;
        fIsOff = fieldConst->IsFieldOff();
    } else if (fType >= 1 && fType <= kMaxFieldMapType) {   // field map
        BmnFieldMap* fieldMap = (BmnFieldMap*)field;
        fBx = fBy = fBz = 0.;
        fXmin = fXmax = fYmin = fYmax = fZmin = fZmax = 0.;
 
        fMapName = field->GetName();
        fPosX = fieldMap->GetPositionX();
        fPosY = fieldMap->GetPositionY();
        fPosZ = fieldMap->GetPositionZ();
        fScale = fieldMap->GetScale();
        fIsOff = fieldMap->IsFieldOff();
    } else {
        cerr << "-W- BmnFieldPar::SetParameters: Unknown field type " << fType << "!" << endl;
        fBx = fBy = fBz = 0.;
        fXmin = fXmax = fYmin = fYmax = fZmin = fZmax = 0.;
        fMapName = "";
        fPosX = fPosY = fPosZ = fScale = 0.;
        fIsOff = kFALSE;
    }
 
    return;
}