BfldMeshVoronoi.cxx

Go to the documentation of this file.

// $Id: BfldMeshVoronoi.cxx,v 1.11 2005/02/03 23:00:29 rhatcher Exp $
// 
// BfldMeshVoronoi 
// 
// This structure provides services for a 2d mesh representing an array of
// points where magnetic field map is "known".
//
// This particular type is specialized to deal with a Voronoi mesh
//
// Author:  R. Hatcher 2000.06.20
//
#include "BField/BfldMeshVoronoi.h"

#include "MessageService/MsgService.h"
CVSID("$Id: BfldMeshVoronoi.cxx,v 1.11 2005/02/03 23:00:29 rhatcher Exp $");

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

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

ClassImp(BfldMeshVoronoi)

//_____________________________________________________________________________
BfldMeshVoronoi::BfldMeshVoronoi() 
   : BfldMesh()
{

}

//_____________________________________________________________________________
BfldMeshVoronoi::BfldMeshVoronoi(BfldGrid::Grid_t grid, Int_t variant) 
   : BfldMesh(grid,variant)
{
   // read in a Mesh ... generator positions currently in Map file

   // map from grid+variant to a file or database record
   switch (grid) {
   case BfldGrid::kNearCoarseV: fFilename = "NearDetMesh.default";
      break;
   default:
      MSG("Bfld",Msg::kWarning)
         << "BfldMeshVoronoi not yet supported for grid " <<
         BfldGrid::AsString(grid) << endl;
      return;
      break;
   }

   Msg::LogLevel_t verbosity = Msg::kDebug;

   MSG("Bfld",verbosity) << "BfldMeshVoronoi::ctor "
                         << fFilename << endl;

   Int_t ngen = 0;
   Int_t i = 0;
   Int_t NodeID = 0;
   Int_t VorCell;
   Int_t VorCellLast = -999;
   Float_t x, y;

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

   DATA >> ngen;
   // allocate a container
   fPositions = new TClonesArray("TVector3",ngen);

   while(!DATA.eof()){
      DATA >> NodeID >> VorCell >> x >> y;
      if (NodeID>=0 && VorCell != VorCellLast) {
         //         fPositions->Add(new TVector3(x,y,VorCell),NodeID); 
         new((*fPositions)[NodeID]) TVector3(x,y,VorCell);
         VorCellLast = VorCell;
      } else {
         MSG("Bfld",verbosity) 
            << "BfldMeshVoronoi::ctor NodeID<0 "
            << NodeID << " ( VorCell " << VorCell << " ) " << endl;
      }
      i++; 
   }
   DATA.close(); 

   MSG("Bfld",verbosity) << "BfldMeshVoronoi file read in " 
                         << i << " nodes " << endl;

   fReadOk = kTRUE;

}

//_____________________________________________________________________________
BfldMeshVoronoi::~BfldMeshVoronoi()
{
   MSG("Bfld",Msg::kDebug) 
     << "~BfldMeshVoroni grid " << BfldGrid::AsString(fGrid) 
     << endl;
}

//_____________________________________________________________________________
TVector3 BfldMeshVoronoi::GetGeneratorPosition(Int_t /* nodeid */)
{
   // return the TVector3 associated with a generator (mesh node)

   Float_t x = 0;
   Float_t y = 0;
   return TVector3(x,y,0.);
}

//_____________________________________________________________________________
Int_t BfldMeshVoronoi::FindANSYSGenerator(Int_t VoronCellID)
{
  Int_t NodeID = 1;
  TIter next(fPositions);
  TVector3 * genVec3;

  while( (genVec3 = (TVector3*)next()) ) {
     if( genVec3->Z() == VoronCellID ) return NodeID;
     NodeID++;  
  }

  return -1;
}

//_____________________________________________________________________________

---