BFLHandler.cxx

Go to the documentation of this file.

#include <iostream>
using std::cerr;

#include <fstream>
#include <string.h>

#ifndef CSTDLIB
#include <cstdlib>
#define CSTDLIB
#endif

#include <TVector3.h>
#include <stdlib.h>

#include "BField/BFLInterpolation.h" 
#include "BField/BFLWingedEdge.h"
#include "BField/BFLHandler.h"
#include "BField/BFLNode.h"

ClassImp(BFLHandler) 

BFLHandler::BFLHandler(Int_t DetectorToLoadDefaults)
{
  Char_t VoronoiFile[] = "NearDetVorD.default";
  Char_t MeshFile[]    = "NearDetMesh.default";
  Char_t BFieldFile[]  = "NearDetBFld.default";
  
  switch(DetectorToLoadDefaults) {
    case (1):
       strcpy(VoronoiFile,"FarDetVorD.default");
       strcpy(MeshFile,   "FarDetMesh.default");
       strcpy(BFieldFile, "FarDetBFld.default");
       break;
    case (0): 
       strcpy(VoronoiFile,"NearDetVorD.default");
       strcpy(MeshFile,   "NearDetMesh.default");
       strcpy(BFieldFile, "NearDetBFld.default");
       break;
    default:      
       strcpy(VoronoiFile,"NearDetVorD.default");
       strcpy(MeshFile,   "NearDetMesh.default");
       strcpy(BFieldFile, "NearDetBFld.default");
       break;
  } 
  
  fBField = new TObjArray(1000);     /* BF Vectors */ 
  fVor = new BFLWingedEdge();        /* A Voronoi diagram object */

  LoadBField(BFieldFile);            /* Fill fBField */ 
  fVor->Load(VoronoiFile);           /* Fill fVor */ 
  fVor->LoadNodalSet(MeshFile);      /* Fill fGenerators */

  fInt = new BFLInterpolation(fBField, fVor); 
}

BFLHandler::BFLHandler(Char_t VoronoiFile[], Char_t MeshFile[], 
                       Char_t BFieldFile[])
{ 
  fBField = new TObjArray(1000);     /* BF Vectors */ 
  fVor = new BFLWingedEdge();        /* A Voronoi diagram object */
  
  LoadBField(BFieldFile);            /* Fill fBField */ 
  fVor->Load(VoronoiFile);           /* Fill fVor */ 
  fVor->LoadNodalSet(MeshFile); 
  
  fInt = new BFLInterpolation(fBField, fVor); 
}

BFLHandler::~BFLHandler()
{
    delete fInt;    fInt=0;
    // already done by BFLInterpolation //delete fVor;    fVor=0;
    // already done by BFLInterpolation //delete fBField; fBField=0;
}

Int_t BFLHandler::LoadBField(Char_t BFieldFile[])
{
  Int_t NodeID = 0;
  Float_t bx, by, bz;
  Int_t NodeIDLast = -999;
  
  ifstream DATA(BFieldFile,std::ios::in);
  if(!DATA) {
     cerr << "No such file: " << BFieldFile << endl;
     exit(1);
  }

  while(!DATA.eof()){
     DATA >> NodeID >> bx >> by >> bz; 
     if (NodeID != NodeIDLast) 
        fBField->AddAtAndExpand(new TVector3(bx,by,bz),NodeID);     
     NodeIDLast = NodeID;
  }
  DATA.close(); 

  return 0;
} 

TVector3 BFLHandler::GetBField(const TVector3& GlobalPosVector,
                               const Bool_t updateSeedPolygon)
{
    return fInt->GetB(GlobalPosVector,updateSeedPolygon);
}

---