ERRTelescopeGeoComponentSingleSi.cxx

/********************************************************************************
 *              Copyright (C) Joint Institute for Nuclear Research              *
 *                                                                              *
 *              This software is distributed under the terms of the             *
 *         GNU Lesser General Public Licence version 3 (LGPL) version 3,        *
 *                  copied verbatim in the file "LICENSE"                       *
 ********************************************************************************/

#include "ERRTelescopeGeoComponentSingleSi.h"

#include "TGeoManager.h"
#include "TGeoMatrix.h"
#include <TDOMParser.h>
#include <TXMLAttr.h>
#include <TXMLNode.h>
#include <TList.h>

#include "FairLogger.h"

#include "ERTelescopeSetup.h"

#include <iostream>
using namespace std;

//--------------------------------------------------------------------------------------------------
void ERRTelescopeGeoComponentSingleSi::ConstructGeometryVolume(void) {
  ParseXmlParameters();
  if (fRMax <= fRMin)
    LOG(FATAL) << "Wrong radial telescope station parameters: r_max <= r_min; r_max = " 
               << fRMax << " r_min = " << fRMin << FairLogger::endl;
  auto* media = CreateMaterial(fMedia);
  media->GetMaterial()->Print();
  const Double_t sphere_r_min = 3000.;
  // Due to a bug in the tubes, the sphere with a big radius is used instead.
  const auto make_tubes = [sphere_r_min](const TString& name, TGeoMedium* material,
                                         const Double_t r_min, const Double_t r_max, const Double_t z, 
                                         const Double_t phi_min, const Double_t phi_max) -> TGeoVolume* {
    const Double_t sphere_r_max = sphere_r_min + z;
    const Double_t sphere_theta_min = TMath::ATan(r_min/sphere_r_min)*TMath::RadToDeg();
    const Double_t sphere_theta_max = TMath::ATan(r_max/sphere_r_max)*TMath::RadToDeg();
    return gGeoManager->MakeSphere(name, material, sphere_r_min, sphere_r_max, 
                                   sphere_theta_min, sphere_theta_max, phi_min, phi_max);
  };
  const auto fullZ = fSensetiveZ + fDeadLayerThicknessFrontSide + fDeadLayerThicknessBackSide;
  fVolume = new TGeoVolumeAssembly(this->GetVolumeName());
  auto* station = make_tubes("r_station", media, fRMin, fRMax, fullZ, 0., 360.);
  fVolume->AddNode(station, 0, new TGeoCombiTrans(0, 0., -sphere_r_min, new TGeoRotation()));
  if (fOrientAroundZ == OrientationAroundZ::Y) {
    const Double_t stripDR = (fSensetiveRMax-fSensetiveRMin) / fStripCount;
    for (Int_t iStrip = 0; iStrip < fStripCount; iStrip++) {
      TGeoVolume* strip = make_tubes("Sensitivestrip", media, fSensetiveRMin + iStrip * stripDR, 
                                     fSensetiveRMin + (iStrip + 1) * stripDR, fSensetiveZ,
                                     0., 360.);
      const Double_t zTranslation = (fDeadLayerThicknessFrontSide  - fDeadLayerThicknessBackSide) / 2.;
      station->AddNode(strip, iStrip, new TGeoCombiTrans(0, 0., zTranslation, new TGeoRotation()));
    }
  } else {
    const Double_t stripPhi = 360. / fStripCount;
    TGeoVolume* strip = make_tubes("Sensitivestrip", media, fSensetiveRMin, fSensetiveRMax, 
                                   fSensetiveZ, -stripPhi/2., stripPhi/2.);
    for (Int_t iStrip = 0; iStrip < fStripCount; iStrip++) {
      auto* rotation = new TGeoRotation();
      rotation->RotateZ(stripPhi * iStrip);
      const Double_t zTranslation = (fDeadLayerThicknessFrontSide  - fDeadLayerThicknessBackSide) / 2.;
      station->AddNode(strip, iStrip, new TGeoCombiTrans(0, 0., zTranslation, rotation));
    }
  }
}
//--------------------------------------------------------------------------------------------------
void ERRTelescopeGeoComponentSingleSi::ParseXmlParameters() {
  TString xmlFile = ERTelescopeSetup::Instance()->GetXMLParametersFile();
  TDOMParser *domParser;
  domParser = new TDOMParser;
  domParser->SetValidate(false); // do not validate with DTD
  Int_t parsecode = domParser->ParseFile(xmlFile);
  if (parsecode < 0) {
     LOG(FATAL) << domParser->GetParseCodeMessage(parsecode) << FairLogger::FairLogger::endl;
    return ;
  }
  TXMLNode *rootNode = domParser->GetXMLDocument()->GetRootNode();
  TXMLNode *detPartNode = rootNode->GetChildren();
  TXMLNode *SiTypeNodes;
  for ( ; detPartNode; detPartNode = detPartNode->GetNextNode()) { // detector's part
    if (!strcasecmp(detPartNode->GetNodeName(), "RSiTypes")) {
      SiTypeNodes = detPartNode->GetChildren();
      for ( ; SiTypeNodes; SiTypeNodes = SiTypeNodes->GetNextNode()) {
        if (!strcasecmp(SiTypeNodes->GetNodeName(), "singleSiTypes")) {
          TXMLNode* curNode = SiTypeNodes->GetChildren()->GetNextNode();
          for(; curNode; curNode = curNode->GetNextNode()) {
            TList *attrList;
            TXMLAttr *attr = 0;
            if (curNode->HasAttributes()){
              attrList = curNode->GetAttributes();
              TIter next(attrList);
              while ((attr=(TXMLAttr*)next())) {
                if (!strcasecmp("id", attr->GetName())) { 
                  break;
                }
              }
            }
            else {
              continue;
            }
            if(!strcasecmp(fComponentId, attr->GetValue())) {
              FillTwoSidedChannelAttribute(curNode->GetAttributes());
              TXMLNode* curNode2 = curNode->GetChildren();
              for(; curNode2; curNode2 = curNode2->GetNextNode()) {
                if(!strcasecmp(curNode2->GetNodeName(), "size")) {
                  attrList = curNode2->GetAttributes();
                  attr = 0;
                  TIter nextSizeAttr(attrList);
                  while ((attr=(TXMLAttr*)nextSizeAttr())) {
                    if (!strcasecmp("r_min", attr->GetName())) {
                      fRMin = atof(attr->GetValue());
                    }
                    if (!strcasecmp("r_max", attr->GetName())) {
                      fRMax = atof(attr->GetValue());
                    }
                  }
                }
                if(!strcasecmp(curNode2->GetNodeName(), "sensetive_size")) {
                  attrList = curNode2->GetAttributes();
                  attr = 0;
                  TIter nextSensSizeAttr(attrList);
                  while ((attr=(TXMLAttr*)nextSensSizeAttr())) {
                    if (!strcasecmp("r_min", attr->GetName())) {
                      fSensetiveRMin = atof(attr->GetValue());
                    }
                    if (!strcasecmp("r_max", attr->GetName())) {
                      fSensetiveRMax = atof(attr->GetValue());
                    }
                    if (!strcasecmp("z", attr->GetName())) {
                      fSensetiveZ = atof(attr->GetValue());
                    }
                  }
                }
                if(!strcasecmp(curNode2->GetNodeName(), "dead_layer_front")) {
                  fDeadLayerThicknessFrontSide = atof(curNode2->GetText());
                }
                if(!strcasecmp(curNode2->GetNodeName(), "dead_layer_back")) {
                  fDeadLayerThicknessBackSide = atof(curNode2->GetText());
                }
                if(!strcasecmp(curNode2->GetNodeName(), "strip_count")) {
                  fStripCount = atof(curNode2->GetText());
                }
                if(!strcasecmp(curNode2->GetNodeName(), "media")) {
                  fMedia = curNode2->GetText();
                }
              }
            }
          }
        }
      }
    }
  }
}
//--------------------------------------------------------------------------------------------------
ClassImp(ERRTelescopeGeoComponentSingleSi)