#include "hvector.h"
#include "hgeotube.h"

typedef HVector<float> FVec;

int HGeoTube::readPoints(HGeoInfo & geoInfo, HString & volName, HGeoParam & ele) {
#ifdef WITHORACLE
  ele.vol=FVec(0.0F,3);
  char* vName=volName;
// not yet implemeted !!!!!!!!
  cerr << "code for reading a TUBE from database is not yet implemeted" << endl;
  return HFAILURE;
#endif
  return HSUCCESS;
}


int HGeoTube::readPoints(HGeoInfo & geoInfo, HGeoParam & ele) {
  ele.vol=FVec(0.0F,3);
  geoInfo.fin >> ele.vol[0];
  geoInfo.fin >> ele.vol[1][0] >> ele.vol[1][1];
  geoInfo.fin >> ele.vol[2];
  return HSUCCESS;
}

void HGeoTube::calcShape(HGeoParam & ele, HGeoShapePar & p) {
  float fac=10.F;
  p.noPar=3;
  p.volPar=FVec(0.0F,3);
  p.volPar[0]=ele.vol[1][0] /fac;
  p.volPar[1]=ele.vol[1][1] /fac;
  p.volPar[2]=(float)ele.vol[2].norm(ele.vol[0]) /fac /2.F;
  return; 
}

int HGeoTube::posShape(HGeoParam & ele, HGeoParam & mo, HGeoShapePar & p) {
  if (mo.shape=="PGON" || mo.shape=="TUBE") {
    posInPgon(ele,mo,p);
    return HSUCCESS;
  }
  cerr << ele.shape << " in " << mo.shape << " is not yet implemeted" << endl;
  return HFAILURE;
}


void HGeoTube::posInPgon(HGeoParam & ele, HGeoParam & mo, HGeoShapePar & p) {
  float fac=10.F;  
  p.nRot=0;
  ele.geaPos[0]=0.0F;
  ele.geaPos[1]=0.0F;
  ele.geaPos[2]=ele.vol[0][2] + fac * p.volPar[2];
  p.pos=(ele.trans - mo.geaPos + ele.geaPos) / fac;
  return;
}