G4tgrSolid.cc

Go to the documentation of this file.

//
// ********************************************************************
// * License and Disclaimer                                           *
// *                                                                  *
// * The  Geant4 software  is  copyright of the Copyright Holders  of *
// * the Geant4 Collaboration.  It is provided  under  the terms  and *
// * conditions of the Geant4 Software License,  included in the file *
// * LICENSE and available at  http://cern.ch/geant4/license .  These *
// * include a list of copyright holders.                             *
// *                                                                  *
// * Neither the authors of this software system, nor their employing *
// * institutes,nor the agencies providing financial support for this *
// * work  make  any representation or  warranty, express or implied, *
// * regarding  this  software system or assume any liability for its *
// * use.  Please see the license in the file  LICENSE  and URL above *
// * for the full disclaimer and the limitation of liability.         *
// *                                                                  *
// * This  code  implementation is the result of  the  scientific and *
// * technical work of the GEANT4 collaboration.                      *
// * By using,  copying,  modifying or  distributing the software (or *
// * any work based  on the software)  you  agree  to acknowledge its *
// * use  in  resulting  scientific  publications,  and indicate your *
// * acceptance of all terms of the Geant4 Software license.          *
// ********************************************************************
//
//
// $Id: G4tgrSolid.cc 66363 2012-12-18 09:12:54Z gcosmo $
//
//
// class G4tgrSolid

// History:
// - Created.                                 P.Arce, CIEMAT (November 2007)
// -------------------------------------------------------------------------

#include <map>
#include <set>

#include "G4tgrSolid.hh"

#include "G4SystemOfUnits.hh"
#include "G4tgrUtils.hh"
#include "G4tgrMessenger.hh"
#include "G4tgrVolumeMgr.hh"

// -------------------------------------------------------------------------
G4tgrSolid::G4tgrSolid()
{
}


// -------------------------------------------------------------------------
G4tgrSolid::~G4tgrSolid()
{
}


// -------------------------------------------------------------------------
G4tgrSolid::G4tgrSolid( const std::vector<G4String>& wl) 
{
  //---------- set name 
  theName = G4tgrUtils::GetString( wl[1] ); 

  //---------- set solid type
  theType = G4tgrUtils::GetString( wl[2] ); 

  //---------- create only vector<double> of theSolidParams
  FillSolidParams( wl );

  G4tgrVolumeMgr::GetInstance()->RegisterMe( this );


#ifdef G4VERBOSE
  if( G4tgrMessenger::GetVerboseLevel() >= 1 )
  {
     G4cout << " Created " << *this << G4endl;
  }
#endif

}


// -------------------------------------------------------------------------
const std::vector< std::vector<G4double>* > G4tgrSolid::GetSolidParams() const
{
  return theSolidParams;
}


// -------------------------------------------------------------------------
const G4String& G4tgrSolid::GetRelativeRotMatName() const
{
  return theName;  // Dummy ...
}


// -------------------------------------------------------------------------
G4ThreeVector G4tgrSolid::GetRelativePlace() const
{
  return G4ThreeVector(0,0,0);  // Dummy...
}


// -------------------------------------------------------------------------
void G4tgrSolid::FillSolidParams( const std::vector<G4String>& wl )
{
  //---- Setting which are angle parameters (for dimensions...)
  std::map< G4String, std::set<G4int> > angleParams;
  std::set<G4int> apar;
  apar.clear(); apar.insert(3); apar.insert(4);
  angleParams["TUBS"] = apar;
  apar.clear(); apar.insert(5); apar.insert(6);
  angleParams["CONS"] = apar;
  apar.clear(); apar.insert(3); apar.insert(4); apar.insert(5);
  angleParams["PARA"] = apar;
  apar.clear(); apar.insert(1); apar.insert(2); apar.insert(6); apar.insert(10);
  angleParams["TRAP"] = apar;
  apar.clear(); apar.insert(2); apar.insert(3); apar.insert(4); apar.insert(5);
  angleParams["SPHERE"] = apar;
  apar.clear(); apar.insert(3); apar.insert(4);
  angleParams["TORUS"] = apar;
  apar.clear(); apar.insert(0); apar.insert(1);
  angleParams["POLYCONE"] = apar;
  apar.clear(); apar.insert(0); apar.insert(1);
  angleParams["POLYHEDRA"] = apar;
  apar.clear(); apar.insert(2); apar.insert(3);
  angleParams["HYPE"] = apar;
  apar.clear(); apar.insert(0);
  angleParams["TWISTED_BOX"] = apar;
  apar.clear(); apar.insert(0); apar.insert(2); apar.insert(3); apar.insert(10);
  angleParams["TWISTED_TRAP"] = apar;
  apar.clear(); apar.insert(5);
  angleParams["TWISTED_TRD"] = apar;
  apar.clear(); apar.insert(0); apar.insert(4);
  angleParams["TWISTED_TUBS"] = apar;

  std::vector<G4double>* vd = new std::vector<G4double>;
  theSolidParams.push_back( vd );
  size_t noParRead = wl.size()-3;

  G4String solidType = wl[2];
  //--- Default unit (mm) if length, deg if angle
  for(size_t ii = 0; ii < noParRead; ii++)
  {
    G4bool isAngle = 0;
    std::map< G4String, std::set<G4int> >::iterator ite
       = angleParams.find(solidType);
    if( ite != angleParams.end() )
    {
      std::set<G4int> apar2 = (*ite).second;
      if( apar2.find(ii) != apar2.end() )
      {
        isAngle = 1;
        vd->push_back( G4tgrUtils::GetDouble( wl[3+ii], deg ));
#ifdef G4VERBOSE
        if( G4tgrMessenger::GetVerboseLevel() >= 3 )
        {
          G4cout << " G4tgrSolid::FillSolidParams() - Angle param found "
                 << solidType << " " << ii << G4endl;
        }
#endif
      }
    } 
    if(!isAngle)
    {
      vd->push_back( G4tgrUtils::GetDouble( wl[3+ii] ) );
    }
  }
}


// -------------------------------------------------------------------------
std::ostream& operator<<(std::ostream& os, const G4tgrSolid& sol)
{
  os << "G4tgrSolid= " << sol.theName
     << " of type " << sol.theType << " PARAMS: ";
  if( sol.theSolidParams.size() != 0 )
  {
    std::vector<G4double> solpar = *(sol.theSolidParams[0]); 
    for( size_t ii = 0; ii < solpar.size(); ii++)
    {
      os << solpar[ii] << " " ;
    }
  }
  os << G4endl;

  return os;
}