NonProjCalGeomParams.cc

Go to the documentation of this file.

#include "NonProjCalGeomParams.hh"
#include "xercesc/dom/deprecated/DOM_Element.hpp"
#include "LCD_Document.hh"
#include "LCD_DocumentUtil.hh"
#include "globals.hh"
#include <cassert>
#include <iostream>
using std::cerr;
using std::endl;

XERCES_CPP_NAMESPACE_USE

// Full constructor
NonProjCalGeomParams::
NonProjCalGeomParams(double Rmin, double Rmax, double Zmax, int nLayers,
                     double preAbsThick, double postAbsThick, double sensThick,
                     double cell_len_Z, double cell_len_Phi )
  : _Rmin(0.), _Rmax(0.), _Zmax(0.)
    , _nLayers(0), _preAbs_thick(0.), _postAbs_thick(0.), _sens_thick(0.)
    , _cell_len_Z(0.), _cell_len_Phi(0.)
{
  _Rmin = Rmin;
  _Rmax = Rmax;
  _Zmax = Zmax;
  _nLayers = nLayers;
  _preAbs_thick = preAbsThick;
  _postAbs_thick = postAbsThick;
  _sens_thick = sensThick;
  _cell_len_Z = cell_len_Z;
  _cell_len_Phi = cell_len_Phi;
}

NonProjCalGeomParams::NonProjCalGeomParams(const DOM_Element& voldb)
  : _Rmin(0.), _Rmax(0.), _Zmax(0.)
    , _nLayers(0), _preAbs_thick(0.), _postAbs_thick(0.), _sens_thick(0.)
    , _cell_len_Z(0.), _cell_len_Phi(0.)
{
  DOM_Element celldb;    /* DOM pointer */

  double offset;         /* offset from 0 to 1 */

  /* retrieve parameters for cell calculations from volume element */
  if( voldb == DOM_Element() ) {
    cerr << "error: XML voldb is invalid" << endl;
    exit(1);
  }

  /* retrieve and convert number of layers from layering element */
  celldb = LCD_DocumentUtil::firstDesc(voldb,"layering");
  LCD_DocumentUtil::convertValue(celldb.getAttribute("n"),&_nLayers);

  /* retrieve and convert inner_r and outer_z from barrel_dimensions */
  celldb = LCD_DocumentUtil::firstDesc(voldb,"barrel_dimensions");
  LCD_DocumentUtil::convertValue(celldb.getAttribute("inner_r"),&_Rmin);
  LCD_DocumentUtil::convertValue(celldb.getAttribute("outer_z"),&_Zmax);
  _Rmin *= cm;
  _Zmax *= cm;

  /* retrieve and convert all variables from cell_info */
  celldb = LCD_DocumentUtil::firstDesc(voldb,"cell_info");
  if ( celldb == DOM_Element() ) {
    G4cerr <<"***** NonProjCalGeomParams Error: could not find cell_info!!!"
           << G4endl;
    exit(1);
  }
  LCD_DocumentUtil::convertValue(celldb.getAttribute("height"),&_cell_len_Z);
  LCD_DocumentUtil::convertValue(celldb.getAttribute("length"),&_cell_len_Phi);
  LCD_DocumentUtil::convertValue(celldb.getAttribute("offset"),&offset);
  LCD_DocumentUtil::convertValue(celldb.getAttribute("outer_r"),&_Rmax);
  _cell_len_Z *= cm;
  _cell_len_Phi *= cm;
  offset *= cm;
  _Rmax *= cm;

  if (_Rmax == 0) {
    printf("_Rmax was not set properly; defaulting to 256 cm\n");
    _Rmax = 256.0*cm;
  }

  /* retrieve and convert slice thickness information
     relies on sensitive slice being 2nd in XML */
  DOM_NodeList sliceList = voldb.getElementsByTagName("slice");
  int listLen = sliceList.getLength();

  // Make sure a multi-layering scheme can be used (just one of the sublayers
  // can be sensitive, otherwise assert(false) will crash)
  double temp;
  bool sens;
  bool sensRead = false;
  for(int iitem=0; iitem<listLen; ++iitem) {
    DOM_Node sliceNode = sliceList.item(iitem);
    celldb = (DOM_Element &)sliceNode;
    LCD_DocumentUtil::convertValue(celldb.getAttribute("width"),&temp);
    LCD_DocumentUtil::convertValue(celldb.getAttribute("sensitive"),&sens);
    if(sens) {
      if(sensRead) {
        G4cerr<< "*** Geometry problem: >1 sensitive layer!!"<< G4endl;
        assert(false);
      }
      sensRead = true;
      _sens_thick = temp * cm;
    }
    else {
      if(sensRead) _postAbs_thick += temp*cm;
      else         _preAbs_thick += temp*cm;
    }
  }

  // cross-validation of geometry
  double test = _Rmin + _nLayers*(_preAbs_thick+_sens_thick+_postAbs_thick);
  test = (test - _Rmax)/_Rmax;
  if( fabs(test) > 0.0001 ) {
    G4cout <<"***** Geometry cross-check error:\n"
           <<"  Please check outer_r for non-projective components"<< G4endl;
    assert( false );
  }
}

---