From bed4c9327b4dcda59da80c95acb2a670b09ee491 Mon Sep 17 00:00:00 2001
From: slebedev <s.lebedev@gsi.de>
Date: Wed, 3 Aug 2022 17:17:24 +0200
Subject: [PATCH] Major update in RICH geometry. Now RICH geometry is created
completely with TGeoManager.
---
macro/rich/Import_GDML_Export_ROOT.c | 2 +-
macro/rich/create_rich_geo.C | 12 +
macro/rich/geotest/run_many.py | 12 +
macro/rich/geotest/run_one.py | 8 +-
macro/rich/run/run_many.py | 1 -
macro/rich/run/run_one.py | 19 +-
reco/detectors/rich/qa/CbmRichGeoTest.cxx | 4 +-
sim/detectors/rich/CMakeLists.txt | 3 +
sim/detectors/rich/CbmRichSimLinkDef.h | 2 +
sim/detectors/rich/geo/RichGeoCreator.cxx | 865 ++++++++++++++++++++++
sim/detectors/rich/geo/RichGeoCreator.h | 142 ++++
sim/detectors/rich/geo/RichGeoMaterials.h | 71 ++
12 files changed, 1124 insertions(+), 17 deletions(-)
create mode 100644 macro/rich/create_rich_geo.C
create mode 100644 sim/detectors/rich/geo/RichGeoCreator.cxx
create mode 100644 sim/detectors/rich/geo/RichGeoCreator.h
create mode 100644 sim/detectors/rich/geo/RichGeoMaterials.h
diff --git a/macro/rich/Import_GDML_Export_ROOT.c b/macro/rich/Import_GDML_Export_ROOT.c
index 7dcfc09372..1504c8a202 100644
--- a/macro/rich/Import_GDML_Export_ROOT.c
+++ b/macro/rich/Import_GDML_Export_ROOT.c
@@ -5,7 +5,7 @@ void addCopyNumbersToNodes(TGeoVolume* volume);
void Import_GDML_Export_ROOT()
{
- TString richGeoFilename = "rich_v21a.gdml";
+ TString richGeoFilename = "ShBox_20220628.gdml";
cout << "Importing '" << richGeoFilename << "'." << endl;
// Load FairRunSim to ensure the correct unit system
diff --git a/macro/rich/create_rich_geo.C b/macro/rich/create_rich_geo.C
new file mode 100644
index 0000000000..dc1a3fb6ca
--- /dev/null
+++ b/macro/rich/create_rich_geo.C
@@ -0,0 +1,12 @@
+/* Copyright (C) 2022-2022 UGiessen/GSI, Giessen/Darmstadt
+ SPDX-License-Identifier: GPL-3.0-only
+ Authors: Semen Lebedev [committer] */
+
+void create_rich_geo()
+{
+ RichGeoCreator* geoCreator = new RichGeoCreator();
+ geoCreator->SetGeoName("rich_v22b");
+ geoCreator->SetAddShieldingBox(false);
+ geoCreator->SetVolumeColors(false);
+ geoCreator->Create();
+}
diff --git a/macro/rich/geotest/run_many.py b/macro/rich/geotest/run_many.py
index c4840d8921..28e544fbf5 100755
--- a/macro/rich/geotest/run_many.py
+++ b/macro/rich/geotest/run_many.py
@@ -7,3 +7,15 @@ os.system(commandXterm1)
commandXterm11 = ('xterm -hold -e python3 run_one.py {} {}&').format(1, "geotest")
os.system(commandXterm11)
+commandXterm1 = ('xterm -hold -e python3 run_one.py {} {}&').format(2, "urqmdtest")
+os.system(commandXterm1)
+
+commandXterm11 = ('xterm -hold -e python3 run_one.py {} {}&').format(2, "geotest")
+os.system(commandXterm11)
+
+commandXterm1 = ('xterm -hold -e python3 run_one.py {} {}&').format(3, "urqmdtest")
+os.system(commandXterm1)
+
+commandXterm11 = ('xterm -hold -e python3 run_one.py {} {}&').format(3, "geotest")
+os.system(commandXterm11)
+
diff --git a/macro/rich/geotest/run_one.py b/macro/rich/geotest/run_one.py
index d638f9d157..982a629da6 100755
--- a/macro/rich/geotest/run_one.py
+++ b/macro/rich/geotest/run_one.py
@@ -9,9 +9,9 @@ def main():
runId = sys.argv[1]
testType = sys.argv[2] # geotest or urqmdtest
- #runId = "0"
- geoSetup = "sis100_electron"
- nEvents = 3000 if testType == "geotest" else 100
+
+ geoSetup = "sis100_electron_" + runId
+ nEvents = 5000 if testType == "geotest" else 250
targetZ = -44.0
energy = "10gev"
@@ -22,7 +22,7 @@ def main():
plutoFile = ""
dataDir = "/Users/slebedev/Development/cbm/data/sim/rich/" + testType
- resultDir = "results_" + testType + "_" + runId + "_"+ energy +"/"
+ resultDir = "results_" + testType + "_"+ energy + "_" + runId + "/"
traFile = dataDir + "/tra." + runId + ".root"
diff --git a/macro/rich/run/run_many.py b/macro/rich/run/run_many.py
index c24eae4d36..be6fb4e5c2 100755
--- a/macro/rich/run/run_many.py
+++ b/macro/rich/run/run_many.py
@@ -2,6 +2,5 @@ import os
for taskId in range(1,4):
commandXterm = ('xterm -hold -e python3 run_one.py {}&').format(taskId)
- print(commandXterm + " " + str(taskId))
os.system(commandXterm)
diff --git a/macro/rich/run/run_one.py b/macro/rich/run/run_one.py
index abfaf2ff48..2068bca28c 100755
--- a/macro/rich/run/run_one.py
+++ b/macro/rich/run/run_one.py
@@ -3,22 +3,23 @@
import os
import sys
-taskId = sys.argv[1]
+runId = sys.argv[1]
cbmrootConfigPath = "/Users/slebedev/Development/cbm/git/build/config.sh"
macroDir = "/Users/slebedev/Development/cbm/git/cbmroot/macro/rich/"
nofEvents = 100
urqmdFile = "/Users/slebedev/Development/cbm/data/urqmd/auau/10gev/centr/urqmd.auau.10gev.centr.00001.root"
dataDir= "/Users/slebedev/Development/cbm/data/sim/rich/reco/"
-traFile = dataDir + "/tra."+ taskId + ".root"
-parFile = dataDir + "/par."+ taskId + ".root"
-digiFile = dataDir + "/digi."+ taskId + ".root"
-recoFile = dataDir + "/reco."+ taskId + ".root"
-qaFile = dataDir + "/qa."+ taskId + ".root"
-geoSimFile = dataDir + "/geosim."+ taskId + ".root"
+traFile = dataDir + "/tra."+ runId + ".root"
+parFile = dataDir + "/par."+ runId + ".root"
+digiFile = dataDir + "/digi."+ runId + ".root"
+recoFile = dataDir + "/reco."+ runId + ".root"
+qaFile = dataDir + "/qa."+ runId + ".root"
+geoSimFile = dataDir + "/geosim."+ runId + ".root"
-geoSetup = "sis100_electron"
-resultDir = "results_recoqa/"
+# geoSetup = "sis100_electron"
+geoSetup = "sis100_electron_" + runId
+resultDir = "results_recoqa_" + runId + "/"
targetZ = -44.0
traCmd=('root -l -b -q {}/run/run_transport.C\(\\"{}\\",\\"{}\\",\\"{}\\",\\"{}\\",{},{},\\"{}\\",\\"{}\\",{},{}\)').format(macroDir, urqmdFile, traFile, parFile, geoSimFile, 5, 5, "", geoSetup, nofEvents, targetZ)
diff --git a/reco/detectors/rich/qa/CbmRichGeoTest.cxx b/reco/detectors/rich/qa/CbmRichGeoTest.cxx
index 98b93c400d..3e1d5d4eda 100644
--- a/reco/detectors/rich/qa/CbmRichGeoTest.cxx
+++ b/reco/detectors/rich/qa/CbmRichGeoTest.cxx
@@ -109,10 +109,10 @@ void CbmRichGeoTest::InitHistograms()
{
double xMin = -120.;
double xMax = 120.;
- int nBinsX = 300;
+ int nBinsX = 240;
double yMin = -210;
double yMax = 210.;
- int nBinsY = 105;
+ int nBinsY = 420;
fHM = new CbmHistManager();
diff --git a/sim/detectors/rich/CMakeLists.txt b/sim/detectors/rich/CMakeLists.txt
index 1fc0b0cd7c..9b758fa7bc 100644
--- a/sim/detectors/rich/CMakeLists.txt
+++ b/sim/detectors/rich/CMakeLists.txt
@@ -1,6 +1,8 @@
set(INCLUDE_DIRECTORIES
${CMAKE_CURRENT_SOURCE_DIR}
+ ${CMAKE_CURRENT_SOURCE_DIR}/geo
+
${CBMDETECTORBASE_DIR}/rich
${CBMROOT_SOURCE_DIR}/sim/transport/base
@@ -32,6 +34,7 @@ link_directories(${LINK_DIRECTORIES})
set(SRCS
CbmRich.cxx
CbmRichDigitizer.cxx
+ geo/RichGeoCreator.cxx
)
IF (SSE_FOUND)
diff --git a/sim/detectors/rich/CbmRichSimLinkDef.h b/sim/detectors/rich/CbmRichSimLinkDef.h
index 01e4238136..8bf710e814 100644
--- a/sim/detectors/rich/CbmRichSimLinkDef.h
+++ b/sim/detectors/rich/CbmRichSimLinkDef.h
@@ -13,4 +13,6 @@
#pragma link C++ class CbmDigitize < CbmRichDigi> + ;
#pragma link C++ class CbmRichDigitizer + ;
+#pragma link C++ class RichGeoCreator + ;
+
#endif
diff --git a/sim/detectors/rich/geo/RichGeoCreator.cxx b/sim/detectors/rich/geo/RichGeoCreator.cxx
new file mode 100644
index 0000000000..d283a67526
--- /dev/null
+++ b/sim/detectors/rich/geo/RichGeoCreator.cxx
@@ -0,0 +1,865 @@
+/* Copyright (C) 2022-2022 UGiessen/GSI, Giessen/Darmstadt
+ SPDX-License-Identifier: GPL-3.0-only
+ Authors: Semen Lebedev [committer] */
+
+#include "RichGeoCreator.h"
+
+#include "FairLogger.h"
+
+#include "TGeoArb8.h"
+#include "TGeoBBox.h"
+#include "TGeoBoolNode.h"
+#include "TGeoCompositeShape.h"
+#include "TGeoCone.h"
+#include "TGeoEltu.h"
+#include "TGeoMatrix.h"
+#include "TGeoSphere.h"
+#include "TGeoTube.h"
+#include "TGeoVolume.h"
+
+#include "RichGeoMaterials.h"
+
+
+using namespace std;
+
+RichGeoCreator::RichGeoCreator() {}
+
+RichGeoCreator::~RichGeoCreator() {}
+
+void RichGeoCreator::Create()
+{
+ LOG(info) << "Create geometry: " << fGeoName;
+ fMaterials.Init();
+
+ CalculateParams();
+
+ // TGeoVolume* vol = CreateMirrorSupportBelts();
+ // vol->Draw("ogl");
+ // vol->Export("rich_volume.root");
+
+ CreateRich();
+}
+
+void RichGeoCreator::CalculateParams()
+{
+ double cosMirRot = cos(ToRad(fMirrorRot));
+ double sinMirRot = sin(ToRad(fMirrorRot));
+ double mirRadCos = fMirrorRadius * cosMirRot;
+ double E = sin(ToRad(90. - fAccAngle - fMirrorRot)) / sin(ToRad(fAccAngle));
+ double A = (4. * E * E + 4. * E * sinMirRot + 1.) / (cosMirRot * cosMirRot);
+ double B = (2. * (2. * E + sinMirRot) / cosMirRot) * (1. - fMirrorPos / mirRadCos);
+ double C = fMirrorPos * (fMirrorPos - 2. * mirRadCos) / (mirRadCos * mirRadCos);
+ double t1 = (-B + sqrt(B * B - 4. * A * C)) / (2. * A);
+ fAlpha = ToDeg(asin(t1)); //rad
+ double sinAlpha = sin(ToRad(fAlpha));
+ double cosAlpha = cos(ToRad(fAlpha));
+ double psi = ToDeg(atan((fMirrorRadius * (sinAlpha * cosMirRot - sinMirRot * (1. - cosAlpha))) / fMirrorPos));
+ double mirRadHalf = fMirrorRadius / 2.;
+
+ fCameraFullRotX = psi + 2. * fMirrorRot - fCameraRot;
+ fCameraOriginY = fMirrorPos * tan(ToRad(psi)) + mirRadHalf * sin(ToRad(psi + 2. * fMirrorRot))
+ + fCameraRadius * sin(ToRad(fCameraFullRotX));
+ fCameraOriginZ =
+ fMirrorPos - mirRadHalf * cos(ToRad((psi + 2. * fMirrorRot))) - fCameraRadius * cos(ToRad(fCameraFullRotX));
+ fMirrorOriginY = fMirrorRadius * sin(ToRad(fMirrorRot + fAlpha));
+ fMirrorOriginZ = fMirrorPos - mirRadCos;
+
+ fCameraHeight = 7. * (2. * fPmtSize + fCameraGap) + 6. * fCameraGap;
+ double camRad1 = fCameraRadius - 1. - fPmtDepth;
+ fCameraTouchRadius = std::sqrt(camRad1 * camRad1 + fCameraTouchWidth * fCameraTouchWidth / 4.);
+ fCameraSegmentAng = ToDeg(2. * std::asin((fCameraTouchWidth / 2.) / fCameraTouchRadius));
+}
+
+TGeoVolume* RichGeoCreator::CreatePmt()
+{
+ double pxSizeHalf = fPixelSize / 2.;
+ double pxThicknessHalf = fPixelThickness / 2.;
+ double pxSizeEdgeHalf = fPixelSizeEdge / 2.;
+ double pmtSizeHalf = fPmtSize / 2.;
+ double pmtDepthHalf = fPmtDepth / 2.;
+ TGeoMedium* medCsI = fMaterials.GetMedium("CsI");
+
+ TGeoBBox* pixelPB = new TGeoBBox("pmt_pixel_PB", pxSizeHalf, pxSizeHalf, pxThicknessHalf);
+ TGeoBBox* pixelCornerPB = new TGeoBBox("pmt_pixel_corner_PB", pxSizeEdgeHalf, pxSizeEdgeHalf, pxThicknessHalf);
+ TGeoBBox* pixelYPB = new TGeoBBox("pmt_pixel_y_PB", pxSizeHalf, pxSizeEdgeHalf, pxThicknessHalf);
+ TGeoBBox* pixelXPB = new TGeoBBox("pmt_pixel_x_PB", pxSizeEdgeHalf, pxSizeHalf, pxThicknessHalf);
+ TGeoBBox* pmtPB = new TGeoBBox("pmt_PB", pmtSizeHalf, pmtSizeHalf, pmtDepthHalf);
+
+ TGeoVolume* pixel = new TGeoVolume("pmt_pixel", pixelPB, medCsI);
+ if (fSetVolumeColors) pixel->SetLineColor(kRed - 9);
+ TGeoVolume* pixelCorner = new TGeoVolume("pmt_pixel_corner", pixelCornerPB, medCsI);
+ if (fSetVolumeColors) pixelCorner->SetLineColor(kRed);
+ TGeoVolume* pixelY = new TGeoVolume("pmt_pixel_y", pixelYPB, medCsI);
+ if (fSetVolumeColors) pixelY->SetLineColor(kRed - 5);
+ TGeoVolume* pixelX = new TGeoVolume("pmt_pixel_x", pixelXPB, medCsI);
+ if (fSetVolumeColors) pixelX->SetLineColor(kRed + 3);
+ TGeoVolume* pmt = new TGeoVolume("pmt", pmtPB, fMaterials.GetMedium("RICHgas_CO2_dis"));
+ if (fSetVolumeColors) pmt->SetLineColor(kRed + 1);
+
+ double edgePxPos = 3. * fPixelSize + pxSizeEdgeHalf;
+ double zC = pmtDepthHalf - pxThicknessHalf;
+ int copyNum = 1;
+ for (int x = 0; x < fNofPixels; x++) {
+ double xC = (x - 3.5) * fPixelSize;
+ if (x == 0) xC = -edgePxPos;
+ if (x == 7) xC = edgePxPos;
+ for (int y = 0; y < fNofPixels; y++) {
+ double yC = (y - 3.5) * fPixelSize;
+ if (y == 0) yC = edgePxPos;
+ if (y == 7) yC = -edgePxPos;
+
+ bool isXEdge = (x == 0 || x == fNofPixels - 1);
+ bool isYEdge = (y == 0 || y == fNofPixels - 1);
+ TGeoVolume* curPmtPixel = pixel;
+ if (isXEdge) curPmtPixel = pixelX;
+ if (isYEdge) curPmtPixel = pixelY;
+ if (isXEdge && isYEdge) curPmtPixel = pixelCorner;
+
+ pmt->AddNode(curPmtPixel, copyNum, MakeTrans(xC, yC, zC));
+ copyNum++;
+ }
+ }
+
+ TGeoBBox* pmtContPB = new TGeoBBox("pmt_cont_PB", pmtSizeHalf, pmtSizeHalf, pmtDepthHalf);
+ TGeoVolume* pmtCont = new TGeoVolume("pmt_cont", pmtContPB, fMaterials.GetMedium("RICHgas_CO2_dis"));
+ pmtCont->AddNode(pmt, 1, MakeTrans(0., 0., 0.));
+ if (fSetVolumeColors) pmtCont->SetLineColor(kRed - 4);
+
+ return pmtCont;
+}
+
+
+TGeoVolume* RichGeoCreator::CreateCameraModule()
+{
+ double pmtSizeGap = fPmtSize + fCameraGap;
+
+ TGeoVolume* pmt = CreatePmt();
+ TGeoBBox* camModulePB = new TGeoBBox("camera_module_PB", (3. * fPmtSize + 2. * fCameraGap) / 2.,
+ (2. * fPmtSize + fCameraGap) / 2., fPmtDepth / 2.);
+
+ TGeoVolume* camModule = new TGeoVolume("camera_module", camModulePB, fMaterials.GetMedium("RICHgas_CO2_dis"));
+ if (fSetVolumeColors) camModule->SetLineColor(kRed - 1);
+ camModule->AddNode(pmt, 1, MakeTrans(-pmtSizeGap, pmtSizeGap / 2., 0.));
+ camModule->AddNode(pmt, 2, MakeTrans(0., pmtSizeGap / 2., 0.));
+ camModule->AddNode(pmt, 3, MakeTrans(pmtSizeGap, pmtSizeGap / 2., 0.));
+ camModule->AddNode(pmt, 4, MakeTrans(-pmtSizeGap, -pmtSizeGap / 2., 0.));
+ camModule->AddNode(pmt, 5, MakeTrans(0., -pmtSizeGap / 2., 0.));
+ camModule->AddNode(pmt, 6, MakeTrans(pmtSizeGap, -pmtSizeGap / 2., 0.));
+
+ return camModule;
+}
+
+TGeoVolume* RichGeoCreator::CreateCameraStrip()
+{
+ double pmtSizeGap = fPmtSize + fCameraGap;
+ double pmtDepthHalf = fPmtDepth / 2.;
+ double aluFrameDepth = 0.5;
+
+ TGeoVolume* camModule = CreateCameraModule();
+
+ double dY = (7. * (2. * fPmtSize + fCameraGap) + 6. * fCameraGap) / 2.;
+ TGeoBBox* aluFramePB =
+ new TGeoBBox("alu_frame_PB", fCameraTouchWidth / 2., (14. * fPmtSize + 13. * fCameraGap) / 2., aluFrameDepth);
+ TGeoBBox* op0 = new TGeoBBox("camera_strip_PB_op0", (3. * fPmtSize + 2. * fCameraGap) / 2., dY, pmtDepthHalf);
+ TGeoBBox* op1 = new TGeoBBox("camera_strip_PB_op1", fCameraTouchWidth / 2., dY, aluFrameDepth);
+ TGeoTranslation* tr = MakeTrans(0., 0., -pmtDepthHalf - aluFrameDepth);
+ TGeoCompositeShape* camStripPB = MakeUnion("camera_strip_PB", op0, op1, nullptr, tr);
+
+ TGeoVolume* aluFrame = new TGeoVolume("alu_frame", aluFramePB, fMaterials.GetMedium("aluminium"));
+ if (fSetVolumeColors) aluFrame->SetLineColor(kBlue - 4);
+ TGeoVolume* camStrip = new TGeoVolume("camera_strip", camStripPB, fMaterials.GetMedium("RICHgas_CO2_dis"));
+ if (fSetVolumeColors) camStrip->SetLineColor(kOrange + 9);
+ int copyNum = 1;
+ for (int i = -6; i <= 6; i += 2) {
+ camStrip->AddNode(camModule, copyNum, MakeTrans(0., static_cast<double>(i) * pmtSizeGap, 0.));
+ copyNum++;
+ }
+ camStrip->AddNode(aluFrame, 1, MakeTrans(0., 0., -pmtDepthHalf - aluFrameDepth));
+
+ return camStrip;
+}
+
+TGeoVolumeAssembly* RichGeoCreator::CreateCameraContainer()
+{
+ TGeoVolume* camStrip = CreateCameraStrip();
+ double len = fCameraRadius - (fPmtDepth / 2.);
+
+ TGeoVolumeAssembly* camCont = new TGeoVolumeAssembly("camera_cont");
+ int copyNum = 1;
+ for (int i = -7; i <= 6; i++) {
+ double ang = (static_cast<double>(i) + 0.5) * fCameraSegmentAng;
+ TGeoMatrix* mat = MakeCombiTrans(len * sin(ToRad(ang)), 0., len * cos(ToRad(ang)), 0., ang, 0.);
+ camCont->AddNode(camStrip, copyNum, mat);
+ copyNum++;
+ }
+ return camCont;
+}
+
+TGeoVolumeAssembly* RichGeoCreator::CreateMirror()
+{
+ double mRadius = fMirrorRadius;
+ double mRadius2 = mRadius + fMirrorThickness;
+ double dTheta = fAlpha / 2. - fMirrorThetaGap * 3. / 4.;
+ double phi1 = -fMirrorPhiSize / 2.;
+ double phi2 = fMirrorPhiSize / 2.;
+ double mRot90 = 90. + fMirrorRot;
+ TGeoMedium* medGlass = fMaterials.GetMedium("RICHglass");
+
+ // Mirror tile type 0
+ double theta0 = mRot90 + fAlpha / 2. + fMirrorThetaGap * 3. / 4.;
+ TGeoSphere* op00 = new TGeoSphere("mirror_tile_type0_PB_op0", mRadius, mRadius2, theta0, theta0 + dTheta, phi1, phi2);
+ TGeoBBox* op01 = new TGeoBBox("mirror_tile_type0_PB_op1", 500., fRichWidth / 2., 10.);
+ TGeoTranslation* tr0 = MakeTrans(0., 0., -fMirrorOriginY - 10.);
+ TGeoCompositeShape* mType0 = MakeSubtraction("mirror_tile_type0_PB", op00, op01, nullptr, tr0);
+ TGeoVolume* mVolType0 = new TGeoVolume("mirror_tile_type0", mType0, medGlass);
+ if (fSetVolumeColors) mVolType0->SetLineColor(kYellow + 1);
+
+ // Mirror tile type 1
+ double theta1 = mRot90 + fMirrorThetaGap / 2.;
+ TGeoSphere* mType1 = new TGeoSphere("mirror_tile_type1_PB", mRadius, mRadius2, theta1, theta1 + dTheta, phi1, phi2);
+ TGeoVolume* mVolType1 = new TGeoVolume("mirror_tile_type1", mType1, medGlass);
+ if (fSetVolumeColors) mVolType1->SetLineColor(kYellow + 2);
+
+ // Mirror tile type 2
+ double theta2 = mRot90 - fAlpha / 2. + fMirrorThetaGap / 4.;
+ TGeoSphere* mType2 = new TGeoSphere("mirror_tile_type2_PB", mRadius, mRadius2, theta2, theta2 + dTheta, phi1, phi2);
+ TGeoVolume* mVolType2 = new TGeoVolume("mirror_tile_type2", mType2, medGlass);
+ if (fSetVolumeColors) mVolType2->SetLineColor(kYellow + 3);
+
+ // Mirror tile type 3
+ double theta3 = mRot90 - fAlpha;
+ TGeoSphere* mType3 = new TGeoSphere("mirror_tile_type3_PB", mRadius, mRadius2, theta3, theta3 + dTheta, phi1, phi2);
+ TGeoVolume* mVolType3 = new TGeoVolume("mirror_tile_type3", mType3, medGlass);
+ if (fSetVolumeColors) mVolType3->SetLineColor(kYellow - 1);
+
+ // Mirror tile type 4
+ double theta4 = mRot90 + fAlpha / 2. + fMirrorThetaGap * 3. / 4.;
+ TGeoSphere* op4_0 =
+ new TGeoSphere("mirror_tile_type4_PB_op0", mRadius, mRadius2, theta4, theta4 + dTheta, phi1, phi2);
+ TGeoBBox* op4_1 = new TGeoBBox("mirror_tile_type4_PB_op1", 500., fRichWidth / 2., 10.);
+ TGeoTranslation* tr4_1 = MakeTrans(0., 0., -fMirrorOriginY - 10.);
+ TGeoCompositeShape* op4_01 = MakeSubtraction("mirror_tile_type4_PB_op01", op4_0, op4_1, nullptr, tr4_1);
+ TGeoTubeSeg* op4_2 = new TGeoTubeSeg("mirror_tile_type4_PB_op2", 0., 80., 320., -90., 45.);
+ TGeoMatrix* mat4 = MakeCombiTrans(0., 0., -mRadius * sin(ToRad(fMirrorRot + fAlpha)), 90., 0., 90.);
+ TGeoCompositeShape* mType4 = MakeSubtraction("mirror_tile_type4_PB", op4_01, op4_2, nullptr, mat4);
+ TGeoVolume* mVolType4 = new TGeoVolume("mirror_tile_type4", mType4, medGlass);
+ if (fSetVolumeColors) mVolType4->SetLineColor(kYellow - 2);
+
+ // Mirror tile type 5
+ double theta5 = mRot90 + fAlpha / 2. + fMirrorThetaGap * 3. / 4.;
+ TGeoSphere* op5_0 =
+ new TGeoSphere("mirror_tile_type5_PB_op0", mRadius, mRadius2, theta5, theta5 + dTheta, phi1, phi2);
+ TGeoBBox* op5_1 = new TGeoBBox("mirror_tile_type5_PB_op1", 500., fRichWidth / 2., 10.);
+ TGeoTranslation* tr5_1 = MakeTrans(0., 0., -fMirrorOriginY - 10.);
+ TGeoCompositeShape* op5_01 = MakeSubtraction("mirror_tile_type5_PB_op01", op5_0, op5_1, nullptr, tr5_1);
+ TGeoTubeSeg* op5_2 = new TGeoTubeSeg("mirror_tile_type5_PB_op2", 0., 80., 320., 135., 270.);
+ TGeoMatrix* mat5 = MakeCombiTrans(0., 0., -mRadius * sin(ToRad(fMirrorRot + fAlpha)), 90., 0., 90.);
+ TGeoCompositeShape* mType5 = MakeSubtraction("mirror_tile_type5_PB", op5_01, op5_2, nullptr, mat5);
+ TGeoVolume* mVolType5 = new TGeoVolume("mirror_tile_type5", mType5, medGlass);
+ if (fSetVolumeColors) mVolType5->SetLineColor(kYellow - 3);
+
+ TGeoVolumeAssembly* mHalf = new TGeoVolumeAssembly("mirror_half");
+
+ //TODO: This is a very strange bug, the first volume which is added to mirror_half is not properly read in simulation
+ // Workaround: create fake volume wth material RichGlass and add it first.
+ // The volume is placed outside of acceptance and has basically 0 volume
+ // You can visualize box, just change dimensions to smth like TGeoBBox(20., 20., 20.)
+ TGeoBBox* fakeBox = new TGeoBBox("fake_box", 0.0, 0.0, 0.0);
+ TGeoVolume* fakeVol = new TGeoVolume("fake_vol", fakeBox, medGlass);
+ mHalf->AddNode(fakeVol, 1, MakeTrans(280., 0., 70.));
+
+ for (int t = 0; t < 4; t++) {
+ TGeoVolume* mVol = nullptr;
+ if (t == 1) mVol = mVolType1;
+ if (t == 2) mVol = mVolType2;
+ if (t == 3) mVol = mVolType3;
+ for (int i = 0; i < 10; i++) {
+ double rotZ = (i - 4.5) * fMirrorPhiStep;
+ int copyNum = i + 1;
+ if (t == 0) {
+ mVol = mVolType0;
+ if (i == 4 || i == 5) { // special case for the 2 middle mirror tiles
+ mVol = (i == 4) ? mVolType5 : mVolType4;
+ copyNum = 1;
+ }
+ }
+ mHalf->AddNode(mVol, copyNum, MakeRot(0., 0., rotZ));
+ }
+ }
+
+ return mHalf;
+}
+
+TGeoVolumeAssembly* RichGeoCreator::CreateMirrorSupportPillars()
+{
+ TGeoMedium* medGas = fMaterials.GetMedium("RICHgas_CO2_dis");
+ TGeoMedium* medAl = fMaterials.GetMedium("aluminium");
+
+ TGeoBBox* pbGas = new TGeoBBox("mirror_support_pillar_gas_PB", 4.8, 176., 4.8);
+ TGeoBBox* pb = new TGeoBBox("mirror_support_pillar_PB", 5., 176., 5.);
+
+ TGeoVolume* volGas = new TGeoVolume("mirror_support_pillar_gas", pbGas, medGas);
+ TGeoVolume* vol = new TGeoVolume("mirror_support_pillar", pb, medAl);
+ if (fSetVolumeColors) volGas->SetLineColor(kAzure + 5);
+ if (fSetVolumeColors) vol->SetLineColor(kAzure + 5);
+ vol->AddNode(volGas, 1);
+
+ double dZ = 41.;
+ TGeoVolumeAssembly* pillars = new TGeoVolumeAssembly("mirror_support_pillar");
+ pillars->AddNode(vol, 1, MakeCombiTrans(210.7759, 0., -29.5298 + dZ, 0., 42., 0.));
+ pillars->AddNode(vol, 2, MakeCombiTrans(134.1222, 0., 21.3989 + dZ, 0., 25.2, 0.));
+ pillars->AddNode(vol, 3, MakeCombiTrans(46.015, 0., 48. + dZ, 0., 8.4, 0.));
+ pillars->AddNode(vol, 4, MakeCombiTrans(-46.015, 0., 48. + dZ, 0., -8.4, 0.));
+ pillars->AddNode(vol, 5, MakeCombiTrans(-134.1222, 0., 21.3989 + dZ, 0., -25.2, 0.));
+ pillars->AddNode(vol, 6, MakeCombiTrans(-210.7759, 0., -29.5298 + dZ, 0., -42., 0.));
+
+ return pillars;
+}
+
+TGeoVolumeAssembly* RichGeoCreator::CreateBeltAssembly()
+{
+ TGeoMedium* medGas = fMaterials.GetMedium("RICHgas_CO2_dis");
+ TGeoMedium* medAl = fMaterials.GetMedium("aluminium");
+
+ TGeoBBox* pbGas1 = new TGeoBBox("belt1_gas_PB", 1.25, 10.367, 1.167);
+ TGeoBBox* pb1 = new TGeoBBox("belt1_PB", 1.5, 10.6, 1.4);
+ TGeoVolume* volGas1 = new TGeoVolume("belt1_gas", pbGas1, medGas);
+ TGeoVolume* vol1 = new TGeoVolume("belt1", pb1, medAl);
+ if (fSetVolumeColors) volGas1->SetLineColor(kViolet + 8);
+ if (fSetVolumeColors) vol1->SetLineColor(kViolet + 8);
+ vol1->AddNode(volGas1, 1, MakeTrans(0., 0., 0.233));
+
+ TGeoBBox* pbGas2 = new TGeoBBox("belt2_gas_PB", 4.7835, 1.25, 1.167);
+ TGeoBBox* pb2 = new TGeoBBox("belt2_PB", 4.9, 1.5, 1.4);
+ TGeoVolume* volGas2 = new TGeoVolume("belt2_gas", pbGas2, medGas);
+ TGeoVolume* vol2 = new TGeoVolume("belt2", pb2, medAl);
+ if (fSetVolumeColors) volGas2->SetLineColor(kViolet + 8);
+ if (fSetVolumeColors) vol2->SetLineColor(kViolet + 8);
+ vol2->AddNode(volGas2, 1, MakeTrans(0.1165, 0., 0.233));
+
+ TGeoBBox* pbGas3 = new TGeoBBox("belt3_gas_PB", 11.6135, 1.25, 1.167);
+ TGeoBBox* pb3 = new TGeoBBox("belt3_PB", 11.73, 1.5, 1.4);
+ TGeoVolume* volGas3 = new TGeoVolume("belt3_gas", pbGas3, medGas);
+ TGeoVolume* vol3 = new TGeoVolume("belt3", pb3, medAl);
+ if (fSetVolumeColors) volGas3->SetLineColor(kViolet + 8);
+ if (fSetVolumeColors) vol3->SetLineColor(kViolet + 8);
+ vol3->AddNode(volGas3, 1, MakeTrans(0.1165, 0., 0.233));
+
+ TGeoTrap* pbGas4 = new TGeoTrap("belt4_gas_PB", 0.7, 16.22, 90., 10.97, 1.3, 1.3, 0., 9.5, 1.3, 1.3, 0.);
+ TGeoTrap* pb4 = new TGeoTrap("belt4_PB", 0.9, 16.22, 90., 11.378, 1.4, 1.4, 0., 9.5, 1.4, 1.4, 0.);
+ TGeoVolume* volGas4 = new TGeoVolume("belt4_gas", pbGas4, medGas);
+ TGeoVolume* vol4 = new TGeoVolume("belt4", pb4, medAl);
+ if (fSetVolumeColors) volGas4->SetLineColor(kViolet + 8);
+ if (fSetVolumeColors) vol4->SetLineColor(kViolet + 8);
+ vol4->AddNode(volGas4, 1, MakeTrans(0.1, 0., 0.));
+
+ TGeoVolumeAssembly* belt = new TGeoVolumeAssembly("belt");
+ belt->AddNode(vol1, 1);
+ belt->AddNode(vol2, 1, MakeTrans(-6.4, 9.1, 0.));
+ belt->AddNode(vol3, 1, MakeTrans(-13.23, -9.1, 0.));
+ belt->AddNode(vol4, 1, MakeCombiTrans(-17.567, 0.75, 0., -127., -90., -90.));
+
+ return belt;
+}
+
+TGeoVolumeAssembly* RichGeoCreator::CreateMirrorSupportBelts()
+{
+ double mRadius2 = fMirrorRadius + fMirrorThickness;
+ double mRadius3 = mRadius2 + fMirrorSupportAdd;
+ double phi1 = -fMirrorPhiSize / 2.;
+ double phi2 = fMirrorPhiSize / 2.;
+ double dTheta = fMirrorThetaSize;
+ double mRot90 = 90. + fMirrorRot;
+ double richWidthHalf = fRichWidth / 2.;
+ TGeoMedium* medGas = fMaterials.GetMedium("RICHgas_CO2_dis");
+
+ TGeoVolumeAssembly* beltAssembly = CreateBeltAssembly();
+
+ // Type 0
+ double theta0 = mRot90 + fAlpha / 2. + fMirrorThetaGap * 3. / 4.;
+ TGeoSphere* sphere0 = new TGeoSphere("belt_type0_PB_op0", mRadius2, mRadius3, theta0, theta0 + dTheta, phi1, phi2);
+ TGeoBBox* box0 = new TGeoBBox("belt_type0_PB_op1", 500., richWidthHalf, 10.);
+ TGeoCompositeShape* pbType0 =
+ MakeSubtraction("belt_type0_PB", sphere0, box0, nullptr, MakeTrans(0., 0., -fMirrorOriginY - 10.));
+ TGeoVolume* volType0 = new TGeoVolume("belt_type0", pbType0, medGas);
+ volType0->AddNode(beltAssembly, 1, MakeCombiTrans(275.5, -13.23, -124.7, 65., 0., -90.));
+ if (fSetVolumeColors) volType0->SetLineColor(kGreen);
+
+ // Type 1
+ double theta1 = mRot90 + fMirrorThetaGap / 2.;
+ TGeoSphere* pbType1 = new TGeoSphere("belt_type1_PB", mRadius2, mRadius3, theta1, theta1 + dTheta, phi1, phi2);
+ TGeoVolume* volType1 = new TGeoVolume("belt_type1", pbType1, medGas);
+ volType1->AddNode(beltAssembly, 1, MakeCombiTrans(290.4, -13.23, -84.0, 73., 0., -90.));
+ if (fSetVolumeColors) volType1->SetLineColor(kGreen - 1);
+
+ // Type 2
+ double theta2 = mRot90 - fAlpha / 2. + fMirrorThetaGap / 4.;
+ TGeoSphere* pbType2 = new TGeoSphere("belt_type2_PB", mRadius2, mRadius3, theta2, theta2 + dTheta, phi1, phi2);
+ TGeoVolume* volType2 = new TGeoVolume("belt_type2", pbType2, medGas);
+ volType2->AddNode(beltAssembly, 1, MakeCombiTrans(299.4, -13.23, -41.5, 82., 0., -90.));
+ if (fSetVolumeColors) volType2->SetLineColor(kGreen - 2);
+
+ // Type 3
+ double theta3 = mRot90 - fAlpha;
+ TGeoSphere* pbType3 = new TGeoSphere("belt_type3_PB", mRadius2, mRadius3, theta3, theta3 + dTheta, phi1, phi2);
+ TGeoVolume* volType3 = new TGeoVolume("belt_type3", pbType3, medGas);
+ volType3->AddNode(beltAssembly, 1, MakeCombiTrans(302.5, -13.23, 1.97, 91., 0., -90.));
+ if (fSetVolumeColors) volType3->SetLineColor(kGreen - 3);
+
+ // Type 4
+ double theta4 = mRot90 + fAlpha / 2. + fMirrorThetaGap * 3. / 4.;
+ TGeoSphere* sphere4 = new TGeoSphere("belt_type4_PB_op0", mRadius2, mRadius3, theta4, theta4 + dTheta, phi1, phi2);
+ TGeoTubeSeg* tube4 = new TGeoTubeSeg("belt_type4_PB_op1", 0., fPipeOuterRadiusExit, 400., 0., 360.);
+ TGeoCompositeShape* sub4_1 =
+ MakeSubtraction("belt_type4_PB_op01", sphere4, tube4, nullptr,
+ MakeCombiTrans(0., 0., -fMirrorOriginY, 90., 0., 90. + fMirrorPhiStep * 0.5));
+ TGeoBBox* box4 = new TGeoBBox("belt_type4_PB_op2", 500., richWidthHalf, 10.);
+ TGeoCompositeShape* pbType4 =
+ MakeSubtraction("belt_type4_PB", sub4_1, box4, nullptr, MakeTrans(0., 0., -fMirrorOriginY - 10.));
+ TGeoVolume* volType4 = new TGeoVolume("belt_type4", pbType4, medGas);
+ volType4->AddNode(beltAssembly, 1, MakeCombiTrans(281.0, 0., -113.0, 0., -66., 0.));
+ if (fSetVolumeColors) volType4->SetLineColor(kGreen + 1);
+
+ // Type 5
+ double theta5 = mRot90 + fAlpha / 2. + fMirrorThetaGap * 3. / 4.;
+ TGeoSphere* sphere5 = new TGeoSphere("belt_type5_PB_op0", mRadius2, mRadius3, theta5, theta5 + dTheta, phi1, phi2);
+ TGeoTubeSeg* tube5 = new TGeoTubeSeg("belt_type5_PB_op1", 0., fPipeOuterRadiusExit, 400., 0., 360.);
+ TGeoCompositeShape* sub5_1 =
+ MakeSubtraction("belt_type5_PB_op01", sphere5, tube5, nullptr,
+ MakeCombiTrans(0., 0., -fMirrorOriginY, 90., 0., 90. - fMirrorPhiStep * 0.5));
+ TGeoBBox* box5 = new TGeoBBox("belt_type5_PB_op2", 500., richWidthHalf, 10.);
+ TGeoCompositeShape* pbType5 =
+ MakeSubtraction("belt_type5_PB", sub5_1, box5, nullptr, MakeTrans(0., 0., -fMirrorOriginY - 10.));
+ TGeoVolume* volType5 = new TGeoVolume("belt_type5", pbType5, medGas);
+ volType5->AddNode(beltAssembly, 1, MakeCombiTrans(281., 0., -113., 0., -114., 180.));
+ if (fSetVolumeColors) volType5->SetLineColor(kGreen + 2);
+
+ TGeoVolumeAssembly* contHalf = new TGeoVolumeAssembly("mirror_support_belts_half");
+ for (int t = 0; t < 4; t++) {
+ TGeoVolume* mVol = nullptr;
+ if (t == 1) mVol = volType1;
+ if (t == 2) mVol = volType2;
+ if (t == 3) mVol = volType3;
+ for (int i = 0; i < 10; i++) {
+ double rotZ = (i - 4.5) * fMirrorPhiStep;
+ int copyNum = i + 1;
+ if (t == 0) {
+ mVol = volType0;
+ if (i == 4 || i == 5) { // special case for the 2 middle belts
+ mVol = (i == 4) ? volType5 : volType4;
+ copyNum = 1;
+ }
+ }
+ contHalf->AddNode(mVol, copyNum, MakeRot(0., 0., rotZ));
+ }
+ }
+ return contHalf;
+}
+
+TGeoVolumeAssembly* RichGeoCreator::CreateMainFrame()
+{
+ TGeoMedium* medGas = fMaterials.GetMedium("RICHgas_CO2_dis");
+ TGeoMedium* medAlum = fMaterials.GetMedium("aluminium");
+
+ double frameSizeHalf = 5.;
+ double frameSizeHalfGas = 4.3;
+ double richHeightHalf = fRichHeight / 2.;
+ double coverThickness = fRichCoveringThickness;
+ double posY = richHeightHalf - coverThickness - frameSizeHalf;
+ double rotY = 31.991651;
+
+ TGeoVolumeAssembly* mainFrame = new TGeoVolumeAssembly("mainframe");
+
+ TGeoBBox* pbGas1 = new TGeoBBox("mainframe1_gas_PB", 139.5, frameSizeHalfGas, frameSizeHalfGas);
+ TGeoBBox* pb1 = new TGeoBBox("mainframe1_PB", 139.5, frameSizeHalf, frameSizeHalf);
+ TGeoVolume* volGas1 = new TGeoVolume("mainframe1_gas", pbGas1, medGas);
+ TGeoVolume* vol1 = new TGeoVolume("mainframe1", pb1, medAlum);
+ if (fSetVolumeColors) vol1->SetLineColor(kBlue);
+ if (fSetVolumeColors) volGas1->SetLineColor(kBlue);
+ vol1->AddNode(volGas1, 1);
+ mainFrame->AddNode(vol1, 1, MakeTrans(0, 100., -128.279439));
+ mainFrame->AddNode(vol1, 2, MakeTrans(0, -100., -128.279439));
+ mainFrame->AddNode(vol1, 3, MakeTrans(0, posY, -147.779439));
+ mainFrame->AddNode(vol1, 4, MakeTrans(0, -posY, -147.779439));
+
+ double posZ2 = -99.4875;
+ double posX2 = 134.5;
+ TGeoBBox* pbGas2 = new TGeoBBox("mainframe2_gas_PB", frameSizeHalfGas, frameSizeHalfGas, 43.2919390);
+ TGeoBBox* pb2 = new TGeoBBox("mainframe2_PB", frameSizeHalf, frameSizeHalf, 43.2919390);
+ TGeoVolume* volGas2 = new TGeoVolume("mainframe2_gas", pbGas2, medGas);
+ TGeoVolume* vol2 = new TGeoVolume("mainframe2", pb2, medAlum);
+ if (fSetVolumeColors) vol2->SetLineColor(kBlue + 1);
+ if (fSetVolumeColors) volGas2->SetLineColor(kBlue + 1);
+ vol2->AddNode(volGas2, 1);
+ mainFrame->AddNode(vol2, 1, MakeTrans(-posX2, posY, posZ2));
+ mainFrame->AddNode(vol2, 2, MakeTrans(-posX2, -posY, posZ2));
+ mainFrame->AddNode(vol2, 3, MakeTrans(posX2, posY, posZ2));
+ mainFrame->AddNode(vol2, 4, MakeTrans(posX2, -posY, posZ2));
+
+ double posX3 = 196.4811997;
+ double posZ3 = -91.7898222;
+ TGeoTrap* pbGas3 = new TGeoTrap("mainframe3_gas_PB", frameSizeHalfGas, -26.016698, 0., frameSizeHalfGas, 71.9697521,
+ 71.9697521, 0., frameSizeHalfGas, 62.4, 62.4, 0.);
+ TGeoTrap* pb3 = new TGeoTrap("mainframe3_PB", frameSizeHalf, -26.016698, 0., frameSizeHalf, 72.7486842, 72.7486842,
+ 0., frameSizeHalf, 61.6210827, 61.6210827, 0.);
+ TGeoVolume* volGas3 = new TGeoVolume("mainframe3_gas", pbGas3, medGas);
+ TGeoVolume* vol3 = new TGeoVolume("mainframe3", pb3, medAlum);
+ if (fSetVolumeColors) vol3->SetLineColor(kBlue + 2);
+ if (fSetVolumeColors) volGas3->SetLineColor(kBlue + 2);
+ vol3->AddNode(volGas3, 1);
+ mainFrame->AddNode(vol3, 1, MakeCombiTrans(posX3, posY, posZ3, 0., -rotY, 0.));
+ mainFrame->AddNode(vol3, 2, MakeCombiTrans(posX3, -posY, posZ3, 0., -rotY, 0.));
+ mainFrame->AddNode(vol3, 3, MakeCombiTrans(-posX3, posY, posZ3, 0., -rotY, -180.0));
+ mainFrame->AddNode(vol3, 4, MakeCombiTrans(-posX3, -posY, posZ3, 0., -rotY, -180.0));
+
+ // long vertical parts
+ double posZ4 = 51.195561;
+ double pos4 = 257.9;
+ TGeoBBox* pbGas4 = new TGeoBBox("mainframe4_gas_PB", frameSizeHalfGas,
+ richHeightHalf - coverThickness - 2. * frameSizeHalf, frameSizeHalfGas);
+ TGeoBBox* pb4 =
+ new TGeoBBox("mainframe4_PB", frameSizeHalf, richHeightHalf - coverThickness - 2. * frameSizeHalf, frameSizeHalf);
+ TGeoVolume* volGas4 = new TGeoVolume("mainframe4_gas", pbGas4, medGas);
+ TGeoVolume* vol4 = new TGeoVolume("mainframe4", pb4, medAlum);
+ if (fSetVolumeColors) vol4->SetLineColor(kBlue + 3);
+ if (fSetVolumeColors) volGas4->SetLineColor(kBlue + 3);
+ vol4->AddNode(volGas4, 1);
+ mainFrame->AddNode(vol4, 1, MakeTrans(-pos4, 0., posZ4));
+ mainFrame->AddNode(vol4, 2, MakeTrans(pos4, 0., posZ4));
+ mainFrame->AddNode(vol4, 3, MakeTrans(-pos4, 0, -posZ4));
+ mainFrame->AddNode(vol4, 4, MakeTrans(pos4, 0, -posZ4));
+ mainFrame->AddNode(vol4, 5, MakeCombiTrans(215.1316481, 0, -80.1395071, 0., -rotY, 0.));
+ mainFrame->AddNode(vol4, 6, MakeCombiTrans(178.111648, 0, -103.2646704, 0., -rotY, 0.));
+ mainFrame->AddNode(vol4, 7, MakeCombiTrans(146.3896048, 0, -123.0803763, 0., -rotY, 0.));
+ mainFrame->AddNode(vol4, 8, MakeCombiTrans(-215.1316481, 0, -80.1395071, 0., -rotY, 0.));
+ mainFrame->AddNode(vol4, 9, MakeCombiTrans(-178.111648, 0, -103.2646704, 0., -rotY, 0.));
+ mainFrame->AddNode(vol4, 10, MakeCombiTrans(-146.3896048, 0, -123.0803763, 0., -rotY, 0.));
+
+ TGeoBBox* pbGas5 = new TGeoBBox("mainframe5_gas_PB", 262.9, frameSizeHalfGas, frameSizeHalfGas);
+ TGeoBBox* pb5 = new TGeoBBox("mainframe5_PB", 262.9, frameSizeHalf, frameSizeHalf);
+ TGeoVolume* volGas5 = new TGeoVolume("mainframe5_gas", pbGas5, medGas);
+ TGeoVolume* vol5 = new TGeoVolume("mainframe5", pb5, medAlum);
+ if (fSetVolumeColors) vol5->SetLineColor(kBlue - 1);
+ if (fSetVolumeColors) volGas5->SetLineColor(kBlue - 1);
+ vol5->AddNode(volGas5, 1);
+ mainFrame->AddNode(vol5, 1, MakeTrans(0., posY, posZ4));
+ mainFrame->AddNode(vol5, 2, MakeTrans(0., posY, -posZ4));
+ mainFrame->AddNode(vol5, 3, MakeTrans(0., -posY, posZ4));
+ mainFrame->AddNode(vol5, 4, MakeTrans(0., -posY, -posZ4));
+
+ // small horizontal parts
+ double lenHalf6 = 46.1955610;
+ TGeoBBox* pbGas6 = new TGeoBBox("mainframe6_gas_PB", frameSizeHalfGas, frameSizeHalfGas, lenHalf6);
+ TGeoBBox* pb6 = new TGeoBBox("mainframe6_PB", frameSizeHalf, frameSizeHalf, lenHalf6);
+ TGeoVolume* volGas6 = new TGeoVolume("mainframe6_gas", pbGas6, medGas);
+ TGeoVolume* vol6 = new TGeoVolume("mainframe6", pb6, medAlum);
+ if (fSetVolumeColors) vol6->SetLineColor(kBlue - 2);
+ if (fSetVolumeColors) volGas6->SetLineColor(kBlue - 2);
+ vol6->AddNode(volGas6, 1);
+ mainFrame->AddNode(vol6, 1, MakeTrans(-pos4, posY, 0.));
+ mainFrame->AddNode(vol6, 2, MakeTrans(-pos4, -posY, 0.));
+ mainFrame->AddNode(vol6, 3, MakeTrans(pos4, posY, 0.));
+ mainFrame->AddNode(vol6, 4, MakeTrans(pos4, -posY, 0.));
+ mainFrame->AddNode(vol6, 5, MakeTrans(-134.5, posY, 0.));
+ mainFrame->AddNode(vol6, 6, MakeTrans(-134.5, -posY, 0.));
+ mainFrame->AddNode(vol6, 7, MakeTrans(134.5, posY, 0.));
+ mainFrame->AddNode(vol6, 8, MakeTrans(134.5, -posY, 0.));
+ mainFrame->AddNode(vol6, 9, MakeTrans(0., posY, 0.));
+ mainFrame->AddNode(vol6, 10, MakeTrans(0., -posY, 0.));
+ mainFrame->AddNode(vol6, 11, MakeTrans(pos4, 87.8, 0.));
+ mainFrame->AddNode(vol6, 12, MakeTrans(pos4, -87.8, 0.));
+ mainFrame->AddNode(vol6, 13, MakeTrans(-pos4, 87.8, 0.));
+ mainFrame->AddNode(vol6, 14, MakeTrans(-pos4, -87.8, 0.));
+ mainFrame->AddNode(vol6, 15, MakeTrans(-192.6, posY, 0.));
+ mainFrame->AddNode(vol6, 16, MakeTrans(-192.6, -posY, 0.));
+ mainFrame->AddNode(vol6, 17, MakeTrans(192.6, posY, 0.));
+ mainFrame->AddNode(vol6, 18, MakeTrans(192.6, -posY, 0.));
+
+ return mainFrame;
+}
+
+TGeoVolume* RichGeoCreator::CreateGas()
+{
+ double richLenHalf = fRichLength / 2.;
+ double richWidthHalf = fRichWidth / 2.;
+ double richHeightHalf = fRichHeight / 2.;
+
+ double x1 = richWidthHalf - fRichCoveringThickness;
+ double y1 = ((2. * richHeightHalf - fRichEntranceY) / 2. - 2. * fRichCoveringThickness) / 2.;
+ double z1 = fRichAddZ / 2. - fRichCoveringThickness / 2. + fRichEntranceThickness / 2.;
+ double y2 = (richHeightHalf + fRichEntranceY / 2.) / 2.;
+ double z2 = -richLenHalf - fRichAddZ / 2. + fRichCoveringThickness;
+
+ TGeoBBox* op0 = new TGeoBBox("rich_gas_PB_op0", x1, richHeightHalf - fRichCoveringThickness,
+ richLenHalf - fRichEntranceThickness / 2. - fRichExitThickness / 2.);
+ TGeoBBox* op1 = new TGeoBBox("rich_gas_PB_op1", x1, y1, z1);
+ TGeoCompositeShape* op01 = MakeUnion("rich_gas_PB_op01", op0, op1, nullptr, MakeTrans(0., y2, z2));
+ TGeoBBox* op2 = new TGeoBBox("rich_gas_PB_op2", x1, y1, z1);
+ TGeoCompositeShape* op012 = MakeUnion("rich_gas_PB_op012", op01, op2, nullptr, MakeTrans(0., -y2, z2));
+
+ TGeoConeSeg* op3 =
+ new TGeoConeSeg("rich_gas_PB_op3", richLenHalf, 0., fPipeOuterRadiusEntrance, 0., fPipeOuterRadiusExit, 0., 360.);
+ TGeoTranslation* tr = MakeTrans(0., 0., -(fRichEntranceThickness - fRichExitThickness) / 2.);
+ TGeoCompositeShape* gasPB = MakeSubtraction("rich_gas_PB", op012, op3, nullptr, tr);
+ TGeoVolume* gas = new TGeoVolume("rich_gas", gasPB, fMaterials.GetMedium("RICHgas_CO2_dis"));
+ if (fSetVolumeColors) gas->SetLineColor(kYellow - 10);
+ return gas;
+}
+
+TGeoVolume* RichGeoCreator::CreateRichContainer()
+{
+ double y1 = (fRichHeight / 2. + fRichEntranceY / 2.) / 2.;
+ double z1 = -fRichLength / 2. - fRichAddZ / 2.;
+ double y2 = (fRichHeight - fRichEntranceY) / 4.;
+ TGeoBBox* op0 = new TGeoBBox("rich_cont_PB_op0", fRichWidth / 2., fRichHeight / 2., fRichLength / 2.);
+ TGeoBBox* op1 = new TGeoBBox("rich_cont_PB_op1", fRichWidth / 2., y2, fRichAddZ / 2.);
+ TGeoCompositeShape* op01 = MakeUnion("rich_cont_PB_op01", op0, op1, nullptr, MakeTrans(0., y1, z1));
+ TGeoBBox* op2 = new TGeoBBox("rich_cont_PB_op2", fRichWidth / 2., y2, fRichAddZ / 2.);
+ TGeoCompositeShape* pb = MakeUnion("rich_cont_PB", op01, op2, nullptr, MakeTrans(0., -y1, z1));
+ double richLenHalf = fRichLength / 2;
+ TGeoConeSeg* pipePB =
+ new TGeoConeSeg("rich_pipe_subtract", richLenHalf, 0, fPipeOuterRadiusEntrance - fPipeThicknessEntrance, 0,
+ fPipeOuterRadiusExit - fPipeThicknessExit, 0., 360.);
+ TGeoCompositeShape* pb2 = MakeSubtraction("rich_cont_PB2", pb, pipePB, nullptr, MakeTrans(0, 0, 0));
+ TGeoVolume* richCont = new TGeoVolume("rich_cont", pb2, fMaterials.GetMedium("aluminium"));
+ if (fSetVolumeColors) richCont->SetLineColor(kBlue - 10);
+ return richCont;
+}
+
+TGeoVolume* RichGeoCreator::CreateSensitivePlane()
+{
+ double richLenHalf = fRichLength / 2.;
+ TGeoBBox* op0 = new TGeoBBox("sens_plane_PB_op0", 170., 170., 0.05);
+ TGeoConeSeg* op1 =
+ new TGeoConeSeg("sens_plane_PB_op1", richLenHalf, 0., fPipeOuterRadiusEntrance, 0., fPipeOuterRadiusExit, 0., 360.);
+ TGeoTranslation* tr = MakeTrans(0., 0., (-fSensPlaneZ + (fRichExitThickness - fRichEntranceThickness) / 2.));
+ TGeoCompositeShape* sensPlanePB = MakeSubtraction("sens_plane_PB", op0, op1, nullptr, tr);
+ TGeoVolume* sensPlane = new TGeoVolume("sens_plane", sensPlanePB, fMaterials.GetMedium("RICHgas_CO2_dis+"));
+ if (fSetVolumeColors) sensPlane->SetLineColor(kRed - 8);
+ return sensPlane;
+}
+
+TGeoVolume* RichGeoCreator::CreateRichEntrance()
+{
+ double richLenHalf = fRichLength / 2.;
+ double richWidthHalf = fRichWidth / 2.;
+ TGeoBBox* op0 = new TGeoBBox("rich_entrance_PB_op0", richWidthHalf, fRichEntranceY / 2., fRichEntranceThickness / 2.);
+ TGeoConeSeg* op1 = new TGeoConeSeg("rich_entrance_PB_op1", richLenHalf, 0., fPipeOuterRadiusEntrance, 0.,
+ fPipeOuterRadiusExit, 0., 360.);
+ TGeoTranslation* tr = MakeTrans(0., 0., richLenHalf - fRichEntranceThickness / 2.);
+ TGeoCompositeShape* pb = MakeSubtraction("rich_entrance_PB", op0, op1, nullptr, tr);
+ TGeoVolume* entrance = new TGeoVolume("rich_entrance", pb, fMaterials.GetMedium("kapton"));
+ if (fSetVolumeColors) entrance->SetLineColor(kGreen - 5);
+ return entrance;
+}
+
+TGeoVolume* RichGeoCreator::CreateRichExit()
+{
+ double richHeightHalf = fRichHeight / 2.;
+ double richWidthHalf = fRichWidth / 2.;
+ double richLenHalf = fRichLength / 2.;
+
+ TGeoBBox* op0 = new TGeoBBox("rich_exit_PB_op0", richWidthHalf, richHeightHalf, fRichExitThickness / 2.);
+ TGeoConeSeg* op1 =
+ new TGeoConeSeg("rich_exit_PB_op1", richLenHalf, 0., fPipeOuterRadiusEntrance, 0., fPipeOuterRadiusExit, 0., 360.);
+ TGeoTranslation* tr = MakeTrans(0., 0., -richLenHalf + fRichExitThickness / 2.);
+ TGeoCompositeShape* pb = MakeSubtraction("rich_exit_PB", op0, op1, nullptr, tr);
+ TGeoVolume* exit = new TGeoVolume("rich_exit", pb, fMaterials.GetMedium("kapton"));
+ if (fSetVolumeColors) exit->SetLineColor(kGreen - 5);
+ return exit;
+}
+
+TGeoVolume* RichGeoCreator::CreatePipe()
+{
+ double richLenHalf = fRichLength / 2;
+ TGeoConeSeg* pipePB = new TGeoConeSeg("rich_pipe_PB", richLenHalf, fPipeOuterRadiusEntrance - fPipeThicknessEntrance,
+ fPipeOuterRadiusEntrance, fPipeOuterRadiusExit - fPipeThicknessExit,
+ fPipeOuterRadiusExit, 0., 360.);
+ TGeoVolume* pipe = new TGeoVolume("rich_pipe", pipePB, fMaterials.GetMedium("aluminium"));
+ if (fSetVolumeColors) pipe->SetLineColor(kSpring + 4);
+ return pipe;
+}
+
+TGeoVolume* RichGeoCreator::CreatePipeCyl()
+{
+ TGeoTubeSeg* pipeCylVacPB = new TGeoTubeSeg(
+ "rich_cyl_pipe_vac_PB", 0., fPipeOuterRadiusEntrance - fPipeThicknessEntrance, fPipeCylPartLen / 2., 0., 360.);
+ TGeoTubeSeg* pipeCylPB =
+ new TGeoTubeSeg("rich_cyl_pipe_PB", 0., fPipeOuterRadiusEntrance, fPipeCylPartLen / 2., 0., 360.);
+ TGeoVolume* pipeCylVac = new TGeoVolume("rich_cyl_pipe_vac", pipeCylVacPB, fMaterials.GetMedium("vacuum"));
+ TGeoVolume* pipeCyl = new TGeoVolume("rich_cyl_pipe", pipeCylPB, fMaterials.GetMedium("aluminium"));
+ if (fSetVolumeColors) pipeCylVac->SetLineColor(kSpring + 5);
+ if (fSetVolumeColors) pipeCyl->SetLineColor(kSpring + 4);
+ pipeCyl->AddNode(pipeCylVac, 1);
+ return pipeCyl;
+}
+
+void RichGeoCreator::CreateRich()
+{
+ double richLenHalf = fRichLength / 2.;
+
+ TGeoVolume* rich_gas = CreateGas();
+
+ TGeoVolumeAssembly* camCont = CreateCameraContainer();
+ double camY = fCameraOriginY + fCameraShiftY;
+ double camZ = -fRichOrigZ - richLenHalf + fCameraOriginZ + fCameraShiftZ;
+ rich_gas->AddNode(camCont, 1, MakeCombiTrans(0., camY, camZ, fCameraFullRotX, 0., 0.));
+ rich_gas->AddNode(camCont, 2, MakeCombiTrans(0., -camY, camZ, -fCameraFullRotX, 0., 0.));
+
+ double mirrorZ = -fRichOrigZ - richLenHalf + fMirrorOriginZ;
+ double beltsGapZ = 0.5; // cm
+
+ TGeoVolumeAssembly* mirrorHalf = CreateMirror();
+ rich_gas->AddNode(mirrorHalf, 1, MakeCombiTrans(0., fMirrorOriginY + fMirrorGapY / 2., mirrorZ, 0., -90., -90.));
+ rich_gas->AddNode(mirrorHalf, 2, MakeCombiTrans(0., -fMirrorOriginY - fMirrorGapY / 2., mirrorZ, 0., -90., 90.));
+
+ TGeoVolumeAssembly* belts_half_container = CreateMirrorSupportBelts();
+
+ rich_gas->AddNode(belts_half_container, 1,
+ MakeCombiTrans(0., fMirrorOriginY + fMirrorGapY, mirrorZ + beltsGapZ, 0., -90., -90.));
+ rich_gas->AddNode(belts_half_container, 2,
+ MakeCombiTrans(0., -fMirrorOriginY - fMirrorGapY, mirrorZ + beltsGapZ, 0., -90., 90.));
+
+ TGeoVolumeAssembly* pillars = CreateMirrorSupportPillars();
+ rich_gas->AddNode(pillars, 1);
+
+ TGeoVolume* sensPlane = CreateSensitivePlane();
+ rich_gas->AddNode(sensPlane, 1, MakeTrans(0., 0., -30.));
+
+ TGeoVolumeAssembly* mainFrame = CreateMainFrame();
+ rich_gas->AddNode(mainFrame, 1, MakeTrans(0., 0., 38.541939));
+
+ if (fAddShieldingBox) {
+ TGeoVolumeAssembly* schieldingBox = CreateSchieldingBox();
+ double schieldingBoxZ = -(fRichEntranceThickness - fRichExitThickness / 2. + fRichOrigZ + richLenHalf);
+ rich_gas->AddNode(schieldingBox, 1, MakeTrans(0., 0., schieldingBoxZ));
+ rich_gas->AddNode(schieldingBox, 2, MakeCombiTrans(0., 0., schieldingBoxZ, 0., 0., 180.));
+ }
+
+ TGeoVolume* richCont = CreateRichContainer();
+ richCont->AddNode(rich_gas, 1, MakeTrans(0., 0., (fRichEntranceThickness - fRichExitThickness) / 2.));
+ richCont->AddNode(CreateRichEntrance(), 1, MakeTrans(0., 0., -richLenHalf + fRichEntranceThickness / 2.));
+ richCont->AddNode(CreateRichExit(), 1, MakeTrans(0., 0., richLenHalf - fRichExitThickness / 2.));
+ richCont->AddNode(CreatePipe(), 1);
+
+ TGeoVolumeAssembly* rich = new TGeoVolumeAssembly(fGeoName.c_str());
+ rich->AddNode(richCont, 1, MakeTrans(0., 0., -40. + fRichOrigZ + richLenHalf));
+ // This pipe segment is removed to be compatible with new beam pipe.
+ // See https://git.cbm.gsi.de/CbmSoft/cbmroot_geometry/-/merge_requests/184
+ //rich->AddNode(CreatePipeCyl(), 1, MakeTrans(0., 0., -40. + fRichOrigZ - fPipeCylPartLen / 2.));
+
+ TGeoVolume* top = new TGeoVolumeAssembly(fGeoName.c_str());
+
+ gGeoManager->SetTopVolume(top);
+ top->AddNode(rich, 1);
+ gGeoManager->CloseGeometry();
+ gGeoManager->CheckOverlaps(0.0001);
+ gGeoManager->PrintOverlaps();
+ gGeoManager->Test();
+
+ string outputFile = fGeoName + ".geo.root";
+ top->Export(outputFile.c_str());
+ LOG(info) << "Geometry file: " << outputFile;
+
+ top->Draw("ogl");
+ gGeoManager->SetVisLevel(7);
+}
+
+TGeoVolumeAssembly* RichGeoCreator::CreateSchieldingBox()
+{
+ // When converting GDML to root please look https://root.cern.ch/root/html520/src/TGDMLParse.cxx.htm
+ // Shielding box is created using a number of boolean operations with shapes
+ // All next shape is marked with increased number in name box1 to box 15
+
+ TGeoEltu* box1 = new TGeoEltu("box1", 183.0, 183.0, 61.3);
+ TGeoBBox* box2 = new TGeoBBox("box2", 229.3 / 2., 73.879 / 2., 113.3 / 2.);
+ TGeoCompositeShape* box3 = MakeIntersection("box3", box1, box2, nullptr, MakeTrans(0., 146.068, -4.65));
+
+ TGeoBBox* box4 = new TGeoBBox("box4", 227.3 / 2., 71.937 / 2., 102.0 / 2.);
+ TGeoCompositeShape* box5 = MakeSubtraction("box5", box3, box4, nullptr, MakeTrans(0., 147.097, 0.));
+
+ TGeoEltu* box6 = new TGeoEltu("box6", 172.0, 172.0, 4.65);
+ TGeoCompositeShape* box7 = MakeSubtraction("box7", box5, box6, nullptr, MakeTrans(0., 0., -56.65));
+
+ TGeoBBox* box8 = new TGeoBBox("box8", 75.0 / 2., 222.422 / 2., 104.0 / 2.);
+ TGeoCompositeShape* box9 = MakeSubtraction("box9", box7, box8, nullptr, MakeTrans(0., -0.082, 0.));
+
+ TGeoBBox* box10 = new TGeoBBox("box10", 145.0 / 2., 10.3 / 2., 11.0 / 2.);
+ TGeoCompositeShape* box11 =
+ MakeSubtraction("box11", box9, box10, nullptr, MakeCombiTrans(61.1525, 177.2736, -56.15, 0., 0., -31.26911));
+
+ TGeoBBox* box12 = new TGeoBBox("box12", 145.0 / 2., 10.3 / 2., 11.0 / 2.);
+ TGeoCompositeShape* box13 =
+ MakeSubtraction("box13", box11, box12, nullptr, MakeCombiTrans(-61.152, 177.274, -56.15, 0., 0., 31.27));
+
+ TGeoBBox* box14 = new TGeoBBox("box14", 145.0 / 2., 10.3 / 2., 16.0 / 2.);
+ TGeoCompositeShape* box15 =
+ MakeSubtraction("box15", box13, box14, nullptr, MakeCombiTrans(0., 181.2639, -59.75, -41.3565, 0., 0.));
+
+ TGeoVolume* shBox = new TGeoVolume("shielding_box", box15, fMaterials.GetMedium("iron"));
+ if (fSetVolumeColors) shBox->SetLineColor(kBlue + 4);
+
+ // Wing part
+ TGeoBBox* wing1 = new TGeoBBox("wing1", 117.0 / 2., 8.9571 / 2., 13.7221 / 2.);
+ TGeoBBox* wing2 = new TGeoBBox("wing2", 117.0 / 2., 8.9571 / 2., 13.7221 / 2.);
+ TGeoCompositeShape* wing3 =
+ MakeSubtraction("wing3", wing1, wing2, nullptr, MakeCombiTrans(49.9989, 0., 3.206, 0., 18.90578, 0.));
+
+ TGeoBBox* wing4 = new TGeoBBox("wing4", 117.0 / 2., 8.9571 / 2., 13.7221 / 2.);
+ TGeoCompositeShape* wing5 =
+ MakeSubtraction("wing5", wing3, wing4, nullptr, MakeCombiTrans(-49.9989, 0., 3.206, 0., -18.90578, 0.));
+
+ TGeoVolume* wing = new TGeoVolume("shieling_wing", wing5, fMaterials.GetMedium("iron"));
+ if (fSetVolumeColors) wing->SetLineColor(kBlue - 1);
+
+ TGeoVolumeAssembly* shBoxFinal = new TGeoVolumeAssembly("shielding_box");
+ shBoxFinal->AddNode(shBox, 1, MakeCombiTrans(0., 240.7246, 79.4162, 68.06551, 0., -180.0));
+ shBoxFinal->AddNode(wing, 1, MakeCombiTrans(0., 118.3074, 225.1717, -26.70901, 0., 0.));
+
+ return shBoxFinal;
+}
+
+TGeoCompositeShape* RichGeoCreator::MakeSubtraction(const string& name, TGeoShape* left, TGeoShape* right,
+ TGeoMatrix* lmat, TGeoMatrix* rmat)
+{
+ return new TGeoCompositeShape(name.c_str(), new TGeoSubtraction(left, right, lmat, rmat));
+}
+
+TGeoCompositeShape* RichGeoCreator::MakeUnion(const string& name, TGeoShape* left, TGeoShape* right, TGeoMatrix* lmat,
+ TGeoMatrix* rmat)
+{
+ return new TGeoCompositeShape(name.c_str(), new TGeoUnion(left, right, lmat, rmat));
+}
+
+TGeoCompositeShape* RichGeoCreator::MakeIntersection(const string& name, TGeoShape* left, TGeoShape* right,
+ TGeoMatrix* lmat, TGeoMatrix* rmat)
+{
+ return new TGeoCompositeShape(name.c_str(), new TGeoIntersection(left, right, lmat, rmat));
+}
+
+
+TGeoCombiTrans* RichGeoCreator::MakeCombiTrans(double dx, double dy, double dz, double rotX, double rotY, double rotZ)
+{
+ // Note that rotX, rotY, rotZ should be provided as in GDML and in deg.
+ // The rotation is done with oposite sign of rot(XYZ), first - Z, second - Y, third - X
+ // Creation of TGeoCombiTrans is done with rot->Inverse()
+ TGeoTranslation* tr = MakeTrans(dx, dy, dz);
+ TGeoRotation* rot = new TGeoRotation();
+ if (rotZ != 0.) rot->RotateZ(-rotZ);
+ if (rotY != 0.) rot->RotateY(-rotY);
+ if (rotX != 0.) rot->RotateX(-rotX);
+ return new TGeoCombiTrans(*tr, rot->Inverse());
+}
+
+TGeoTranslation* RichGeoCreator::MakeTrans(double dx, double dy, double dz) { return new TGeoTranslation(dx, dy, dz); }
+
+TGeoRotation* RichGeoCreator::MakeRot(double rotX, double rotY, double rotZ)
+{
+ TGeoRotation* rot = new TGeoRotation();
+ if (rotZ != 0.) rot->RotateZ(-rotZ);
+ if (rotY != 0.) rot->RotateY(-rotY);
+ if (rotX != 0.) rot->RotateX(-rotX);
+ return rot;
+}
+
+ClassImp(RichGeoCreator)
diff --git a/sim/detectors/rich/geo/RichGeoCreator.h b/sim/detectors/rich/geo/RichGeoCreator.h
new file mode 100644
index 0000000000..fa99f838e1
--- /dev/null
+++ b/sim/detectors/rich/geo/RichGeoCreator.h
@@ -0,0 +1,142 @@
+/* Copyright (C) 2022-2022 UGiessen/GSI, Giessen/Darmstadt
+ SPDX-License-Identifier: GPL-3.0-only
+ Authors: Semen Lebedev [committer] */
+
+#ifndef RICH_GEO_CREATOR
+#define RICH_GEO_CREATOR
+
+#include "TObject.h"
+
+#include <string>
+
+#include <cmath>
+
+#include "RichGeoMaterials.h"
+
+class TGeoVolumeAssembly;
+class TGeoVolume;
+class TGeoShape;
+class TGeoMatrix;
+class TGeoCompositeShape;
+class TGeoRotation;
+
+class RichGeoCreator : public TObject {
+public:
+ RichGeoCreator();
+
+ virtual ~RichGeoCreator();
+
+ void Create();
+
+ void SetGeoName(const std::string& geoName) { fGeoName = geoName; }
+ void SetAddShieldingBox(bool addSchildingBox) { fAddShieldingBox = addSchildingBox; }
+ void SetVolumeColors(bool setVolumeColors) { fSetVolumeColors = setVolumeColors; }
+
+private:
+ RichGeoMaterials fMaterials;
+
+ std::string fGeoName = "rich";
+ bool fAddShieldingBox = true;
+ bool fSetVolumeColors = false;
+
+ // General RICH parameters
+ double fDegToRad = 1.74532925199433E-02;
+ double fRichOrigZ = 180.; //cm
+ double fAccAngle = 25.; //deg
+ double fRichLength = 190.; //cm
+ double fRichHeight = 500.; //cm
+ double fRichWidth = 526.8; //cm
+ double fRichEntranceY = 180.; //cm
+ double fRichCoveringThickness = 0.5; //cm
+ double fRichEntranceThickness = 0.025; //cm
+ double fRichExitThickness = 0.5; //cm
+ double fRichAddZ = 20.; //cm
+
+ // Pipe parameters
+ double fPipeOuterRadiusEntrance = 7.42; //cm
+ double fPipeThicknessEntrance = 0.148; //cm
+ double fPipeOuterRadiusExit = 16.16; //cm
+ double fPipeThicknessExit = 0.32; //cm
+ double fPipeCylPartLen = 10.0; //cm
+
+ // Mirror parameters
+ double fMirrorRadius = 300.; //cm
+ double fMirrorThickness = 0.6; //cm
+ double fMirrorPos = 350.; //cm
+ double fMirrorRot = 12.; //deg
+ double fMirrorPhiSize = 8.425; //deg
+ double fMirrorPhiStep = 8.5; //deg
+ double fMirrorThetaSize = 8.088; //deg
+ double fMirrorThetaStep = 8.1625; //deg
+ double fMirrorThetaGap = 0.0745; //deg
+ double fMirrorSupportAdd = 4.5; //cm
+ double fMirrorGapY = 0.1; // cm
+
+ // PMT parameters
+ int fNofPixels = 8;
+ double fPixelSize = 0.6; //cm
+ double fPixelSizeEdge = 0.625; //cm
+ double fPixelThickness = 0.05; //cm
+ double fPmtSize = 5.2; //cm
+ double fPmtDepth = 3.871; //cm
+ double fCameraRot = 15.; //deg
+ double fCameraShiftY = 2.5; //cm
+ double fCameraShiftZ = 5.; //cm
+ double fCameraRadius = 169.7; //cm
+ double fCameraGap = 0.1; //cm
+ double fCameraTouchWidth = 15.75; //cm
+
+ // Sensitive plane parameters
+ double fSensPlaneZ = -70.; //cm
+
+ // Calculated parameters. These parameters are calculated in CalculateParams()
+ double fCameraFullRotX = 0.;
+ double fCameraOriginY = 0.;
+ double fCameraOriginZ = 0.;
+ double fMirrorOriginY = 0.;
+ double fMirrorOriginZ = 0.;
+ double fAlpha = 0.;
+ double fCameraHeight = 0.;
+ double fCameraTouchRadius = 0.;
+ double fCameraSegmentAng = 0.;
+
+ void CalculateParams();
+
+ TGeoVolume* CreatePmt();
+ TGeoVolume* CreateCameraModule();
+ TGeoVolume* CreateCameraStrip();
+ TGeoVolumeAssembly* CreateCameraContainer();
+ TGeoVolumeAssembly* CreateMirror();
+ TGeoVolumeAssembly* CreateMainFrame();
+ TGeoVolume* CreateGas();
+ TGeoVolume* CreateRichContainer();
+ TGeoVolume* CreateSensitivePlane();
+ TGeoVolume* CreateRichEntrance();
+ TGeoVolume* CreateRichExit();
+ TGeoVolume* CreatePipe();
+ TGeoVolume* CreatePipeCyl();
+ TGeoVolumeAssembly* CreateSchieldingBox();
+ TGeoVolumeAssembly* CreateBeltAssembly();
+ TGeoVolumeAssembly* CreateMirrorSupportBelts();
+ TGeoVolumeAssembly* CreateMirrorSupportPillars();
+
+ void CreateRich();
+
+ double ToDeg(double rad) { return rad / fDegToRad; }
+ double ToRad(double deg) { return deg * fDegToRad; }
+
+ TGeoCompositeShape* MakeSubtraction(const std::string& name, TGeoShape* left, TGeoShape* right,
+ TGeoMatrix* lmat = nullptr, TGeoMatrix* rmat = nullptr);
+ TGeoCompositeShape* MakeUnion(const std::string& name, TGeoShape* left, TGeoShape* right, TGeoMatrix* lmat = nullptr,
+ TGeoMatrix* rmat = nullptr);
+ TGeoCompositeShape* MakeIntersection(const std::string& name, TGeoShape* left, TGeoShape* right,
+ TGeoMatrix* lmat = nullptr, TGeoMatrix* rmat = nullptr);
+
+ TGeoCombiTrans* MakeCombiTrans(double dx, double dy, double dz, double rotX, double rotY, double rotZ);
+ TGeoTranslation* MakeTrans(double dx, double dy, double dz);
+ TGeoRotation* MakeRot(double rotX, double rotY, double rotZ);
+
+ ClassDef(RichGeoCreator, 1)
+};
+
+#endif
diff --git a/sim/detectors/rich/geo/RichGeoMaterials.h b/sim/detectors/rich/geo/RichGeoMaterials.h
new file mode 100644
index 0000000000..25f91a32bf
--- /dev/null
+++ b/sim/detectors/rich/geo/RichGeoMaterials.h
@@ -0,0 +1,71 @@
+/* Copyright (C) 2022-2022 UGiessen/GSI, Giessen/Darmstadt
+ SPDX-License-Identifier: GPL-3.0-only
+ Authors: Semen Lebedev [committer] */
+
+#ifndef RICH_GEO_MATERIALS
+#define RICH_GEO_MATERIALS
+
+#include "FairGeoBuilder.h"
+#include "FairGeoInterface.h"
+#include "FairGeoLoader.h"
+#include "FairGeoMedia.h"
+#include "FairGeoMedium.h"
+#include "FairLogger.h"
+
+#include "TGeoManager.h"
+#include "TGeoMedium.h"
+#include "TSystem.h"
+
+#include <map>
+#include <string>
+
+
+class RichGeoMaterials {
+public:
+ void Init()
+ {
+ LOG(info) << "RichGeoMaterials init materials";
+
+ FairGeoLoader* geoLoad = new FairGeoLoader("TGeo", "FairGeoLoader");
+ FairGeoInterface* geoFace = geoLoad->getGeoInterface();
+ std::string mediaPath = std::string(gSystem->Getenv("VMCWORKDIR")) + "/geometry/media.geo";
+ geoFace->setMediaFile(mediaPath.c_str());
+ geoFace->readMedia();
+
+ fGeoMedia = geoFace->getMedia();
+ fGeoBuild = geoLoad->getGeoBuilder();
+
+ fMediums["aluminium"] = InitMedium("aluminium");
+ fMediums["CsI"] = InitMedium("CsI");
+ fMediums["RICHgas_CO2_dis"] = InitMedium("RICHgas_CO2_dis");
+ fMediums["vacuum"] = InitMedium("vacuum");
+ fMediums["RICHglass"] = InitMedium("RICHglass");
+ fMediums["kapton"] = InitMedium("kapton");
+ fMediums["iron"] = InitMedium("iron");
+ fMediums["RICHgas_CO2_dis+"] = InitMedium("RICHgas_CO2_dis+");
+ }
+
+ TGeoMedium* GetMedium(const std::string& name)
+ {
+ if (fMediums.find(name) == fMediums.end()) LOG(fatal) << "RichGeoMaterials::GetMedium " << name << " not found";
+ return fMediums[name];
+ }
+
+private:
+ std::map<std::string, TGeoMedium*> fMediums;
+
+ FairGeoMedia* fGeoMedia = nullptr;
+ FairGeoBuilder* fGeoBuild = nullptr;
+
+ TGeoMedium* InitMedium(const std::string& name)
+ {
+ FairGeoMedium* fairMedium = fGeoMedia->getMedium(name.c_str());
+ if (fairMedium == nullptr) LOG(fatal) << "RichGeoMaterials FairMedium " << name << " not found";
+ fGeoBuild->createMedium(fairMedium);
+ TGeoMedium* rootMedium = gGeoManager->GetMedium(name.c_str());
+ if (rootMedium == nullptr) LOG(fatal) << "RichGeoMaterials Medium " << name << " not found";
+ return rootMedium;
+ }
+};
+
+#endif
--
GitLab