Arb8.cxx

Go to the documentation of this file.

// $Id$
 
// -----------------------------------------------------------------------
// The RootGM package of the Virtual Geometry Model
// Copyright (C) 2007, Ivana Hrivnacova
// All rights reserved.
//
// For the licensing terms see vgm/LICENSE.
// Contact: ivana@ipno.in2p3.fr
// -----------------------------------------------------------------------
 
//
// Class RootGM::Arb8
// --------------------
// VGM implementation for Root Arb8 solid.
//
// Author: Ivana Hrivnacova; IPN Orsay
 
#include "RootGM/solids/Arb8.h"
#include "RootGM/common/Units.h"
#include "RootGM/solids/SolidMap.h"
 
#include "TGeoArb8.h"
#include "TMath.h"
 
#include <cstdlib>
#include <iostream>
#include <math.h>
 
const int RootGM::Arb8::fgkNofVertices = 8;
 
//_____________________________________________________________________________
RootGM::Arb8::Arb8(
  const std::string& name, double hz, std::vector<VGM::TwoVector> vertices)
  : VGM::ISolid(), VGM::IArb8(), BaseVGM::VArb8(), fArb8(0)
{
 
  fArb8 = new TGeoArb8(hz / RootGM::Units::Length());
  fArb8->SetName(name.c_str());
 
  // set vertices
  for (unsigned int i = 0; i < vertices.size(); ++i) {
    fArb8->SetVertex(i, vertices[i].first / RootGM::Units::Length(),
      vertices[i].second / RootGM::Units::Length());
  }
 
  RootGM::SolidMap::Instance()->AddSolid(this, fArb8);
}
 
//_____________________________________________________________________________
RootGM::Arb8::Arb8(TGeoArb8* arb8)
  : VGM::ISolid(), VGM::IArb8(), BaseVGM::VArb8(), fArb8(arb8)
{
 
  RootGM::SolidMap::Instance()->AddSolid(this, fArb8);
}
 
//_____________________________________________________________________________
RootGM::Arb8::Arb8() : VGM::ISolid(), VGM::IArb8(), BaseVGM::VArb8(), fArb8(0)
{
}
 
//_____________________________________________________________________________
RootGM::Arb8::Arb8(const Arb8& rhs)
  : VGM::ISolid(rhs), VGM::IArb8(rhs), BaseVGM::VArb8(rhs), fArb8(0)
{
}
 
//_____________________________________________________________________________
RootGM::Arb8::~Arb8()
{
  //
}
 
//_____________________________________________________________________________
std::string RootGM::Arb8::Name() const { return fArb8->GetName(); }
 
//_____________________________________________________________________________
int RootGM::Arb8::NofVertices() const { return fgkNofVertices; }
 
//_____________________________________________________________________________
VGM::TwoVector RootGM::Arb8::Vertex(int index) const
{
  if (index < 0 || index >= fgkNofVertices) {
    std::cerr << "+++ Error  +++" << std::endl;
    std::cerr << "    Wrong vertex index: " << index << std::endl;
    exit(1);
  }
 
  Double_t* xy = fArb8->GetVertices();
  /*
    for ( int i=0; i< 2*fgkNofVertices; i++ )
      std::cout << xy[i]  <<  ", ";
    std::cout << std::endl;
  */
  // CHECK !!!
  return VGM::TwoVector(xy[2 * index] * RootGM::Units::Length(),
    xy[2 * index + 1] * RootGM::Units::Length());
}
 
//_____________________________________________________________________________
double RootGM::Arb8::TwistAngle(int index) const
{
  if (index < 0 || index >= 4) {
    std::cerr << "+++ Error  +++" << std::endl;
    std::cerr << "    Wrong twist angle index: " << index << std::endl;
    exit(1);
  }
 
  Double_t twistAngle = TMath::ATan(fArb8->GetTwist(index)) * TMath::RadToDeg();
 
  return twistAngle * RootGM::Units::Angle();
}
 
//_____________________________________________________________________________
double RootGM::Arb8::ZHalfLength() const
{
  return fArb8->GetDZ() * RootGM::Units::Length();
}