BmnFieldMapData.cxx

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                    BmnFieldMapData source file                -----
// -----                V. Friese                14/02/06              -----
// -------------------------------------------------------------------------
#include "BmnFieldMapData.h"
 
#include "BmnFieldMap.h"
#include "TArrayF.h"
 
// -------------   Default constructor  ----------------------------------
BmnFieldMapData::BmnFieldMapData()
    : TNamed()
    , fType(1)
    , fXmin(0.)
    , fXmax(0.)
    , fYmin(0.)
    , fYmax(0.)
    , fZmin(0.)
    , fZmax(0.)
    , fNx(0)
    , fNy(0)
    , fNz(0)
    , fBx(NULL)
    , fBy(NULL)
    , fBz(NULL)
{}
// ------------------------------------------------------------------------
 
// -------------   Standard constructor   ---------------------------------
BmnFieldMapData::BmnFieldMapData(const char* mapName)
    : TNamed(mapName, "CBM Field Map Data")
    , fType(1)
    , fXmin(0.)
    , fXmax(0.)
    , fYmin(0.)
    , fYmax(0.)
    , fZmin(0.)
    , fZmax(0.)
    , fNx(0)
    , fNy(0)
    , fNz(0)
    , fBx(NULL)
    , fBy(NULL)
    , fBz(NULL)
{}
// ------------------------------------------------------------------------
 
// -----   Constructor from BmnFieldMap   ------------------------------
BmnFieldMapData::BmnFieldMapData(const char* name, const BmnFieldMap& map)
    : TNamed(name, "CBM Field Map Data")
    , fType(map.GetType())
    , fXmin(map.GetXmin())
    , fXmax(map.GetXmax())
    , fYmin(map.GetYmin())
    , fYmax(map.GetYmax())
    , fZmin(map.GetZmin())
    , fZmax(map.GetZmax())
    , fNx(map.GetNx())
    , fNy(map.GetNy())
    , fNz(map.GetNz())
    , fBx(new TArrayF(*(map.GetBx())))
    , fBy(new TArrayF(*(map.GetBy())))
    , fBz(new TArrayF(*(map.GetBz())))
{
    // Take out scaling factor and convert from kG to T
    if (map.GetScale() == 0.) {
        Fatal("BmnFieldMap", "Zero scaling factor");
    }
    Double_t factor = map.GetScale() * 10.;
    Int_t index = 0;
    for (Int_t ix = 0; ix < fNx; ix++) {
        for (Int_t iy = 0; iy < fNy; iy++) {
            for (Int_t iz = 0; iz < fNz; iz++) {
                index = ix * fNy * fNz + iy * fNz + iz;
                if (fBx)
                    (*fBx)[index] = (*fBx)[index] / factor;
                if (fBy)
                    (*fBy)[index] = (*fBy)[index] / factor;
                if (fBz)
                    (*fBz)[index] = (*fBz)[index] / factor;
            }   // z loop
        }   // y loop
    }   // x loop
}
// ------------------------------------------------------------------------
 
// ------------   Destructor   --------------------------------------------
BmnFieldMapData::~BmnFieldMapData()
{
    if (fBx)
        delete fBx;
    if (fBy)
        delete fBy;
    if (fBz)
        delete fBz;
}
// ------------------------------------------------------------------------