From d6be89bc7913aec54866c3ea1a00efc14f7d5fb2 Mon Sep 17 00:00:00 2001
From: Mehulkumar Shiroya <m.shiroya@gsi.de>
Date: Fri, 19 Nov 2021 16:20:11 +0100
Subject: [PATCH] STS geometries version v21f and v21g
v21f: STS geometry is similar to the version v21b. The target position considered from -4 cm.
The geometry can be used with the beampipe with diameter 40 mm upstream side and 104 mm downstream side.
The geometry has beampipe flange on the back-side of the wall.
v21g: STS geometry is similar to the version v21e. The target position is at -40 cm considering global origin
from the center of magnet. The STS local coordinate is from the last station.
The geometry can be used with the beampipe with diameter 40 mm upstream side and 104 mm downstream side.
In both geometry, the STS box and the beam pipe flange is made of Carbon fiber: The material has 70% carbon fiber with
density 1.75 g/cm3 and epoxy resin (30%) with density 1.2 g/cm3. The material is defined in media.geo file as
STS_CarbonFiber_21.
---
macro/sts/geometry/create_stsgeo_v21f.C | 2188 ++++++++++++++++++++++
macro/sts/geometry/create_stsgeo_v21g.C | 2197 +++++++++++++++++++++++
2 files changed, 4385 insertions(+)
create mode 100644 macro/sts/geometry/create_stsgeo_v21f.C
create mode 100644 macro/sts/geometry/create_stsgeo_v21g.C
diff --git a/macro/sts/geometry/create_stsgeo_v21f.C b/macro/sts/geometry/create_stsgeo_v21f.C
new file mode 100644
index 0000000000..b97c1a89af
--- /dev/null
+++ b/macro/sts/geometry/create_stsgeo_v21f.C
@@ -0,0 +1,2188 @@
+/* Copyright (C) 2012-2021 Goethe-Universitaet Frankfurt, Frankfurt
+ SPDX-License-Identifier: GPL-3.0-only
+ Authors: Mehul Shiroya [committer], Volker Friese, Evgeny Lavrik*/
+
+/******************************************************************************
+ ** Creation of STS geometry in ROOT format (TGeo).
+ **
+ ** @file create_stsgeo_v21f.C
+ ** @author Volker Friese <v.friese@gsi.de>
+ ** @since 15 June 2012
+ ** @date 09.05.2014
+ ** @author Tomas Balog <T.Balog@gsi.de>
+ **
+ ** v21f: Similar to the STS geometry version v21b
+ * [ Airex foam in the front and back side wall (Solid).
+ * New beam pipe R = 2.0 cm upstream-side with thickness 1 mm.
+ * Central ladders adjusted 5mm away from the beampipe.
+ * Nomenclature for the Halfladders has been changed from the upside(u) to top (t)
+ * and from downside(d) to bottom (b). ]
+ * The only difference is that it includes beampipe flange in the back side of the wall,
+ * which is considered to be made of carbon fiber(70%) and epoxy resin(30%).
+ *
+ ** v19r: bugfix of v19q - align all halfladders
+ ** v19q: based on v19l - align ladders to virtual plane in station center, closing the gaps in z
+ ** v19p: based on v19k - parameters : delta Z prime = 0.70 cm - delta Z pitch = 0.20 cm
+ ** v19o: based on v19k - parameters : delta Z prime = 0.30 cm - delta Z pitch = 0.20 cm
+ ** v19n: based on v19k - parameters : delta Z prime = 0.50 cm - delta Z pitch = 0.20 cm (bug fix of v19m)
+ ** v19m: based on v19k - parameters : delta Z prime = 0.55 cm - delta Z pitch = 0.20 cm (bug)
+ ** v19l: based on v19k - parameters : delta Z prime = 0.50 cm - delta Z pitch = 0.15 cm
+ ** v19k: ladders on upstream side of units get upper half ladders installed first,
+ ** ladders on downstream side of units get lower half ladders installed first,
+ ** this saves 1.5 mm space in z per station, 12 mm in total (LadderType went from 3 to 4 digits)
+ ** parameters : delta Z prime = 1.00 cm - delta Z pitch = 0.15 cm
+ ** v19j: use overlap and distance parameters from CAD model
+ ** v19h: put STS stations from v19d at z-positions = 260; 365; 470; 575; 680; 785; 890; 995 mm
+ ** v19g: place a box with services around v19e
+ ** v19f: place a box with services around v19d
+ ** v19e: increase spacing between stations by +10 mm from 100 mm
+ ** v19d: increase spacing between stations by + 5 mm from 100 mm
+ ** v19c: drop station 8 and increase spacing between remaining 7 stations from 10 cm to 12 c
+ ** v19b: introduce FEB orientation in ladder numbering (LadderType went from 2 to 3 digits)
+ ** v19a: import passive materials from gdml file
+ ** extend CF ladder structures and cables towards FEE plane
+ ** change CF ladder frame shape
+ ** v18d: increases thickness of sensors within a 7.5 degree cone from 300 mu to 400 mu (based on v18b)
+ ** v18c: fixed cut-out windows in cooling plates, improve the box shape/materials
+ ** v18b: increases thickness of sensors within a 7.5 degree cone from 300 mu to 400 mu
+ ** v18a: adds 9 cooling/holding plates and a box around the setup
+ ** v16g: v16g is the new standard geometry from November 2017
+ ** v16g: switch from stations to units - left / right ("Unit01L", "Unit01R")
+ ** v16f: switch from stations to units
+ ** - split in upstream / downstream and left / right parts
+ ** - named Unit0xUR, Unit0xUL, Unit0xDR, Unit0xDL
+ ** v16e: switch from stations to units - upstream / downstream ("Unit01U", "Unit01D")
+ ** v16d: skip keeping volumes of sts and stations
+ ** v16c: like v16b, but senors of ladders beampipe next to beampipe
+ ** shifted closer to the pipe, like in the CAD model
+ ** v16b: like v16a, but yellow sensors removed
+ ** v16a: derived from v15c (no cones), but with sensor types renamed:
+ ** 2 -> 1, 3 -> 2, 4 -> 3, 5 -> 4, 1 -> 5
+ **
+ ** v15c: as v15b without cones
+ ** v15b: introduce modified carbon ladders from v13z
+ ** v15a: with flipped ladder orientation for stations 0,2,4,6 to match CAD design
+ **
+ ** TODO:
+ **
+ ** DONE:
+ ** v15b - use carbon macaroni as ladder support
+ ** v15b - introduce a small gap between lowest sensor and carbon ladder
+ ** v15b - build small cones for the first 2 stations
+ ** v15b - within a station the ladders of adjacent units should not touch eachother - set gkLadderGapZ to 10 mm
+ ** v15b - for all ladders set an even number of ladder elements
+ ** v15b - z offset of cones to ladders should not be 0.3 by default, but 0.26
+ ** v15b - within a station the ladders should be aligned in z, defined either by the unit or the ladder with most sensors
+ ** v15b - get rid of cone overlap in stations 7 and 8 - done by adapting rHole size
+ **
+ ** The geometry hierarachy is:
+ **
+ ** 1. Sensors (see function CreateSensors)
+ ** The sensors are the active volumes and the lowest geometry level.
+ ** They are built as TGeoVolumes, shape box, material silicon.
+ ** x size is determined by strip pitch 58 mu and 1024 strips
+ ** plus guard ring of 1.3 mm at each border -> 6.1992 cm.
+ ** Sensor type 1 is half of that (3.0792 cm).
+ ** y size is determined by strip length (2.2 / 4.2 / 6.3 cm) plus
+ ** guard ring of 1.3 mm at top and bottom -> 2.46 / 4.46 / 6.46 cm.
+ ** z size is a parameter, to be set by gkSensorThickness.
+ **
+ ** 2. Sectors (see function CreateSectors)
+ ** Sectors consist of several chained sensors. These are arranged
+ ** vertically on top of each other with a gap to be set by
+ ** gkChainGapY. Sectors are constructed as TGeoVolumeAssembly.
+ ** The sectors are auxiliary volumes used for proper placement
+ ** of the sensor(s) in the module. They do not show up in the
+ ** final geometry.
+ **
+ ** 3. Modules (see function ConstructModule)
+ ** A module is a readout unit, consisting of one sensor or
+ ** a chain of sensors (see sector) and a cable.
+ ** The cable extends from the top of the sector vertically to the
+ ** top of the halfladder the module is placed in. The cable and module
+ ** volume thus depend on the vertical position of the sector in
+ ** the halfladder. The cables consist of silicon with a thickness to be
+ ** set by gkCableThickness.
+ ** Modules are constructed as TGeoVolume, shape box, medium gStsMedium.
+ ** The module construction can be switched off (gkConstructCables)
+ ** to reproduce older geometries.
+ **
+ ** 4. Halfladders (see function ConstructHalfLadder)
+ ** A halfladder is a vertical assembly of several modules. The modules
+ ** are placed vertically such that their sectors overlap by
+ ** gkSectorOverlapY. They are displaced in z direction to allow for the
+ ** overlap in y by gkSectorGapZ.
+ ** The horizontal placement of modules in the halfladder can be choosen
+ ** to left aligned or right aligned, which only matters if sensors of
+ ** different x size are involved.
+ ** Halfladders are constructed as TGeoVolumeAssembly.
+ **
+ ** 5. Ladders (see function CreateLadders and ConstructLadder)
+ ** A ladder is a vertical assembly of two halfladders, and is such the
+ ** vertical building block of a station. The second (bottom) half ladder
+ ** is rotated upside down. The vertical arrangement is such that the
+ ** inner sectors of the two halfladders have the overlap gkSectorOverlapY
+ ** (function CreateLadder) or that there is a vertical gap for the beam
+ ** hole (function CreateLadderWithGap).
+ ** Ladders are constructed as TGeoVolumeAssembly.
+ **
+ ** 6. Stations (see function ConstructStation)
+ ** A station represents one layer of the STS geometry: one measurement
+ ** at (approximately) a given z position. It consist of several ladders
+ ** arranged horizontally to cover the acceptance.
+ ** The ladders are arranged such that there is a horizontal overlap
+ ** between neighbouring ladders (gkLadderOverLapX) and a vertical gap
+ ** to allow for this overlap (gkLadderGapZ). Each second ladder is
+ ** rotated around its y axis to face away from or into the beam.
+ ** Stations are constructed as TGeoVolumes, shape box minus tube (for
+ ** the beam hole), material gStsMedium.
+ **
+ ** 7. STS
+ ** The STS is a volume hosting the entire detectors system. It consists
+ ** of several stations located at different z positions.
+ ** The STS is constructed as TGeoVolume, shape box minus cone (for the
+ ** beam pipe), material gStsMedium. The size of the box is computed to
+ ** enclose all stations.
+ *****************************************************************************/
+
+
+// Remark: With the proper steering variables, this should exactly reproduce
+// the geometry version v11b of A. Kotynia's described in the ASCII format.
+// The only exception is a minimal difference in the z position of the
+// sectors/sensors. This is because of ladder types 2 and 4 containing the half
+// sensors around the beam hole (stations 1,2 and 3). In v11b, the two ladders
+// covering the beam hole cannot be transformed into each other by rotations,
+// but only by a reflection. This means they are constructionally different.
+// To avoid introducing another two ladder types, the difference in z position
+// was accepted.
+
+
+// Differences to v12:
+// gkChainGap reduced from 1 mm to 0
+// gkCableThickness increased from 100 mum to 200 mum (2 cables per module)
+// gkSectorOverlapY reduced from 3 mm to 2.4 mm
+// New sensor types 05 and 06
+// New sector types 07 and 08
+// Re-definiton of ladders (17 types instead of 8)
+// Re-definiton of station from new ladders
+
+
+#include "TGeoCompositeShape.h"
+#include "TGeoCone.h"
+#include "TGeoManager.h"
+#include "TGeoPara.h"
+#include "TGeoPhysicalNode.h"
+#include "TGeoTrd2.h"
+#include "TGeoTube.h"
+#include "TGeoXtru.h"
+
+#include <iomanip>
+#include <iostream>
+
+// forward declarations
+Int_t CreateSensors();
+Int_t CreateSectors();
+Int_t CreateLadders();
+TGeoVolume* ConstructModule(const char* name, TGeoVolume* sector, Double_t cableLength);
+TGeoVolume* ConstructHalfLadder(const TString& name, Int_t nSectors, Int_t* sectorTypes, char align,
+ Double_t ladderLength, Double_t offsetY);
+TGeoVolume* ConstructLadder(Int_t LadderIndex, TGeoVolume* halfLadderU, TGeoVolume* halfLadderD, Double_t gapY,
+ Double_t shiftZ, Double_t pitchZ, Int_t nSectors);
+TGeoVolume* ConstructUnit(Int_t iSide, Int_t iUnit, Int_t nLadders, Int_t* ladderTypes, Int_t iStation);
+void ImportPassive(TGeoVolume* stsVolume, TString geoTag, fstream& infoFile);
+void PostProcessGdml(TGeoVolume* gdmlTop);
+void CheckVolume(TGeoVolume* volume);
+void CheckVolume(TGeoVolume* volume, fstream& file, Bool_t listChildren = kTRUE);
+Double_t BeamPipeRadius(Double_t z);
+TGeoVolume* ConstructFrameElement(const TString& name, TGeoVolume* frameBoxVol, Double_t x);
+TGeoVolume* ConstructSmallCone(Double_t coneDz);
+TGeoVolume* ConstructBigCone(Double_t coneDz);
+
+// ------------- Version highlight -----------------------------------
+
+const std::string gVersionHighlight = R"(
+Summary:
+ This version adds passive materials imported from GDML model to the STS geometry:
+ * Taken from and largely correspond to mechanical CAD drawings of the detector
+ * Thermal insulation box:
+ - made out of carbon sandwitch panel (2mm carbon fiber sheet + layer of carbon foam + 2mm carbon fiber sheet)
+ - front window of complex shape with interface to MVD / target chamber
+ - back window with large aperture (2000 x 1200 mm) square cut into carbon foam
+ * Structural units:
+ - made of 2 complex shape aluminum C-Frames, 15mm thick
+ - placed at 25, 35, ... ,105 cm absolute Z
+ - contain front-end and power distribution boxes with equivalent X_0 values
+
+ Scripted geometry tweaks:
+ * Ladders and cables are extended towards the read-out planes having same lengths in respective rows
+ * Adjusted form and shape of carbon ladder structures from L-type to X-type
+ * Reduced verbosity of this file
+
+ Sensor arrangement is the same as in version v16g
+
+ !! Important for this version is the discrepancy from the mechanical CAD w.r.t. front wall.
+ The square window was replaced by a round one to avoid overlaps with present beam pipe designs, e.g. pipe_v16b_1e
+)";
+
+// ------------- Steering variables -----------------------------------
+
+// ---> Horizontal width of sensors [cm]
+const Double_t gkSensorSizeX = 6.2; // was 6.2092; // 6.2 - Oleg CAD 15/05/2020
+
+// ---> Thickness of sensors [cm]
+const Double_t gkSensorThickness = 0.03;
+
+// ---> Vertical gap between chained sensors [cm]
+const Double_t gkChainGapY = 0.00;
+
+// ---> Thickness of cables [cm]
+const Double_t gkCableThickness = 0.02;
+
+// ---> Horizontal overlap of neighbouring ladders [cm]
+const Double_t gkLadderOverlapX = 0.25; // delta X - Oleg CAD 14/05/2020
+
+// ---> Vertical overlap of neighbouring sectors in a ladder [cm]
+const Double_t gkSectorOverlapY = 0.46; // delta Y - Oleg CAD 14/05/2020
+
+// ---> Gap in z between neighbouring sectors in a ladder [cm]
+const Double_t gkSectorGapZ = 0.12; // gap + thickness = pitch // delta Z pitch = 0.15 - Oleg CAD 14/05/2020
+
+// ---> Gap in z between neighbouring ladders [cm]
+const Double_t gkLadderGapZ = 0.50 - 0.15; // for asym // 0.5 for sym // delta Z prime
+
+// ---> Gap in z between lowest sector to carbon support structure [cm]
+const Double_t gkSectorGapZFrame =
+ 0.280
+ - 0.025; // Oleg CAD 05/05/2020 // there is a 2.8 mm gap between the bottom side of the sensor and the top ledge of the carbon ladder
+
+// ---> Switch to construct / not to construct readout cables
+const Bool_t gkConstructCables = kTRUE;
+
+// ---> Switch to construct / not to construct frames
+const Bool_t gkConstructCones = kFALSE; // kTRUE; // switch this false by default for v15c and v16x
+const Bool_t gkConstructFrames = kTRUE; // kFALSE; // switch this true by default for v15c and v16x
+const Bool_t gkConstructSmallFrames = kTRUE; // kFALSE;
+const Bool_t gkCylindricalFrames = kTRUE; // kFALSE;
+
+// ---> Size of the frame
+const Double_t gkFrameThickness = 0.2;
+const Double_t gkThinFrameThickness = 0.05;
+const Double_t gkFrameStep = 4.0; // size of frame cell along y direction
+const Double_t gkCylinderDiaInner =
+ 0.07; // properties of cylindrical carbon supports, see CBM-STS Integration Meeting (10 Jul 2015)
+const Double_t gkCylinderDiaOuter =
+ 0.15; // properties of cylindrical carbon supports, see CBM-STS Integration Meeting (10 Jul 2015)
+
+// ---> Switch to import / not to import the Passive materials from GDML file
+const Bool_t gkImportPassive = kTRUE;
+
+// ----------------------------------------------------------------------------
+
+
+// -------------- Parameters of beam pipe in the STS region --------------
+// ---> Needed to compute stations and STS such as to avoid overlaps
+const Double_t gkPipeZ1 = 18.0;
+const Double_t gkPipeR1 = 2.0;
+const Double_t gkPipeZ2 = 125.0;
+const Double_t gkPipeR2 = 5.0;
+//const Double_t gkPipeZ3 = 125.0;
+//const Double_t gkPipeR3 = 5.5;
+
+//const Double_t gkPipeZ2 = 50.0;
+//const Double_t gkPipeR2 = 1.8;
+
+//DE const Double_t gkPipeZ1 = 27.0;
+//DE const Double_t gkPipeR1 = 1.05;
+//DE const Double_t gkPipeZ2 = 160.0;
+//DE const Double_t gkPipeR2 = 3.25;
+// ----------------------------------------------------------------------------
+
+//TString unitName[16] = // names of units for v16e
+// { "Unit00D",
+// "Unit01U", "Unit01D",
+// "Unit02U", "Unit02D",
+// "Unit03U", "Unit03D",
+// "Unit04U", "Unit04D",
+// "Unit05U", "Unit05D",
+// "Unit06U", "Unit06D",
+// "Unit07U", "Unit07D",
+// "Unit08U" };
+
+//TString unitName[32] = // names of units for v16f
+// { "Unit00DR", "Unit00DL",
+// "Unit01UR", "Unit01UL", "Unit01DR", "Unit01DL",
+// "Unit02UR", "Unit02UL", "Unit02DR", "Unit02DL",
+// "Unit03UR", "Unit03UL", "Unit03DR", "Unit03DL",
+// "Unit04UR", "Unit04UL", "Unit04DR", "Unit04DL",
+// "Unit05UR", "Unit05UL", "Unit05DR", "Unit05DL",
+// "Unit06UR", "Unit06UL", "Unit06DR", "Unit06DL",
+// "Unit07UR", "Unit07UL", "Unit07DR", "Unit07DL",
+// "Unit08UR", "Unit08UL" };
+
+TString unitName[32] = // names of units for v16g - while merging D and U parts
+ {"Unit00R", "Unit00L", "Unit01R", "Unit01L", "Unit01R", "Unit01L", "Unit02R", "Unit02L",
+ "Unit02R", "Unit02L", "Unit03R", "Unit03L", "Unit03R", "Unit03L", "Unit04R", "Unit04L",
+ "Unit04R", "Unit04L", "Unit05R", "Unit05L", "Unit05R", "Unit05L", "Unit06R", "Unit06L",
+ "Unit06R", "Unit06L", "Unit07R", "Unit07L", "Unit07R", "Unit07L", "Unit08R", "Unit08L"};
+
+TString unitName18[18] = // names of units for v16g
+ {"Unit00R", "Unit00L", "Unit01R", "Unit01L", "Unit02R", "Unit02L", "Unit03R", "Unit03L", "Unit04R",
+ "Unit04L", "Unit05R", "Unit05L", "Unit06R", "Unit06L", "Unit07R", "Unit07L", "Unit08R", "Unit08L"};
+
+// ------------- Other global variables -----------------------------------
+// ---> STS medium (for every volume except silicon)
+TGeoMedium* gStsMedium = NULL; // will be set later
+// ---> TGeoManager (too lazy to write out 'Manager' all the time
+TGeoManager* gGeoMan = NULL; // will be set later
+// ----------------------------------------------------------------------------
+
+
+// ============================================================================
+// ====== Main function =====
+// ============================================================================
+
+void create_stsgeo_v21f(const char* geoTag = "v21f")
+{
+
+ // ------- Geometry file name (output) ----------------------------------
+ TString geoFileName = "sts_";
+ geoFileName = geoFileName + geoTag + ".geo.root";
+ // --------------------------------------------------------------------------
+
+
+ // ------- Open info file -----------------------------------------------
+ TString infoFileName = geoFileName;
+ infoFileName.ReplaceAll("root", "info");
+ fstream infoFile;
+ infoFile.open(infoFileName.Data(), fstream::out);
+ infoFile << "STS geometry created with create_stsgeo_v21f.C" << endl;
+ infoFile << gVersionHighlight << endl;
+ infoFile << "Global variables: " << endl;
+ infoFile << "Sensor thickness = " << gkSensorThickness << " cm" << endl;
+ infoFile << "Vertical gap in sensor chain = " << gkChainGapY << " cm" << endl;
+ infoFile << "Vertical overlap of sensors = " << gkSectorOverlapY << " cm" << endl;
+ infoFile << "Gap in z between neighbour sensors = " << gkSectorGapZ << " cm" << endl;
+ infoFile << "Horizontal overlap of sensors = " << gkLadderOverlapX << " cm" << endl;
+ infoFile << "Gap in z between neighbour ladders = " << gkLadderGapZ << " cm" << endl;
+ if (gkConstructCables) infoFile << "Cable thickness = " << gkCableThickness << " cm" << endl;
+ else
+ infoFile << "No cables" << endl;
+ infoFile << endl;
+ infoFile << "Beam pipe: R1 = " << gkPipeR1 << " cm at z = " << gkPipeZ1 << " cm" << endl;
+ infoFile << "Beam pipe: R2 = " << gkPipeR2 << " cm at z = " << gkPipeZ2 << " cm" << endl;
+ //infoFile << "Beam pipe: R3 = " << gkPipeR3 << " cm at z = " << gkPipeZ3
+ // << " cm" << endl;
+ // --------------------------------------------------------------------------
+
+
+ // ------- Load media from media file -----------------------------------
+ FairGeoLoader* geoLoad = new FairGeoLoader("TGeo", "FairGeoLoader");
+ FairGeoInterface* geoFace = geoLoad->getGeoInterface();
+ TString geoPath = gSystem->Getenv("VMCWORKDIR");
+ TString medFile = geoPath + "/geometry/media.geo";
+ geoFace->setMediaFile(medFile);
+ geoFace->readMedia();
+ gGeoMan = gGeoManager;
+ gSystem->Load("libGeom");
+ // --------------------------------------------------------------------------
+
+
+ // ----------------- Get and create the required media -----------------
+ FairGeoMedia* geoMedia = geoFace->getMedia();
+ FairGeoBuilder* geoBuild = geoLoad->getGeoBuilder();
+
+ // ---> air
+ FairGeoMedium* mAir = geoMedia->getMedium("air");
+ if (!mAir) Fatal("Main", "FairMedium air not found");
+ geoBuild->createMedium(mAir);
+ TGeoMedium* air = gGeoMan->GetMedium("air");
+ if (!air) Fatal("Main", "Medium air not found");
+
+ // ---> air
+ FairGeoMedium* mvacuum = geoMedia->getMedium("vacuum");
+ if (!mvacuum) Fatal("Main", "FairMedium vacuum not found");
+ geoBuild->createMedium(mvacuum);
+ TGeoMedium* Vacuum = gGeoMan->GetMedium("vacuum");
+ if (!Vacuum) Fatal("Main", "Medium vacuum not found");
+
+
+ // ---> silicon
+ FairGeoMedium* mSilicon = geoMedia->getMedium("silicon");
+ if (!mSilicon) Fatal("Main", "FairMedium silicon not found");
+ geoBuild->createMedium(mSilicon);
+ TGeoMedium* silicon = gGeoMan->GetMedium("silicon");
+ if (!silicon) Fatal("Main", "Medium silicon not found");
+
+ // ---> carbon
+ FairGeoMedium* mCarbon = geoMedia->getMedium("carbon");
+ if (!mCarbon) Fatal("Main", "FairMedium carbon not found");
+ geoBuild->createMedium(mCarbon);
+ TGeoMedium* carbon = gGeoMan->GetMedium("carbon");
+ if (!carbon) Fatal("Main", "Medium carbon not found");
+
+ // ---> STSBoxCarbonFoam
+ FairGeoMedium* mSTSBoxCarbonFoam = geoMedia->getMedium("STSBoxCarbonFoam");
+ if (!mSTSBoxCarbonFoam) Fatal("Main", "FairMedium STSBoxCarbonFoam not found");
+ geoBuild->createMedium(mSTSBoxCarbonFoam);
+ TGeoMedium* STSBoxCarbonFoam = gGeoMan->GetMedium("STSBoxCarbonFoam");
+ if (!STSBoxCarbonFoam) Fatal("Main", "Medium STSBoxCarbonFoam not found");
+
+ // ---> STSBoxCarbonFibre
+ FairGeoMedium* mSTSBoxCarbonFibre = geoMedia->getMedium("STSBoxCarbonFibre");
+ if (!mSTSBoxCarbonFibre) Fatal("Main", "FairMedium STSBoxCarbonFibre not found");
+ geoBuild->createMedium(mSTSBoxCarbonFibre);
+ TGeoMedium* STSBoxCarbonFibre = gGeoMan->GetMedium("STSBoxCarbonFibre");
+ if (!STSBoxCarbonFibre) Fatal("Main", "Medium STSBoxCarbonFibre not found");
+
+ // ---> STScable
+ FairGeoMedium* mSTScable = geoMedia->getMedium("STScable");
+ if (!mSTScable) Fatal("Main", "FairMedium STScable not found");
+ geoBuild->createMedium(mSTScable);
+ TGeoMedium* STScable = gGeoMan->GetMedium("STScable");
+ if (!STScable) Fatal("Main", "Medium STScable not found");
+
+ // ---> Aluminium
+ FairGeoMedium* mAluminium = geoMedia->getMedium("aluminium");
+ if (!mAluminium) Fatal("Main", "FairMedium aluminium not found");
+ geoBuild->createMedium(mAluminium);
+ TGeoMedium* aluminium = gGeoMan->GetMedium("aluminium");
+ if (!aluminium) Fatal("Main", "Medium aluminium not found");
+
+ // ---
+ //gStsMedium = air;
+ gStsMedium = gGeoMan->GetMedium("air");
+ if (!gStsMedium) Fatal("Main", "medium sts_air not found");
+ // --------------------------------------------------------------------------
+ // --------------------------------------------------------------------------
+
+
+ // -------------- Create geometry and top volume -------------------------
+ gGeoMan = (TGeoManager*) gROOT->FindObject("FAIRGeom");
+ //gGeoMan->SetName("STSgeom");
+ TGeoVolume* top = new TGeoVolumeAssembly("top");
+ // TGeoBBox* topbox= new TGeoBBox("", 120., 120., 120.);
+ // TGeoVolume* top = new TGeoVolume("top", topbox, gGeoMan->GetMedium("air"));
+ gGeoMan->SetTopVolume(top);
+
+ // --------------------------------------------------------------------------
+
+
+ // -------------- Create media ------------------------------------------
+ /*
+ cout << endl;
+ cout << "===> Creating media....";
+ cout << CreateMedia();
+ cout << " media created" << endl;
+ TList* media = gGeoMan->GetListOfMedia();
+ for (Int_t iMedium = 0; iMedium < media->GetSize(); iMedium++ ) {
+ cout << "Medium " << iMedium << ": "
+ << ((TGeoMedium*) media->At(iMedium))->GetName() << endl;
+ }
+
+ gStsMedium = gGeoMan->GetMedium("air");
+ if ( ! gStsMedium ) Fatal("Main", "medium sts_air not found");
+ */
+ // --------------------------------------------------------------------------
+
+
+ // --------------- Create sensors ---------------------------------------
+ cout << endl << endl;
+ cout << "===> Creating sensors...." << endl << endl;
+ infoFile << endl << "Sensors: " << endl;
+ Int_t nSensors = CreateSensors();
+ for (Int_t iSensor = 1; iSensor <= nSensors; iSensor++) {
+ TString name = Form("Sensor%02d", iSensor);
+ TGeoVolume* sensor = gGeoMan->GetVolume(name);
+
+ // add color to sensors
+ if (iSensor == 1) sensor->SetLineColor(kRed);
+ if (iSensor == 2) sensor->SetLineColor(kGreen + 3);
+ if (iSensor == 3) sensor->SetLineColor(kBlue + 3);
+ if (iSensor == 4) sensor->SetLineColor(kAzure - 7);
+ if (iSensor == 5) sensor->SetLineColor(kYellow);
+ if (iSensor == 6) sensor->SetLineColor(kYellow);
+ if (iSensor == 7) sensor->SetLineColor(kYellow);
+
+ CheckVolume(sensor);
+ CheckVolume(sensor, infoFile);
+ }
+ // --------------------------------------------------------------------------
+
+
+ // ---------------- Create sectors --------------------------------------
+ cout << endl << endl;
+ cout << "===> Creating sectors...." << endl;
+ // infoFile << endl << "Sectors: " << endl;
+ Int_t nSectors = CreateSectors();
+ for (Int_t iSector = 1; iSector <= nSectors; iSector++) {
+ // cout << endl;
+ TString name = Form("Sector%02d", iSector);
+ TGeoVolume* sector = gGeoMan->GetVolume(name);
+ CheckVolume(sector);
+ // CheckVolume(sector, infoFile);
+ }
+ // --------------------------------------------------------------------------
+
+
+ // ---------------- Create ladders --------------------------------------
+ cout << endl << endl;
+ cout << "===> Creating ladders...." << endl;
+ infoFile << endl << "Ladders:" << endl;
+
+ TString name = "";
+ TGeoVolume* ladder;
+
+
+ Int_t nLadders = CreateLadders();
+
+ for (Int_t iLadder = 1; iLadder <= nLadders; iLadder++) {
+ cout << endl;
+
+ // name = Form("LadderType%02d", iLadder); // v19a
+ // name = Form("LadderType0%02d", iLadder); // v19b
+ name = Form("LadderType00%02d", iLadder); // v19k
+ ladder = gGeoMan->GetVolume(name);
+ CheckVolume(ladder);
+ CheckVolume(ladder, infoFile, kFALSE);
+
+ cout << "DF1: ladder name: " << name << endl << endl;
+
+ // name = Form("LadderType%02d", iLadder); // v19a
+ // name = Form("LadderType1%02d", iLadder); // v19b
+ name = Form("LadderType01%02d", iLadder); // v19k
+ ladder = gGeoMan->GetVolume(name);
+ CheckVolume(ladder);
+ CheckVolume(ladder, infoFile, kFALSE);
+
+ cout << "DF2: ladder name: " << name << endl << endl;
+
+ // name = Form("LadderType%02d", iLadder); // v19a
+ name = Form("LadderType10%02d", iLadder); // v19k
+ ladder = gGeoMan->GetVolume(name);
+ CheckVolume(ladder);
+ CheckVolume(ladder, infoFile, kFALSE);
+
+ cout << "DF3: ladder name: " << name << endl << endl;
+
+ // name = Form("LadderType%02d", iLadder); // v19a
+ name = Form("LadderType11%02d", iLadder); // v19k
+ ladder = gGeoMan->GetVolume(name);
+ CheckVolume(ladder);
+ CheckVolume(ladder, infoFile, kFALSE);
+
+ cout << "DF4: ladder name: " << name << endl << endl;
+ }
+ // --------------------------------------------------------------------------
+
+
+ // ---------------- Create cones ----------------------------------------
+ Double_t coneDz = 1.64;
+ TGeoVolume* coneSmallVolum = ConstructSmallCone(coneDz);
+ if (!coneSmallVolum) Fatal("ConstructSmallCone", "Volume Cone not found");
+ TGeoVolume* coneBigVolum = ConstructBigCone(coneDz);
+ if (!coneBigVolum) Fatal("ConstructBigCone", "Volume Cone not found");
+ // --------------------------------------------------------------------------
+
+ // ---------------- Create stations -------------------------------------
+ cout << endl << endl;
+ cout << "===> Creating stations...." << endl;
+ infoFile << endl << "Stations: " << endl;
+ Int_t angle = 0;
+ nLadders = 0;
+ Int_t ladderTypes[16]; // there are max 16 ladders in one layer
+ TGeoTranslation* statTrans = NULL;
+
+ TGeoVolume* myunit[32]; // units
+
+ // Int_t statPos[16] = { 30, 30, 40, 40, 50, 50, 60, 60,
+ // 70, 70, 80, 80, 90, 90, 100, 100 }; // z positions of units
+ // Int_t statPos18[18] = { 30, 30, // expanded for placement of Unit00
+ // 30, 30, 40, 40, 50, 50, 60, 60,
+ // 70, 70, 80, 80, 90, 90, 100, 100 }; // z positions of units
+
+ Double_t statPos[16] = {30.0, 30.0, 40.5, 40.5, 51.0, 51.0, 61.5, 61.5,
+ 72.0, 72.0, 82.5, 82.5, 93.0, 93.0, 103.5, 103.5}; // z positions of units
+ Double_t statPos18[18] = {30.0, 30.0, // expanded for placement of Unit00
+ 30.0, 30.0, 40.5, 40.5, 51.0, 51.0, 61.5, 61.5,
+ 72.0, 72.0, 82.5, 82.5, 93.0, 93.0, 103.5, 103.5}; // z positions of unit
+
+ //Double_t statPos[16] = {26.0, 26.0, 36.5, 36.5, 47.0, 47.0, 57.5, 57.5,
+ // 68.0, 68.0, 78.5, 78.5, 89.0, 89.0, 99.5, 99.5}; // z positions of units
+ //Double_t statPos18[18] = {26.0, 26.0, // expanded for placement of Unit00
+ // 26.0, 26.0, 36.5, 36.5, 47.0, 47.0, 57.5, 57.5,
+ // 68.0, 68.0, 78.5, 78.5, 89.0, 89.0, 99.5, 99.5}; // z positions of unit
+
+
+ ////Double_t rHole[8] = { 2.0, 2.0, 2.0, 2.9 , 3.7 , 3.7 , 4.2 , 4.2 }; // size of cutouts in stations
+ // Double_t rHole[8] = { 2.0, 2.0, 2.0, 2.43, 3.04, 3.35, 3.96, 4.2 }; // size of cutouts in stations, derived from gapXYZ[x][1]/2
+
+ Int_t cone_size[8] = {0, 0, 0, 1, 1, 1, 1, 1}; // size of cones: 0 = small, 1 = large
+
+ Double_t cone_offset[2] = {0.305, 0.285};
+
+ //==============================================================================================
+
+ // explanation: type xyzz
+ // where x = carbon ladder orientation
+ // where y = FEB box orientation
+ // where zz = sensor arrangement on ladder
+ // with FEB orientation - v19b
+ Int_t allUnitTypes[16][16] = {
+ {-1, -1, -1, -1, 10, 0, 9, 0, 101, 0, 109, 0, -1, -1, -1, -1}, // unit00D Station01 00
+ {-1, -1, -1, -1, 0, 1109, 0, 1101, 0, 1009, 0, 1010, -1, -1, -1, -1}, // unit01U Station01 01
+
+ {-1, -1, 0, 10, 0, 9, 0, 2, 0, 109, 0, 110, 0, 111, -1, -1}, // unit01D Station02 02
+ {-1, -1, 1130, 0, 1129, 0, 1128, 0, 1002, 0, 1028, 0, 1029, 0, -1, -1}, // unit02U Station02 03
+
+ {-1, -1, 14, 0, 12, 0, 12, 0, 103, 0, 112, 0, 113, 0, -1, -1}, // unit02D Station03 04
+ {-1, -1, 0, 1113, 0, 1112, 0, 1103, 0, 1012, 0, 1012, 0, 1014, -1, -1}, // unit03U Station03 05
+
+ {-1, 15, 0, 13, 0, 12, 0, 4, 0, 112, 0, 112, 0, 114, 0, -1}, // unit03D Station04 06
+ {-1, 0, 1133, 0, 1131, 0, 1131, 0, 1004, 0, 1031, 0, 1032, 0, 1034, -1}, // unit04U Station04 07
+
+ {-1, 0, 18, 0, 17, 0, 16, 0, 105, 0, 116, 0, 117, 0, 119, -1}, // unit04D Station05 08
+ {-1, 1119, 0, 1117, 0, 1116, 0, 1105, 0, 1016, 0, 1017, 0, 1018, 0, -1}, // unit05U Station05 09
+
+ {-1, 19, 0, 17, 0, 16, 0, 6, 0, 116, 0, 117, 0, 118, 0, -1}, // unit05D Station06 10
+ {-1, 0, 1137, 0, 1136, 0, 1135, 0, 1006, 0, 1035, 0, 1036, 0, 1038, -1}, // unit06U Station06 11
+
+ {21, 0, 25, 0, 20, 0, 20, 0, 107, 0, 120, 0, 120, 0, 127, 0}, // unit06D Station07 12
+ {0, 1127, 0, 1120, 0, 1120, 0, 1107, 0, 1020, 0, 1020, 0, 1025, 0, 1021}, // unit07U Station07 13
+
+ {0, 24, 0, 22, 0, 22, 0, 8, 0, 122, 0, 122, 0, 123, 0, 126}, // unit07D Station08 14
+ {1126, 0, 1123, 0, 1122, 0, 1122, 0, 1008, 0, 1022, 0, 1022, 0, 1024, 0}}; // unit08U Station08 15
+
+
+ //==============================================================================================
+
+ // without FEB orientation - v19a
+ // v19a Int_t allUnitTypes[16][16]=
+ // v19a { { -1, -1, -1, -1, 10, 0, 9, 0, 1, 0, 9, 0, -1, -1, -1, -1 }, // unit00D Station01 00
+ // v19a { -1, -1, -1, -1, 0, 109, 0, 101, 0, 109, 0, 110, -1, -1, -1, -1 }, // unit01U Station01 01
+ // v19a { -1, -1, 0, 10, 0, 9, 0, 2, 0, 9, 0, 10, 0, 11, -1, -1 }, // unit01D Station02 02
+ // v19a { -1, -1, 111, 0, 110, 0, 109, 0, 102, 0, 109, 0, 110, 0, -1, -1 }, // unit02U Station02 03
+ // v19a { -1, -1, 14, 0, 12, 0, 12, 0, 3, 0, 12, 0, 13, 0, -1, -1 }, // unit02D Station03 04
+ // v19a { -1, -1, 0, 113, 0, 112, 0, 103, 0, 112, 0, 112, 0, 114, -1, -1 }, // unit03U Station03 05
+ // v19a { -1, 15, 0, 13, 0, 12, 0, 4, 0, 12, 0, 12, 0, 14, 0, -1 }, // unit03D Station04 06
+ // v19a { -1, 0, 114, 0, 112, 0, 112, 0, 104, 0, 112, 0, 113, 0, 115, -1 }, // unit04U Station04 07
+ // v19a { -1, 0, 18, 0, 17, 0, 16, 0, 5, 0, 16, 0, 17, 0, 19, -1 }, // unit04D Station05 08
+ // v19a { -1, 119, 0, 117, 0, 116, 0, 105, 0, 116, 0, 117, 0, 118, 0, -1 }, // unit05U Station05 09
+ // v19a { -1, 19, 0, 17, 0, 16, 0, 6, 0, 16, 0, 17, 0, 18, 0, -1 }, // unit05D Station06 10
+ // v19a { -1, 0, 118, 0, 117, 0, 116, 0, 106, 0, 116, 0, 117, 0, 119, -1 }, // unit06U Station06 11
+ // v19a { 21, 0, 25, 0, 20, 0, 20, 0, 7, 0, 20, 0, 20, 0, 27, 0 }, // unit06D Station07 12
+ // v19a { 0, 127, 0, 120, 0, 120, 0, 107, 0, 120, 0, 120, 0, 125, 0, 121 }, // unit07U Station07 13
+ // v19a { 0, 24, 0, 22, 0, 22, 0, 8, 0, 22, 0, 22, 0, 23, 0, 26 }, // unit07D Station08 14
+ // v19a { 126, 0, 123, 0, 122, 0, 122, 0, 108, 0, 122, 0, 122, 0, 124, 0 } }; // unit08U Station08 15
+
+ // // generate unit
+ // for (Int_t iUnit = 0; iUnit < 16; iUnit++)
+ // for (Int_t iLadder = 0; iLadder < 16; iLadder++)
+ // {
+ // allUnitTypes[iUnit][iLadder] = 0;
+ // if ((iUnit % 2 == 0) && (allLadderTypes[iUnit/2][iLadder] < 100)) // if carbon structure is oriented upstream
+ // allUnitTypes[iUnit][iLadder] = allLadderTypes[iUnit/2][iLadder];
+ // if ((iUnit % 2 == 1) && (allLadderTypes[iUnit/2][iLadder] >= 100)) // if carbon structure is oriented downstream
+ // allUnitTypes[iUnit][iLadder] = allLadderTypes[iUnit/2][iLadder];
+ // }
+
+
+ // dump unit
+ for (Int_t iUnit = 0; iUnit < 16; iUnit++) {
+ cout << "DE unitTypes[" << iUnit << "] = { ";
+ for (Int_t iLadder = 0; iLadder < 16; iLadder++) {
+ cout << allUnitTypes[iUnit][iLadder];
+ if (iLadder < 15) cout << ", ";
+ else
+ cout << " };";
+ }
+ cout << endl;
+ }
+
+
+ // --- Units 01 - 16
+ for (Int_t iUnit = 0; iUnit < 16; iUnit++) {
+ cout << endl;
+
+ nLadders = 0;
+ //int ladderTypes[nLadders];
+ for (Int_t iLadder = 0; iLadder < 16; iLadder++)
+ if (allUnitTypes[iUnit][iLadder] >= 0) {
+ ladderTypes[nLadders] = allUnitTypes[iUnit][iLadder];
+ cout << "DE ladderTypes[" << nLadders << "] = " << allUnitTypes[iUnit][iLadder] << ";" << endl;
+ nLadders++;
+ }
+ myunit[iUnit * 2 + 0] = ConstructUnit(0, iUnit * 2 + 0, nLadders, ladderTypes, iUnit / 2 + 1);
+ myunit[iUnit * 2 + 1] = ConstructUnit(1, iUnit * 2 + 1, nLadders, ladderTypes, iUnit / 2 + 1);
+
+ // if (gkConstructCones) {
+ // if (iUnit%2 == 0)
+ // angle = 90;
+ // else
+ // angle = -90;
+ //
+ // // upstream
+ // TGeoRotation* coneRot11 = new TGeoRotation;
+ // coneRot11->RotateZ(angle);
+ // coneRot11->RotateY(180);
+ // TGeoCombiTrans* conePosRot11 = new TGeoCombiTrans(name+"conePosRot2", 0., 0., -coneDz-cone_offset[cone_size[iUnit]]-gkLadderGapZ/2., coneRot11);
+ // if (cone_size[iUnit] == 0)
+ // myunit[iUnit]->AddNode(coneSmallVolum, 1, conePosRot11);
+ // else
+ // myunit[iUnit]->AddNode(coneBigVolum, 1, conePosRot11);
+ //
+ // // downstream
+ // TGeoRotation* coneRot12 = new TGeoRotation;
+ // coneRot12->RotateZ(angle);
+ // TGeoCombiTrans* conePosRot12 = new TGeoCombiTrans(name+"conePosRot1", 0., 0., coneDz+cone_offset[cone_size[iUnit]]+gkLadderGapZ/2., coneRot12);
+ // if (cone_size[iUnit] == 0)
+ // myunit[iUnit]->AddNode(coneSmallVolum, 2, conePosRot12);
+ // else
+ // myunit[iUnit]->AddNode(coneBigVolum, 2, conePosRot12);
+ //
+ // myunit[iUnit]->GetShape()->ComputeBBox();
+ // }
+
+ // CheckVolume(myunit[iUnit]);
+ // CheckVolume(myunit[iUnit], infoFile);
+ if ((iUnit % 2 == 0) || (iUnit == 15)) {
+ CheckVolume(myunit[iUnit * 2 + 0]);
+ CheckVolume(myunit[iUnit * 2 + 0], infoFile);
+ CheckVolume(myunit[iUnit * 2 + 1]);
+ CheckVolume(myunit[iUnit * 2 + 1], infoFile);
+ }
+ infoFile << "Position z = " << statPos[iUnit] << endl;
+ }
+ // --------------------------------------------------------------------------
+
+
+ // --------------- Create STS volume ------------------------------------
+ cout << endl << endl;
+ cout << "===> Creating STS...." << endl;
+
+ // --- Create STS volume
+ TString stsName = "sts_";
+ stsName += geoTag;
+
+ // TGeoShape* stsShape = new TGeoCompositeShape("stsShape",
+ // "stsBox-stsCone1:trans1-stsCone2:trans2");
+ // TGeoVolume* sts = new TGeoVolume(stsName.Data(), stsShape, gStsMedium);
+
+ Double_t stsBorder = 2 * 5.;
+
+ TGeoVolume* sts = new TGeoVolumeAssembly(stsName.Data());
+
+ // --- Place stations in the STS
+ Double_t stsPosZ = 0.5 * (statPos[15] + statPos[0]); // todo units: update statPos[7]
+ // cout << "stsPosZ " << stsPosZ << " " << statPos[15] << " " << statPos[0] << "*****" << endl;
+
+ // for (Int_t iUnit = 0; iUnit < 16; iUnit++) {
+ for (Int_t iUnit = 0; iUnit < 18; iUnit++) {
+ // for (Int_t iUnit = 0; iUnit < 32; iUnit++) {
+ TGeoVolume* station = gGeoMan->GetVolume(unitName18[iUnit]);
+ // Double_t posZ = statPos[iUnit] - stsPosZ;
+ Double_t posZ = statPos18[iUnit] - stsPosZ;
+ // Double_t posZ = statPos[iUnit/2] - stsPosZ;
+ TGeoTranslation* trans = new TGeoTranslation(0., 0., posZ);
+ sts->AddNode(station, iUnit + 1, trans);
+ sts->GetShape()->ComputeBBox();
+ }
+
+ // --- Import passive elements from GDML file
+ if (gkImportPassive) { ImportPassive(sts, geoTag, infoFile); }
+
+ // --------------------------------------------------------------------------
+ // --------------------------------------------------------------------------
+
+
+ cout << endl;
+ CheckVolume(sts);
+
+
+ // --------------- Finish -----------------------------------------------
+ TGeoTranslation* stsTrans = new TGeoTranslation(0., 0., stsPosZ - 4.0);
+ top->AddNode(sts, 1, stsTrans);
+ top->GetShape()->ComputeBBox();
+ cout << endl << endl;
+
+ CheckVolume(top);
+ cout << endl << endl;
+ gGeoMan->CloseGeometry();
+ gGeoMan->CheckOverlaps(0.0001);
+ gGeoMan->PrintOverlaps();
+ gGeoMan->CheckOverlaps(0.0001, "s");
+ gGeoMan->PrintOverlaps();
+ gGeoMan->Test();
+
+ //TFile* geoFile = new TFile(geoFileName, "RECREATE");
+ //top->Write();
+ sts->Export(geoFileName); // an another way of writing the stsvolume placement
+
+ TFile* geoFile = new TFile(geoFileName, "UPDATE");
+ TGeoTranslation* sts_placement = new TGeoTranslation("sts_trans", 0., 0., stsPosZ - 4.0);
+ sts_placement->Write();
+ //*/
+ cout << endl;
+ cout << "Geometry " << top->GetName() << " written to " << geoFileName << endl;
+ geoFile->Close();
+
+ TString geoFileName_ = "sts_";
+ geoFileName_ = geoFileName_ + geoTag + "_geo.root";
+
+ geoFile = new TFile(geoFileName_, "RECREATE");
+ gGeoMan->Write(); // use this is you want GeoManager format in the output
+ geoFile->Close();
+
+ TString geoFileName__ = "sts_";
+ geoFileName_ = geoFileName__ + geoTag + "-geo.root";
+ sts->Export(geoFileName_);
+
+ geoFile = new TFile(geoFileName_, "UPDATE");
+ stsTrans->Write();
+ geoFile->Close();
+
+
+ // gGeoManager->FindVolumeFast("LadderType10_CarbonElement")->Draw("ogl");
+ top->Draw("ogl");
+ //sts->Draw("ogl");
+ gGeoManager->SetVisLevel(9);
+
+ infoFile.close();
+}
+// ============================================================================
+// ====== End of main function =====
+// ============================================================================
+
+
+// ****************************************************************************
+// ***** Definition of media, sensors, sectors and ladders *****
+// ***** *****
+// ***** Decoupled from main function for better readability *****
+// ****************************************************************************
+
+/** ===========================================================================
+ ** Create media
+ **
+ ** Currently created: air, active silicon, passive silion
+ **
+ ** Not used for the time being
+ **/
+Int_t CreateMedia()
+{
+
+ Int_t nMedia = 0;
+ Double_t density = 0.;
+
+ // --- Material air
+ density = 1.205e-3; // [g/cm^3]
+ TGeoMixture* matAir = new TGeoMixture("sts_air", 3, density);
+ matAir->AddElement(14.0067, 7, 0.755); // Nitrogen
+ matAir->AddElement(15.999, 8, 0.231); // Oxygen
+ matAir->AddElement(39.948, 18, 0.014); // Argon
+
+ // --- Material silicon
+ density = 2.33; // [g/cm^3]
+ TGeoElement* elSi = gGeoMan->GetElementTable()->GetElement(14);
+ TGeoMaterial* matSi = new TGeoMaterial("matSi", elSi, density);
+
+
+ // --- Air (passive)
+ TGeoMedium* medAir = new TGeoMedium("air", nMedia++, matAir);
+ medAir->SetParam(0, 0.); // is passive
+ medAir->SetParam(1, 1.); // is in magnetic field
+ medAir->SetParam(2, 20.); // max. field [kG]
+ medAir->SetParam(6, 0.001); // boundary crossing precision [cm]
+
+
+ // --- Active silicon for sensors
+ TGeoMedium* medSiAct = new TGeoMedium("silicon", nMedia++, matSi);
+ medSiAct->SetParam(0, 1.); // is active
+ medSiAct->SetParam(1, 1.); // is in magnetic field
+ medSiAct->SetParam(2, 20.); // max. field [kG]
+ medSiAct->SetParam(6, 0.001); // boundary crossing precisison [cm]
+
+ // --- Passive silicon for cables
+ TGeoMedium* medSiPas = new TGeoMedium("carbon", nMedia++, matSi);
+ medSiPas->SetParam(0, 0.); // is passive
+ medSiPas->SetParam(1, 1.); // is in magnetic field
+ medSiPas->SetParam(2, 20.); // max. field [kG]
+ medSiPas->SetParam(6, 0.001); // boundary crossing precisison [cm]
+
+ return nMedia;
+}
+/** ======================================================================= **/
+
+
+/** ===========================================================================
+ ** Create sensors
+ **
+ ** Sensors are created as volumes with box shape and active silicon as medium.
+ ** Four kinds of sensors: 3.2x2.2, 6.2x2.2, 6.2x4.2, 6.2x6.2
+ **/
+Int_t CreateSensors()
+{
+
+ Int_t nSensors = 0;
+
+ Double_t xSize = 0.;
+ Double_t ySize = 0.;
+ Double_t zSize = gkSensorThickness;
+ TGeoMedium* silicon = gGeoMan->GetMedium("silicon");
+
+
+ // --- Sensor type 01: Small sensor (6.2 cm x 2.2 cm)
+ xSize = gkSensorSizeX;
+ ySize = 2.2;
+ TGeoBBox* shape_sensor01 = new TGeoBBox("sensor01", xSize / 2., ySize / 2., zSize / 2.);
+ new TGeoVolume("Sensor01", shape_sensor01, silicon);
+ nSensors++;
+
+
+ // --- Sensor type 02: Medium sensor (6.2 cm x 4.2 cm)
+ xSize = gkSensorSizeX;
+ ySize = 4.2;
+ TGeoBBox* shape_sensor02 = new TGeoBBox("sensor02", xSize / 2., ySize / 2., zSize / 2.);
+ new TGeoVolume("Sensor02", shape_sensor02, silicon);
+ nSensors++;
+
+
+ // --- Sensor type 03: Big sensor (6.2 cm x 6.2 cm)
+ xSize = gkSensorSizeX;
+ ySize = 6.2;
+ TGeoBBox* shape_sensor03 = new TGeoBBox("sensor03", xSize / 2., ySize / 2., zSize / 2.);
+ new TGeoVolume("Sensor03", shape_sensor03, silicon);
+ nSensors++;
+
+
+ // --- Sensor type 04: Big sensor (6.2 cm x 12.4 cm)
+ xSize = gkSensorSizeX;
+ ySize = 12.4;
+ TGeoBBox* shape_sensor04 = new TGeoBBox("sensor04", xSize / 2., ySize / 2., zSize / 2.);
+ new TGeoVolume("Sensor04", shape_sensor04, silicon);
+ nSensors++;
+
+
+ // below are extra small sensors, those are not available in the CAD model
+
+ // --- Sensor Type 05: Half small sensor (4 cm x 2.5 cm)
+ xSize = 4.0;
+ ySize = 2.5;
+ TGeoBBox* shape_sensor05 = new TGeoBBox("sensor05", xSize / 2., ySize / 2., zSize / 2.);
+ new TGeoVolume("Sensor05", shape_sensor05, silicon);
+ nSensors++;
+
+
+ // --- Sensor type 06: Additional "in hole" sensor (3.1 cm x 4.2 cm)
+ xSize = 3.1;
+ ySize = 4.2;
+ TGeoBBox* shape_sensor06 = new TGeoBBox("sensor06", xSize / 2., ySize / 2., zSize / 2.);
+ new TGeoVolume("Sensor06", shape_sensor06, silicon);
+ nSensors++;
+
+
+ // --- Sensor type 07: Mini Medium sensor (1.5 cm x 4.2 cm)
+ xSize = 1.5;
+ ySize = 4.2;
+ TGeoBBox* shape_sensor07 = new TGeoBBox("sensor07", xSize / 2., ySize / 2., zSize / 2.);
+ new TGeoVolume("Sensor07", shape_sensor07, silicon);
+ nSensors++;
+
+
+ return nSensors;
+}
+/** ======================================================================= **/
+
+
+/** ===========================================================================
+ ** Create sectors
+ **
+ ** A sector is either a single sensor or several chained sensors.
+ ** It is implemented as TGeoVolumeAssembly.
+ ** Currently available:
+ ** - single sensors of type 1 - 4
+ ** - two chained sensors of type 4
+ ** - three chained sensors of type 4
+ **/
+Int_t CreateSectors()
+{
+
+ Int_t nSectors = 0;
+
+ TGeoVolume* sensor01 = gGeoMan->GetVolume("Sensor01");
+ TGeoVolume* sensor02 = gGeoMan->GetVolume("Sensor02");
+ TGeoVolume* sensor03 = gGeoMan->GetVolume("Sensor03");
+ TGeoVolume* sensor04 = gGeoMan->GetVolume("Sensor04");
+ TGeoVolume* sensor05 = gGeoMan->GetVolume("Sensor05");
+ TGeoVolume* sensor06 = gGeoMan->GetVolume("Sensor06");
+ TGeoVolume* sensor07 = gGeoMan->GetVolume("Sensor07");
+ // TGeoBBox* box4 = (TGeoBBox*) sensor04->GetShape();
+
+ // --- Sector type 1: single sensor of type 1
+ TGeoVolumeAssembly* sector01 = new TGeoVolumeAssembly("Sector01");
+ sector01->AddNode(sensor01, 1);
+ sector01->GetShape()->ComputeBBox();
+ nSectors++;
+
+ // --- Sector type 2: single sensor of type 2
+ TGeoVolumeAssembly* sector02 = new TGeoVolumeAssembly("Sector02");
+ sector02->AddNode(sensor02, 1);
+ sector02->GetShape()->ComputeBBox();
+ nSectors++;
+
+ // --- Sector type 3: single sensor of type 3
+ TGeoVolumeAssembly* sector03 = new TGeoVolumeAssembly("Sector03");
+ sector03->AddNode(sensor03, 1);
+ sector03->GetShape()->ComputeBBox();
+ nSectors++;
+
+ // --- Sector type 4: single sensor of type 4
+ TGeoVolumeAssembly* sector04 = new TGeoVolumeAssembly("Sector04");
+ sector04->AddNode(sensor04, 1);
+ sector04->GetShape()->ComputeBBox();
+ nSectors++;
+
+ // --- Sector type 5: single sensor of type 5
+ TGeoVolumeAssembly* sector05 = new TGeoVolumeAssembly("Sector05");
+ sector05->AddNode(sensor05, 1);
+ sector05->GetShape()->ComputeBBox();
+ nSectors++;
+
+ // --- Sector type 6: single sensor of type 6
+ TGeoVolumeAssembly* sector06 = new TGeoVolumeAssembly("Sector06");
+ sector06->AddNode(sensor06, 1);
+ sector06->GetShape()->ComputeBBox();
+ nSectors++;
+
+ // --- Sector type 7: single sensor of type 7
+ TGeoVolumeAssembly* sector07 = new TGeoVolumeAssembly("Sector07");
+ sector07->AddNode(sensor07, 1);
+ sector07->GetShape()->ComputeBBox();
+ nSectors++;
+
+ // // --- Sector type 5: two sensors of type 4
+ // TGeoVolumeAssembly* sector05 = new TGeoVolumeAssembly("Sector05");
+ // Double_t shift5 = 0.5 * gkChainGapY + box4->GetDY();
+ // TGeoTranslation* transD5 =
+ // new TGeoTranslation("td", 0., -1. * shift5, 0.);
+ // TGeoTranslation* transU5 =
+ // new TGeoTranslation("tu", 0., shift5, 0.);
+ // sector05->AddNode(sensor04, 1, transD5);
+ // sector05->AddNode(sensor04, 2, transU5);
+ // sector05->GetShape()->ComputeBBox();
+ // nSectors++;
+
+ return nSectors;
+}
+/** ======================================================================= **/
+
+
+/** ===========================================================================
+ ** Create ladders
+ **
+ ** Ladders are the building blocks of the stations. They contain
+ ** several modules placed one after the other along the z axis
+ ** such that the sectors are arranged vertically (with overlap).
+ **
+ ** A ladder is constructed out of two half ladders, the second of which
+ ** is rotated in the x-y plane by 180 degrees and displaced
+ ** in z direction.
+ **/
+Int_t CreateLadders()
+{
+
+ Int_t nLadders = 0;
+
+ // --- Some variables
+ Int_t nSectors = 0;
+ Int_t sectorTypes[10];
+ TGeoBBox* shape = NULL;
+ TString s0name;
+ TString hlname;
+ char align;
+ TGeoVolume* s0vol = NULL;
+ TGeoVolume* halfLadderU = NULL;
+ TGeoVolume* halfLadderD = NULL;
+
+ // --- Ladders 01-23
+ // --- Ladders 01-27
+ Int_t allSectorTypes[38][6] = {
+ {1, 2, 3, 3, 0, -1}, // ladder 01 - 5 - last column defines alignment of small sensors
+ {1, 2, 3, 3, 0, 0}, // ladder 02 - 5 - last column defines alignment of small sensors
+ {2, 2, 3, 4, 0, -1}, // ladder 03 - 6 - last column defines alignment of small sensors
+ {2, 2, 3, 4, 0, 0}, // ladder 04 - 6 - last column defines alignment of small sensors
+ {2, 3, 4, 4, 0, -1}, // ladder 05 - 7 - last column defines alignment of small sensors
+ {2, 3, 4, 4, 0, 0}, // ladder 06 - 7 - last column defines alignment of small sensors
+ {2, 2, 3, 4, 4, 0}, // ladder 07 - last column defines alignment of small sensors
+ {3, 4, 4, 4, 0, 0}, // ladder 08 - last column defines alignment of small sensors
+ {1, 1, 2, 3, 3, 0}, // ladder 09 - last column defines alignment of small sensors
+ {1, 1, 2, 2, 3, 0}, // ladder 10 - last column defines alignment of small sensors
+ {2, 2, 0, 0, 0, 0}, // ladder 11 - last column defines alignment of small sensors
+ {2, 2, 2, 3, 4, 0}, // ladder 12 - last column defines alignment of small sensors
+ {2, 2, 3, 4, 0, 0}, // ladder 13 - last column defines alignment of small sensors
+ {2, 3, 4, 0, 0, 0}, // ladder 14 - last column defines alignment of small sensors
+ {3, 3, 0, 0, 0, 0}, // ladder 15 - last column defines alignment of small sensors
+ {2, 2, 3, 4, 4, 0}, // ladder 16 - last column defines alignment of small sensors
+ {2, 3, 4, 4, 0, 0}, // ladder 17 - last column defines alignment of small sensors
+ {3, 4, 4, 0, 0, 0}, // ladder 18 - last column defines alignment of small sensors
+ {4, 4, 0, 0, 0, 0}, // ladder 19 - last column defines alignment of small sensors
+ {1, 2, 4, 4, 4, 0}, // ladder 20 - last column defines alignment of small sensors
+ {4, 0, 0, 0, 0, 0}, // ladder 21 - last column defines alignment of small sensors
+ {2, 3, 4, 4, 4, 0}, // ladder 22 - last column defines alignment of small sensors
+ {2, 3, 3, 4, 4, 0}, // ladder 23 - last column defines alignment of small sensors
+ {2, 3, 4, 4, 0, 0}, // ladder 24 - copy of 17 with different total length
+ {3, 4, 4, 0, 0, 0}, // ladder 25 - copy of 18 with different total length
+ {4, 4, 0, 0, 0, 0}, // ladder 26 - copy of 19 with different total length
+ {4, 4, 0, 0, 0, 0}, // ladder 27 - copy of 19 with different total length
+
+ {1, 1, 2, 3, 3, 0}, // ladder 28 - copy of 09 with different total length
+ {1, 1, 2, 2, 3, 0}, // ladder 29 - copy of 10 with different total length
+ {2, 2, 0, 0, 0, 0}, // ladder 30 - copy of 11 with different total length
+
+ {2, 2, 2, 3, 4, 0}, // ladder 31 - copy of 12 with different total length
+ {2, 2, 3, 4, 0, 0}, // ladder 32 - copy of 13 with different total length
+ {2, 3, 4, 0, 0, 0}, // ladder 33 - copy of 14 with different total length
+ {3, 3, 0, 0, 0, 0}, // ladder 34 - copy of 15 with different total length
+
+ {2, 2, 3, 4, 4, 0}, // ladder 35 - copy of 16 with different total length
+ {2, 3, 4, 4, 0, 0}, // ladder 36 - copy of 17 with different total length
+ {3, 4, 4, 0, 0, 0}, // ladder 37 - copy of 18 with different total length
+ {4, 4, 0, 0, 0, 0}, // ladder 38 - copy of 19 with different total length
+
+ }; // 8 full - 19 partial ladders
+
+ // Issue #405
+ // Counting from the most upstream ladder, the gaps between sensors are as follows:
+ // 01 (most upstream): 41.3mm
+ // 02: 41.3mm
+ // 03: 42.0mm
+ // 04: 48.6mm
+ // 05: 60.8mm
+ // 06: 67.0mm
+ // 07: 79.2mm
+ // 08 (most downstream): 88.0mm
+
+ Double_t pitchZ = 0;
+ Double_t gapXYZ[38][3] = {
+ {0., 5.43, 0.}, // ladder 01
+ {0., 6.00, 0.}, // ladder 02
+ {0., 6.70, 0.}, // ladder 03
+ {0., 7.35, 0.}, // ladder 04
+ {0., 8.02, 0.}, // ladder 05
+ {0., 8.70, 0.}, // ladder 06
+ {0., 9.22, 0.}, // ladder 07
+ {0., 10.00, 0.}, // ladder 08
+ {0., -gkSectorOverlapY, 0.}, // ladder 09
+ {0., -gkSectorOverlapY, 0.}, // ladder 10
+ {0., -gkSectorOverlapY, 0.}, // ladder 11
+ {0., -gkSectorOverlapY, 0.}, // ladder 12
+ {0., -gkSectorOverlapY, 0.}, // ladder 13
+ {0., -gkSectorOverlapY, 0.}, // ladder 14
+ {0., -gkSectorOverlapY, 0.}, // ladder 15
+ {0., -gkSectorOverlapY, 0.}, // ladder 16
+ {0., -gkSectorOverlapY, 0.}, // ladder 17
+ {0., -gkSectorOverlapY, 0.}, // ladder 18
+ {0., -gkSectorOverlapY, 0.}, // ladder 19
+ {0., -gkSectorOverlapY, 0.}, // ladder 20
+ {0., -gkSectorOverlapY, 0.}, // ladder 21
+ {0., -gkSectorOverlapY, 0.}, // ladder 22
+ {0., -gkSectorOverlapY, 0.}, // ladder 23
+
+ {0., -gkSectorOverlapY, 0.}, // ladder 24 - copy of 17 with different total length
+ {0., -gkSectorOverlapY, 0.}, // ladder 25 - copy of 18 with different total length
+ {0., -gkSectorOverlapY, 0.}, // ladder 26 - copy of 19 with different total length
+ {0., -gkSectorOverlapY, 0.}, // ladder 27 - copy of 19 with different total length
+
+ {0., -gkSectorOverlapY, 0.}, // ladder 28 - copy of 09 with different total length
+ {0., -gkSectorOverlapY, 0.}, // ladder 29 - copy of 10 with different total length
+ {0., -gkSectorOverlapY, 0.}, // ladder 30 - copy of 11 with different total length
+
+ {0., -gkSectorOverlapY, 0.}, // ladder 31 - copy of 12 with different total length
+ {0., -gkSectorOverlapY, 0.}, // ladder 32 - copy of 13 with different total length
+ {0., -gkSectorOverlapY, 0.}, // ladder 33 - copy of 14 with different total length
+ {0., -gkSectorOverlapY, 0.}, // ladder 34 - copy of 15 with different total length
+
+ {0., -gkSectorOverlapY, 0.}, // ladder 35 - copy of 16 with different total length
+ {0., -gkSectorOverlapY, 0.}, // ladder 36 - copy of 17 with different total length
+ {0., -gkSectorOverlapY, 0.}, // ladder 37 - copy of 18 with different total length
+ {0., -gkSectorOverlapY, 0.}, // ladder 38 - copy of 19 with different total length
+
+ };
+
+ Double_t ladderLength[38] = {
+ 48.0, // ladder 01
+ 48.0, // ladder 02
+ 66.8, // ladder 03
+ 66.8, // ladder 04
+ 81.8, // ladder 05
+ 81.8, // ladder 06
+ 89.2, // ladder 07
+ 96.7, // ladder 08
+ 48.0, // ladder 09
+ 48.0, // ladder 10
+ 48.0, // ladder 11
+ 66.8, // ladder 12
+ 66.8, // ladder 13
+ 66.8, // ladder 14
+ 66.8, // ladder 15
+ 81.8, // ladder 16
+ 81.8, // ladder 17
+ 81.8, // ladder 18
+ 81.8, // ladder 19
+ 89.2, // ladder 20
+ 89.2, // ladder 21
+ 96.7, // ladder 22
+ 96.7, // ladder 23
+
+ 96.7, // ladder 24 - copy of 17 with different total length
+ 89.2, // ladder 25 - copy of 18 with different total length
+ 96.7, // ladder 26 - copy of 19 with different total length
+ 89.2, // ladder 27 - copy of 19 with different total length
+
+ 66.8, // ladder 28 - copy of 09 with different total length
+ 66.8, // ladder 29 - copy of 10 with different total length
+ 66.8, // ladder 30 - copy of 11 with different total length
+
+ 81.8, // ladder 31 - copy of 12 with different total length
+ 81.8, // ladder 32 - copy of 13 with different total length
+ 81.8, // ladder 33 - copy of 14 with different total length
+ 81.8, // ladder 34 - copy of 15 with different total length
+
+ 89.2, // ladder 35 - copy of 16 with different total length
+ 89.2, // ladder 36 - copy of 17 with different total length
+ 89.2, // ladder 37 - copy of 18 with different total length
+ 89.2, // ladder 38 - copy of 19 with different total length
+
+ };
+ // ========================================================================
+
+ // calculate Z shift for ladders with and without gaps in the center
+ s0name = Form("Sector%02d", allSectorTypes[0][0]);
+ s0vol = gGeoMan->GetVolume(s0name);
+ shape = (TGeoBBox*) s0vol->GetShape();
+
+ // ========================================================================
+
+ for (Int_t iLadder = 0; iLadder < 38; iLadder++) {
+ cout << endl;
+ nSectors = 0;
+ for (Int_t i = 0; i < 5; i++)
+ if (allSectorTypes[iLadder][i] != 0) {
+ sectorTypes[nSectors] = allSectorTypes[iLadder][i]; // copy sectors for this ladder
+ cout << "DE iLadder " << iLadder + 1 << " sectorTypes[" << nSectors << "] = " << allSectorTypes[iLadder][i]
+ << ";" << endl;
+ nSectors++; // count how many sectors are in this ladder
+ }
+
+ // always set displacement in z between upper and lower half ladder
+ gapXYZ[iLadder][2] = 2. * shape->GetDZ() + gkSectorGapZ;
+
+ // define additional offset to carbon ladder for half ladders with less than 5 sensors
+ pitchZ = 2. * shape->GetDZ() + gkSectorGapZ;
+
+ if (allSectorTypes[iLadder][5] == 0) align = 'l';
+ else
+ align = 'r';
+ hlname = Form("HalfLadder%02dt", iLadder + 1);
+ // build upper half ladder
+ halfLadderU = ConstructHalfLadder(hlname, nSectors, sectorTypes, align, ladderLength[iLadder] / 2.,
+ gapXYZ[iLadder][1] / 2.); // mirrored
+
+ if (allSectorTypes[iLadder][5] == 0) align = 'r';
+ else
+ align = 'l';
+ hlname = Form("HalfLadder%02db", iLadder + 1);
+ // build lower half ladder
+ halfLadderD = ConstructHalfLadder(hlname, nSectors, sectorTypes, align, ladderLength[iLadder] / 2.,
+ gapXYZ[iLadder][1] / 2.); // mirrored
+
+ // at this point half ladders are constructed
+
+ // build all 4 possible ladders types for this sensor arrangement
+ ConstructLadder(iLadder + 1, halfLadderU, halfLadderD, gapXYZ[iLadder][1], gapXYZ[iLadder][2], pitchZ,
+ nSectors); // v19a
+ ConstructLadder(100 + iLadder + 1, halfLadderU, halfLadderD, gapXYZ[iLadder][1], gapXYZ[iLadder][2], pitchZ,
+ nSectors); // v19b
+
+ ConstructLadder(1000 + iLadder + 1, halfLadderU, halfLadderD, gapXYZ[iLadder][1], gapXYZ[iLadder][2], pitchZ,
+ nSectors); // v19k
+ ConstructLadder(1100 + iLadder + 1, halfLadderU, halfLadderD, gapXYZ[iLadder][1], gapXYZ[iLadder][2], pitchZ,
+ nSectors); // v19k
+
+ nLadders++;
+ }
+
+ return nLadders;
+}
+/** ======================================================================= **/
+
+
+// ****************************************************************************
+// ***** *****
+// ***** Generic functions for the construction of STS elements *****
+// ***** *****
+// ***** module: volume (made of a sector and a cable) *****
+// ***** haf ladder: assembly (made of modules) *****
+// ***** ladder: assembly (made of two half ladders) *****
+// ***** station: volume (made of ladders) *****
+// ***** *****
+// ****************************************************************************
+
+
+/** ===========================================================================
+ ** Construct a module
+ **
+ ** A module is a sector plus the readout cable extending from the
+ ** top of the sector. The cable is made from passive silicon.
+ ** The cable has the same x size as the sector.
+ ** Its thickness is given by the global variable gkCableThickness.
+ ** The cable length is a parameter.
+ ** The sensor(s) of the sector is/are placed directly in the module;
+ ** the sector is just auxiliary for the proper placement.
+ **
+ ** Arguments:
+ ** name volume name
+ ** sector pointer to sector volume
+ ** cableLength length of cable
+ **/
+TGeoVolume* ConstructModule(const char* name, TGeoVolume* sector, Double_t cableLength)
+{
+
+ // --- Check sector volume
+ if (!sector) Fatal("CreateModule", "Sector volume not found!");
+
+ // --- Get size of sector
+ TGeoBBox* box = (TGeoBBox*) sector->GetShape();
+ Double_t sectorX = 2. * box->GetDX();
+ Double_t sectorY = 2. * box->GetDY();
+ Double_t sectorZ = 2. * box->GetDZ();
+
+ // --- Get size of cable
+ Double_t cableX = sectorX;
+ Double_t cableY = cableLength;
+ Double_t cableZ = gkCableThickness;
+
+ // --- Create module volume
+ Double_t moduleX = TMath::Max(sectorX, cableX);
+ Double_t moduleY = sectorY + cableLength;
+
+ Double_t moduleZ = TMath::Max(sectorZ, cableZ);
+
+ TGeoVolume* module = gGeoManager->MakeBox(name, gStsMedium, moduleX / 2., moduleY / 2., moduleZ / 2.);
+
+ // --- Position of sector in module
+ // --- Sector is centred in x and z and aligned to the bottom
+ Double_t sectorXpos = 0.;
+ Double_t sectorYpos = 0.5 * (sectorY - moduleY);
+ Double_t sectorZpos = 0.;
+
+
+ // --- Get sensor(s) from sector
+ Int_t nSensors = sector->GetNdaughters();
+ for (Int_t iSensor = 0; iSensor < nSensors; iSensor++) {
+ TGeoNode* sensor = sector->GetNode(iSensor);
+
+ // --- Calculate position of sensor in module
+ const Double_t* xSensTrans = sensor->GetMatrix()->GetTranslation();
+ Double_t sensorXpos = 0.;
+ Double_t sensorYpos = sectorYpos + xSensTrans[1];
+ Double_t sensorZpos = 0.;
+ TGeoTranslation* sensTrans = new TGeoTranslation("sensTrans", sensorXpos, sensorYpos, sensorZpos);
+
+ // --- Add sensor volume to module
+ TGeoVolume* sensVol = sensor->GetVolume();
+ module->AddNode(sensor->GetVolume(), iSensor + 1, sensTrans);
+ module->GetShape()->ComputeBBox();
+ }
+
+
+ // --- Create cable volume, if necessary, and place it in module
+ // --- Cable is centred in x and z and aligned to the top
+ if (gkConstructCables && cableLength > 0.0001) {
+ TString cableName = TString(name) + "_cable";
+ TGeoMedium* cableMedium = gGeoMan->GetMedium("STScable");
+ if (!cableMedium) Fatal("CreateModule", "Medium STScable not found!");
+ TGeoVolume* cable = gGeoManager->MakeBox(cableName.Data(), cableMedium, cableX / 2., cableY / 2., cableZ / 2.);
+ // add color to cables
+ cable->SetLineColor(kOrange);
+ cable->SetTransparency(60);
+ Double_t cableXpos = 0.;
+ Double_t cableYpos = sectorY + 0.5 * cableY - 0.5 * moduleY;
+ Double_t cableZpos = 0.;
+ TGeoTranslation* cableTrans = new TGeoTranslation("cableTrans", cableXpos, cableYpos, cableZpos);
+ module->AddNode(cable, 1, cableTrans);
+ module->GetShape()->ComputeBBox();
+ }
+
+ return module;
+}
+/** ======================================================================= **/
+
+
+/** ===========================================================================
+ ** Construct a half ladder
+ **
+ ** A half ladder is a virtual volume (TGeoVolumeAssembly) consisting
+ ** of several modules arranged on top of each other. The modules
+ ** have a given overlap in y and a displacement in z to allow for the
+ ** overlap.
+ **
+ ** The typ of sectors / modules to be placed must be specified:
+ ** 1 = sensor01
+ ** 2 = sensor02
+ ** 3 = sensor03
+ ** 4 = sensor04
+ ** 5 = 2 x sensor04 (chained)
+ ** 6 = 3 x sensor04 (chained)
+ ** The cable is added automatically from the top of each sensor to
+ ** the top of the half ladder.
+ ** The alignment can be left (l) or right (r), which matters in the
+ ** case of different x sizes of sensors (e.g. SensorType01).
+ **
+ ** Arguments:
+ ** name volume name
+ ** nSectors number of sectors
+ ** sectorTypes array with sector types
+ ** align horizontal alignment of sectors
+ * ladderLength full length of the ladder towards FEE
+ * offsetY gap in the beam-pipe region
+ **/
+TGeoVolume* ConstructHalfLadder(const TString& name, Int_t nSectors, Int_t* sectorTypes, char align,
+ Double_t ladderLength, Double_t offsetY)
+{
+
+ // --- Create half ladder volume assembly
+ TGeoVolumeAssembly* halfLadder = new TGeoVolumeAssembly(name);
+
+ // --- Determine size of ladder
+ Double_t ladderX = 0.;
+ Double_t ladderY = 0.;
+ Double_t ladderZ = 0.;
+ for (Int_t iSector = 0; iSector < nSectors; iSector++) {
+ TString sectorName = Form("Sector%02d", sectorTypes[iSector]);
+ TGeoVolume* sector = gGeoMan->GetVolume(sectorName);
+ if (!sector) Fatal("ConstructHalfLadder", Form("Volume %s not found", sectorName.Data()));
+ TGeoBBox* box = (TGeoBBox*) sector->GetShape();
+ // --- Ladder x size equals largest sector x size
+ ladderX = TMath::Max(ladderX, 2. * box->GetDX());
+ // --- Ladder y size is sum of sector ysizes
+ ladderY += 2. * box->GetDY();
+ // --- Ladder z size is sum of sector z sizes
+ ladderZ += 2. * box->GetDZ();
+ }
+ // --- Subtract overlaps in y
+ ladderY -= Double_t(nSectors - 1) * gkSectorOverlapY;
+ // --- Add gaps in z direction
+ ladderZ += Double_t(nSectors - 1) * gkSectorGapZ;
+
+ ladderY = TMath::Max(ladderLength - offsetY, ladderY);
+
+ // --- Create and place modules
+ Double_t yPosSect = -0.5 * ladderY;
+ Double_t zPosMod = -0.5 * ladderZ;
+ for (Int_t iSector = 0; iSector < nSectors; iSector++) {
+ TString sectorName = Form("Sector%02d", sectorTypes[iSector]);
+ TGeoVolume* sector = gGeoMan->GetVolume(sectorName);
+ TGeoBBox* box = (TGeoBBox*) sector->GetShape();
+ Double_t sectorX = 2. * box->GetDX();
+ Double_t sectorY = 2. * box->GetDY();
+ Double_t sectorZ = 2. * box->GetDZ();
+ yPosSect += 0.5 * sectorY; // Position of sector in ladder
+ Double_t cableLength = 0.5 * ladderY - yPosSect - 0.5 * sectorY;
+ TString moduleName = name + "_" + Form("Module%02d", sectorTypes[iSector]);
+ TGeoVolume* module = ConstructModule(moduleName.Data(), sector, cableLength);
+
+ TGeoBBox* shapeMod = (TGeoBBox*) module->GetShape();
+ Double_t moduleX = 2. * shapeMod->GetDX();
+ Double_t moduleY = 2. * shapeMod->GetDY();
+ Double_t moduleZ = 2. * shapeMod->GetDZ();
+ Double_t xPosMod = 0.;
+ if (align == 'l') xPosMod = 0.5 * (moduleX - ladderX); // left aligned
+ else if (align == 'r')
+ xPosMod = 0.5 * (ladderX - moduleX); // right aligned
+ else
+ xPosMod = 0.; // centred in x
+ Double_t yPosMod = 0.5 * (ladderY - moduleY); // top aligned
+ zPosMod += 0.5 * moduleZ;
+ TGeoTranslation* trans = new TGeoTranslation("t", xPosMod, yPosMod, zPosMod);
+ halfLadder->AddNode(module, iSector + 1, trans);
+ halfLadder->GetShape()->ComputeBBox();
+ yPosSect += 0.5 * sectorY - gkSectorOverlapY;
+ zPosMod += 0.5 * moduleZ + gkSectorGapZ;
+ }
+
+ CheckVolume(halfLadder);
+ cout << endl;
+
+ return halfLadder;
+}
+/** ======================================================================= **/
+
+
+/** ===========================================================================
+ ** Add a carbon support to a ladder
+ **
+ ** Arguments:
+ ** LadderIndex ladder number
+ ** ladder pointer to ladder
+ ** xu size of halfladder
+ ** ladderY height of ladder along y
+ ** ladderZ thickness of ladder along z
+ **/
+void AddCarbonLadder(Int_t LadderIndex, TGeoVolume* ladder, Double_t xu, Double_t ladderY, Double_t ladderZ)
+{
+
+ // --- Some variables
+ TString name = Form("LadderType%04d", LadderIndex); // v19k
+ Int_t i;
+ Double_t j;
+
+ Int_t YnumOfFrameBoxes = round(ladderY / gkFrameStep);
+
+ // cout << "DEXZ: lad " << LadderIndex << " inum " << YnumOfFrameBoxes << endl;
+
+ Double_t ladderDZ = (xu / 2. + sqrt(2.) * gkFrameThickness / 2. + gkSectorGapZFrame * 2) / 2.;
+ cout << "DEFR: frame Z size " << 2 * ladderDZ << " cm" << endl;
+
+ TGeoBBox* fullFrameShp = new TGeoBBox(name + "_CarbonElement_shp", xu / 2., gkFrameStep / 2., ladderDZ);
+ TGeoVolume* fullFrameBoxVol = new TGeoVolume(name + "_CarbonElement", fullFrameShp, gStsMedium);
+
+ ConstructFrameElement("CarbonElement", fullFrameBoxVol, xu / 2.);
+ TGeoRotation* fullFrameRot = new TGeoRotation;
+ fullFrameRot->RotateY(180);
+
+ Int_t inum = YnumOfFrameBoxes - 3.5; // - 3.5 was done to short frame box to remove overlaps
+ for (i = 1; i <= inum; i++) {
+ j = -(inum - 1) / 2. + (i - 1);
+ // -(10-1)/2. +0 +10-1 -> -4.5 .. +4.5 -> -0.5, +0.5 (= 2)
+ // -(11-1)/2. +0 +11-1 -> -5.0 .. +5.0 -> -1, 0, 1 (= 3)
+ // cout << "DE: i " << i << " j " << j << endl;
+
+ if (LadderIndex % 100 <= 3) // central ladders in stations 1 to 3
+ {
+ if ((j >= -1) && (j <= 1)) // keep the inner 2 (even) or 3 (odd) elements free for the cone
+ continue;
+ }
+ else if (LadderIndex % 100 <= 8) // central ladders in stations 4 to 8
+ {
+ if ((j >= -2) && (j <= 2)) // keep the inner 4 elements free for the cone
+ continue;
+ }
+
+ cout << "DELZ: ladderDZ " << ladderDZ << " cm " << -ladderZ / 2. - ladderDZ << " cm " << endl;
+ ladder->AddNode(fullFrameBoxVol, i,
+ new TGeoCombiTrans(name + "_CarbonElement_posrot", 0., j * gkFrameStep, -(ladderZ / 2. + ladderDZ),
+ fullFrameRot));
+ }
+
+ ladder->GetShape()->ComputeBBox();
+}
+/** ======================================================================= **/
+
+
+/** ===========================================================================
+ ** Construct a ladder out of two half ladders with vertical gap
+ **
+ ** The second half ladder will be rotated by 180 degrees
+ ** in the x-y plane. The two half ladders will be put on top of each
+ ** other with a vertical gap.
+ **
+ ** Arguments:
+ ** name volume name
+ ** halfLadderU pointer to upper half ladder
+ ** halfLadderD pointer to lower half ladder
+ ** gapY vertical gap
+ ** shiftZ relative displacement along the z axis
+ **/
+
+TGeoVolume* ConstructLadder(Int_t LadderIndex, TGeoVolume* halfLadderU, TGeoVolume* halfLadderD, Double_t gapY,
+ Double_t shiftZ, Double_t pitchZ, Int_t nSectors)
+{
+
+ // --- Some variables
+ TGeoBBox* shape = NULL;
+
+ // define additional offset to carbon ladder for half ladders with less than 5 sensors
+ Double_t offsetZ = (5. - nSectors) * pitchZ;
+
+ // --- Dimensions of half ladders
+ shape = (TGeoBBox*) halfLadderU->GetShape(); // up
+ Double_t xu = 2. * shape->GetDX();
+ Double_t yu = 2. * shape->GetDY();
+ Double_t zu = 2. * shape->GetDZ();
+
+ shape = (TGeoBBox*) halfLadderD->GetShape(); // down
+ Double_t xd = 2. * shape->GetDX();
+ Double_t yd = 2. * shape->GetDY();
+ Double_t zd = 2. * shape->GetDZ();
+
+ // --- Create ladder volume assembly
+ TString name = Form("LadderType%04d", LadderIndex); // v19k
+ TGeoVolumeAssembly* ladder = new TGeoVolumeAssembly(name);
+ Double_t ladderX = TMath::Max(xu, xd);
+ Double_t ladderY = yu + yd + gapY;
+ Double_t ladderZ = TMath::Max(zu,
+ zd + shiftZ + offsetZ); // there are 6 slots - 5 x 1.5 mm + 0.3 mm = 7.8 mm
+ // Double_t ladderZ = TMath::Max(zu, zd + shiftZ);
+
+ cout << "DERR iladder " << LadderIndex << " nSec " << nSectors << " zu " << zu << " zd " << zd << " zd+shi "
+ << zd + shiftZ << " ladderZ " << ladderZ << " offsetZ " << offsetZ << endl;
+
+ // --- Place half ladders
+ Double_t xPosU = 0.0; // centred in x
+ Double_t yPosU = 0.5 * (ladderY - yu); // top aligned
+ Double_t zPosU = 0;
+ zPosU = 0.5 * (ladderZ - zu); // front aligned
+ if (LadderIndex >= 1000) zPosU = -0.5 * (ladderZ - zu) + offsetZ; // back aligned with possible offset
+ //zPosU = -zPosU;
+
+ TGeoTranslation* tu = new TGeoTranslation("tu", xPosU, yPosU, zPosU);
+ ladder->AddNode(halfLadderU, 1, tu);
+
+ Double_t xPosD = 0.0; // centred in x
+ Double_t yPosD = 0.5 * -(ladderY - yd); // bottom aligned
+ Double_t zPosD = 0;
+ zPosD = 0.5 * -(ladderZ - zd) + offsetZ; // back aligned with possible offset
+ if (LadderIndex >= 1000) zPosD = -0.5 * -(ladderZ - zd); // front aligned
+ //zPosD = -zPosD;
+
+ TGeoRotation* rd = new TGeoRotation();
+ rd->RotateZ(180.);
+ TGeoCombiTrans* cd = new TGeoCombiTrans(xPosD, yPosD, zPosD, rd);
+ ladder->AddNode(halfLadderD, 2, cd);
+
+ ladder->GetShape()->ComputeBBox();
+
+ shape = (TGeoBBox*) ladder->GetShape();
+ cout << "DDDD0ladder" << LadderIndex << " ladderX " << 2. * shape->GetDX() << " ladderY " << 2. * shape->GetDY()
+ << " ladderZ " << 2. * shape->GetDZ() << endl;
+
+ cout << "DDDD ladder" << LadderIndex << endl;
+ cout << "DDDD1ladder" << LadderIndex << " ladderX " << ladderX << " ladderY " << ladderY << " ladderZ " << ladderZ
+ << endl;
+
+ // ---------------- Create and place frame boxes ------------------------
+
+ if (gkConstructFrames) AddCarbonLadder(LadderIndex, ladder, ladderX, ladderY, ladderZ);
+
+ // --------------------------------------------------------------------------
+
+ shape = (TGeoBBox*) ladder->GetShape();
+ cout << "DDDD2ladder" << LadderIndex << " ladderX " << 2. * shape->GetDX() << " ladderY " << 2. * shape->GetDY()
+ << " ladderZ " << 2. * shape->GetDZ() << endl;
+
+ return ladder;
+}
+/** ======================================================================= **/
+
+
+/** ===========================================================================
+ ** Construct a unit
+ **
+ ** The unit volume is the minimal box comprising all ladders
+ ** minus a tube accomodating the beam pipe.
+ **
+ ** The ladders are arranged horizontally from left to right with
+ ** a given overlap in x.
+ ** Every second ladder is slightly displaced upstream from the centre
+ ** z plane and facing downstream, the others are slightly displaced
+ ** downstream and facing upstream (rotated around the y axis).
+ **
+ ** Arguments:
+ ** name volume name
+ ** nLadders number of ladders
+ ** ladderTypes array of ladder types
+ **/
+
+TGeoVolume* ConstructUnit(Int_t iSide, Int_t iUnit, Int_t nLadders, Int_t* ladderTypes, Int_t iStation)
+{
+
+ Bool_t isFirstPartOfHalfUnit = kFALSE;
+
+ // TString name = Form("Unit%02d", iUnit); // 0,1,2,3,4,5,6,7 - Unit00 missing in output
+ // TString name = Form("Unit%02d", iUnit+1); // 1,2,3,4,5,6,7,8
+
+ TGeoVolume* unit = gGeoMan->GetVolume(unitName[iUnit]);
+ if (!unit) // if it does not yet exist, create a new one
+ {
+ unit = new TGeoVolumeAssembly(unitName[iUnit]);
+ isFirstPartOfHalfUnit = kTRUE;
+ }
+
+ // --- Some local variables
+ TGeoBBox* ladderShape = NULL;
+ TGeoVolume* ladder = NULL;
+ TString ladderName;
+ Double_t subtractedVal;
+
+ // --- Determine size of unit from ladders
+ Double_t statX = 0.;
+ // Double_t statY = 0.;
+
+ for (Int_t iLadder = 0; iLadder < nLadders; iLadder++) {
+ Int_t ladderType = ladderTypes[iLadder]; // v19k
+
+ // if (iSide == 0) cout << "DWER " << ladderTypes[iLadder] << " " << ladderType << endl;
+
+ if (ladderType > 0) {
+ ladderName = Form("LadderType%04d", ladderType); // v19k
+
+ ladder = gGeoManager->GetVolume(ladderName);
+ if (!ladder) Fatal("ConstructUnit", Form("Volume %s not found", ladderName.Data()));
+ ladderShape = (TGeoBBox*) ladder->GetShape();
+ statX += 2. * ladderShape->GetDX();
+ }
+ else
+ statX += gkSensorSizeX; // empty ladder in unit
+ }
+ statX -= Double_t(nLadders - 1) * gkLadderOverlapX;
+
+ // if (iSide == 0) cout << "DWER -" << endl;
+
+ // --- Place ladders in unit
+ cout << "xPos0: " << statX << endl;
+ Double_t xPos = -0.5 * statX;
+ cout << "xPos1: " << xPos << endl;
+ Double_t yPos = 0.;
+ Double_t zPos = 0.;
+
+ Double_t maxdz = 0.;
+ for (Int_t iLadder = 0; iLadder < nLadders; iLadder++) // find maximum dz in this unit
+ {
+ Int_t ladderType = ladderTypes[iLadder]; // v19k
+
+ if (ladderType > 0) {
+ ladderName = Form("LadderType%04d", ladderType); // v19k
+
+ ladder = gGeoManager->GetVolume(ladderName);
+ ladderShape = (TGeoBBox*) ladder->GetShape();
+ if (maxdz < ladderShape->GetDZ()) maxdz = ladderShape->GetDZ();
+ }
+ }
+
+ for (Int_t iLadder = 0; iLadder < nLadders; iLadder++) {
+ Int_t ladderType = ladderTypes[iLadder]; // v19k
+
+ if (ladderType > 0) {
+ ladderName = Form("LadderType%04d", ladderType); // v19k
+
+ ladder = gGeoManager->GetVolume(ladderName);
+ ladderShape = (TGeoBBox*) ladder->GetShape();
+ xPos += ladderShape->GetDX();
+ cout << "xPos2: " << xPos << endl;
+ yPos = 0.; // vertically centred
+ TGeoRotation* rot = new TGeoRotation();
+
+ if (gkConstructFrames)
+ subtractedVal = ladderShape->GetDX() + sqrt(2.) * gkFrameThickness / 2. + gkSectorGapZFrame * 2;
+ else
+ subtractedVal = 0.;
+
+ zPos = 0.5 * gkLadderGapZ + (2 * maxdz - ladderShape->GetDZ() - subtractedVal / 2.); // z-aligned ladders
+
+ // v19k
+ cout << "DE ladder" << ladderTypes[iLadder] << " dx: " << ladderShape->GetDX()
+ << " dy: " << ladderShape->GetDY() << " dz: " << ladderShape->GetDZ() << " max dz: " << maxdz << endl;
+
+ // v19k
+ cout << "DE ladder" << ladderTypes[iLadder] << " fra: " << gkFrameThickness / 2. << " sub: " << subtractedVal
+ << " zpo: " << zPos << endl
+ << endl;
+
+ // v19k
+ if (ladderTypes[iLadder] / 1000 == 1) // flip some of the ladders to reproduce the CAD layout
+ rot->RotateY(180.);
+ else
+ zPos = -zPos;
+
+ if (!isFirstPartOfHalfUnit) zPos += 10.5; // v19d
+ // zPos += 10.0; // initial version
+
+ TGeoCombiTrans* trans = new TGeoCombiTrans(xPos, yPos, zPos, rot);
+ // start
+ // cout << "DEEE** iLadder " << iLadder << " " << nLadders/2 << " " << nLadders << endl;
+
+ if (iSide == 0) {
+ if (iLadder < nLadders / 2) // right side - only half unit -x
+ {
+ unit->AddNode(ladder, iStation * 100 + iLadder + 1, trans);
+ // calculate Station number to encode the ladder copy number
+ cout << "DEFG** iLadder: " << iLadder << " iStation: " << iStation << endl;
+ }
+ }
+ else {
+ if (iLadder >= nLadders / 2) // left side - only half unit +x
+ {
+ unit->AddNode(ladder, iStation * 100 + iLadder + 1, trans);
+ // calculate Station number to encode the ladder copy number
+ cout << "DEFG** iLadder: " << iLadder << " iStation: " << iStation << endl;
+ }
+ }
+ unit->GetShape()->ComputeBBox();
+ // stop
+ xPos += ladderShape->GetDX() - gkLadderOverlapX;
+ cout << "xPos3: " << xPos << endl;
+ }
+ else
+ xPos += gkSensorSizeX - gkLadderOverlapX;
+ }
+
+ return unit;
+}
+/** ======================================================================= **/
+
+
+/** ===========================================================================
+ ** Import and add the passive materials to the STS volume
+ **/
+void ImportPassive(TGeoVolume* stsVolume, TString geoTag, fstream& infoFile)
+{
+ TString passiveName = TString("sts_passive_") + geoTag;
+ TString basePath = gSystem->Getenv("VMCWORKDIR");
+ TString relPath = "/geometry/sts/passive/" + passiveName + ".gdml";
+ TString passiveFileName = basePath + relPath;
+ infoFile << std::endl << std::endl;
+ infoFile << "Importing STS passive materials from GDML file '" << relPath << "'." << std::endl;
+
+ TGDMLParse parser;
+ TGeoVolume* gdmlVolume = parser.GDMLReadFile(passiveFileName);
+ PostProcessGdml(gdmlVolume);
+ gdmlVolume->SetName(passiveName);
+
+ TGeoTranslation* passiveTrans = new TGeoTranslation(0., 0., 4.68 - 2.);
+ infoFile << "Passive assembly is translated for Z=2.68 cm downstream with "
+ "respect to parent volume"
+ << std::endl
+ << std::endl;
+
+ gdmlVolume->GetShape()->ComputeBBox();
+ CheckVolume(gdmlVolume, infoFile);
+
+ infoFile << std::endl;
+ for (Int_t iNode = 0; iNode < gdmlVolume->GetNdaughters(); iNode++) {
+ CheckVolume(gdmlVolume->GetNode(iNode)->GetVolume(), infoFile, kFALSE);
+ }
+
+ stsVolume->AddNode(gdmlVolume, stsVolume->GetNdaughters(), passiveTrans, "");
+}
+
+/** ===========================================================================
+ ** Assign visual properties to the imported gdml volumes
+ **/
+void PostProcessGdml(TGeoVolume* gdmlVolume)
+{
+ const UInt_t kPOBColor = kRed - 6;
+ const UInt_t kPOBTransparency = 0; // 5;
+
+ const UInt_t kFEBColor = kOrange - 6;
+ const UInt_t kFEBTransparency = 0; // 5;
+
+ const UInt_t kUnitColor = kCyan - 10;
+ const UInt_t kUnitTransparency = 0; // 5;
+
+ const UInt_t kCfColor = kCyan - 1;
+ const UInt_t kCfTransparency = 0; // 10;
+
+ // name <Color, Transparency>
+ std::map<std::string, std::tuple<UInt_t, UInt_t>> props {
+ {"passive_POB", std::tuple<UInt_t, UInt_t> {kPOBColor, kPOBTransparency}},
+ {"passive_FEB", std::tuple<UInt_t, UInt_t> {kFEBColor, kFEBTransparency}},
+ {"passive_unit", std::tuple<UInt_t, UInt_t> {kUnitColor, kUnitTransparency}},
+ {"passive_Box_Wall", std::tuple<UInt_t, UInt_t> {kCfColor, kCfTransparency}},
+ {"passive_Box_Wall_Front_CF2", std::tuple<UInt_t, UInt_t> {kCfColor + 3, kCfTransparency}},
+ {"passive_Box_Wall_Back_CF2", std::tuple<UInt_t, UInt_t> {kCfColor + 3, kCfTransparency}},
+ {"passive_Box_Wall_Back_Light_Foam", std::tuple<UInt_t, UInt_t> {kCfColor - 8, kCfTransparency}},
+ {"passive_Box_Wall_Front_Light_Foam", std::tuple<UInt_t, UInt_t> {kCfColor - 8, kCfTransparency}},
+ };
+
+ // Match volume name and apply visual properties passive_Box_Wall_Front_Light_Foam
+ const TObjArray* volumes = gGeoManager->GetListOfVolumes();
+ for (auto& entry : props) {
+ TIter next(volumes);
+ TGeoVolume* vol = nullptr;
+ while ((vol = (TGeoVolume*) next())) {
+ if (TString(vol->GetName()).Contains(entry.first.c_str())) {
+ vol->SetLineColor(std::get<0>(entry.second));
+ vol->SetTransparency(std::get<1>(entry.second));
+ }
+ }
+ }
+}
+
+/** ===========================================================================
+ ** Volume information for debugging
+ **/
+void CheckVolume(TGeoVolume* volume)
+{
+
+ TGeoBBox* shape = (TGeoBBox*) volume->GetShape();
+ cout << volume->GetName() << ": size " << fixed << setprecision(4) << setw(7) << 2. * shape->GetDX() << " x "
+ << setw(7) << 2. * shape->GetDY() << " x " << setw(7) << 2. * shape->GetDZ();
+ if (volume->IsAssembly()) cout << ", assembly";
+ else {
+ if (volume->GetMedium()) cout << ", medium " << volume->GetMedium()->GetName();
+ else
+ cout << ", "
+ << "\033[31m"
+ << " no medium"
+ << "\033[0m";
+ }
+ cout << endl;
+ if (volume->GetNdaughters()) {
+ cout << "Daughters: " << endl;
+ for (Int_t iNode = 0; iNode < volume->GetNdaughters(); iNode++) {
+ TGeoNode* node = volume->GetNode(iNode);
+ TGeoBBox* shape = (TGeoBBox*) node->GetVolume()->GetShape();
+ cout << setw(15) << node->GetName() << ", size " << fixed << setprecision(3) << setw(6) << 2. * shape->GetDX()
+ << " x " << setw(6) << 2. * shape->GetDY() << " x " << setw(6) << 2. * shape->GetDZ() << ", position ( ";
+ TGeoMatrix* matrix = node->GetMatrix();
+ const Double_t* pos = matrix->GetTranslation();
+ cout << setfill(' ');
+ cout << fixed << setw(8) << pos[0] << ", " << setw(8) << pos[1] << ", " << setw(8) << pos[2] << " )" << endl;
+ }
+ }
+}
+/** ======================================================================= **/
+
+/** ===========================================================================
+ ** Volume information for output to file
+ **/
+void CheckVolume(TGeoVolume* volume, fstream& file, Bool_t listChildren)
+{
+ if (!file) return;
+
+ TGeoBBox* shape = (TGeoBBox*) volume->GetShape();
+ file << volume->GetName() << ": size " << fixed << setprecision(4) << setw(7) << 2. * shape->GetDX() << " x "
+ << setw(7) << 2. * shape->GetDY() << " x " << setw(7) << 2. * shape->GetDZ();
+ if (volume->IsAssembly()) file << ", assembly";
+ else {
+ if (volume->GetMedium()) file << ", medium " << volume->GetMedium()->GetName();
+ else
+ file << ", "
+ << "\033[31m"
+ << " no medium"
+ << "\033[0m";
+ }
+ file << endl;
+ if (volume->GetNdaughters() && listChildren) {
+ file << "Contains: ";
+ for (Int_t iNode = 0; iNode < volume->GetNdaughters(); iNode++)
+ file << volume->GetNode(iNode)->GetVolume()->GetName() << " ";
+ file << endl;
+ }
+}
+
+/** ======================================================================= **/
+
+
+/** ===========================================================================
+ ** Calculate beam pipe outer radius for a given z
+ **/
+Double_t BeamPipeRadius(Double_t z)
+{
+ if (z < gkPipeZ2) return gkPipeR1;
+ Double_t slope = (gkPipeR2 - gkPipeR1) / (gkPipeZ2 - gkPipeZ1);
+ return gkPipeR2 + slope * (z - gkPipeZ2);
+}
+/** ======================================================================= **/
+
+
+/** ======================================================================= **/
+TGeoVolume* ConstructFrameElement(const TString& name, TGeoVolume* frameBoxVol, Double_t x)
+{
+ // --- Material of the frames
+ TGeoMedium* framesMaterial = gGeoMan->GetMedium("carbon");
+
+ TGeoBBox* frameVertPillarShp;
+
+ Double_t t = gkFrameThickness / 2.;
+
+ // --- Main vertical pillars
+ // TGeoBBox* frameVertPillarShp = new TGeoBBox(name + "_vertpillar_shape", t, gkFrameStep/2., t); // square crossection, along y
+ // TGeoVolume* frameVertPillarVol = new TGeoVolume(name + "_vertpillar", frameVertPillarShp, framesMaterial);
+ // frameVertPillarVol->SetLineColor(kGreen);
+ // frameBoxVol->AddNode(frameVertPillarVol, 1, new TGeoTranslation(name + "_vertpillar_pos_1", x-t, 0., -(x+sqrt(2.)*t-2.*t)/2.));
+ // frameBoxVol->AddNode(frameVertPillarVol, 2, new TGeoTranslation(name + "_vertpillar_pos_2", -(x-t), 0., -(x+sqrt(2.)*t-2.*t)/2.));
+
+ if (gkCylindricalFrames)
+ // TGeoBBox* frameVertPillarShp = new TGeoTube(name + "_vertpillar_shape", 0, t, gkFrameStep/2.); // circle crossection, along z
+ frameVertPillarShp = new TGeoTube(name + "_vertpillar_shape", gkCylinderDiaInner / 2., gkCylinderDiaOuter / 2.,
+ gkFrameStep / 2.); // circle crossection, along z
+ else
+ frameVertPillarShp = new TGeoBBox(name + "_vertpillar_shape", t, t,
+ gkFrameStep / 2.); // square crossection, along z
+ TGeoVolume* frameVertPillarVol = new TGeoVolume(name + "_vertpillar", frameVertPillarShp, framesMaterial);
+ frameVertPillarVol->SetLineColor(kGreen);
+
+ TGeoRotation* xRot90 = new TGeoRotation;
+ xRot90->RotateX(90.);
+ frameBoxVol->AddNode(
+ frameVertPillarVol, 1,
+ new TGeoCombiTrans(name + "_vertpillar_pos_1", x - t, 0., -(x + sqrt(2.) * t - 2. * t) / 2., xRot90));
+ frameBoxVol->AddNode(
+ frameVertPillarVol, 2,
+ new TGeoCombiTrans(name + "_vertpillar_pos_2", -(x - t), 0., -(x + sqrt(2.) * t - 2. * t) / 2., xRot90));
+
+ // TGeoRotation* vertRot = new TGeoRotation(name + "_vertpillar_rot_1", 90., 45., -90.);
+ TGeoRotation* vertRot = new TGeoRotation;
+ vertRot->RotateX(90.);
+ vertRot->RotateY(45.);
+ frameBoxVol->AddNode(frameVertPillarVol, 3,
+ new TGeoCombiTrans(name + "_vertpillar_pos_3", 0., 0., (x - sqrt(2.) * t) / 2., vertRot));
+
+ // --- Small horizontal pillar
+ // TGeoBBox* frameHorPillarShp = new TGeoBBox(name + "_horpillar_shape", x-2.*t, gkThinFrameThickness/2., gkThinFrameThickness/2.);
+ // TGeoVolume* frameHorPillarVol = new TGeoVolume(name + "_horpillar", frameHorPillarShp, framesMaterial);
+ // frameHorPillarVol->SetLineColor(kCyan);
+ // frameBoxVol->AddNode(frameHorPillarVol, 1, new TGeoTranslation(name + "_horpillar_pos_1", 0., -gkFrameStep/2.+gkThinFrameThickness/2., -(x+sqrt(2.)*t-2.*t)/2.));
+
+ if (gkConstructSmallFrames) {
+
+ // --- Small sloping pillars
+ TGeoPara* frameSlopePillarShp =
+ new TGeoPara(name + "_slopepillar_shape", (x - 2. * t) / TMath::Cos(31.4 / 180. * TMath::Pi()),
+ gkThinFrameThickness / 2., gkThinFrameThickness / 2., 31.4, 0., 90.);
+ TGeoVolume* frameSlopePillarVol = new TGeoVolume(name + "_slopepillar", frameSlopePillarShp, framesMaterial);
+ frameSlopePillarVol->SetLineColor(kCyan);
+ TGeoRotation* slopeRot = new TGeoRotation(name + "_slopepillar_rot_1", 0., 0., 31.4);
+ TGeoRotation* slopeRot2 = new TGeoRotation(name + "_slopepillar_rot_2", 0., 0., -31.4);
+ TGeoCombiTrans* slopeTrRot =
+ new TGeoCombiTrans(name + "_slopepillar_posrot_1", 0., 0., -(x + sqrt(2.) * t - 2. * t) / 2., slopeRot);
+ TGeoCombiTrans* slopeTrRot2 =
+ new TGeoCombiTrans(name + "_slopepillar_posrot_2", 0., 0., -(x + sqrt(2.) * t - 2. * t) / 2., slopeRot2);
+
+ frameBoxVol->AddNode(frameSlopePillarVol, 1, slopeTrRot);
+ frameBoxVol->AddNodeOverlap(frameSlopePillarVol, 2, slopeTrRot2);
+
+
+ Double_t angl = 23.;
+ // --- Small sub pillar
+ TGeoPara* frameSubPillarShp =
+ new TGeoPara(name + "_subpillar_shape", (sqrt(2) * (x / 2. - t) - t / 2.) / TMath::Cos(angl / 180. * TMath::Pi()),
+ gkThinFrameThickness / 2., gkThinFrameThickness / 2., angl, 0., 90.);
+ TGeoVolume* frameSubPillarVol = new TGeoVolume(name + "_subpillar", frameSubPillarShp, framesMaterial);
+ frameSubPillarVol->SetLineColor(kMagenta);
+
+ Double_t posZ = t * (1. - 3. / (2. * sqrt(2.)));
+
+ // one side of X direction
+ TGeoRotation* subRot1 = new TGeoRotation(name + "_subpillar_rot_1", 90., 45., -90. + angl);
+ TGeoCombiTrans* subTrRot1 =
+ new TGeoCombiTrans(name + "_subpillar_posrot_1", -(-x / 2. + t - t / (2. * sqrt(2.))), 1., posZ, subRot1);
+
+ TGeoRotation* subRot2 = new TGeoRotation(name + "_subpillar_rot_2", 90., -90. - 45., -90. + angl);
+ TGeoCombiTrans* subTrRot2 =
+ new TGeoCombiTrans(name + "_subpillar_posrot_2", -(-x / 2. + t - t / (2. * sqrt(2.))), -1., posZ, subRot2);
+
+ // other side of X direction
+ TGeoRotation* subRot3 = new TGeoRotation(name + "_subpillar_rot_3", 90., 90. + 45., -90. + angl);
+ TGeoCombiTrans* subTrRot3 =
+ new TGeoCombiTrans(name + "_subpillar_posrot_3", -x / 2. + t - t / (2. * sqrt(2.)), 1., posZ, subRot3);
+
+ TGeoRotation* subRot4 = new TGeoRotation(name + "_subpillar_rot_4", 90., -45., -90. + angl);
+ TGeoCombiTrans* subTrRot4 =
+ new TGeoCombiTrans(name + "_subpillar_posrot_4", -x / 2. + t - t / (2. * sqrt(2.)), -1., posZ, subRot4);
+
+ frameBoxVol->AddNode(frameSubPillarVol, 1, subTrRot1);
+ frameBoxVol->AddNode(frameSubPillarVol, 2, subTrRot2);
+ frameBoxVol->AddNode(frameSubPillarVol, 3, subTrRot3);
+ frameBoxVol->AddNode(frameSubPillarVol, 4, subTrRot4);
+ // frameBoxVol->GetShape()->ComputeBBox();
+ }
+
+ return frameBoxVol;
+}
+/** ======================================================================= **/
+
+/** ======================================================================= **/
+TGeoVolume* ConstructSmallCone(Double_t coneDz)
+{
+ // --- Material of the frames
+ TGeoMedium* framesMaterial = gGeoMan->GetMedium("carbon");
+
+ // --- Outer cone
+ // TGeoConeSeg* A = new TGeoConeSeg ("A", coneDz, 6., 7.6, 6., 6.04, 0., 180.);
+ // TGeoBBox* B = new TGeoBBox ("B", 8., 6., 10.);
+
+ Double_t radius = 3.0;
+ Double_t thickness = 0.04; // 0.4 mm
+ // TGeoConeSeg* A = new TGeoConeSeg ("A", coneDz, 3., 3.2, 3., 3.2, 0., 180.);
+ TGeoConeSeg* A = new TGeoConeSeg("A", coneDz, radius, radius + thickness, radius, radius + thickness, 0., 180.);
+ TGeoBBox* B = new TGeoBBox("B", 8., 6., 10.);
+
+ TGeoCombiTrans* M = new TGeoCombiTrans("M");
+ M->RotateX(45.);
+ M->SetDy(-5.575);
+ M->SetDz(6.935);
+ M->RegisterYourself();
+
+ TGeoShape* coneShp = new TGeoCompositeShape("Cone_shp", "A-B:M");
+ TGeoVolume* coneVol = new TGeoVolume("Cone", coneShp, framesMaterial);
+ coneVol->SetLineColor(kGreen);
+ // coneVol->RegisterYourself();
+
+ // // --- Inner cone
+ // Double_t thickness = 0.02;
+ // Double_t thickness2 = 0.022;
+ // // TGeoConeSeg* A2 = new TGeoConeSeg ("A2", coneDz-thickness, 6.+thickness, 7.6-thickness2, 5.99+thickness, 6.05-thickness2, 0., 180.);
+ // TGeoConeSeg* A2 = new TGeoConeSeg ("A2", coneDz-thickness, 3.+thickness, 4.6-thickness2, 2.99+thickness, 3.05-thickness2, 0., 180.);
+ //
+ // TGeoCombiTrans* M2 = new TGeoCombiTrans ("M2");
+ // M2->RotateX (45.);
+ // M2->SetDy (-5.575+thickness*sqrt(2.));
+ // M2->SetDz (6.935);
+ // M2->RegisterYourself();
+ //
+ // TGeoShape* coneShp2 = new TGeoCompositeShape ("Cone2_shp", "A2-B:M2");
+ // TGeoVolume* coneVol2 = new TGeoVolume ("Cone2", coneShp2, gStsMedium);
+ // coneVol2->SetLineColor(kGreen);
+ //// coneVol2->RegisterYourself();
+ //
+ // coneVol->AddNode(coneVol2, 1);
+
+ return coneVol;
+}
+/** ======================================================================= **/
+
+/** ======================================================================= **/
+TGeoVolume* ConstructBigCone(Double_t coneDz)
+{
+ // --- Material of the frames
+ TGeoMedium* framesMaterial = gGeoMan->GetMedium("carbon");
+
+ // --- Outer cone
+ TGeoConeSeg* bA = new TGeoConeSeg("bA", coneDz, 6., 7.6, 6., 6.04, 0., 180.);
+ TGeoBBox* bB = new TGeoBBox("bB", 8., 6., 10.);
+
+ TGeoCombiTrans* bM = new TGeoCombiTrans("bM");
+ bM->RotateX(45.);
+ bM->SetDy(-5.575);
+ bM->SetDz(6.935);
+ bM->RegisterYourself();
+
+ TGeoShape* coneBigShp = new TGeoCompositeShape("ConeBig_shp", "bA-bB:bM");
+ TGeoVolume* coneBigVol = new TGeoVolume("ConeBig", coneBigShp, framesMaterial);
+ coneBigVol->SetLineColor(kGreen);
+ // coneBigVol->RegisterYourself();
+
+ // --- Inner cone
+ Double_t thickness = 0.02;
+ Double_t thickness2 = 0.022;
+ TGeoConeSeg* bA2 = new TGeoConeSeg("bA2", coneDz - thickness, 6. + thickness, 7.6 - thickness2, 5.99 + thickness,
+ 6.05 - thickness2, 0., 180.);
+
+ TGeoCombiTrans* bM2 = new TGeoCombiTrans("bM2");
+ bM2->RotateX(45.);
+ bM2->SetDy(-5.575 + thickness * sqrt(2.));
+ bM2->SetDz(6.935);
+ bM2->RegisterYourself();
+
+ TGeoShape* coneBigShp2 = new TGeoCompositeShape("ConeBig2_shp", "bA2-bB:bM2");
+ TGeoVolume* coneBigVol2 = new TGeoVolume("ConeBig2", coneBigShp2, gStsMedium);
+ coneBigVol2->SetLineColor(kGreen);
+ // coneBigVol2->RegisterYourself();
+
+ coneBigVol->AddNode(coneBigVol2, 1);
+
+ return coneBigVol;
+}
+
+/** ======================================================================= **/
diff --git a/macro/sts/geometry/create_stsgeo_v21g.C b/macro/sts/geometry/create_stsgeo_v21g.C
new file mode 100644
index 0000000000..73e85ea83a
--- /dev/null
+++ b/macro/sts/geometry/create_stsgeo_v21g.C
@@ -0,0 +1,2197 @@
+/* Copyright (C) 2012-2021 Goethe-Universitaet Frankfurt, Frankfurt
+ SPDX-License-Identifier: GPL-3.0-only
+ Authors: Mehul Shiroya [committer], Volker Friese, Evgeny Lavrik*/
+
+/******************************************************************************
+ ** Creation of STS geometry in ROOT format (TGeo).
+ **
+ ** @file create_stsgeo_v21g.C
+ ** @author Volker Friese <v.friese@gsi.de>
+ ** @since 15 June 2012
+ ** @date 09.05.2014
+ ** @author Tomas Balog <T.Balog@gsi.de>
+ **
+ ** v21g: Similar to the geometry version v21e
+ * [ Airex foam in the front and back side wall (Solid).
+ * New beam pipe R = 2.0 cm upstream-side with thickness 1mm.
+ * Central ladders adjusted 5mm away from the beampipe. Global origin is taken from center of magnet.
+ * Origin of the STS is taken from last station which is located at(0,0,0) position.
+ * Global origin is considered from the center of the magnet.
+ * Nomenclature for the Halfladders has been changed from the upside(u) to top (t)
+ * and from downside(d) to bottom (b). ]
+ * The only difference is that it includes beampipe flange in the back side of the wall,
+ * which is considered to be made of carbon fiber(70%) and epoxy resin(30%).
+ **
+ ** v19r: bugfix of v19q - align all halfladders
+ ** v19q: based on v19l - align ladders to virtual plane in station center, closing the gaps in z
+ ** v19p: based on v19k - parameters : delta Z prime = 0.70 cm - delta Z pitch = 0.20 cm
+ ** v19o: based on v19k - parameters : delta Z prime = 0.30 cm - delta Z pitch = 0.20 cm
+ ** v19n: based on v19k - parameters : delta Z prime = 0.50 cm - delta Z pitch = 0.20 cm (bug fix of v19m)
+ ** v19m: based on v19k - parameters : delta Z prime = 0.55 cm - delta Z pitch = 0.20 cm (bug)
+ ** v19l: based on v19k - parameters : delta Z prime = 0.50 cm - delta Z pitch = 0.15 cm
+ ** v19k: ladders on upstream side of units get upper half ladders installed first,
+ ** ladders on downstream side of units get lower half ladders installed first,
+ ** this saves 1.5 mm space in z per station, 12 mm in total (LadderType went from 3 to 4 digits)
+ ** parameters : delta Z prime = 1.00 cm - delta Z pitch = 0.15 cm
+ ** v19j: use overlap and distance parameters from CAD model
+ ** v19h: put STS stations from v19d at z-positions = 260; 365; 470; 575; 680; 785; 890; 995 mm
+ ** v19g: place a box with services around v19e
+ ** v19f: place a box with services around v19d
+ ** v19e: increase spacing between stations by +10 mm from 100 mm
+ ** v19d: increase spacing between stations by + 5 mm from 100 mm
+ ** v19c: drop station 8 and increase spacing between remaining 7 stations from 10 cm to 12 c
+ ** v19b: introduce FEB orientation in ladder numbering (LadderType went from 2 to 3 digits)
+ ** v19a: import passive materials from gdml file
+ ** extend CF ladder structures and cables towards FEE plane
+ ** change CF ladder frame shape
+ ** v18d: increases thickness of sensors within a 7.5 degree cone from 300 mu to 400 mu (based on v18b)
+ ** v18c: fixed cut-out windows in cooling plates, improve the box shape/materials
+ ** v18b: increases thickness of sensors within a 7.5 degree cone from 300 mu to 400 mu
+ ** v18a: adds 9 cooling/holding plates and a box around the setup
+ ** v16g: v16g is the new standard geometry from November 2017
+ ** v16g: switch from stations to units - left / right ("Unit01L", "Unit01R")
+ ** v16f: switch from stations to units
+ ** - split in upstream / downstream and left / right parts
+ ** - named Unit0xUR, Unit0xUL, Unit0xDR, Unit0xDL
+ ** v16e: switch from stations to units - upstream / downstream ("Unit01U", "Unit01D")
+ ** v16d: skip keeping volumes of sts and stations
+ ** v16c: like v16b, but senors of ladders beampipe next to beampipe
+ ** shifted closer to the pipe, like in the CAD model
+ ** v16b: like v16a, but yellow sensors removed
+ ** v16a: derived from v15c (no cones), but with sensor types renamed:
+ ** 2 -> 1, 3 -> 2, 4 -> 3, 5 -> 4, 1 -> 5
+ **
+ ** v15c: as v15b without cones
+ ** v15b: introduce modified carbon ladders from v13z
+ ** v15a: with flipped ladder orientation for stations 0,2,4,6 to match CAD design
+ **
+ ** TODO:
+ **
+ ** DONE:
+ ** v15b - use carbon macaroni as ladder support
+ ** v15b - introduce a small gap between lowest sensor and carbon ladder
+ ** v15b - build small cones for the first 2 stations
+ ** v15b - within a station the ladders of adjacent units should not touch eachother - set gkLadderGapZ to 10 mm
+ ** v15b - for all ladders set an even number of ladder elements
+ ** v15b - z offset of cones to ladders should not be 0.3 by default, but 0.26
+ ** v15b - within a station the ladders should be aligned in z, defined either by the unit or the ladder with most sensors
+ ** v15b - get rid of cone overlap in stations 7 and 8 - done by adapting rHole size
+ **
+ ** The geometry hierarachy is:
+ **
+ ** 1. Sensors (see function CreateSensors)
+ ** The sensors are the active volumes and the lowest geometry level.
+ ** They are built as TGeoVolumes, shape box, material silicon.
+ ** x size is determined by strip pitch 58 mu and 1024 strips
+ ** plus guard ring of 1.3 mm at each border -> 6.1992 cm.
+ ** Sensor type 1 is half of that (3.0792 cm).
+ ** y size is determined by strip length (2.2 / 4.2 / 6.3 cm) plus
+ ** guard ring of 1.3 mm at top and bottom -> 2.46 / 4.46 / 6.46 cm.
+ ** z size is a parameter, to be set by gkSensorThickness.
+ **
+ ** 2. Sectors (see function CreateSectors)
+ ** Sectors consist of several chained sensors. These are arranged
+ ** vertically on top of each other with a gap to be set by
+ ** gkChainGapY. Sectors are constructed as TGeoVolumeAssembly.
+ ** The sectors are auxiliary volumes used for proper placement
+ ** of the sensor(s) in the module. They do not show up in the
+ ** final geometry.
+ **
+ ** 3. Modules (see function ConstructModule)
+ ** A module is a readout unit, consisting of one sensor or
+ ** a chain of sensors (see sector) and a cable.
+ ** The cable extends from the top of the sector vertically to the
+ ** top of the halfladder the module is placed in. The cable and module
+ ** volume thus depend on the vertical position of the sector in
+ ** the halfladder. The cables consist of silicon with a thickness to be
+ ** set by gkCableThickness.
+ ** Modules are constructed as TGeoVolume, shape box, medium gStsMedium.
+ ** The module construction can be switched off (gkConstructCables)
+ ** to reproduce older geometries.
+ **
+ ** 4. Halfladders (see function ConstructHalfLadder)
+ ** A halfladder is a vertical assembly of several modules. The modules
+ ** are placed vertically such that their sectors overlap by
+ ** gkSectorOverlapY. They are displaced in z direction to allow for the
+ ** overlap in y by gkSectorGapZ.
+ ** The horizontal placement of modules in the halfladder can be choosen
+ ** to left aligned or right aligned, which only matters if sensors of
+ ** different x size are involved.
+ ** Halfladders are constructed as TGeoVolumeAssembly.
+ **
+ ** 5. Ladders (see function CreateLadders and ConstructLadder)
+ ** A ladder is a vertical assembly of two halfladders, and is such the
+ ** vertical building block of a station. The second (bottom) half ladder
+ ** is rotated upside down. The vertical arrangement is such that the
+ ** inner sectors of the two halfladders have the overlap gkSectorOverlapY
+ ** (function CreateLadder) or that there is a vertical gap for the beam
+ ** hole (function CreateLadderWithGap).
+ ** Ladders are constructed as TGeoVolumeAssembly.
+ **
+ ** 6. Stations (see function ConstructStation)
+ ** A station represents one layer of the STS geometry: one measurement
+ ** at (approximately) a given z position. It consist of several ladders
+ ** arranged horizontally to cover the acceptance.
+ ** The ladders are arranged such that there is a horizontal overlap
+ ** between neighbouring ladders (gkLadderOverLapX) and a vertical gap
+ ** to allow for this overlap (gkLadderGapZ). Each second ladder is
+ ** rotated around its y axis to face away from or into the beam.
+ ** Stations are constructed as TGeoVolumes, shape box minus tube (for
+ ** the beam hole), material gStsMedium.
+ **
+ ** 7. STS
+ ** The STS is a volume hosting the entire detectors system. It consists
+ ** of several stations located at different z positions.
+ ** The STS is constructed as TGeoVolume, shape box minus cone (for the
+ ** beam pipe), material gStsMedium. The size of the box is computed to
+ ** enclose all stations.
+ *****************************************************************************/
+
+
+// Remark: With the proper steering variables, this should exactly reproduce
+// the geometry version v11b of A. Kotynia's described in the ASCII format.
+// The only exception is a minimal difference in the z position of the
+// sectors/sensors. This is because of ladder types 2 and 4 containing the half
+// sensors around the beam hole (stations 1,2 and 3). In v11b, the two ladders
+// covering the beam hole cannot be transformed into each other by rotations,
+// but only by a reflection. This means they are constructionally different.
+// To avoid introducing another two ladder types, the difference in z position
+// was accepted.
+
+
+// Differences to v12:
+// gkChainGap reduced from 1 mm to 0
+// gkCableThickness increased from 100 mum to 200 mum (2 cables per module)
+// gkSectorOverlapY reduced from 3 mm to 2.4 mm
+// New sensor types 05 and 06
+// New sector types 07 and 08
+// Re-definiton of ladders (17 types instead of 8)
+// Re-definiton of station from new ladders
+
+
+#include "TGeoCompositeShape.h"
+#include "TGeoCone.h"
+#include "TGeoManager.h"
+#include "TGeoPara.h"
+#include "TGeoPhysicalNode.h"
+#include "TGeoTrd2.h"
+#include "TGeoTube.h"
+#include "TGeoXtru.h"
+
+#include <iomanip>
+#include <iostream>
+
+// forward declarations
+Int_t CreateSensors();
+Int_t CreateSectors();
+Int_t CreateLadders();
+TGeoVolume* ConstructModule(const char* name, TGeoVolume* sector, Double_t cableLength);
+TGeoVolume* ConstructHalfLadder(const TString& name, Int_t nSectors, Int_t* sectorTypes, char align,
+ Double_t ladderLength, Double_t offsetY);
+TGeoVolume* ConstructLadder(Int_t LadderIndex, TGeoVolume* halfLadderU, TGeoVolume* halfLadderD, Double_t gapY,
+ Double_t shiftZ, Double_t pitchZ, Int_t nSectors);
+TGeoVolume* ConstructUnit(Int_t iSide, Int_t iUnit, Int_t nLadders, Int_t* ladderTypes, Int_t iStation);
+void ImportPassive(TGeoVolume* stsVolume, TString geoTag, fstream& infoFile);
+void PostProcessGdml(TGeoVolume* gdmlTop);
+void CheckVolume(TGeoVolume* volume);
+void CheckVolume(TGeoVolume* volume, fstream& file, Bool_t listChildren = kTRUE);
+Double_t BeamPipeRadius(Double_t z);
+TGeoVolume* ConstructFrameElement(const TString& name, TGeoVolume* frameBoxVol, Double_t x);
+TGeoVolume* ConstructSmallCone(Double_t coneDz);
+TGeoVolume* ConstructBigCone(Double_t coneDz);
+
+// ------------- Version highlight -----------------------------------
+
+const std::string gVersionHighlight = R"(
+Summary:
+ This version adds passive materials imported from GDML model to the STS geometry:
+ * Taken from and largely correspond to mechanical CAD drawings of the detector
+ * Thermal insulation box:
+ - made out of carbon sandwitch panel (2mm carbon fiber sheet + layer of carbon foam + 2mm carbon fiber sheet)
+ - front window of complex shape with interface to MVD / target chamber
+ - back window with large aperture (2000 x 1200 mm) square cut into carbon foam
+ * Structural units:
+ - made of 2 complex shape aluminum C-Frames, 15mm thick
+ - placed at 25, 35, ... ,105 cm absolute Z
+ - contain front-end and power distribution boxes with equivalent X_0 values
+
+ Scripted geometry tweaks:
+ * Ladders and cables are extended towards the read-out planes having same lengths in respective rows
+ * Adjusted form and shape of carbon ladder structures from L-type to X-type
+ * Reduced verbosity of this file
+
+ Sensor arrangement is the same as in version v16g
+
+ !! Important for this version is the discrepancy from the mechanical CAD w.r.t. front wall.
+ The square window was replaced by a round one to avoid overlaps with present beam pipe designs, e.g. pipe_v16b_1e
+)";
+
+// ------------- Steering variables -----------------------------------
+
+// ---> Horizontal width of sensors [cm]
+const Double_t gkSensorSizeX = 6.2; // was 6.2092; // 6.2 - Oleg CAD 15/05/2020
+
+// ---> Thickness of sensors [cm]
+const Double_t gkSensorThickness = 0.03;
+
+// ---> Vertical gap between chained sensors [cm]
+const Double_t gkChainGapY = 0.00;
+
+// ---> Thickness of cables [cm]
+const Double_t gkCableThickness = 0.02;
+
+// ---> Horizontal overlap of neighbouring ladders [cm]
+const Double_t gkLadderOverlapX = 0.25; // delta X - Oleg CAD 14/05/2020
+
+// ---> Vertical overlap of neighbouring sectors in a ladder [cm]
+const Double_t gkSectorOverlapY = 0.46; // delta Y - Oleg CAD 14/05/2020
+
+// ---> Gap in z between neighbouring sectors in a ladder [cm]
+const Double_t gkSectorGapZ = 0.12; // gap + thickness = pitch // delta Z pitch = 0.15 - Oleg CAD 14/05/2020
+
+// ---> Gap in z between neighbouring ladders [cm]
+const Double_t gkLadderGapZ = 0.50 - 0.15; // for asym // 0.5 for sym // delta Z prime
+
+// ---> Gap in z between lowest sector to carbon support structure [cm]
+const Double_t gkSectorGapZFrame =
+ 0.280
+ - 0.025; // Oleg CAD 05/05/2020 // there is a 2.8 mm gap between the bottom side of the sensor and the top ledge of the carbon ladder
+
+// ---> Switch to construct / not to construct readout cables
+const Bool_t gkConstructCables = kTRUE;
+
+// ---> Switch to construct / not to construct frames
+const Bool_t gkConstructCones = kFALSE; // kTRUE; // switch this false by default for v15c and v16x
+const Bool_t gkConstructFrames = kTRUE; // kFALSE; // switch this true by default for v15c and v16x
+const Bool_t gkConstructSmallFrames = kTRUE; // kFALSE;
+const Bool_t gkCylindricalFrames = kTRUE; // kFALSE;
+
+// ---> Size of the frame
+const Double_t gkFrameThickness = 0.2;
+const Double_t gkThinFrameThickness = 0.05;
+const Double_t gkFrameStep = 4.0; // size of frame cell along y direction
+const Double_t gkCylinderDiaInner =
+ 0.07; // properties of cylindrical carbon supports, see CBM-STS Integration Meeting (10 Jul 2015)
+const Double_t gkCylinderDiaOuter =
+ 0.15; // properties of cylindrical carbon supports, see CBM-STS Integration Meeting (10 Jul 2015)
+
+// ---> Switch to import / not to import the Passive materials from GDML file
+const Bool_t gkImportPassive = kTRUE;
+
+// ----------------------------------------------------------------------------
+
+
+// -------------- Parameters of beam pipe in the STS region --------------
+// ---> Needed to compute stations and STS such as to avoid overlaps
+const Double_t gkPipeZ1 = 18.0;
+const Double_t gkPipeR1 = 2.0;
+const Double_t gkPipeZ2 = 125.0;
+const Double_t gkPipeR2 = 5.0;
+//const Double_t gkPipeZ3 = 125.0;
+//const Double_t gkPipeR3 = 5.5;
+
+//const Double_t gkPipeZ2 = 50.0;
+//const Double_t gkPipeR2 = 1.8;
+
+//DE const Double_t gkPipeZ1 = 27.0;
+//DE const Double_t gkPipeR1 = 1.05;
+//DE const Double_t gkPipeZ2 = 160.0;
+//DE const Double_t gkPipeR2 = 3.25;
+// ----------------------------------------------------------------------------
+
+//TString unitName[16] = // names of units for v16e
+// { "Unit00D",
+// "Unit01U", "Unit01D",
+// "Unit02U", "Unit02D",
+// "Unit03U", "Unit03D",
+// "Unit04U", "Unit04D",
+// "Unit05U", "Unit05D",
+// "Unit06U", "Unit06D",
+// "Unit07U", "Unit07D",
+// "Unit08U" };
+
+//TString unitName[32] = // names of units for v16f
+// { "Unit00DR", "Unit00DL",
+// "Unit01UR", "Unit01UL", "Unit01DR", "Unit01DL",
+// "Unit02UR", "Unit02UL", "Unit02DR", "Unit02DL",
+// "Unit03UR", "Unit03UL", "Unit03DR", "Unit03DL",
+// "Unit04UR", "Unit04UL", "Unit04DR", "Unit04DL",
+// "Unit05UR", "Unit05UL", "Unit05DR", "Unit05DL",
+// "Unit06UR", "Unit06UL", "Unit06DR", "Unit06DL",
+// "Unit07UR", "Unit07UL", "Unit07DR", "Unit07DL",
+// "Unit08UR", "Unit08UL" };
+
+TString unitName[32] = // names of units for v16g - while merging D and U parts
+ {"Unit00R", "Unit00L", "Unit01R", "Unit01L", "Unit01R", "Unit01L", "Unit02R", "Unit02L",
+ "Unit02R", "Unit02L", "Unit03R", "Unit03L", "Unit03R", "Unit03L", "Unit04R", "Unit04L",
+ "Unit04R", "Unit04L", "Unit05R", "Unit05L", "Unit05R", "Unit05L", "Unit06R", "Unit06L",
+ "Unit06R", "Unit06L", "Unit07R", "Unit07L", "Unit07R", "Unit07L", "Unit08R", "Unit08L"};
+
+TString unitName18[18] = // names of units for v16g
+ {"Unit00R", "Unit00L", "Unit01R", "Unit01L", "Unit02R", "Unit02L", "Unit03R", "Unit03L", "Unit04R",
+ "Unit04L", "Unit05R", "Unit05L", "Unit06R", "Unit06L", "Unit07R", "Unit07L", "Unit08R", "Unit08L"};
+
+// ------------- Other global variables -----------------------------------
+// ---> STS medium (for every volume except silicon)
+TGeoMedium* gStsMedium = NULL; // will be set later
+// ---> TGeoManager (too lazy to write out 'Manager' all the time
+TGeoManager* gGeoMan = NULL; // will be set later
+// ----------------------------------------------------------------------------
+
+
+// ============================================================================
+// ====== Main function =====
+// ============================================================================
+
+void create_stsgeo_v21g(const char* geoTag = "v21g")
+{
+
+ // ------- Geometry file name (output) ----------------------------------
+ TString geoFileName = "sts_";
+ geoFileName = geoFileName + geoTag + ".geo.root";
+ // --------------------------------------------------------------------------
+
+
+ // ------- Open info file -----------------------------------------------
+ TString infoFileName = geoFileName;
+ infoFileName.ReplaceAll("root", "info");
+ fstream infoFile;
+ infoFile.open(infoFileName.Data(), fstream::out);
+ infoFile << "STS geometry created with create_stsgeo_v21g.C" << endl;
+ infoFile << gVersionHighlight << endl;
+ infoFile << "Global variables: " << endl;
+ infoFile << "Sensor thickness = " << gkSensorThickness << " cm" << endl;
+ infoFile << "Vertical gap in sensor chain = " << gkChainGapY << " cm" << endl;
+ infoFile << "Vertical overlap of sensors = " << gkSectorOverlapY << " cm" << endl;
+ infoFile << "Gap in z between neighbour sensors = " << gkSectorGapZ << " cm" << endl;
+ infoFile << "Horizontal overlap of sensors = " << gkLadderOverlapX << " cm" << endl;
+ infoFile << "Gap in z between neighbour ladders = " << gkLadderGapZ << " cm" << endl;
+ if (gkConstructCables) infoFile << "Cable thickness = " << gkCableThickness << " cm" << endl;
+ else
+ infoFile << "No cables" << endl;
+ infoFile << endl;
+ infoFile << "Beam pipe: R1 = " << gkPipeR1 << " cm at z = " << gkPipeZ1 << " cm" << endl;
+ infoFile << "Beam pipe: R2 = " << gkPipeR2 << " cm at z = " << gkPipeZ2 << " cm" << endl;
+ //infoFile << "Beam pipe: R3 = " << gkPipeR3 << " cm at z = " << gkPipeZ3
+ // << " cm" << endl;
+ // --------------------------------------------------------------------------
+
+
+ // ------- Load media from media file -----------------------------------
+ FairGeoLoader* geoLoad = new FairGeoLoader("TGeo", "FairGeoLoader");
+ FairGeoInterface* geoFace = geoLoad->getGeoInterface();
+ TString geoPath = gSystem->Getenv("VMCWORKDIR");
+ TString medFile = geoPath + "/geometry/media.geo";
+ geoFace->setMediaFile(medFile);
+ geoFace->readMedia();
+ gGeoMan = gGeoManager;
+ gSystem->Load("libGeom");
+ // --------------------------------------------------------------------------
+
+
+ // ----------------- Get and create the required media -----------------
+ FairGeoMedia* geoMedia = geoFace->getMedia();
+ FairGeoBuilder* geoBuild = geoLoad->getGeoBuilder();
+
+ // ---> air
+ FairGeoMedium* mAir = geoMedia->getMedium("air");
+ if (!mAir) Fatal("Main", "FairMedium air not found");
+ geoBuild->createMedium(mAir);
+ TGeoMedium* air = gGeoMan->GetMedium("air");
+ if (!air) Fatal("Main", "Medium air not found");
+
+ // ---> air
+ FairGeoMedium* mvacuum = geoMedia->getMedium("vacuum");
+ if (!mvacuum) Fatal("Main", "FairMedium vacuum not found");
+ geoBuild->createMedium(mvacuum);
+ TGeoMedium* Vacuum = gGeoMan->GetMedium("vacuum");
+ if (!Vacuum) Fatal("Main", "Medium vacuum not found");
+
+ // ---> silicon
+ FairGeoMedium* mSilicon = geoMedia->getMedium("silicon");
+ if (!mSilicon) Fatal("Main", "FairMedium silicon not found");
+ geoBuild->createMedium(mSilicon);
+ TGeoMedium* silicon = gGeoMan->GetMedium("silicon");
+ if (!silicon) Fatal("Main", "Medium silicon not found");
+
+ // ---> carbon
+ FairGeoMedium* mCarbon = geoMedia->getMedium("carbon");
+ if (!mCarbon) Fatal("Main", "FairMedium carbon not found");
+ geoBuild->createMedium(mCarbon);
+ TGeoMedium* carbon = gGeoMan->GetMedium("carbon");
+ if (!carbon) Fatal("Main", "Medium carbon not found");
+
+ // ---> STSBoxCarbonFoam
+ FairGeoMedium* mSTSBoxCarbonFoam = geoMedia->getMedium("STSBoxCarbonFoam");
+ if (!mSTSBoxCarbonFoam) Fatal("Main", "FairMedium STSBoxCarbonFoam not found");
+ geoBuild->createMedium(mSTSBoxCarbonFoam);
+ TGeoMedium* STSBoxCarbonFoam = gGeoMan->GetMedium("STSBoxCarbonFoam");
+ if (!STSBoxCarbonFoam) Fatal("Main", "Medium STSBoxCarbonFoam not found");
+
+ // ---> STSBoxCarbonFibre
+ FairGeoMedium* mSTSBoxCarbonFibre = geoMedia->getMedium("STSBoxCarbonFibre");
+ if (!mSTSBoxCarbonFibre) Fatal("Main", "FairMedium STSBoxCarbonFibre not found");
+ geoBuild->createMedium(mSTSBoxCarbonFibre);
+ TGeoMedium* STSBoxCarbonFibre = gGeoMan->GetMedium("STSBoxCarbonFibre");
+ if (!STSBoxCarbonFibre) Fatal("Main", "Medium STSBoxCarbonFibre not found");
+
+ // ---> STScable
+ FairGeoMedium* mSTScable = geoMedia->getMedium("STScable");
+ if (!mSTScable) Fatal("Main", "FairMedium STScable not found");
+ geoBuild->createMedium(mSTScable);
+ TGeoMedium* STScable = gGeoMan->GetMedium("STScable");
+ if (!STScable) Fatal("Main", "Medium STScable not found");
+
+ // ---> Aluminium
+ FairGeoMedium* mAluminium = geoMedia->getMedium("aluminium");
+ if (!mAluminium) Fatal("Main", "FairMedium aluminium not found");
+ geoBuild->createMedium(mAluminium);
+ TGeoMedium* aluminium = gGeoMan->GetMedium("aluminium");
+ if (!aluminium) Fatal("Main", "Medium aluminium not found");
+
+ // ---
+ //gStsMedium = air;
+ gStsMedium = gGeoMan->GetMedium("air");
+ if (!gStsMedium) Fatal("Main", "medium sts_air not found");
+ // --------------------------------------------------------------------------
+ // --------------------------------------------------------------------------
+
+ // -------------- Create geometry and top volume -------------------------
+ gGeoMan = (TGeoManager*) gROOT->FindObject("FAIRGeom");
+ //gGeoMan->SetName("STSgeom");
+ TGeoVolume* top = new TGeoVolumeAssembly("top");
+ // TGeoBBox* topbox= new TGeoBBox("", 120., 120., 120.);
+ // TGeoVolume* top = new TGeoVolume("top", topbox, gGeoMan->GetMedium("air"));
+ gGeoMan->SetTopVolume(top);
+
+ // --------------------------------------------------------------------------
+
+
+ // -------------- Create media ------------------------------------------
+ /*
+ cout << endl;
+ cout << "===> Creating media....";
+ cout << CreateMedia();
+ cout << " media created" << endl;
+ TList* media = gGeoMan->GetListOfMedia();
+ for (Int_t iMedium = 0; iMedium < media->GetSize(); iMedium++ ) {
+ cout << "Medium " << iMedium << ": "
+ << ((TGeoMedium*) media->At(iMedium))->GetName() << endl;
+ }
+ gStsMedium = gGeoMan->GetMedium("air");
+ if ( ! gStsMedium ) Fatal("Main", "medium sts_air not found");
+ */
+ // --------------------------------------------------------------------------
+
+
+ // --------------- Create sensors ---------------------------------------
+ cout << endl << endl;
+ cout << "===> Creating sensors...." << endl << endl;
+ infoFile << endl << "Sensors: " << endl;
+ Int_t nSensors = CreateSensors();
+ for (Int_t iSensor = 1; iSensor <= nSensors; iSensor++) {
+ TString name = Form("Sensor%02d", iSensor);
+ TGeoVolume* sensor = gGeoMan->GetVolume(name);
+
+ // add color to sensors
+ if (iSensor == 1) sensor->SetLineColor(kRed);
+ if (iSensor == 2) sensor->SetLineColor(kGreen + 3);
+ if (iSensor == 3) sensor->SetLineColor(kBlue + 3);
+ if (iSensor == 4) sensor->SetLineColor(kAzure - 7);
+ if (iSensor == 5) sensor->SetLineColor(kYellow);
+ if (iSensor == 6) sensor->SetLineColor(kYellow);
+ if (iSensor == 7) sensor->SetLineColor(kYellow);
+
+ CheckVolume(sensor);
+ CheckVolume(sensor, infoFile);
+ }
+ // --------------------------------------------------------------------------
+
+
+ // ---------------- Create sectors --------------------------------------
+ cout << endl << endl;
+ cout << "===> Creating sectors...." << endl;
+ // infoFile << endl << "Sectors: " << endl;
+ Int_t nSectors = CreateSectors();
+ for (Int_t iSector = 1; iSector <= nSectors; iSector++) {
+ // cout << endl;
+ TString name = Form("Sector%02d", iSector);
+ TGeoVolume* sector = gGeoMan->GetVolume(name);
+ CheckVolume(sector);
+ // CheckVolume(sector, infoFile);
+ }
+ // --------------------------------------------------------------------------
+
+
+ // ---------------- Create ladders --------------------------------------
+ cout << endl << endl;
+ cout << "===> Creating ladders...." << endl;
+ infoFile << endl << "Ladders:" << endl;
+
+ TString name = "";
+ TGeoVolume* ladder;
+
+
+ Int_t nLadders = CreateLadders();
+
+ for (Int_t iLadder = 1; iLadder <= nLadders; iLadder++) {
+ cout << endl;
+
+ // name = Form("LadderType%02d", iLadder); // v19a
+ // name = Form("LadderType0%02d", iLadder); // v19b
+ name = Form("LadderType00%02d", iLadder); // v19k
+ ladder = gGeoMan->GetVolume(name);
+ CheckVolume(ladder);
+ CheckVolume(ladder, infoFile, kFALSE);
+
+ cout << "DF1: ladder name: " << name << endl << endl;
+
+ // name = Form("LadderType%02d", iLadder); // v19a
+ // name = Form("LadderType1%02d", iLadder); // v19b
+ name = Form("LadderType01%02d", iLadder); // v19k
+ ladder = gGeoMan->GetVolume(name);
+ CheckVolume(ladder);
+ CheckVolume(ladder, infoFile, kFALSE);
+
+ cout << "DF2: ladder name: " << name << endl << endl;
+
+ // name = Form("LadderType%02d", iLadder); // v19a
+ name = Form("LadderType10%02d", iLadder); // v19k
+ ladder = gGeoMan->GetVolume(name);
+ CheckVolume(ladder);
+ CheckVolume(ladder, infoFile, kFALSE);
+
+ cout << "DF3: ladder name: " << name << endl << endl;
+
+ // name = Form("LadderType%02d", iLadder); // v19a
+ name = Form("LadderType11%02d", iLadder); // v19k
+ ladder = gGeoMan->GetVolume(name);
+ CheckVolume(ladder);
+ CheckVolume(ladder, infoFile, kFALSE);
+
+ cout << "DF4: ladder name: " << name << endl << endl;
+ }
+ // --------------------------------------------------------------------------
+
+
+ // ---------------- Create cones ----------------------------------------
+ Double_t coneDz = 1.64;
+ TGeoVolume* coneSmallVolum = ConstructSmallCone(coneDz);
+ if (!coneSmallVolum) Fatal("ConstructSmallCone", "Volume Cone not found");
+ TGeoVolume* coneBigVolum = ConstructBigCone(coneDz);
+ if (!coneBigVolum) Fatal("ConstructBigCone", "Volume Cone not found");
+ // --------------------------------------------------------------------------
+
+ // ---------------- Create stations -------------------------------------
+ cout << endl << endl;
+ cout << "===> Creating stations...." << endl;
+ infoFile << endl << "Stations: " << endl;
+ Int_t angle = 0;
+ nLadders = 0;
+ Int_t ladderTypes[16]; // there are max 16 ladders in one layer
+ TGeoTranslation* statTrans = NULL;
+
+ TGeoVolume* myunit[32]; // units
+
+ // Int_t statPos[16] = { 30, 30, 40, 40, 50, 50, 60, 60,
+ // 70, 70, 80, 80, 90, 90, 100, 100 }; // z positions of units
+ // Int_t statPos18[18] = { 30, 30, // expanded for placement of Unit00
+ // 30, 30, 40, 40, 50, 50, 60, 60,
+ // 70, 70, 80, 80, 90, 90, 100, 100 }; // z positions of units
+ // v19h
+
+ Double_t statPos[16] = {-73.5, -73.5, -63.0, -63.0, -52.5, -52.5, -42.0, -42.0,
+ -31.5, -31.5, -21.0, -21.0, -10.5, -10.5, 0, 0};
+
+ Double_t statPos18[18] = {-73.5, -73.5, -73.5, -73.5, -63.0, -63.0, -52.5, -52.5, -42.0,
+ -42.0, -31.5, -31.5, -21.0, -21.0, -10.5, -10.5, 0, 0}; // z positions of unit
+
+ // // v19d
+ // Double_t statPos[16] = { 30.0, 30.0, 40.5, 40.5, 51.0, 51.0, 61.5, 61.5,
+ // 72.0, 72.0, 82.5, 82.5, 93.0, 93.0, 103.5, 103.5 }; // z positions of units
+ // Double_t statPos18[18] = { 30.0, 30.0, // expanded for placement of Unit00
+ // 30.0, 30.0, 40.5, 40.5, 51.0, 51.0, 61.5, 61.5,
+ // 72.0, 72.0, 82.5, 82.5, 93.0, 93.0, 103.5, 103.5 }; // z positions of units
+
+ ////Double_t rHole[8] = { 2.0, 2.0, 2.0, 2.9 , 3.7 , 3.7 , 4.2 , 4.2 }; // size of cutouts in stations
+ // Double_t rHole[8] = { 2.0, 2.0, 2.0, 2.43, 3.04, 3.35, 3.96, 4.2 }; // size of cutouts in stations, derived from gapXYZ[x][1]/2
+
+ Int_t cone_size[8] = {0, 0, 0, 1, 1, 1, 1, 1}; // size of cones: 0 = small, 1 = large
+
+ Double_t cone_offset[2] = {0.305, 0.285};
+
+ //==============================================================================================
+
+ // explanation: type xyzz
+ // where x = carbon ladder orientation
+ // where y = FEB box orientation
+ // where zz = sensor arrangement on ladder
+ // with FEB orientation - v19b
+ Int_t allUnitTypes[16][16] = {
+ {-1, -1, -1, -1, 10, 0, 9, 0, 101, 0, 109, 0, -1, -1, -1, -1}, // unit00D Station01 00
+ {-1, -1, -1, -1, 0, 1109, 0, 1101, 0, 1009, 0, 1010, -1, -1, -1, -1}, // unit01U Station01 01
+
+ {-1, -1, 0, 10, 0, 9, 0, 2, 0, 109, 0, 110, 0, 111, -1, -1}, // unit01D Station02 02
+ {-1, -1, 1130, 0, 1129, 0, 1128, 0, 1002, 0, 1028, 0, 1029, 0, -1, -1}, // unit02U Station02 03
+
+ {-1, -1, 14, 0, 12, 0, 12, 0, 103, 0, 112, 0, 113, 0, -1, -1}, // unit02D Station03 04
+ {-1, -1, 0, 1113, 0, 1112, 0, 1103, 0, 1012, 0, 1012, 0, 1014, -1, -1}, // unit03U Station03 05
+
+ {-1, 15, 0, 13, 0, 12, 0, 4, 0, 112, 0, 112, 0, 114, 0, -1}, // unit03D Station04 06
+ {-1, 0, 1133, 0, 1131, 0, 1131, 0, 1004, 0, 1031, 0, 1032, 0, 1034, -1}, // unit04U Station04 07
+
+ {-1, 0, 18, 0, 17, 0, 16, 0, 105, 0, 116, 0, 117, 0, 119, -1}, // unit04D Station05 08
+ {-1, 1119, 0, 1117, 0, 1116, 0, 1105, 0, 1016, 0, 1017, 0, 1018, 0, -1}, // unit05U Station05 09
+
+ {-1, 19, 0, 17, 0, 16, 0, 6, 0, 116, 0, 117, 0, 118, 0, -1}, // unit05D Station06 10
+ {-1, 0, 1137, 0, 1136, 0, 1135, 0, 1006, 0, 1035, 0, 1036, 0, 1038, -1}, // unit06U Station06 11
+
+ {21, 0, 25, 0, 20, 0, 20, 0, 107, 0, 120, 0, 120, 0, 127, 0}, // unit06D Station07 12
+ {0, 1127, 0, 1120, 0, 1120, 0, 1107, 0, 1020, 0, 1020, 0, 1025, 0, 1021}, // unit07U Station07 13
+
+ {0, 24, 0, 22, 0, 22, 0, 8, 0, 122, 0, 122, 0, 123, 0, 126}, // unit07D Station08 14
+ {1126, 0, 1123, 0, 1122, 0, 1122, 0, 1008, 0, 1022, 0, 1022, 0, 1024, 0}}; // unit08U Station08 15
+
+
+ //==============================================================================================
+
+ // without FEB orientation - v19a
+ // v19a Int_t allUnitTypes[16][16]=
+ // v19a { { -1, -1, -1, -1, 10, 0, 9, 0, 1, 0, 9, 0, -1, -1, -1, -1 }, // unit00D Station01 00
+ // v19a { -1, -1, -1, -1, 0, 109, 0, 101, 0, 109, 0, 110, -1, -1, -1, -1 }, // unit01U Station01 01
+ // v19a { -1, -1, 0, 10, 0, 9, 0, 2, 0, 9, 0, 10, 0, 11, -1, -1 }, // unit01D Station02 02
+ // v19a { -1, -1, 111, 0, 110, 0, 109, 0, 102, 0, 109, 0, 110, 0, -1, -1 }, // unit02U Station02 03
+ // v19a { -1, -1, 14, 0, 12, 0, 12, 0, 3, 0, 12, 0, 13, 0, -1, -1 }, // unit02D Station03 04
+ // v19a { -1, -1, 0, 113, 0, 112, 0, 103, 0, 112, 0, 112, 0, 114, -1, -1 }, // unit03U Station03 05
+ // v19a { -1, 15, 0, 13, 0, 12, 0, 4, 0, 12, 0, 12, 0, 14, 0, -1 }, // unit03D Station04 06
+ // v19a { -1, 0, 114, 0, 112, 0, 112, 0, 104, 0, 112, 0, 113, 0, 115, -1 }, // unit04U Station04 07
+ // v19a { -1, 0, 18, 0, 17, 0, 16, 0, 5, 0, 16, 0, 17, 0, 19, -1 }, // unit04D Station05 08
+ // v19a { -1, 119, 0, 117, 0, 116, 0, 105, 0, 116, 0, 117, 0, 118, 0, -1 }, // unit05U Station05 09
+ // v19a { -1, 19, 0, 17, 0, 16, 0, 6, 0, 16, 0, 17, 0, 18, 0, -1 }, // unit05D Station06 10
+ // v19a { -1, 0, 118, 0, 117, 0, 116, 0, 106, 0, 116, 0, 117, 0, 119, -1 }, // unit06U Station06 11
+ // v19a { 21, 0, 25, 0, 20, 0, 20, 0, 7, 0, 20, 0, 20, 0, 27, 0 }, // unit06D Station07 12
+ // v19a { 0, 127, 0, 120, 0, 120, 0, 107, 0, 120, 0, 120, 0, 125, 0, 121 }, // unit07U Station07 13
+ // v19a { 0, 24, 0, 22, 0, 22, 0, 8, 0, 22, 0, 22, 0, 23, 0, 26 }, // unit07D Station08 14
+ // v19a { 126, 0, 123, 0, 122, 0, 122, 0, 108, 0, 122, 0, 122, 0, 124, 0 } }; // unit08U Station08 15
+
+ // // generate unit
+ // for (Int_t iUnit = 0; iUnit < 16; iUnit++)
+ // for (Int_t iLadder = 0; iLadder < 16; iLadder++)
+ // {
+ // allUnitTypes[iUnit][iLadder] = 0;
+ // if ((iUnit % 2 == 0) && (allLadderTypes[iUnit/2][iLadder] < 100)) // if carbon structure is oriented upstream
+ // allUnitTypes[iUnit][iLadder] = allLadderTypes[iUnit/2][iLadder];
+ // if ((iUnit % 2 == 1) && (allLadderTypes[iUnit/2][iLadder] >= 100)) // if carbon structure is oriented downstream
+ // allUnitTypes[iUnit][iLadder] = allLadderTypes[iUnit/2][iLadder];
+ // }
+
+
+ // dump unit
+ for (Int_t iUnit = 0; iUnit < 16; iUnit++) {
+ cout << "DE unitTypes[" << iUnit << "] = { ";
+ for (Int_t iLadder = 0; iLadder < 16; iLadder++) {
+ cout << allUnitTypes[iUnit][iLadder];
+ if (iLadder < 15) cout << ", ";
+ else
+ cout << " };";
+ }
+ cout << endl;
+ }
+
+
+ // --- Units 01 - 16
+ for (Int_t iUnit = 0; iUnit < 16; iUnit++) {
+ cout << endl;
+
+ nLadders = 0;
+ //int ladderTypes[nLadders];
+ for (Int_t iLadder = 0; iLadder < 16; iLadder++)
+ if (allUnitTypes[iUnit][iLadder] >= 0) {
+ ladderTypes[nLadders] = allUnitTypes[iUnit][iLadder];
+ cout << "DE ladderTypes[" << nLadders << "] = " << allUnitTypes[iUnit][iLadder] << ";" << endl;
+ nLadders++;
+ }
+ myunit[iUnit * 2 + 0] = ConstructUnit(0, iUnit * 2 + 0, nLadders, ladderTypes, iUnit / 2 + 1);
+ myunit[iUnit * 2 + 1] = ConstructUnit(1, iUnit * 2 + 1, nLadders, ladderTypes, iUnit / 2 + 1);
+
+ // if (gkConstructCones) {
+ // if (iUnit%2 == 0)
+ // angle = 90;
+ // else
+ // angle = -90;
+ //
+ // // upstream
+ // TGeoRotation* coneRot11 = new TGeoRotation;
+ // coneRot11->RotateZ(angle);
+ // coneRot11->RotateY(180);
+ // TGeoCombiTrans* conePosRot11 = new TGeoCombiTrans(name+"conePosRot2", 0., 0., -coneDz-cone_offset[cone_size[iUnit]]-gkLadderGapZ/2., coneRot11);
+ // if (cone_size[iUnit] == 0)
+ // myunit[iUnit]->AddNode(coneSmallVolum, 1, conePosRot11);
+ // else
+ // myunit[iUnit]->AddNode(coneBigVolum, 1, conePosRot11);
+ //
+ // // downstream
+ // TGeoRotation* coneRot12 = new TGeoRotation;
+ // coneRot12->RotateZ(angle);
+ // TGeoCombiTrans* conePosRot12 = new TGeoCombiTrans(name+"conePosRot1", 0., 0., coneDz+cone_offset[cone_size[iUnit]]+gkLadderGapZ/2., coneRot12);
+ // if (cone_size[iUnit] == 0)
+ // myunit[iUnit]->AddNode(coneSmallVolum, 2, conePosRot12);
+ // else
+ // myunit[iUnit]->AddNode(coneBigVolum, 2, conePosRot12);
+ //
+ // myunit[iUnit]->GetShape()->ComputeBBox();
+ // }
+
+ // CheckVolume(myunit[iUnit]);
+ // CheckVolume(myunit[iUnit], infoFile);
+ if ((iUnit % 2 == 0) || (iUnit == 15)) {
+ CheckVolume(myunit[iUnit * 2 + 0]);
+ CheckVolume(myunit[iUnit * 2 + 0], infoFile);
+ CheckVolume(myunit[iUnit * 2 + 1]);
+ CheckVolume(myunit[iUnit * 2 + 1], infoFile);
+ }
+ infoFile << "Position z = " << statPos[iUnit] << endl;
+ }
+ // --------------------------------------------------------------------------
+
+
+ // --------------- Create STS volume ------------------------------------
+ cout << endl << endl;
+ cout << "===> Creating STS...." << endl;
+
+ // --- Create STS volume
+ TString stsName = "sts_";
+ stsName += geoTag;
+
+ // TGeoShape* stsShape = new TGeoCompositeShape("stsShape",
+ // "stsBox-stsCone1:trans1-stsCone2:trans2");
+ // TGeoVolume* sts = new TGeoVolume(stsName.Data(), stsShape, gStsMedium);
+
+ Double_t stsBorder = 2 * 5.;
+
+ TGeoVolume* sts = new TGeoVolumeAssembly(stsName.Data());
+
+ // --- Place stations in the STS
+ Double_t stsPosZ = 0.5 * (statPos[15] + statPos[0]); // todo units: update statPos[7]
+ //cout << "stsPosZ " << stsPosZ << " " << statPos[15] << " " << statPos[0] << "*****" << endl;
+
+ // for (Int_t iUnit = 0; iUnit < 16; iUnit++) {
+ for (Int_t iUnit = 0; iUnit < 18; iUnit++) {
+ // for (Int_t iUnit = 0; iUnit < 32; iUnit++) {
+ TGeoVolume* station = gGeoMan->GetVolume(unitName18[iUnit]);
+ // Double_t posZ = statPos[iUnit] - stsPosZ;
+ //Double_t posZ = statPos18[iUnit] - stsPosZ;
+ Double_t posZ =
+ statPos18[iUnit]; // change has been implemented because local sts origin is taken from last station of sts.
+ cout << "Here is position_Z of stations : " << posZ << endl;
+
+ // Double_t posZ = statPos[iUnit/2] - stsPosZ;
+ TGeoTranslation* trans = new TGeoTranslation(0., 0., posZ);
+ sts->AddNode(station, iUnit + 1, trans);
+ sts->GetShape()->ComputeBBox();
+ }
+
+ // --- Import passive elements from GDML file
+ if (gkImportPassive) { ImportPassive(sts, geoTag, infoFile); }
+
+ // --------------------------------------------------------------------------
+ // --------------------------------------------------------------------------
+
+ cout << endl;
+ CheckVolume(sts);
+
+
+ // --------------- Finish -----------------------------------------------
+ cout << "Here is sts_position_Z : " << stsPosZ << endl;
+ // cout << "Here is position_Z of units : " << posZ << endl;
+
+ // TGeoTranslation* stsTrans = new TGeoTranslation(0., 0., stsPosZ);
+
+ // to make translation to the center of magnet (Reason: Origin point{x=0, y=0, z=0} is taken from the last sts station.)
+ TGeoTranslation* stsTrans = new TGeoTranslation(0., 0., 59.5);
+ top->AddNode(sts, 1, stsTrans);
+ top->GetShape()->ComputeBBox();
+ cout << endl << endl;
+
+ CheckVolume(top);
+ cout << endl << endl;
+ gGeoMan->CloseGeometry();
+ gGeoMan->CheckOverlaps(0.0001);
+ gGeoMan->PrintOverlaps();
+ gGeoMan->CheckOverlaps(0.0001, "s");
+ gGeoMan->PrintOverlaps();
+ gGeoMan->Test();
+
+ //TFile* geoFile = new TFile(geoFileName, "RECREATE");
+ //top->Write();
+ sts->Export(geoFileName); // an another way of writing the stsvolume placement
+
+ TFile* geoFile = new TFile(geoFileName, "UPDATE");
+ TGeoTranslation* sts_volume_placement = new TGeoTranslation("sts_trans", 0., 0., 59.5);
+ sts_volume_placement->Write();
+ //*/
+ cout << endl;
+ cout << "Geometry " << top->GetName() << " written to " << geoFileName << endl;
+ geoFile->Close();
+
+ TString geoFileName_ = "sts_";
+ geoFileName_ = geoFileName_ + geoTag + "_geo.root";
+
+ geoFile = new TFile(geoFileName_, "RECREATE");
+ gGeoMan->Write(); // use this is you want GeoManager format in the output
+ geoFile->Close();
+
+ TString geoFileName__ = "sts_";
+ geoFileName_ = geoFileName__ + geoTag + "-geo.root";
+ sts->Export(geoFileName_);
+
+ geoFile = new TFile(geoFileName_, "UPDATE");
+ stsTrans->Write();
+ geoFile->Close();
+
+
+ // gGeoManager->FindVolumeFast("LadderType10_CarbonElement")->Draw("ogl");
+ top->Draw("ogl");
+ //sts->Draw("ogl");
+ gGeoManager->SetVisLevel(9);
+
+ infoFile.close();
+}
+// ============================================================================
+// ====== End of main function =====
+// ============================================================================
+
+
+// ****************************************************************************
+// ***** Definition of media, sensors, sectors and ladders *****
+// ***** *****
+// ***** Decoupled from main function for better readability *****
+// ****************************************************************************
+
+/** ===========================================================================
+ ** Create media
+ **
+ ** Currently created: air, active silicon, passive silion
+ **
+ ** Not used for the time being
+ **/
+Int_t CreateMedia()
+{
+
+ Int_t nMedia = 0;
+ Double_t density = 0.;
+
+ // --- Material air
+ density = 1.205e-3; // [g/cm^3]
+ TGeoMixture* matAir = new TGeoMixture("sts_air", 3, density);
+ matAir->AddElement(14.0067, 7, 0.755); // Nitrogen
+ matAir->AddElement(15.999, 8, 0.231); // Oxygen
+ matAir->AddElement(39.948, 18, 0.014); // Argon
+
+ // --- Material silicon
+ density = 2.33; // [g/cm^3]
+ TGeoElement* elSi = gGeoMan->GetElementTable()->GetElement(14);
+ TGeoMaterial* matSi = new TGeoMaterial("matSi", elSi, density);
+
+
+ // --- Air (passive)
+ TGeoMedium* medAir = new TGeoMedium("air", nMedia++, matAir);
+ medAir->SetParam(0, 0.); // is passive
+ medAir->SetParam(1, 1.); // is in magnetic field
+ medAir->SetParam(2, 20.); // max. field [kG]
+ medAir->SetParam(6, 0.001); // boundary crossing precision [cm]
+
+
+ // --- Active silicon for sensors
+ TGeoMedium* medSiAct = new TGeoMedium("silicon", nMedia++, matSi);
+ medSiAct->SetParam(0, 1.); // is active
+ medSiAct->SetParam(1, 1.); // is in magnetic field
+ medSiAct->SetParam(2, 20.); // max. field [kG]
+ medSiAct->SetParam(6, 0.001); // boundary crossing precisison [cm]
+
+ // --- Passive silicon for cables
+ TGeoMedium* medSiPas = new TGeoMedium("carbon", nMedia++, matSi);
+ medSiPas->SetParam(0, 0.); // is passive
+ medSiPas->SetParam(1, 1.); // is in magnetic field
+ medSiPas->SetParam(2, 20.); // max. field [kG]
+ medSiPas->SetParam(6, 0.001); // boundary crossing precisison [cm]
+
+ return nMedia;
+}
+/** ======================================================================= **/
+
+
+/** ===========================================================================
+ ** Create sensors
+ **
+ ** Sensors are created as volumes with box shape and active silicon as medium.
+ ** Four kinds of sensors: 3.2x2.2, 6.2x2.2, 6.2x4.2, 6.2x6.2
+ **/
+Int_t CreateSensors()
+{
+
+ Int_t nSensors = 0;
+
+ Double_t xSize = 0.;
+ Double_t ySize = 0.;
+ Double_t zSize = gkSensorThickness;
+ TGeoMedium* silicon = gGeoMan->GetMedium("silicon");
+
+
+ // --- Sensor type 01: Small sensor (6.2 cm x 2.2 cm)
+ xSize = gkSensorSizeX;
+ ySize = 2.2;
+ TGeoBBox* shape_sensor01 = new TGeoBBox("sensor01", xSize / 2., ySize / 2., zSize / 2.);
+ new TGeoVolume("Sensor01", shape_sensor01, silicon);
+ nSensors++;
+
+
+ // --- Sensor type 02: Medium sensor (6.2 cm x 4.2 cm)
+ xSize = gkSensorSizeX;
+ ySize = 4.2;
+ TGeoBBox* shape_sensor02 = new TGeoBBox("sensor02", xSize / 2., ySize / 2., zSize / 2.);
+ new TGeoVolume("Sensor02", shape_sensor02, silicon);
+ nSensors++;
+
+
+ // --- Sensor type 03: Big sensor (6.2 cm x 6.2 cm)
+ xSize = gkSensorSizeX;
+ ySize = 6.2;
+ TGeoBBox* shape_sensor03 = new TGeoBBox("sensor03", xSize / 2., ySize / 2., zSize / 2.);
+ new TGeoVolume("Sensor03", shape_sensor03, silicon);
+ nSensors++;
+
+
+ // --- Sensor type 04: Big sensor (6.2 cm x 12.4 cm)
+ xSize = gkSensorSizeX;
+ ySize = 12.4;
+ TGeoBBox* shape_sensor04 = new TGeoBBox("sensor04", xSize / 2., ySize / 2., zSize / 2.);
+ new TGeoVolume("Sensor04", shape_sensor04, silicon);
+ nSensors++;
+
+
+ // below are extra small sensors, those are not available in the CAD model
+
+ // --- Sensor Type 05: Half small sensor (4 cm x 2.5 cm)
+ xSize = 4.0;
+ ySize = 2.5;
+ TGeoBBox* shape_sensor05 = new TGeoBBox("sensor05", xSize / 2., ySize / 2., zSize / 2.);
+ new TGeoVolume("Sensor05", shape_sensor05, silicon);
+ nSensors++;
+
+
+ // --- Sensor type 06: Additional "in hole" sensor (3.1 cm x 4.2 cm)
+ xSize = 3.1;
+ ySize = 4.2;
+ TGeoBBox* shape_sensor06 = new TGeoBBox("sensor06", xSize / 2., ySize / 2., zSize / 2.);
+ new TGeoVolume("Sensor06", shape_sensor06, silicon);
+ nSensors++;
+
+
+ // --- Sensor type 07: Mini Medium sensor (1.5 cm x 4.2 cm)
+ xSize = 1.5;
+ ySize = 4.2;
+ TGeoBBox* shape_sensor07 = new TGeoBBox("sensor07", xSize / 2., ySize / 2., zSize / 2.);
+ new TGeoVolume("Sensor07", shape_sensor07, silicon);
+ nSensors++;
+
+
+ return nSensors;
+}
+/** ======================================================================= **/
+
+
+/** ===========================================================================
+ ** Create sectors
+ **
+ ** A sector is either a single sensor or several chained sensors.
+ ** It is implemented as TGeoVolumeAssembly.
+ ** Currently available:
+ ** - single sensors of type 1 - 4
+ ** - two chained sensors of type 4
+ ** - three chained sensors of type 4
+ **/
+Int_t CreateSectors()
+{
+
+ Int_t nSectors = 0;
+
+ TGeoVolume* sensor01 = gGeoMan->GetVolume("Sensor01");
+ TGeoVolume* sensor02 = gGeoMan->GetVolume("Sensor02");
+ TGeoVolume* sensor03 = gGeoMan->GetVolume("Sensor03");
+ TGeoVolume* sensor04 = gGeoMan->GetVolume("Sensor04");
+ TGeoVolume* sensor05 = gGeoMan->GetVolume("Sensor05");
+ TGeoVolume* sensor06 = gGeoMan->GetVolume("Sensor06");
+ TGeoVolume* sensor07 = gGeoMan->GetVolume("Sensor07");
+ // TGeoBBox* box4 = (TGeoBBox*) sensor04->GetShape();
+
+ // --- Sector type 1: single sensor of type 1
+ TGeoVolumeAssembly* sector01 = new TGeoVolumeAssembly("Sector01");
+ sector01->AddNode(sensor01, 1);
+ sector01->GetShape()->ComputeBBox();
+ nSectors++;
+
+ // --- Sector type 2: single sensor of type 2
+ TGeoVolumeAssembly* sector02 = new TGeoVolumeAssembly("Sector02");
+ sector02->AddNode(sensor02, 1);
+ sector02->GetShape()->ComputeBBox();
+ nSectors++;
+
+ // --- Sector type 3: single sensor of type 3
+ TGeoVolumeAssembly* sector03 = new TGeoVolumeAssembly("Sector03");
+ sector03->AddNode(sensor03, 1);
+ sector03->GetShape()->ComputeBBox();
+ nSectors++;
+
+ // --- Sector type 4: single sensor of type 4
+ TGeoVolumeAssembly* sector04 = new TGeoVolumeAssembly("Sector04");
+ sector04->AddNode(sensor04, 1);
+ sector04->GetShape()->ComputeBBox();
+ nSectors++;
+
+ // --- Sector type 5: single sensor of type 5
+ TGeoVolumeAssembly* sector05 = new TGeoVolumeAssembly("Sector05");
+ sector05->AddNode(sensor05, 1);
+ sector05->GetShape()->ComputeBBox();
+ nSectors++;
+
+ // --- Sector type 6: single sensor of type 6
+ TGeoVolumeAssembly* sector06 = new TGeoVolumeAssembly("Sector06");
+ sector06->AddNode(sensor06, 1);
+ sector06->GetShape()->ComputeBBox();
+ nSectors++;
+
+ // --- Sector type 7: single sensor of type 7
+ TGeoVolumeAssembly* sector07 = new TGeoVolumeAssembly("Sector07");
+ sector07->AddNode(sensor07, 1);
+ sector07->GetShape()->ComputeBBox();
+ nSectors++;
+
+ // // --- Sector type 5: two sensors of type 4
+ // TGeoVolumeAssembly* sector05 = new TGeoVolumeAssembly("Sector05");
+ // Double_t shift5 = 0.5 * gkChainGapY + box4->GetDY();
+ // TGeoTranslation* transD5 =
+ // new TGeoTranslation("td", 0., -1. * shift5, 0.);
+ // TGeoTranslation* transU5 =
+ // new TGeoTranslation("tu", 0., shift5, 0.);
+ // sector05->AddNode(sensor04, 1, transD5);
+ // sector05->AddNode(sensor04, 2, transU5);
+ // sector05->GetShape()->ComputeBBox();
+ // nSectors++;
+
+ return nSectors;
+}
+/** ======================================================================= **/
+
+
+/** ===========================================================================
+ ** Create ladders
+ **
+ ** Ladders are the building blocks of the stations. They contain
+ ** several modules placed one after the other along the z axis
+ ** such that the sectors are arranged vertically (with overlap).
+ **
+ ** A ladder is constructed out of two half ladders, the second of which
+ ** is rotated in the x-y plane by 180 degrees and displaced
+ ** in z direction.
+ **/
+Int_t CreateLadders()
+{
+
+ Int_t nLadders = 0;
+
+ // --- Some variables
+ Int_t nSectors = 0;
+ Int_t sectorTypes[10];
+ TGeoBBox* shape = NULL;
+ TString s0name;
+ TString hlname;
+ char align;
+ TGeoVolume* s0vol = NULL;
+ TGeoVolume* halfLadderU = NULL;
+ TGeoVolume* halfLadderD = NULL;
+
+ // --- Ladders 01-23
+ // --- Ladders 01-27
+ Int_t allSectorTypes[38][6] = {
+ {1, 2, 3, 3, 0, -1}, // ladder 01 - 5 - last column defines alignment of small sensors
+ {1, 2, 3, 3, 0, 0}, // ladder 02 - 5 - last column defines alignment of small sensors
+ {2, 2, 3, 4, 0, -1}, // ladder 03 - 6 - last column defines alignment of small sensors
+ {2, 2, 3, 4, 0, 0}, // ladder 04 - 6 - last column defines alignment of small sensors
+ {2, 3, 4, 4, 0, -1}, // ladder 05 - 7 - last column defines alignment of small sensors
+ {2, 3, 4, 4, 0, 0}, // ladder 06 - 7 - last column defines alignment of small sensors
+ {2, 2, 3, 4, 4, 0}, // ladder 07 - last column defines alignment of small sensors
+ {3, 4, 4, 4, 0, 0}, // ladder 08 - last column defines alignment of small sensors
+ {1, 1, 2, 3, 3, 0}, // ladder 09 - last column defines alignment of small sensors
+ {1, 1, 2, 2, 3, 0}, // ladder 10 - last column defines alignment of small sensors
+ {2, 2, 0, 0, 0, 0}, // ladder 11 - last column defines alignment of small sensors
+ {2, 2, 2, 3, 4, 0}, // ladder 12 - last column defines alignment of small sensors
+ {2, 2, 3, 4, 0, 0}, // ladder 13 - last column defines alignment of small sensors
+ {2, 3, 4, 0, 0, 0}, // ladder 14 - last column defines alignment of small sensors
+ {3, 3, 0, 0, 0, 0}, // ladder 15 - last column defines alignment of small sensors
+ {2, 2, 3, 4, 4, 0}, // ladder 16 - last column defines alignment of small sensors
+ {2, 3, 4, 4, 0, 0}, // ladder 17 - last column defines alignment of small sensors
+ {3, 4, 4, 0, 0, 0}, // ladder 18 - last column defines alignment of small sensors
+ {4, 4, 0, 0, 0, 0}, // ladder 19 - last column defines alignment of small sensors
+ {1, 2, 4, 4, 4, 0}, // ladder 20 - last column defines alignment of small sensors
+ {4, 0, 0, 0, 0, 0}, // ladder 21 - last column defines alignment of small sensors
+ {2, 3, 4, 4, 4, 0}, // ladder 22 - last column defines alignment of small sensors
+ {2, 3, 3, 4, 4, 0}, // ladder 23 - last column defines alignment of small sensors
+ {2, 3, 4, 4, 0, 0}, // ladder 24 - copy of 17 with different total length
+ {3, 4, 4, 0, 0, 0}, // ladder 25 - copy of 18 with different total length
+ {4, 4, 0, 0, 0, 0}, // ladder 26 - copy of 19 with different total length
+ {4, 4, 0, 0, 0, 0}, // ladder 27 - copy of 19 with different total length
+
+ {1, 1, 2, 3, 3, 0}, // ladder 28 - copy of 09 with different total length
+ {1, 1, 2, 2, 3, 0}, // ladder 29 - copy of 10 with different total length
+ {2, 2, 0, 0, 0, 0}, // ladder 30 - copy of 11 with different total length
+
+ {2, 2, 2, 3, 4, 0}, // ladder 31 - copy of 12 with different total length
+ {2, 2, 3, 4, 0, 0}, // ladder 32 - copy of 13 with different total length
+ {2, 3, 4, 0, 0, 0}, // ladder 33 - copy of 14 with different total length
+ {3, 3, 0, 0, 0, 0}, // ladder 34 - copy of 15 with different total length
+
+ {2, 2, 3, 4, 4, 0}, // ladder 35 - copy of 16 with different total length
+ {2, 3, 4, 4, 0, 0}, // ladder 36 - copy of 17 with different total length
+ {3, 4, 4, 0, 0, 0}, // ladder 37 - copy of 18 with different total length
+ {4, 4, 0, 0, 0, 0}, // ladder 38 - copy of 19 with different total length
+
+ }; // 8 full - 19 partial ladders
+
+ // Issue #405
+ // Counting from the most upstream ladder, the gaps between sensors are as follows:
+ // 01 (most upstream): 41.3mm
+ // 02: 41.3mm
+ // 03: 42.0mm
+ // 04: 48.6mm
+ // 05: 60.8mm
+ // 06: 67.0mm
+ // 07: 79.2mm
+ // 08 (most downstream): 88.0mm
+
+ Double_t pitchZ = 0;
+ Double_t gapXYZ[38][3] = {
+ {0., 5.43, 0.}, // ladder 01
+ {0., 6.00, 0.}, // ladder 02
+ {0., 6.70, 0.}, // ladder 03
+ {0., 7.35, 0.}, // ladder 04
+ {0., 8.02, 0.}, // ladder 05
+ {0., 8.70, 0.}, // ladder 06
+ {0., 9.22, 0.}, // ladder 07
+ {0., 10.00, 0.}, // ladder 08
+ {0., -gkSectorOverlapY, 0.}, // ladder 09
+ {0., -gkSectorOverlapY, 0.}, // ladder 10
+ {0., -gkSectorOverlapY, 0.}, // ladder 11
+ {0., -gkSectorOverlapY, 0.}, // ladder 12
+ {0., -gkSectorOverlapY, 0.}, // ladder 13
+ {0., -gkSectorOverlapY, 0.}, // ladder 14
+ {0., -gkSectorOverlapY, 0.}, // ladder 15
+ {0., -gkSectorOverlapY, 0.}, // ladder 16
+ {0., -gkSectorOverlapY, 0.}, // ladder 17
+ {0., -gkSectorOverlapY, 0.}, // ladder 18
+ {0., -gkSectorOverlapY, 0.}, // ladder 19
+ {0., -gkSectorOverlapY, 0.}, // ladder 20
+ {0., -gkSectorOverlapY, 0.}, // ladder 21
+ {0., -gkSectorOverlapY, 0.}, // ladder 22
+ {0., -gkSectorOverlapY, 0.}, // ladder 23
+
+ {0., -gkSectorOverlapY, 0.}, // ladder 24 - copy of 17 with different total length
+ {0., -gkSectorOverlapY, 0.}, // ladder 25 - copy of 18 with different total length
+ {0., -gkSectorOverlapY, 0.}, // ladder 26 - copy of 19 with different total length
+ {0., -gkSectorOverlapY, 0.}, // ladder 27 - copy of 19 with different total length
+
+ {0., -gkSectorOverlapY, 0.}, // ladder 28 - copy of 09 with different total length
+ {0., -gkSectorOverlapY, 0.}, // ladder 29 - copy of 10 with different total length
+ {0., -gkSectorOverlapY, 0.}, // ladder 30 - copy of 11 with different total length
+
+ {0., -gkSectorOverlapY, 0.}, // ladder 31 - copy of 12 with different total length
+ {0., -gkSectorOverlapY, 0.}, // ladder 32 - copy of 13 with different total length
+ {0., -gkSectorOverlapY, 0.}, // ladder 33 - copy of 14 with different total length
+ {0., -gkSectorOverlapY, 0.}, // ladder 34 - copy of 15 with different total length
+
+ {0., -gkSectorOverlapY, 0.}, // ladder 35 - copy of 16 with different total length
+ {0., -gkSectorOverlapY, 0.}, // ladder 36 - copy of 17 with different total length
+ {0., -gkSectorOverlapY, 0.}, // ladder 37 - copy of 18 with different total length
+ {0., -gkSectorOverlapY, 0.}, // ladder 38 - copy of 19 with different total length
+
+ };
+
+ Double_t ladderLength[38] = {
+ 48.0, // ladder 01
+ 48.0, // ladder 02
+ 66.8, // ladder 03
+ 66.8, // ladder 04
+ 81.8, // ladder 05
+ 81.8, // ladder 06
+ 89.2, // ladder 07
+ 96.7, // ladder 08
+ 48.0, // ladder 09
+ 48.0, // ladder 10
+ 48.0, // ladder 11
+ 66.8, // ladder 12
+ 66.8, // ladder 13
+ 66.8, // ladder 14
+ 66.8, // ladder 15
+ 81.8, // ladder 16
+ 81.8, // ladder 17
+ 81.8, // ladder 18
+ 81.8, // ladder 19
+ 89.2, // ladder 20
+ 89.2, // ladder 21
+ 96.7, // ladder 22
+ 96.7, // ladder 23
+
+ 96.7, // ladder 24 - copy of 17 with different total length
+ 89.2, // ladder 25 - copy of 18 with different total length
+ 96.7, // ladder 26 - copy of 19 with different total length
+ 89.2, // ladder 27 - copy of 19 with different total length
+
+ 66.8, // ladder 28 - copy of 09 with different total length
+ 66.8, // ladder 29 - copy of 10 with different total length
+ 66.8, // ladder 30 - copy of 11 with different total length
+
+ 81.8, // ladder 31 - copy of 12 with different total length
+ 81.8, // ladder 32 - copy of 13 with different total length
+ 81.8, // ladder 33 - copy of 14 with different total length
+ 81.8, // ladder 34 - copy of 15 with different total length
+
+ 89.2, // ladder 35 - copy of 16 with different total length
+ 89.2, // ladder 36 - copy of 17 with different total length
+ 89.2, // ladder 37 - copy of 18 with different total length
+ 89.2, // ladder 38 - copy of 19 with different total length
+
+ };
+ // ========================================================================
+
+ // calculate Z shift for ladders with and without gaps in the center
+ s0name = Form("Sector%02d", allSectorTypes[0][0]);
+ s0vol = gGeoMan->GetVolume(s0name);
+ shape = (TGeoBBox*) s0vol->GetShape();
+
+ // ========================================================================
+
+ for (Int_t iLadder = 0; iLadder < 38; iLadder++) {
+ cout << endl;
+ nSectors = 0;
+ for (Int_t i = 0; i < 5; i++)
+ if (allSectorTypes[iLadder][i] != 0) {
+ sectorTypes[nSectors] = allSectorTypes[iLadder][i]; // copy sectors for this ladder
+ cout << "DE iLadder " << iLadder + 1 << " sectorTypes[" << nSectors << "] = " << allSectorTypes[iLadder][i]
+ << ";" << endl;
+ nSectors++; // count how many sectors are in this ladder
+ }
+
+ // always set displacement in z between upper and lower half ladder
+ gapXYZ[iLadder][2] = 2. * shape->GetDZ() + gkSectorGapZ;
+
+ // define additional offset to carbon ladder for half ladders with less than 5 sensors
+ pitchZ = 2. * shape->GetDZ() + gkSectorGapZ;
+
+ if (allSectorTypes[iLadder][5] == 0) align = 'l';
+ else
+ align = 'r';
+ hlname = Form("HalfLadder%02dt", iLadder + 1);
+ // build upper half ladder
+ halfLadderU = ConstructHalfLadder(hlname, nSectors, sectorTypes, align, ladderLength[iLadder] / 2.,
+ gapXYZ[iLadder][1] / 2.); // mirrored
+
+ if (allSectorTypes[iLadder][5] == 0) align = 'r';
+ else
+ align = 'l';
+ hlname = Form("HalfLadder%02db", iLadder + 1);
+ // build lower half ladder
+ halfLadderD = ConstructHalfLadder(hlname, nSectors, sectorTypes, align, ladderLength[iLadder] / 2.,
+ gapXYZ[iLadder][1] / 2.); // mirrored
+
+ // at this point half ladders are constructed
+
+ // build all 4 possible ladders types for this sensor arrangement
+ ConstructLadder(iLadder + 1, halfLadderU, halfLadderD, gapXYZ[iLadder][1], gapXYZ[iLadder][2], pitchZ,
+ nSectors); // v19a
+ ConstructLadder(100 + iLadder + 1, halfLadderU, halfLadderD, gapXYZ[iLadder][1], gapXYZ[iLadder][2], pitchZ,
+ nSectors); // v19b
+
+ ConstructLadder(1000 + iLadder + 1, halfLadderU, halfLadderD, gapXYZ[iLadder][1], gapXYZ[iLadder][2], pitchZ,
+ nSectors); // v19k
+ ConstructLadder(1100 + iLadder + 1, halfLadderU, halfLadderD, gapXYZ[iLadder][1], gapXYZ[iLadder][2], pitchZ,
+ nSectors); // v19k
+
+ nLadders++;
+ }
+
+ return nLadders;
+}
+/** ======================================================================= **/
+
+
+// ****************************************************************************
+// ***** *****
+// ***** Generic functions for the construction of STS elements *****
+// ***** *****
+// ***** module: volume (made of a sector and a cable) *****
+// ***** haf ladder: assembly (made of modules) *****
+// ***** ladder: assembly (made of two half ladders) *****
+// ***** station: volume (made of ladders) *****
+// ***** *****
+// ****************************************************************************
+
+
+/** ===========================================================================
+ ** Construct a module
+ **
+ ** A module is a sector plus the readout cable extending from the
+ ** top of the sector. The cable is made from passive silicon.
+ ** The cable has the same x size as the sector.
+ ** Its thickness is given by the global variable gkCableThickness.
+ ** The cable length is a parameter.
+ ** The sensor(s) of the sector is/are placed directly in the module;
+ ** the sector is just auxiliary for the proper placement.
+ **
+ ** Arguments:
+ ** name volume name
+ ** sector pointer to sector volume
+ ** cableLength length of cable
+ **/
+TGeoVolume* ConstructModule(const char* name, TGeoVolume* sector, Double_t cableLength)
+{
+
+ // --- Check sector volume
+ if (!sector) Fatal("CreateModule", "Sector volume not found!");
+
+ // --- Get size of sector
+ TGeoBBox* box = (TGeoBBox*) sector->GetShape();
+ Double_t sectorX = 2. * box->GetDX();
+ Double_t sectorY = 2. * box->GetDY();
+ Double_t sectorZ = 2. * box->GetDZ();
+
+ // --- Get size of cable
+ Double_t cableX = sectorX;
+ Double_t cableY = cableLength;
+ Double_t cableZ = gkCableThickness;
+
+ // --- Create module volume
+ Double_t moduleX = TMath::Max(sectorX, cableX);
+ Double_t moduleY = sectorY + cableLength;
+
+ Double_t moduleZ = TMath::Max(sectorZ, cableZ);
+
+ TGeoVolume* module = gGeoManager->MakeBox(name, gStsMedium, moduleX / 2., moduleY / 2., moduleZ / 2.);
+
+ // --- Position of sector in module
+ // --- Sector is centred in x and z and aligned to the bottom
+ Double_t sectorXpos = 0.;
+ Double_t sectorYpos = 0.5 * (sectorY - moduleY);
+ Double_t sectorZpos = 0.;
+
+
+ // --- Get sensor(s) from sector
+ Int_t nSensors = sector->GetNdaughters();
+ for (Int_t iSensor = 0; iSensor < nSensors; iSensor++) {
+ TGeoNode* sensor = sector->GetNode(iSensor);
+
+ // --- Calculate position of sensor in module
+ const Double_t* xSensTrans = sensor->GetMatrix()->GetTranslation();
+ Double_t sensorXpos = 0.;
+ Double_t sensorYpos = sectorYpos + xSensTrans[1];
+ Double_t sensorZpos = 0.;
+ TGeoTranslation* sensTrans = new TGeoTranslation("sensTrans", sensorXpos, sensorYpos, sensorZpos);
+
+ // --- Add sensor volume to module
+ TGeoVolume* sensVol = sensor->GetVolume();
+ module->AddNode(sensor->GetVolume(), iSensor + 1, sensTrans);
+ module->GetShape()->ComputeBBox();
+ }
+
+
+ // --- Create cable volume, if necessary, and place it in module
+ // --- Cable is centred in x and z and aligned to the top
+ if (gkConstructCables && cableLength > 0.0001) {
+ TString cableName = TString(name) + "_cable";
+ TGeoMedium* cableMedium = gGeoMan->GetMedium("STScable");
+ if (!cableMedium) Fatal("CreateModule", "Medium STScable not found!");
+ TGeoVolume* cable = gGeoManager->MakeBox(cableName.Data(), cableMedium, cableX / 2., cableY / 2., cableZ / 2.);
+ // add color to cables
+ cable->SetLineColor(kOrange);
+ cable->SetTransparency(60);
+ Double_t cableXpos = 0.;
+ Double_t cableYpos = sectorY + 0.5 * cableY - 0.5 * moduleY;
+ Double_t cableZpos = 0.;
+ TGeoTranslation* cableTrans = new TGeoTranslation("cableTrans", cableXpos, cableYpos, cableZpos);
+ module->AddNode(cable, 1, cableTrans);
+ module->GetShape()->ComputeBBox();
+ }
+
+ return module;
+}
+/** ======================================================================= **/
+
+
+/** ===========================================================================
+ ** Construct a half ladder
+ **
+ ** A half ladder is a virtual volume (TGeoVolumeAssembly) consisting
+ ** of several modules arranged on top of each other. The modules
+ ** have a given overlap in y and a displacement in z to allow for the
+ ** overlap.
+ **
+ ** The typ of sectors / modules to be placed must be specified:
+ ** 1 = sensor01
+ ** 2 = sensor02
+ ** 3 = sensor03
+ ** 4 = sensor04
+ ** 5 = 2 x sensor04 (chained)
+ ** 6 = 3 x sensor04 (chained)
+ ** The cable is added automatically from the top of each sensor to
+ ** the top of the half ladder.
+ ** The alignment can be left (l) or right (r), which matters in the
+ ** case of different x sizes of sensors (e.g. SensorType01).
+ **
+ ** Arguments:
+ ** name volume name
+ ** nSectors number of sectors
+ ** sectorTypes array with sector types
+ ** align horizontal alignment of sectors
+ * ladderLength full length of the ladder towards FEE
+ * offsetY gap in the beam-pipe region
+ **/
+TGeoVolume* ConstructHalfLadder(const TString& name, Int_t nSectors, Int_t* sectorTypes, char align,
+ Double_t ladderLength, Double_t offsetY)
+{
+
+ // --- Create half ladder volume assembly
+ TGeoVolumeAssembly* halfLadder = new TGeoVolumeAssembly(name);
+
+ // --- Determine size of ladder
+ Double_t ladderX = 0.;
+ Double_t ladderY = 0.;
+ Double_t ladderZ = 0.;
+ for (Int_t iSector = 0; iSector < nSectors; iSector++) {
+ TString sectorName = Form("Sector%02d", sectorTypes[iSector]);
+ TGeoVolume* sector = gGeoMan->GetVolume(sectorName);
+ if (!sector) Fatal("ConstructHalfLadder", Form("Volume %s not found", sectorName.Data()));
+ TGeoBBox* box = (TGeoBBox*) sector->GetShape();
+ // --- Ladder x size equals largest sector x size
+ ladderX = TMath::Max(ladderX, 2. * box->GetDX());
+ // --- Ladder y size is sum of sector ysizes
+ ladderY += 2. * box->GetDY();
+ // --- Ladder z size is sum of sector z sizes
+ ladderZ += 2. * box->GetDZ();
+ }
+ // --- Subtract overlaps in y
+ ladderY -= Double_t(nSectors - 1) * gkSectorOverlapY;
+ // --- Add gaps in z direction
+ ladderZ += Double_t(nSectors - 1) * gkSectorGapZ;
+
+ ladderY = TMath::Max(ladderLength - offsetY, ladderY);
+
+ // --- Create and place modules
+ Double_t yPosSect = -0.5 * ladderY;
+ Double_t zPosMod = -0.5 * ladderZ;
+ for (Int_t iSector = 0; iSector < nSectors; iSector++) {
+ TString sectorName = Form("Sector%02d", sectorTypes[iSector]);
+ TGeoVolume* sector = gGeoMan->GetVolume(sectorName);
+ TGeoBBox* box = (TGeoBBox*) sector->GetShape();
+ Double_t sectorX = 2. * box->GetDX();
+ Double_t sectorY = 2. * box->GetDY();
+ Double_t sectorZ = 2. * box->GetDZ();
+ yPosSect += 0.5 * sectorY; // Position of sector in ladder
+ Double_t cableLength = 0.5 * ladderY - yPosSect - 0.5 * sectorY;
+ TString moduleName = name + "_" + Form("Module%02d", sectorTypes[iSector]);
+ TGeoVolume* module = ConstructModule(moduleName.Data(), sector, cableLength);
+
+ TGeoBBox* shapeMod = (TGeoBBox*) module->GetShape();
+ Double_t moduleX = 2. * shapeMod->GetDX();
+ Double_t moduleY = 2. * shapeMod->GetDY();
+ Double_t moduleZ = 2. * shapeMod->GetDZ();
+ Double_t xPosMod = 0.;
+ if (align == 'l') xPosMod = 0.5 * (moduleX - ladderX); // left aligned
+ else if (align == 'r')
+ xPosMod = 0.5 * (ladderX - moduleX); // right aligned
+ else
+ xPosMod = 0.; // centred in x
+ Double_t yPosMod = 0.5 * (ladderY - moduleY); // top aligned
+ zPosMod += 0.5 * moduleZ;
+ TGeoTranslation* trans = new TGeoTranslation("t", xPosMod, yPosMod, zPosMod);
+ halfLadder->AddNode(module, iSector + 1, trans);
+ halfLadder->GetShape()->ComputeBBox();
+ yPosSect += 0.5 * sectorY - gkSectorOverlapY;
+ zPosMod += 0.5 * moduleZ + gkSectorGapZ;
+ }
+
+ CheckVolume(halfLadder);
+ cout << endl;
+
+ return halfLadder;
+}
+/** ======================================================================= **/
+
+
+/** ===========================================================================
+ ** Add a carbon support to a ladder
+ **
+ ** Arguments:
+ ** LadderIndex ladder number
+ ** ladder pointer to ladder
+ ** xu size of halfladder
+ ** ladderY height of ladder along y
+ ** ladderZ thickness of ladder along z
+ **/
+void AddCarbonLadder(Int_t LadderIndex, TGeoVolume* ladder, Double_t xu, Double_t ladderY, Double_t ladderZ)
+{
+
+ // --- Some variables
+ TString name = Form("LadderType%04d", LadderIndex); // v19k
+ Int_t i;
+ Double_t j;
+
+ Int_t YnumOfFrameBoxes = round(ladderY / gkFrameStep);
+
+ // cout << "DEXZ: lad " << LadderIndex << " inum " << YnumOfFrameBoxes << endl;
+
+ Double_t ladderDZ = (xu / 2. + sqrt(2.) * gkFrameThickness / 2. + gkSectorGapZFrame * 2) / 2.;
+ cout << "DEFR: frame Z size " << 2 * ladderDZ << " cm" << endl;
+
+ TGeoBBox* fullFrameShp = new TGeoBBox(name + "_CarbonElement_shp", xu / 2., gkFrameStep / 2., ladderDZ);
+ TGeoVolume* fullFrameBoxVol = new TGeoVolume(name + "_CarbonElement", fullFrameShp, gStsMedium);
+
+ ConstructFrameElement("CarbonElement", fullFrameBoxVol, xu / 2.);
+ TGeoRotation* fullFrameRot = new TGeoRotation;
+ fullFrameRot->RotateY(180);
+
+ Int_t inum = YnumOfFrameBoxes - 3.5; // - 3.5 was done to short frame box to remove overlaps
+ for (i = 1; i <= inum; i++) {
+ j = -(inum - 1) / 2. + (i - 1);
+ // -(10-1)/2. +0 +10-1 -> -4.5 .. +4.5 -> -0.5, +0.5 (= 2)
+ // -(11-1)/2. +0 +11-1 -> -5.0 .. +5.0 -> -1, 0, 1 (= 3)
+ // cout << "DE: i " << i << " j " << j << endl;
+
+ if (LadderIndex % 100 <= 3) // central ladders in stations 1 to 3
+ {
+ if ((j >= -1) && (j <= 1)) // keep the inner 2 (even) or 3 (odd) elements free for the cone
+ continue;
+ }
+ else if (LadderIndex % 100 <= 8) // central ladders in stations 4 to 8
+ {
+ if ((j >= -2) && (j <= 2)) // keep the inner 4 elements free for the cone
+ continue;
+ }
+
+ cout << "DELZ: ladderDZ " << ladderDZ << " cm " << -ladderZ / 2. - ladderDZ << " cm " << endl;
+ ladder->AddNode(fullFrameBoxVol, i,
+ new TGeoCombiTrans(name + "_CarbonElement_posrot", 0., j * gkFrameStep, -(ladderZ / 2. + ladderDZ),
+ fullFrameRot));
+ }
+
+ ladder->GetShape()->ComputeBBox();
+}
+/** ======================================================================= **/
+
+
+/** ===========================================================================
+ ** Construct a ladder out of two half ladders with vertical gap
+ **
+ ** The second half ladder will be rotated by 180 degrees
+ ** in the x-y plane. The two half ladders will be put on top of each
+ ** other with a vertical gap.
+ **
+ ** Arguments:
+ ** name volume name
+ ** halfLadderU pointer to upper half ladder
+ ** halfLadderD pointer to lower half ladder
+ ** gapY vertical gap
+ ** shiftZ relative displacement along the z axis
+ **/
+
+TGeoVolume* ConstructLadder(Int_t LadderIndex, TGeoVolume* halfLadderU, TGeoVolume* halfLadderD, Double_t gapY,
+ Double_t shiftZ, Double_t pitchZ, Int_t nSectors)
+{
+
+ // --- Some variables
+ TGeoBBox* shape = NULL;
+
+ // define additional offset to carbon ladder for half ladders with less than 5 sensors
+ Double_t offsetZ = (5. - nSectors) * pitchZ;
+
+ // --- Dimensions of half ladders
+ shape = (TGeoBBox*) halfLadderU->GetShape(); // up
+ Double_t xu = 2. * shape->GetDX();
+ Double_t yu = 2. * shape->GetDY();
+ Double_t zu = 2. * shape->GetDZ();
+
+ shape = (TGeoBBox*) halfLadderD->GetShape(); // down
+ Double_t xd = 2. * shape->GetDX();
+ Double_t yd = 2. * shape->GetDY();
+ Double_t zd = 2. * shape->GetDZ();
+
+ // --- Create ladder volume assembly
+ TString name = Form("LadderType%04d", LadderIndex); // v19k
+ TGeoVolumeAssembly* ladder = new TGeoVolumeAssembly(name);
+ Double_t ladderX = TMath::Max(xu, xd);
+ Double_t ladderY = yu + yd + gapY;
+ Double_t ladderZ = TMath::Max(zu,
+ zd + shiftZ + offsetZ); // there are 6 slots - 5 x 1.5 mm + 0.3 mm = 7.8 mm
+ // Double_t ladderZ = TMath::Max(zu, zd + shiftZ);
+
+ cout << "DERR iladder " << LadderIndex << " nSec " << nSectors << " zu " << zu << " zd " << zd << " zd+shi "
+ << zd + shiftZ << " ladderZ " << ladderZ << " offsetZ " << offsetZ << endl;
+
+ // --- Place half ladders
+ Double_t xPosU = 0.0; // centred in x
+ Double_t yPosU = 0.5 * (ladderY - yu); // top aligned
+ Double_t zPosU = 0;
+ zPosU = 0.5 * (ladderZ - zu); // front aligned
+ if (LadderIndex >= 1000) zPosU = -0.5 * (ladderZ - zu) + offsetZ; // back aligned with possible offset
+ //zPosU = -zPosU;
+
+ TGeoTranslation* tu = new TGeoTranslation("tu", xPosU, yPosU, zPosU);
+ ladder->AddNode(halfLadderU, 1, tu);
+
+ Double_t xPosD = 0.0; // centred in x
+ Double_t yPosD = 0.5 * -(ladderY - yd); // bottom aligned
+ Double_t zPosD = 0;
+ zPosD = 0.5 * -(ladderZ - zd) + offsetZ; // back aligned with possible offset
+ if (LadderIndex >= 1000) zPosD = -0.5 * -(ladderZ - zd); // front aligned
+ //zPosD = -zPosD;
+
+ TGeoRotation* rd = new TGeoRotation();
+ rd->RotateZ(180.);
+ TGeoCombiTrans* cd = new TGeoCombiTrans(xPosD, yPosD, zPosD, rd);
+ ladder->AddNode(halfLadderD, 2, cd);
+
+ ladder->GetShape()->ComputeBBox();
+
+ shape = (TGeoBBox*) ladder->GetShape();
+ cout << "DDDD0ladder" << LadderIndex << " ladderX " << 2. * shape->GetDX() << " ladderY " << 2. * shape->GetDY()
+ << " ladderZ " << 2. * shape->GetDZ() << endl;
+
+ cout << "DDDD ladder" << LadderIndex << endl;
+ cout << "DDDD1ladder" << LadderIndex << " ladderX " << ladderX << " ladderY " << ladderY << " ladderZ " << ladderZ
+ << endl;
+
+ // ---------------- Create and place frame boxes ------------------------
+
+ if (gkConstructFrames) AddCarbonLadder(LadderIndex, ladder, ladderX, ladderY, ladderZ);
+
+ // --------------------------------------------------------------------------
+
+ shape = (TGeoBBox*) ladder->GetShape();
+ cout << "DDDD2ladder" << LadderIndex << " ladderX " << 2. * shape->GetDX() << " ladderY " << 2. * shape->GetDY()
+ << " ladderZ " << 2. * shape->GetDZ() << endl;
+
+ return ladder;
+}
+/** ======================================================================= **/
+
+
+/** ===========================================================================
+ ** Construct a unit
+ **
+ ** The unit volume is the minimal box comprising all ladders
+ ** minus a tube accomodating the beam pipe.
+ **
+ ** The ladders are arranged horizontally from left to right with
+ ** a given overlap in x.
+ ** Every second ladder is slightly displaced upstream from the centre
+ ** z plane and facing downstream, the others are slightly displaced
+ ** downstream and facing upstream (rotated around the y axis).
+ **
+ ** Arguments:
+ ** name volume name
+ ** nLadders number of ladders
+ ** ladderTypes array of ladder types
+ **/
+
+TGeoVolume* ConstructUnit(Int_t iSide, Int_t iUnit, Int_t nLadders, Int_t* ladderTypes, Int_t iStation)
+{
+
+ Bool_t isFirstPartOfHalfUnit = kFALSE;
+
+ // TString name = Form("Unit%02d", iUnit); // 0,1,2,3,4,5,6,7 - Unit00 missing in output
+ // TString name = Form("Unit%02d", iUnit+1); // 1,2,3,4,5,6,7,8
+
+ TGeoVolume* unit = gGeoMan->GetVolume(unitName[iUnit]);
+ if (!unit) // if it does not yet exist, create a new one
+ {
+ unit = new TGeoVolumeAssembly(unitName[iUnit]);
+ isFirstPartOfHalfUnit = kTRUE;
+ }
+
+ // --- Some local variables
+ TGeoBBox* ladderShape = NULL;
+ TGeoVolume* ladder = NULL;
+ TString ladderName;
+ Double_t subtractedVal;
+
+ // --- Determine size of unit from ladders
+ Double_t statX = 0.;
+ // Double_t statY = 0.;
+
+ for (Int_t iLadder = 0; iLadder < nLadders; iLadder++) {
+ Int_t ladderType = ladderTypes[iLadder]; // v19k
+
+ // if (iSide == 0) cout << "DWER " << ladderTypes[iLadder] << " " << ladderType << endl;
+
+ if (ladderType > 0) {
+ ladderName = Form("LadderType%04d", ladderType); // v19k
+
+ ladder = gGeoManager->GetVolume(ladderName);
+ if (!ladder) Fatal("ConstructUnit", Form("Volume %s not found", ladderName.Data()));
+ ladderShape = (TGeoBBox*) ladder->GetShape();
+ statX += 2. * ladderShape->GetDX();
+ }
+ else
+ statX += gkSensorSizeX; // empty ladder in unit
+ }
+ statX -= Double_t(nLadders - 1) * gkLadderOverlapX;
+
+ // if (iSide == 0) cout << "DWER -" << endl;
+
+ // --- Place ladders in unit
+ cout << "xPos0: " << statX << endl;
+ Double_t xPos = -0.5 * statX;
+ cout << "xPos1: " << xPos << endl;
+ Double_t yPos = 0.;
+ Double_t zPos = 0.;
+
+ Double_t maxdz = 0.;
+ for (Int_t iLadder = 0; iLadder < nLadders; iLadder++) // find maximum dz in this unit
+ {
+ Int_t ladderType = ladderTypes[iLadder]; // v19k
+
+ if (ladderType > 0) {
+ ladderName = Form("LadderType%04d", ladderType); // v19k
+
+ ladder = gGeoManager->GetVolume(ladderName);
+ ladderShape = (TGeoBBox*) ladder->GetShape();
+ if (maxdz < ladderShape->GetDZ()) maxdz = ladderShape->GetDZ();
+ }
+ }
+
+ for (Int_t iLadder = 0; iLadder < nLadders; iLadder++) {
+ Int_t ladderType = ladderTypes[iLadder]; // v19k
+
+ if (ladderType > 0) {
+ ladderName = Form("LadderType%04d", ladderType); // v19k
+
+ ladder = gGeoManager->GetVolume(ladderName);
+ ladderShape = (TGeoBBox*) ladder->GetShape();
+ xPos += ladderShape->GetDX();
+ cout << "xPos2: " << xPos << endl;
+ yPos = 0.; // vertically centred
+ TGeoRotation* rot = new TGeoRotation();
+
+ if (gkConstructFrames)
+ subtractedVal = ladderShape->GetDX() + sqrt(2.) * gkFrameThickness / 2. + gkSectorGapZFrame * 2;
+ else
+ subtractedVal = 0.;
+
+ zPos = 0.5 * gkLadderGapZ + (2 * maxdz - ladderShape->GetDZ() - subtractedVal / 2.); // z-aligned ladders
+
+ // v19k
+ cout << "DE ladder" << ladderTypes[iLadder] << " dx: " << ladderShape->GetDX()
+ << " dy: " << ladderShape->GetDY() << " dz: " << ladderShape->GetDZ() << " max dz: " << maxdz << endl;
+
+ // v19k
+ cout << "DE ladder" << ladderTypes[iLadder] << " fra: " << gkFrameThickness / 2. << " sub: " << subtractedVal
+ << " zpo: " << zPos << endl
+ << endl;
+
+ // v19k
+ if (ladderTypes[iLadder] / 1000 == 1) // flip some of the ladders to reproduce the CAD layout
+ rot->RotateY(180.);
+ else
+ zPos = -zPos;
+
+ if (!isFirstPartOfHalfUnit) zPos += 10.5; // v19d
+ // zPos += 10.0; // initial version
+
+ TGeoCombiTrans* trans = new TGeoCombiTrans(xPos, yPos, zPos, rot);
+ // start
+ // cout << "DEEE** iLadder " << iLadder << " " << nLadders/2 << " " << nLadders << endl;
+
+ if (iSide == 0) {
+ if (iLadder < nLadders / 2) // right side - only half unit -x
+ {
+ unit->AddNode(ladder, iStation * 100 + iLadder + 1, trans);
+ // calculate Station number to encode the ladder copy number
+ cout << "DEFG** iLadder: " << iLadder << " iStation: " << iStation << endl;
+ }
+ }
+ else {
+ if (iLadder >= nLadders / 2) // left side - only half unit +x
+ {
+ unit->AddNode(ladder, iStation * 100 + iLadder + 1, trans);
+ // calculate Station number to encode the ladder copy number
+ cout << "DEFG** iLadder: " << iLadder << " iStation: " << iStation << endl;
+ }
+ }
+ unit->GetShape()->ComputeBBox();
+ // stop
+ xPos += ladderShape->GetDX() - gkLadderOverlapX;
+ cout << "xPos3: " << xPos << endl;
+ }
+ else
+ xPos += gkSensorSizeX - gkLadderOverlapX;
+ }
+
+ return unit;
+}
+/** ======================================================================= **/
+
+
+/** ===========================================================================
+ ** Import and add the passive materials to the STS volume
+ **/
+void ImportPassive(TGeoVolume* stsVolume, TString geoTag, fstream& infoFile)
+{
+ TString passiveName = TString("sts_passive_") + geoTag;
+ TString basePath = gSystem->Getenv("VMCWORKDIR");
+ TString relPath = "/geometry/sts/passive/" + passiveName + ".gdml";
+ TString passiveFileName = basePath + relPath;
+ infoFile << std::endl << std::endl;
+ infoFile << "Importing STS passive materials from GDML file '" << relPath << "'." << std::endl;
+
+ TGDMLParse parser;
+ TGeoVolume* gdmlVolume = parser.GDMLReadFile(passiveFileName);
+ PostProcessGdml(gdmlVolume);
+ gdmlVolume->SetName(passiveName);
+
+ TGeoTranslation* passiveTrans = new TGeoTranslation(0., 0., 4.68 - 2.);
+ infoFile << "Passive assembly is translated for Z=2.68 cm downstream with "
+ "respect to parent volume"
+ << std::endl
+ << std::endl;
+
+ gdmlVolume->GetShape()->ComputeBBox();
+ CheckVolume(gdmlVolume, infoFile);
+
+ infoFile << std::endl;
+ for (Int_t iNode = 0; iNode < gdmlVolume->GetNdaughters(); iNode++) {
+ CheckVolume(gdmlVolume->GetNode(iNode)->GetVolume(), infoFile, kFALSE);
+ }
+
+ stsVolume->AddNode(gdmlVolume, stsVolume->GetNdaughters(), passiveTrans, "");
+}
+
+/** ===========================================================================
+ ** Assign visual properties to the imported gdml volumes
+ **/
+void PostProcessGdml(TGeoVolume* gdmlVolume)
+{
+ const UInt_t kPOBColor = kRed - 6;
+ const UInt_t kPOBTransparency = 0; // 5;
+
+ const UInt_t kFEBColor = kOrange - 6;
+ const UInt_t kFEBTransparency = 0; // 5;
+
+ const UInt_t kUnitColor = kCyan - 10;
+ const UInt_t kUnitTransparency = 0; // 5;
+
+ const UInt_t kCfColor = kCyan - 1;
+ const UInt_t kCfTransparency = 0; // 10;
+
+ // name <Color, Transparency>
+ std::map<std::string, std::tuple<UInt_t, UInt_t>> props {
+ {"passive_POB", std::tuple<UInt_t, UInt_t> {kPOBColor, kPOBTransparency}},
+ {"passive_FEB", std::tuple<UInt_t, UInt_t> {kFEBColor, kFEBTransparency}},
+ {"passive_unit", std::tuple<UInt_t, UInt_t> {kUnitColor, kUnitTransparency}},
+ {"passive_Box_Wall", std::tuple<UInt_t, UInt_t> {kCfColor, kCfTransparency}},
+ {"passive_Box_Wall_Front_CF2", std::tuple<UInt_t, UInt_t> {kCfColor + 3, kCfTransparency}},
+ {"passive_Box_Wall_Back_CF2", std::tuple<UInt_t, UInt_t> {kCfColor + 3, kCfTransparency}},
+ {"passive_Box_Wall_Back_Light_Foam", std::tuple<UInt_t, UInt_t> {kCfColor - 8, kCfTransparency}},
+ {"passive_Box_Wall_Front_Light_Foam", std::tuple<UInt_t, UInt_t> {kCfColor - 8, kCfTransparency}},
+ };
+
+ // Match volume name and apply visual properties passive_Box_Wall_Front_Light_Foam
+ const TObjArray* volumes = gGeoManager->GetListOfVolumes();
+ for (auto& entry : props) {
+ TIter next(volumes);
+ TGeoVolume* vol = nullptr;
+ while ((vol = (TGeoVolume*) next())) {
+ if (TString(vol->GetName()).Contains(entry.first.c_str())) {
+ vol->SetLineColor(std::get<0>(entry.second));
+ vol->SetTransparency(std::get<1>(entry.second));
+ }
+ }
+ }
+}
+
+/** ===========================================================================
+ ** Volume information for debugging
+ **/
+void CheckVolume(TGeoVolume* volume)
+{
+
+ TGeoBBox* shape = (TGeoBBox*) volume->GetShape();
+ cout << volume->GetName() << ": size " << fixed << setprecision(4) << setw(7) << 2. * shape->GetDX() << " x "
+ << setw(7) << 2. * shape->GetDY() << " x " << setw(7) << 2. * shape->GetDZ();
+ if (volume->IsAssembly()) cout << ", assembly";
+ else {
+ if (volume->GetMedium()) cout << ", medium " << volume->GetMedium()->GetName();
+ else
+ cout << ", "
+ << "\033[31m"
+ << " no medium"
+ << "\033[0m";
+ }
+ cout << endl;
+ if (volume->GetNdaughters()) {
+ cout << "Daughters: " << endl;
+ for (Int_t iNode = 0; iNode < volume->GetNdaughters(); iNode++) {
+ TGeoNode* node = volume->GetNode(iNode);
+ TGeoBBox* shape = (TGeoBBox*) node->GetVolume()->GetShape();
+ cout << setw(15) << node->GetName() << ", size " << fixed << setprecision(3) << setw(6) << 2. * shape->GetDX()
+ << " x " << setw(6) << 2. * shape->GetDY() << " x " << setw(6) << 2. * shape->GetDZ() << ", position ( ";
+ TGeoMatrix* matrix = node->GetMatrix();
+ const Double_t* pos = matrix->GetTranslation();
+ cout << setfill(' ');
+ cout << fixed << setw(8) << pos[0] << ", " << setw(8) << pos[1] << ", " << setw(8) << pos[2] << " )" << endl;
+ }
+ }
+}
+/** ======================================================================= **/
+
+/** ===========================================================================
+ ** Volume information for output to file
+ **/
+void CheckVolume(TGeoVolume* volume, fstream& file, Bool_t listChildren)
+{
+ if (!file) return;
+
+ TGeoBBox* shape = (TGeoBBox*) volume->GetShape();
+ file << volume->GetName() << ": size " << fixed << setprecision(4) << setw(7) << 2. * shape->GetDX() << " x "
+ << setw(7) << 2. * shape->GetDY() << " x " << setw(7) << 2. * shape->GetDZ();
+ if (volume->IsAssembly()) file << ", assembly";
+ else {
+ if (volume->GetMedium()) file << ", medium " << volume->GetMedium()->GetName();
+ else
+ file << ", "
+ << "\033[31m"
+ << " no medium"
+ << "\033[0m";
+ }
+ file << endl;
+ if (volume->GetNdaughters() && listChildren) {
+ file << "Contains: ";
+ for (Int_t iNode = 0; iNode < volume->GetNdaughters(); iNode++)
+ file << volume->GetNode(iNode)->GetVolume()->GetName() << " ";
+ file << endl;
+ }
+}
+
+/** ======================================================================= **/
+
+
+/** ===========================================================================
+ ** Calculate beam pipe outer radius for a given z
+ **/
+Double_t BeamPipeRadius(Double_t z)
+{
+ if (z < gkPipeZ2) return gkPipeR1;
+ Double_t slope = (gkPipeR2 - gkPipeR1) / (gkPipeZ2 - gkPipeZ1);
+ return gkPipeR2 + slope * (z - gkPipeZ2);
+}
+/** ======================================================================= **/
+
+
+/** ======================================================================= **/
+TGeoVolume* ConstructFrameElement(const TString& name, TGeoVolume* frameBoxVol, Double_t x)
+{
+ // --- Material of the frames
+ TGeoMedium* framesMaterial = gGeoMan->GetMedium("carbon");
+
+ TGeoBBox* frameVertPillarShp;
+
+ Double_t t = gkFrameThickness / 2.;
+
+ // --- Main vertical pillars
+ // TGeoBBox* frameVertPillarShp = new TGeoBBox(name + "_vertpillar_shape", t, gkFrameStep/2., t); // square crossection, along y
+ // TGeoVolume* frameVertPillarVol = new TGeoVolume(name + "_vertpillar", frameVertPillarShp, framesMaterial);
+ // frameVertPillarVol->SetLineColor(kGreen);
+ // frameBoxVol->AddNode(frameVertPillarVol, 1, new TGeoTranslation(name + "_vertpillar_pos_1", x-t, 0., -(x+sqrt(2.)*t-2.*t)/2.));
+ // frameBoxVol->AddNode(frameVertPillarVol, 2, new TGeoTranslation(name + "_vertpillar_pos_2", -(x-t), 0., -(x+sqrt(2.)*t-2.*t)/2.));
+
+ if (gkCylindricalFrames)
+ // TGeoBBox* frameVertPillarShp = new TGeoTube(name + "_vertpillar_shape", 0, t, gkFrameStep/2.); // circle crossection, along z
+ frameVertPillarShp = new TGeoTube(name + "_vertpillar_shape", gkCylinderDiaInner / 2., gkCylinderDiaOuter / 2.,
+ gkFrameStep / 2.); // circle crossection, along z
+ else
+ frameVertPillarShp = new TGeoBBox(name + "_vertpillar_shape", t, t,
+ gkFrameStep / 2.); // square crossection, along z
+ TGeoVolume* frameVertPillarVol = new TGeoVolume(name + "_vertpillar", frameVertPillarShp, framesMaterial);
+ frameVertPillarVol->SetLineColor(kGreen);
+
+ TGeoRotation* xRot90 = new TGeoRotation;
+ xRot90->RotateX(90.);
+ frameBoxVol->AddNode(
+ frameVertPillarVol, 1,
+ new TGeoCombiTrans(name + "_vertpillar_pos_1", x - t, 0., -(x + sqrt(2.) * t - 2. * t) / 2., xRot90));
+ frameBoxVol->AddNode(
+ frameVertPillarVol, 2,
+ new TGeoCombiTrans(name + "_vertpillar_pos_2", -(x - t), 0., -(x + sqrt(2.) * t - 2. * t) / 2., xRot90));
+
+ // TGeoRotation* vertRot = new TGeoRotation(name + "_vertpillar_rot_1", 90., 45., -90.);
+ TGeoRotation* vertRot = new TGeoRotation;
+ vertRot->RotateX(90.);
+ vertRot->RotateY(45.);
+ frameBoxVol->AddNode(frameVertPillarVol, 3,
+ new TGeoCombiTrans(name + "_vertpillar_pos_3", 0., 0., (x - sqrt(2.) * t) / 2., vertRot));
+
+ // --- Small horizontal pillar
+ // TGeoBBox* frameHorPillarShp = new TGeoBBox(name + "_horpillar_shape", x-2.*t, gkThinFrameThickness/2., gkThinFrameThickness/2.);
+ // TGeoVolume* frameHorPillarVol = new TGeoVolume(name + "_horpillar", frameHorPillarShp, framesMaterial);
+ // frameHorPillarVol->SetLineColor(kCyan);
+ // frameBoxVol->AddNode(frameHorPillarVol, 1, new TGeoTranslation(name + "_horpillar_pos_1", 0., -gkFrameStep/2.+gkThinFrameThickness/2., -(x+sqrt(2.)*t-2.*t)/2.));
+
+ if (gkConstructSmallFrames) {
+
+ // --- Small sloping pillars
+ TGeoPara* frameSlopePillarShp =
+ new TGeoPara(name + "_slopepillar_shape", (x - 2. * t) / TMath::Cos(31.4 / 180. * TMath::Pi()),
+ gkThinFrameThickness / 2., gkThinFrameThickness / 2., 31.4, 0., 90.);
+ TGeoVolume* frameSlopePillarVol = new TGeoVolume(name + "_slopepillar", frameSlopePillarShp, framesMaterial);
+ frameSlopePillarVol->SetLineColor(kCyan);
+ TGeoRotation* slopeRot = new TGeoRotation(name + "_slopepillar_rot_1", 0., 0., 31.4);
+ TGeoRotation* slopeRot2 = new TGeoRotation(name + "_slopepillar_rot_2", 0., 0., -31.4);
+ TGeoCombiTrans* slopeTrRot =
+ new TGeoCombiTrans(name + "_slopepillar_posrot_1", 0., 0., -(x + sqrt(2.) * t - 2. * t) / 2., slopeRot);
+ TGeoCombiTrans* slopeTrRot2 =
+ new TGeoCombiTrans(name + "_slopepillar_posrot_2", 0., 0., -(x + sqrt(2.) * t - 2. * t) / 2., slopeRot2);
+
+ frameBoxVol->AddNode(frameSlopePillarVol, 1, slopeTrRot);
+ frameBoxVol->AddNodeOverlap(frameSlopePillarVol, 2, slopeTrRot2);
+
+
+ Double_t angl = 23.;
+ // --- Small sub pillar
+ TGeoPara* frameSubPillarShp =
+ new TGeoPara(name + "_subpillar_shape", (sqrt(2) * (x / 2. - t) - t / 2.) / TMath::Cos(angl / 180. * TMath::Pi()),
+ gkThinFrameThickness / 2., gkThinFrameThickness / 2., angl, 0., 90.);
+ TGeoVolume* frameSubPillarVol = new TGeoVolume(name + "_subpillar", frameSubPillarShp, framesMaterial);
+ frameSubPillarVol->SetLineColor(kMagenta);
+
+ Double_t posZ = t * (1. - 3. / (2. * sqrt(2.)));
+
+ // one side of X direction
+ TGeoRotation* subRot1 = new TGeoRotation(name + "_subpillar_rot_1", 90., 45., -90. + angl);
+ TGeoCombiTrans* subTrRot1 =
+ new TGeoCombiTrans(name + "_subpillar_posrot_1", -(-x / 2. + t - t / (2. * sqrt(2.))), 1., posZ, subRot1);
+
+ TGeoRotation* subRot2 = new TGeoRotation(name + "_subpillar_rot_2", 90., -90. - 45., -90. + angl);
+ TGeoCombiTrans* subTrRot2 =
+ new TGeoCombiTrans(name + "_subpillar_posrot_2", -(-x / 2. + t - t / (2. * sqrt(2.))), -1., posZ, subRot2);
+
+ // other side of X direction
+ TGeoRotation* subRot3 = new TGeoRotation(name + "_subpillar_rot_3", 90., 90. + 45., -90. + angl);
+ TGeoCombiTrans* subTrRot3 =
+ new TGeoCombiTrans(name + "_subpillar_posrot_3", -x / 2. + t - t / (2. * sqrt(2.)), 1., posZ, subRot3);
+
+ TGeoRotation* subRot4 = new TGeoRotation(name + "_subpillar_rot_4", 90., -45., -90. + angl);
+ TGeoCombiTrans* subTrRot4 =
+ new TGeoCombiTrans(name + "_subpillar_posrot_4", -x / 2. + t - t / (2. * sqrt(2.)), -1., posZ, subRot4);
+
+ frameBoxVol->AddNode(frameSubPillarVol, 1, subTrRot1);
+ frameBoxVol->AddNode(frameSubPillarVol, 2, subTrRot2);
+ frameBoxVol->AddNode(frameSubPillarVol, 3, subTrRot3);
+ frameBoxVol->AddNode(frameSubPillarVol, 4, subTrRot4);
+ // frameBoxVol->GetShape()->ComputeBBox();
+ }
+
+ return frameBoxVol;
+}
+/** ======================================================================= **/
+
+/** ======================================================================= **/
+TGeoVolume* ConstructSmallCone(Double_t coneDz)
+{
+ // --- Material of the frames
+ TGeoMedium* framesMaterial = gGeoMan->GetMedium("carbon");
+
+ // --- Outer cone
+ // TGeoConeSeg* A = new TGeoConeSeg ("A", coneDz, 6., 7.6, 6., 6.04, 0., 180.);
+ // TGeoBBox* B = new TGeoBBox ("B", 8., 6., 10.);
+
+ Double_t radius = 3.0;
+ Double_t thickness = 0.04; // 0.4 mm
+ // TGeoConeSeg* A = new TGeoConeSeg ("A", coneDz, 3., 3.2, 3., 3.2, 0., 180.);
+ TGeoConeSeg* A = new TGeoConeSeg("A", coneDz, radius, radius + thickness, radius, radius + thickness, 0., 180.);
+ TGeoBBox* B = new TGeoBBox("B", 8., 6., 10.);
+
+ TGeoCombiTrans* M = new TGeoCombiTrans("M");
+ M->RotateX(45.);
+ M->SetDy(-5.575);
+ M->SetDz(6.935);
+ M->RegisterYourself();
+
+ TGeoShape* coneShp = new TGeoCompositeShape("Cone_shp", "A-B:M");
+ TGeoVolume* coneVol = new TGeoVolume("Cone", coneShp, framesMaterial);
+ coneVol->SetLineColor(kGreen);
+ // coneVol->RegisterYourself();
+
+ // // --- Inner cone
+ // Double_t thickness = 0.02;
+ // Double_t thickness2 = 0.022;
+ // // TGeoConeSeg* A2 = new TGeoConeSeg ("A2", coneDz-thickness, 6.+thickness, 7.6-thickness2, 5.99+thickness, 6.05-thickness2, 0., 180.);
+ // TGeoConeSeg* A2 = new TGeoConeSeg ("A2", coneDz-thickness, 3.+thickness, 4.6-thickness2, 2.99+thickness, 3.05-thickness2, 0., 180.);
+ //
+ // TGeoCombiTrans* M2 = new TGeoCombiTrans ("M2");
+ // M2->RotateX (45.);
+ // M2->SetDy (-5.575+thickness*sqrt(2.));
+ // M2->SetDz (6.935);
+ // M2->RegisterYourself();
+ //
+ // TGeoShape* coneShp2 = new TGeoCompositeShape ("Cone2_shp", "A2-B:M2");
+ // TGeoVolume* coneVol2 = new TGeoVolume ("Cone2", coneShp2, gStsMedium);
+ // coneVol2->SetLineColor(kGreen);
+ //// coneVol2->RegisterYourself();
+ //
+ // coneVol->AddNode(coneVol2, 1);
+
+ return coneVol;
+}
+/** ======================================================================= **/
+
+/** ======================================================================= **/
+TGeoVolume* ConstructBigCone(Double_t coneDz)
+{
+ // --- Material of the frames
+ TGeoMedium* framesMaterial = gGeoMan->GetMedium("carbon");
+
+ // --- Outer cone
+ TGeoConeSeg* bA = new TGeoConeSeg("bA", coneDz, 6., 7.6, 6., 6.04, 0., 180.);
+ TGeoBBox* bB = new TGeoBBox("bB", 8., 6., 10.);
+
+ TGeoCombiTrans* bM = new TGeoCombiTrans("bM");
+ bM->RotateX(45.);
+ bM->SetDy(-5.575);
+ bM->SetDz(6.935);
+ bM->RegisterYourself();
+
+ TGeoShape* coneBigShp = new TGeoCompositeShape("ConeBig_shp", "bA-bB:bM");
+ TGeoVolume* coneBigVol = new TGeoVolume("ConeBig", coneBigShp, framesMaterial);
+ coneBigVol->SetLineColor(kGreen);
+ // coneBigVol->RegisterYourself();
+
+ // --- Inner cone
+ Double_t thickness = 0.02;
+ Double_t thickness2 = 0.022;
+ TGeoConeSeg* bA2 = new TGeoConeSeg("bA2", coneDz - thickness, 6. + thickness, 7.6 - thickness2, 5.99 + thickness,
+ 6.05 - thickness2, 0., 180.);
+
+ TGeoCombiTrans* bM2 = new TGeoCombiTrans("bM2");
+ bM2->RotateX(45.);
+ bM2->SetDy(-5.575 + thickness * sqrt(2.));
+ bM2->SetDz(6.935);
+ bM2->RegisterYourself();
+
+ TGeoShape* coneBigShp2 = new TGeoCompositeShape("ConeBig2_shp", "bA2-bB:bM2");
+ TGeoVolume* coneBigVol2 = new TGeoVolume("ConeBig2", coneBigShp2, gStsMedium);
+ coneBigVol2->SetLineColor(kGreen);
+ // coneBigVol2->RegisterYourself();
+
+ coneBigVol->AddNode(coneBigVol2, 1);
+
+ return coneBigVol;
+}
+
+/** ======================================================================= **/
--
GitLab