BfldMapVoronoi.cxx

Go to the documentation of this file.

// $Id: BfldMapVoronoi.cxx,v 1.9 2005/02/09 04:38:11 rhatcher Exp $
// 
// BfldMapVoronoi 
// 
// This structure provides services for holding the magnetic field components
// at a series of points (cf. BfldMesh).  It encapsulates knowledge 
// about the field values themselves, an possibly how to adjust them
// for external effects (chemistry, coil current, etc.) but doesn't
// explicitly hold any knowledge about the relative spatial placement
// of the field values; that it must glean from BfldMesh.
//
// This class is specialized to handle maps associated with Voronoi meshes
//
// Author:  R. Hatcher 2000.06.20
//
#include "BField/BfldMapVoronoi.h"

#include "BField/BfldMeshVoronoi.h"

#include "MessageService/MsgService.h"
CVSID("$Id: BfldMapVoronoi.cxx,v 1.9 2005/02/09 04:38:11 rhatcher Exp $");

// necessary so that one can check "InheritsFrom"
#include "TClass.h"
// for expanding filename
#include "TSystem.h"
#include "TMath.h"

#include <iostream>
#include <fstream>
#include <iomanip>

ClassImp(BfldMapVoronoi)

//_____________________________________________________________________________
BfldMapVoronoi::BfldMapVoronoi() 
   : BfldMap()
{

}

//_____________________________________________________________________________
BfldMapVoronoi::BfldMapVoronoi(BfldGrid::Grid_t grid, Int_t variant) 
   : BfldMap(grid,variant)
{
   // read in a Map

   // map from grid+variant to a file or database record
   switch (grid) {
   case BfldGrid::kNearCoarseV: 
      fDataSource = "NearDetBFld.default";
      fGeneratedCoilCurrent = 1.0; // should be part of the file
      break;
   default:
      MSG("Bfld",Msg::kWarning)
         << "BfldMapVoronoi not yet supported for grid " <<
         BfldGrid::AsString(grid) << " variant " << variant << endl;
      return;
      break;
   }

   Msg::LogLevel_t verbosity = Msg::kDebug;

   MSG("Bfld",verbosity) << "BfldMapVoronoi::ctor "
                         << fDataSource << endl;

   Int_t i = 0;
   Int_t NodeID = 0;
   Float_t bx, by, bz;
   Int_t NodeIDLast = -999;

   // allocate a container, don't know the size needed yet so 
   // take a guess ... actually since there's just but the one
   // file (see above) 
   Int_t initSize = 4092;
   ResizeVectors(initSize);

   ifstream DATA(fDataSource.c_str(),ios::in);
   if(!DATA) {
      MSG("Bfld",Msg::kWarning) << "BfldMapVoronoi No such file: " 
                                << fDataSource << endl;
      return;
   }

   while(!DATA.eof()){
      DATA >> NodeID >> bx >> by >> bz; 
      if (NodeID != NodeIDLast) AddGenerator(NodeID,bx,by,bz);
      NodeIDLast = NodeID;
      i++;
   }
   DATA.close(); 

   MSG("Bfld",verbosity) << "BfldMapVoronoi file read in " << i 
                         << " vectors " <<endl;

   fFilledOkay = kTRUE;
}

//_____________________________________________________________________________
BfldMapVoronoi::~BfldMapVoronoi()
{
   // dtor deletes owned components

   MSG("Bfld",Msg::kInfo) << "~BfldMapVoronoi variant " << fVariant << endl;

   fGrid    = BfldGrid::kUndefined;
   fVariant = -1;

}

//_____________________________________________________________________________

---