create_rootgeom_GEMS_Run8_detailed_aligned.C

Go to the documentation of this file.

/*
 * Baranov D.
 * 14.08.2020 (created)
 * 18.01.2021 (updated)
 * 02.03.2021 (updated): radius of hole in active zone (4.0 -> 6.0 cm)
 * 12.12.2022 (updated): small GEM added
 * Use this macro to create ROOT geometry for the FutureConfig2020 GEM configuration
 * with passive volumes (frames, al. cable supports, ...)
 *
 */
 
// Set Parameters of GEMS -------------------------------------------------------
const Int_t NStations = 8;     // stations in the detector
const Int_t NMaxModules = 4;   // max. number of modules in a station
 
//(X-Y-Z)Positions of stations (sensitive volumes)
const Double_t XStationPositions[NStations] = {+1.474 /*0*/, +2.614 /*1*/, +3.393 /*2*/, +4.699 /*3*/,
                                               +6.337 /*4*/, +8.189 /*5*/, +9.893 /*6*/, +50.494 /*7*/};
const Double_t YStationPositions[NStations] = {-0.086 /*0*/, -0.107 /*1*/,  +0.003 /*2*/, -0.097 /*3*/,
                                               -0.003 /*4*/, +0.0245 /*5*/, +0.022 /*6*/, -0.186 /*7*/};
const Double_t ZStationPositions[NStations] = {+60.76 /*0*/,   +92.242 /*1*/,  +123.512 /*2*/, +154.692 /*3*/,
                                               +185.989 /*4*/, +217.287 /*5*/, +248.401 /*6*/, +792.532 /*7*/};
 
//(X-Y-Z)Shifts of modules in each station
Double_t XModuleShifts[NStations][NMaxModules] = {
    {163.2 * 0.25, -163.2 * 0.25, 163.2 * 0.25, -163.2 * 0.25},   // st. 0
    {163.2 * 0.25, -163.2 * 0.25, 163.2 * 0.25, -163.2 * 0.25},   // st. 1
    {163.2 * 0.25, -163.2 * 0.25, 163.2 * 0.25, -163.2 * 0.25},   // st. 2
    {163.2 * 0.25, -163.2 * 0.25, 163.2 * 0.25, -163.2 * 0.25},   // st. 3
    {163.2 * 0.25, -163.2 * 0.25, 163.2 * 0.25, -163.2 * 0.25},   // st. 4
    {163.2 * 0.25, -163.2 * 0.25, 163.2 * 0.25, -163.2 * 0.25},   // st. 5
    {163.2 * 0.25, -163.2 * 0.25, 163.2 * 0.25, -163.2 * 0.25},   // st. 6
    {0.0}                                                         // st. 7
};
Double_t YModuleShifts[NStations][NMaxModules] = {
    {+22.5, +22.5, -19.5, -19.5},   // st. 0
    {+22.5, +22.5, -19.5, -19.5},   // st. 1
    {+22.5, +22.5, -19.5, -19.5},   // st. 2
    {+22.5, +22.5, -19.5, -19.5},   // st. 3
    {+22.5, +22.5, -19.5, -19.5},   // st. 4
    {+22.5, +22.5, -19.5, -19.5},   // st. 5
    {+22.5, +22.5, -19.5, -19.5},   // st. 6
    {0.0}                           // st. 7
};
Double_t ZModuleShifts[NStations][NMaxModules] = {
    {+4.1, +4.1, 0.0, 0.0},   // st. 0
    {0.0, 0.0, +4.1, +4.1},   // st. 1
    {+4.1, +4.1, 0.0, 0.0},   // st. 2
    {0.0, 0.0, +4.1, +4.1},   // st. 3
    {+4.1, +4.1, 0.0, 0.0},   // st. 4
    {0.0, 0.0, +4.1, +4.1},   // st. 5
    {+4.1, +4.1, 0.0, 0.0},   // st. 6
    {0.0}                     // st. 7
};
 
// alignment corrections --------------------------------------------------------
const Double_t StsDX[NStations][8] = {
    {2.28613e-01, 2.34606e-01, 3.70449e-01, 4.39267e-01, 2.67664e-01, 2.84859e-01, 3.41005e-01, 4.16753e-01},
    {3.17030e-01, 2.45618e-01, 4.29878e-01, 4.29169e-01, 2.60443e-01, 1.79925e-01, 3.76252e-01, 3.71674e-01},
    {1.25543e-01, 1.37935e-01, 2.77096e-01, 3.51704e-01, 2.16012e-01, 2.12139e-01, 2.88912e-01, 3.57699e-01},
    {3.17647e-01, 2.36295e-01, 3.75431e-01, 3.80575e-01, 1.68884e-01, 8.80524e-02, 2.51805e-01, 2.56513e-01},
    {-3.01883e-02, -2.87880e-02, 1.09744e-01, 1.90429e-01, 1.78162e-02, 1.32671e-02, 9.23760e-02, 1.67779e-01},
    {1.43966e-01, 5.53303e-02, 2.48555e-01, 2.49193e-01, 3.88544e-02, -3.31393e-02, 7.96085e-02, 7.41273e-02},
    {0., -5.63645e-04, 0.05, 8.74027e-04, 0., 3.60763e-03, 0.074, 1.49797e-01}};
const Double_t StsDY[NStations][8] = {
    {-2.40830e-01, -2.27117e-01, -4.41393e-01, -2.98662e-01, -1.47376e-01, -5.64919e-01, 3.60250e-02, -2.15290e-01},
    {-4.83443e-01, -3.63780e-01, -3.52015e-01, -3.36903e-01, -1.14008e-01, -3.49136e-01, -2.14579e-01, -6.23524e-01},
    {-1.91062e-01, -2.04515e-01, -3.79972e-01, -2.57092e-01, -1.01400e-01, -5.18137e-01, 5.39497e-02, -1.85679e-01},
    {-4.55871e-01, -3.48696e-01, -2.33191e-01, -2.40044e-01, -5.37329e-02, -3.00596e-01, -1.98740e-01, -6.04983e-01},
    {-2.31062e-01, -2.58276e-01, -4.03759e-01, -2.70196e-01, -8.53213e-02, -5.14891e-01, 3.15495e-02, -2.08464e-01},
    {-2.58821e-01, -1.58621e-01, -9.84743e-02, -1.05756e-01, 1.54258e-01, -1.04935e-01, 1.83280e-02, -3.86622e-01},
    {0., -5.86937e-03, -0.23, -1.45342e-01, 0., -3.98117e-01, 0.15, -1.04839e-01}};
 
const Double_t XAlignmentCorrections[NStations][NMaxModules] = {
    {StsDX[0][0], StsDX[0][2], StsDX[0][4], StsDX[0][6]},   // st. 0
    {StsDX[1][0], StsDX[1][2], StsDX[1][4], StsDX[1][6]},   // st. 1
    {StsDX[2][0], StsDX[2][2], StsDX[2][4], StsDX[2][6]},   // st. 2
    {StsDX[3][0], StsDX[3][2], StsDX[3][4], StsDX[3][6]},   // st. 3
    {StsDX[4][0], StsDX[4][2], StsDX[4][4], StsDX[4][6]},   // st. 4
    {StsDX[5][0], StsDX[5][2], StsDX[5][4], StsDX[5][6]},   // st. 5
    {StsDX[6][0], StsDX[6][2], StsDX[6][4], StsDX[6][6]},   // st. 6
    {0.0}                                                   // st. 7
};
 
const Double_t YAlignmentCorrections[NStations][NMaxModules] = {
    {StsDY[0][0], StsDY[0][2], StsDY[0][4], StsDY[0][6]},   // st. 0
    {StsDY[1][0], StsDY[1][2], StsDY[1][4], StsDY[1][6]},   // st. 1
    {StsDY[2][0], StsDY[2][2], StsDY[2][4], StsDY[2][6]},   // st. 2
    {StsDY[3][0], StsDY[3][2], StsDY[3][4], StsDY[3][6]},   // st. 3
    {StsDY[4][0], StsDY[4][2], StsDY[4][4], StsDY[4][6]},   // st. 4
    {StsDY[5][0], StsDY[5][2], StsDY[5][4], StsDY[5][6]},   // st. 5
    {StsDY[6][0], StsDY[6][2], StsDY[6][4], StsDY[6][6]},   // st. 6
    {0.0}                                                   // st. 7
};
 
const Double_t ZAlignmentCorrections[NStations][NMaxModules] = {
    {0.243654, 0.250775, 0.173749, 0.161366},           // st. 0
    {-0.143876, -0.150758, 0.109959, 0.156433},         // st. 1
    {-0.110145, -0.165981, -0.213889, -0.229543},       // st. 2
    {0.103517, 0.036342, 0.086135, 0.0941369},          // st. 3
    {-0.037648, -0.0941807, -0.180758, -0.176327},      // st. 4
    {-0.0674491, -0.0367844, -0.0972473, -0.0584354},   // st. 5
    {0, 0, 0, 0},                                       // st. 6
    {0.0}                                               // st. 7
};
//------------------------------------------------------------------------------
 
/*
 * Basic station:
 *
 *  -------------------- --------------------
 * |                    |                    |
 * |  module 0          |  module 1          |    ^ beam direction (Z)
 * |                    |                    |   /
 * |        ------------|--------------      |  /
 * |       /            |             /      | /
 * |      /             |            /       |/
 *  -----------------------------------------
 * |      \             |            \       |
 * |       \            |             \      |
 * |        \___________|______________\     |
 * |                    |                    |
 * |                    |                    |
 * |  module 3          |  module 4          |
 *  -----------------------------------------
 *
 */
 
// rotations of stations around y-axis by 180 deg.
const Bool_t YStationRotations[NStations] = {true /*0*/, false /*1*/, true /*2*/, false /*3*/,
                                             true /*4*/, false /*5*/, true /*6*/, false /*7*/};
 
// rotations of stations around x-axis by 180 deg.
const Bool_t XStationRotations[NStations] = {false /*0*/, false /*1*/, false /*2*/, false /*3*/,
                                             false /*4*/, false /*5*/, false /*6*/, false /*7*/};
 
// GEM plane sizes (66x41 type) -------------------------------------------------
const Double_t XModuleSize_Station66x41 = 66.0;
const Double_t YModuleSize_Station66x41 = 41.0;
const Double_t ZModuleSize_Station66x41 = 0.9;
 
const Double_t dXFrame_Station66x41 = 2.0;
const Double_t dYFrame_Station66x41 = 2.0;
const Double_t dZFrame_Station66x41 = 4.1;
//------------------------------------------------------------------------------
 
// GEM plane sizes (163x45 type) ------------------------------------------------
const Double_t XModuleSize_Station163x45 = 163.2 * 0.5;   // 81.6
const Double_t YModuleSize_Station163x45 = 45.0;
const Double_t ZModuleSize_Station163x45 = 0.9;   // sensitive volume
const Double_t ZGEMSize_Station163x45 = 4.1;      // common thickness of GEM (including sens. vol)
 
const Double_t dXFrame_Station163x45 = 2.0;
const Double_t dYFrame_Station163x45 = 2.0;
const Double_t dZFrame_Station163x45 = 4.1;
//------------------------------------------------------------------------------
 
// GEM plane sizes (163x39 type) ------------------------------------------------
const Double_t XModuleSize_Station163x39 = 163.2 * 0.5;   // 81.6
const Double_t YModuleSize_Station163x39 = 39.0;
const Double_t ZModuleSize_Station163x39 = 0.9;   // sensitive volume
const Double_t ZGEMSize_Station163x39 = 4.1;      // common thickness of GEM (including sens. vol)
 
const Double_t dXFrame_Station163x39 = 2.0;
const Double_t dYFrame_Station163x39 = 2.0;
const Double_t dZFrame_Station163x39 = 4.1;
//------------------------------------------------------------------------------
 
// GEM plane sizes (10x10 type) -------------------------------------------------
const Double_t XModuleSize_Station10x10 = 10.0;   // 81.6
const Double_t YModuleSize_Station10x10 = 10.0;
const Double_t ZModuleSize_Station10x10 = 0.9;   // sensitive volume
const Double_t ZGEMSize_Station10x10 = 4.1;      // common thickness of GEM (including sens. vol)
 
const Double_t dXFrame_Station10x10 = 2.0;
const Double_t dYFrame_Station10x10 = 2.0;
const Double_t dZFrame_Station10x10 = 4.1;
//------------------------------------------------------------------------------
 
// GeoManager
TGeoManager* gGeoMan = NULL;
 
// media
TGeoMedium* pMedAir = 0;
TGeoMedium* pMedCarbon = 0;
TGeoMedium* pMedFiberGlass = 0;
TGeoMedium* pMedCopper = 0;
TGeoMedium* pMedArCO27030 = 0;
TGeoMedium* pMedArC4H10_80_20 = 0;
TGeoMedium* pMedAluminium = 0;
TGeoMedium* pMedEpoxideCompound = 0;
TGeoMedium* pMedGlue = 0;
TGeoMedium* pMedBrass = 0;
TGeoMedium* pMedSteel = 0;
TGeoMedium* pMedHoneyComb = 0;
TGeoMedium* pMedPolyacetal = 0;
TGeoMedium* pMedPolyamide = 0;
 
class FairGeoMedia;
class FairGeoBuilder;
 
TGeoVolume* CreateStation(TString station_name);
TGeoVolume* CreateModule_Station66x41(TString module_name, Double_t xsize, Double_t ysize, Double_t zsize);
 
TGeoVolume* CreateFrameForModule_Station66x41(TString frame_name, Double_t dx, Double_t dy, Double_t dz);
TGeoVolume* CreateModule_Station163x45(TString module_name,
                                       Double_t xsize,
                                       Double_t ysize,
                                       Double_t zsize,
                                       Double_t hole_radius);
TGeoVolume* CreateLayersForModule_Station163x45(TString layer_name);
TGeoVolume* CreateFrameForModule_Station163x45(TString frame_name,
                                               Double_t dx,
                                               Double_t dy,
                                               Double_t dz,
                                               Double_t hole_radius);
TGeoVolume* CreateFrameForHalfPlane_Station163x45(TString frame_name, Int_t istation);
 
TGeoVolume* CreateModule_Station163x39(TString module_name,
                                       Double_t xsize,
                                       Double_t ysize,
                                       Double_t zsize,
                                       Double_t hole_radius);
TGeoVolume* CreateLayersForModule_Station163x39(TString layer_name);
TGeoVolume* CreateFrameForModule_Station163x39(TString frame_name,
                                               Double_t dx,
                                               Double_t dy,
                                               Double_t dz,
                                               Double_t hole_radius);
TGeoVolume* CreateFrameForHalfPlane_Station163x39(TString frame_name, Int_t istation);
 
TGeoVolume* CreateModule_Station10x10(TString module_name, Double_t xsize, Double_t ysize, Double_t zsize);
 
void DefineRequiredMedia(FairGeoMedia* geoMedia, FairGeoBuilder* geoBuild)
{
 
    // air medium
    FairGeoMedium* mAir = geoMedia->getMedium("air");
    if (!mAir)
        Fatal("Main", "FairMedium air not found");
    geoBuild->createMedium(mAir);
    pMedAir = gGeoManager->GetMedium("air");
    if (!pMedAir)
        Fatal("Main", "Medium air not found");
 
    // carbon medium
    FairGeoMedium* mCarbon = geoMedia->getMedium("carbon");
    if (!mCarbon)
        Fatal("Main", "FairMedium carbon not found");
    geoBuild->createMedium(mCarbon);
    pMedCarbon = gGeoManager->GetMedium("carbon");
    if (!pMedCarbon)
        Fatal("Main", "Medium carbon not found");
 
    // fiberglass medium
    FairGeoMedium* mFiberGlass = geoMedia->getMedium("fiberglass");
    if (!mFiberGlass)
        Fatal("Main", "FairMedium fiberglass not found");
    geoBuild->createMedium(mFiberGlass);
    pMedFiberGlass = gGeoManager->GetMedium("fiberglass");
    if (!pMedFiberGlass)
        Fatal("Main", "Medium fiberglass not found");
 
    // copper medium
    FairGeoMedium* mCopper = geoMedia->getMedium("copper");
    if (!mCopper)
        Fatal("Main", "FairMedium copper not found");
    geoBuild->createMedium(mCopper);
    pMedCopper = gGeoManager->GetMedium("copper");
    if (!pMedCopper)
        Fatal("Main", "Medium copper not found");
 
    // arco27020 medium
    FairGeoMedium* mArCO27030 = geoMedia->getMedium("arco27030");
    if (!mArCO27030)
        Fatal("Main", "FairMedium arco27030 not found");
    geoBuild->createMedium(mArCO27030);
    pMedArCO27030 = gGeoManager->GetMedium("arco27030");
    if (!pMedArCO27030)
        Fatal("Main", "Medium arco27030 not found");
 
    // ArC4H10_80_20 medium
    FairGeoMedium* mArC4H10_80_20 = geoMedia->getMedium("ArC4H10_80_20");
    if (!mArC4H10_80_20)
        Fatal("Main", "FairMedium ArC4H10_80_20 not found");
    geoBuild->createMedium(mArC4H10_80_20);
    pMedArC4H10_80_20 = gGeoManager->GetMedium("ArC4H10_80_20");
    if (!pMedArC4H10_80_20)
        Fatal("Main", "Medium ArC4H10_80_20 not found");
 
    // aluminium medium
    FairGeoMedium* mAluminium = geoMedia->getMedium("aluminium");
    if (!mAluminium)
        Fatal("Main", "FairMedium aluminium not found");
    geoBuild->createMedium(mAluminium);
    pMedAluminium = gGeoManager->GetMedium("aluminium");
    if (!pMedAluminium)
        Fatal("Main", "Medium aluminium not found");
 
    // epoxide compound (for GEM frames)
    FairGeoMedium* mEpoxideCompound = geoMedia->getMedium("polyurethane_hd");
    if (!mEpoxideCompound)
        Fatal("Main", "FairMedium epoxideCompound not found");
    geoBuild->createMedium(mEpoxideCompound);
    pMedEpoxideCompound = gGeoManager->GetMedium("polyurethane_hd");
    if (!pMedEpoxideCompound)
        Fatal("Main", "Medium epoxideCompound not found");
 
    // glue
    FairGeoMedium* mGlue = geoMedia->getMedium("acrylic_glue");
    if (!mGlue)
        Fatal("Main", "FairMedium acrylic_glue not found");
    geoBuild->createMedium(mGlue);
    pMedGlue = gGeoManager->GetMedium("acrylic_glue");
    if (!pMedGlue)
        Fatal("Main", "Medium acrylic_glue not found");
 
    // brass medium
    FairGeoMedium* mBrass = geoMedia->getMedium("brass");
    if (!mBrass)
        Fatal("Main", "FairMedium brass not found");
    geoBuild->createMedium(mBrass);
    pMedBrass = gGeoManager->GetMedium("brass");
    if (!pMedBrass)
        Fatal("Main", "Medium brass not found");
 
    // steel medium
    FairGeoMedium* mSteel = geoMedia->getMedium("steel");
    if (!mSteel)
        Fatal("Main", "FairMedium steel not found");
    geoBuild->createMedium(mSteel);
    pMedSteel = gGeoManager->GetMedium("steel");
    if (!pMedSteel)
        Fatal("Main", "Medium steel not found");
 
    // honeycomb medium
    FairGeoMedium* mHoneyComb = geoMedia->getMedium("AramidHoneyComb");
    if (!mHoneyComb)
        Fatal("Main", "FairMedium AramidHoneyComb not found");
    geoBuild->createMedium(mHoneyComb);
    pMedHoneyComb = gGeoManager->GetMedium("AramidHoneyComb");
    if (!pMedHoneyComb)
        Fatal("Main", "Medium AramidHoneyComb not found");
 
    // polyacetal medium ((polyoxymethylene))
    FairGeoMedium* mPolyacetal = geoMedia->getMedium("polyacetal");
    if (!mPolyacetal)
        Fatal("Main", "FairMedium polyacetal not found");
    geoBuild->createMedium(mPolyacetal);
    pMedPolyacetal = gGeoManager->GetMedium("polyacetal");
    if (!pMedPolyacetal)
        Fatal("Main", "Medium polyacetal not found");
 
    // polyamide medium
    FairGeoMedium* mPolyamide = geoMedia->getMedium("polyamide");
    if (!mPolyamide)
        Fatal("Main", "FairMedium polyamide not found");
    geoBuild->createMedium(mPolyamide);
    pMedPolyamide = gGeoManager->GetMedium("polyamide");
    if (!pMedPolyamide)
        Fatal("Main", "Medium polyamide not found");
}
 
void create_rootgeom_GEMS_Run8_detailed_aligned()
{
 
    // ----  set working directory  --------------------------------------------
    TString gPath = gSystem->Getenv("VMCWORKDIR");
 
    // -------   Geometry file name (output)   ---------------------------------
    const TString geoDetectorName = "GEMS";
    const TString geoDetectorVersion = "Run8_detailed_aligned";
    const TString geoFileName = gPath + "/geometry/" + geoDetectorName + "_" + geoDetectorVersion + ".root";
 
    // ----  global geometry parameters  ---------------------------------------
    FairGeoLoader* geoLoad = new FairGeoLoader("TGeo", "FairGeoLoader");
    FairGeoInterface* geoFace = geoLoad->getGeoInterface();
 
    // -------   Load media from media file   ----------------------------------
    TString medFile = gPath + "/geometry/media.geo";
    geoFace->setMediaFile(medFile);
    geoFace->readMedia();
    FairGeoMedia* geoMedia = geoFace->getMedia();
    FairGeoBuilder* geoBuild = geoLoad->getGeoBuilder();
 
    DefineRequiredMedia(geoMedia, geoBuild);
    // -------------------------------------------------------------------------
 
    // --------------   Create geometry and top volume  ------------------------
    gGeoMan = (TGeoManager*)gROOT->FindObject("FAIRGeom");
    gGeoMan->SetName(geoDetectorName + "_geom");
    TGeoVolume* top = new TGeoVolumeAssembly("TOP");
    top->SetMedium(pMedAir);
    gGeoMan->SetTopVolume(top);
    // gGeoMan->SetTopVisible(1);
    //  -------------------------------------------------------------------------
 
    // Define TOP Geometry
    TGeoVolume* GEMS = new TGeoVolumeAssembly(geoDetectorName);
    GEMS->SetMedium(pMedAir);
 
    // for(Int_t istation = 0; istation < NStations; ++istation) {
    for (Int_t istation = 0; istation < 8; ++istation) {
 
        Int_t stationNum = istation;
        TGeoVolume* station = CreateStation(TString("station") + TString::Itoa(stationNum, 10));
 
        // Applying alignment corrections to module shifts !!!!!!!!!!
        for (Int_t imodule = 0; imodule < NMaxModules; ++imodule) {
            XModuleShifts[istation][imodule] -= XAlignmentCorrections[istation][imodule];
            YModuleShifts[istation][imodule] -= YAlignmentCorrections[istation][imodule];
            ZModuleShifts[istation][imodule] += ZAlignmentCorrections[istation][imodule];
        }
 
        if (istation < 7) {   // big GEMs
            // upper part
            TGeoVolume* module0 =
                CreateModule_Station163x45(TString("Sensor_module0_") + station->GetName(), XModuleSize_Station163x45,
                                           YModuleSize_Station163x45, ZModuleSize_Station163x45, 4.0 + 2.0);
            TGeoVolume* module1 =
                CreateModule_Station163x45(TString("Sensor_module1_") + station->GetName(), XModuleSize_Station163x45,
                                           YModuleSize_Station163x45, ZModuleSize_Station163x45, 4.0 + 2.0);
 
            TGeoVolume* layers0 = CreateLayersForModule_Station163x45(TString("layer0_") + station->GetName());
            TGeoVolume* layers1 = CreateLayersForModule_Station163x45(TString("layer1_") + station->GetName());
 
            TGeoVolume* frame_upper =
                CreateFrameForHalfPlane_Station163x45(TString("frame_upper_") + station->GetName(), istation);
 
            // lower part
            TGeoVolume* module2 =
                CreateModule_Station163x39(TString("Sensor_module2_") + station->GetName(), XModuleSize_Station163x39,
                                           YModuleSize_Station163x39, ZModuleSize_Station163x39, 4.0 + 2.0);
            TGeoVolume* module3 =
                CreateModule_Station163x39(TString("Sensor_module3_") + station->GetName(), XModuleSize_Station163x39,
                                           YModuleSize_Station163x39, ZModuleSize_Station163x39, 4.0 + 2.0);
 
            TGeoVolume* layers2 = CreateLayersForModule_Station163x39(TString("layer3_") + station->GetName());
            TGeoVolume* layers3 = CreateLayersForModule_Station163x39(TString("layer4_") + station->GetName());
 
            TGeoVolume* frame_lower =
                CreateFrameForHalfPlane_Station163x39(TString("frame_lower_") + station->GetName(), istation);
 
            TGeoCombiTrans* module0_transform = new TGeoCombiTrans();
            module0_transform->SetTranslation(XModuleShifts[stationNum][0], YModuleShifts[stationNum][0],
                                              ZModuleShifts[stationNum][0] + ZModuleSize_Station163x45 * 0.5);
 
            TGeoCombiTrans* module1_transform = new TGeoCombiTrans();
            module1_transform->ReflectX(true);
            module1_transform->SetTranslation(XModuleShifts[stationNum][1], YModuleShifts[stationNum][1],
                                              ZModuleShifts[stationNum][1] + ZModuleSize_Station163x45 * 0.5);
 
            TGeoCombiTrans* module2_transform = new TGeoCombiTrans();
            module2_transform->ReflectY(true);
            module2_transform->SetTranslation(XModuleShifts[stationNum][2], YModuleShifts[stationNum][2],
                                              ZModuleShifts[stationNum][2] + ZModuleSize_Station163x39 * 0.5);
 
            TGeoCombiTrans* module3_transform = new TGeoCombiTrans();
            module3_transform->ReflectY(true);
            module3_transform->ReflectX(true);
            module3_transform->SetTranslation(XModuleShifts[stationNum][3], YModuleShifts[stationNum][3],
                                              ZModuleShifts[stationNum][3] + ZModuleSize_Station163x39 * 0.5);
 
            TGeoCombiTrans* layers0_transform = new TGeoCombiTrans();
            if (YStationRotations[stationNum] == true)
                layers0_transform->ReflectZ(true);
            if (XStationRotations[stationNum] == true)
                layers0_transform->RotateX(180.0);
            layers0_transform->SetTranslation(XModuleShifts[stationNum][0], YModuleShifts[stationNum][0],
                                              ZModuleShifts[stationNum][0] + ZModuleSize_Station163x45 * 0.5);
 
            TGeoCombiTrans* layers1_transform = new TGeoCombiTrans();
            if (YStationRotations[stationNum] == true)
                layers1_transform->ReflectZ(true);
            if (XStationRotations[stationNum] == true)
                layers1_transform->RotateX(180.0);
            layers1_transform->ReflectX(true);
            layers1_transform->SetTranslation(XModuleShifts[stationNum][1], YModuleShifts[stationNum][1],
                                              ZModuleShifts[stationNum][1] + ZModuleSize_Station163x45 * 0.5);
 
            TGeoCombiTrans* layers2_transform = new TGeoCombiTrans();
            layers2_transform->ReflectY(true);
            if (YStationRotations[stationNum] == false)
                layers2_transform->ReflectZ(true);
            if (XStationRotations[stationNum] == true)
                layers2_transform->RotateX(180.0);
            layers2_transform->SetTranslation(XModuleShifts[stationNum][2], YModuleShifts[stationNum][2],
                                              ZModuleShifts[stationNum][2] + ZModuleSize_Station163x45 * 0.5);
 
            TGeoCombiTrans* layers3_transform = new TGeoCombiTrans();
            layers3_transform->ReflectY(true);
            layers3_transform->ReflectX(true);
            if (YStationRotations[stationNum] == false)
                layers3_transform->ReflectZ(true);
            if (XStationRotations[stationNum] == true)
                layers3_transform->RotateX(180.0);
            layers3_transform->SetTranslation(XModuleShifts[stationNum][3], YModuleShifts[stationNum][3],
                                              ZModuleShifts[stationNum][3] + ZModuleSize_Station163x45 * 0.5);
 
            TGeoCombiTrans* frame_upper_transform = new TGeoCombiTrans();
            if (YStationRotations[stationNum] == true)
                frame_upper_transform->RotateY(180.0);
            if (XStationRotations[stationNum] == true)
                frame_upper_transform->RotateX(180.0);
            frame_upper_transform->SetTranslation((XModuleShifts[stationNum][0] + XModuleShifts[stationNum][1]),
                                                  YModuleShifts[stationNum][0],
                                                  ZModuleShifts[stationNum][0] + ZModuleSize_Station163x45 * 0.5);
 
            TGeoCombiTrans* frame_lower_transform = new TGeoCombiTrans();
            frame_lower_transform->ReflectY(true);
            if (YStationRotations[stationNum] == false)
                frame_lower_transform->RotateY(180.0);
            if (XStationRotations[stationNum] == true)
                frame_lower_transform->RotateX(180.0);
            frame_lower_transform->SetTranslation((XModuleShifts[stationNum][2] + XModuleShifts[stationNum][3]),
                                                  YModuleShifts[stationNum][2],
                                                  ZModuleShifts[stationNum][2] + ZModuleSize_Station163x45 * 0.5);
 
            station->AddNode(module0, 0, new TGeoCombiTrans(*module0_transform));   // module
            station->AddNode(module1, 0, new TGeoCombiTrans(*module1_transform));   // module
            station->AddNode(module2, 0, new TGeoCombiTrans(*module2_transform));   // module
            station->AddNode(module3, 0, new TGeoCombiTrans(*module3_transform));   // module
 
            station->AddNode(layers0, 0, new TGeoCombiTrans(*layers0_transform));   // layers
            station->AddNode(layers1, 0, new TGeoCombiTrans(*layers1_transform));   // layers
            station->AddNode(layers2, 0, new TGeoCombiTrans(*layers2_transform));   // layers
            station->AddNode(layers3, 0, new TGeoCombiTrans(*layers3_transform));   // layers
 
            station->AddNode(frame_upper, 0, frame_upper_transform);
            station->AddNode(frame_lower, 0, frame_lower_transform);
        }
 
        if (istation == 7) {   // small GEM
            TGeoVolume* module0 =
                CreateModule_Station10x10(TString("Sensor_module0_") + station->GetName(), XModuleSize_Station10x10,
                                          YModuleSize_Station10x10, ZModuleSize_Station10x10);
 
            TGeoCombiTrans* module0_transform = new TGeoCombiTrans();
            module0_transform->SetTranslation(XModuleShifts[stationNum][0], YModuleShifts[stationNum][0],
                                              ZModuleShifts[stationNum][0] + ZModuleSize_Station10x10 * 0.5);
 
            station->AddNode(module0, 0, new TGeoCombiTrans(*module0_transform));   // module
        }
 
        TGeoCombiTrans* station_transform = new TGeoCombiTrans();
        station_transform->SetTranslation(XStationPositions[stationNum], YStationPositions[stationNum],
                                          ZStationPositions[stationNum]);
 
        GEMS->AddNode(station, 0, station_transform);
    }
 
    top->AddNode(GEMS, 0);
    top->SetVisContainers(kTRUE);
 
    // ---------------   Finish   -----------------------------------------------
    gGeoMan->CloseGeometry();
    // gGeoMan->CheckOverlaps(0.0001);
    // gGeoMan->PrintOverlaps();
    gGeoMan->Test();
    // gGeoMan->RestoreMasterVolume();
 
    TFile* geoFile = new TFile(geoFileName, "RECREATE");
    top->Write();
    geoFile->Close();
    top->Draw("ogl");
}
 
//------------------------------------------------------------------------------
 
TGeoVolume* CreateStation(TString station_name)
{
 
    TGeoVolume* stationA = new TGeoVolumeAssembly(station_name);
    stationA->SetMedium(pMedAir);
 
    return stationA;
}
 
TGeoVolume* CreateModule_Station66x41(TString module_name, Double_t xsize, Double_t ysize, Double_t zsize)
{
 
    // module shape
    TGeoShape* moduleS = new TGeoBBox("moduleS", xsize * 0.5, ysize * 0.5, zsize * 0.5);
 
    // module volume
    TGeoVolume* moduleV = new TGeoVolume(module_name, moduleS);
 
    // medium
    if (pMedArC4H10_80_20)
        moduleV->SetMedium(pMedArC4H10_80_20);
    else
        Fatal("Main", "Invalid medium for sensitive modules!");
 
    // visual parameters
    moduleV->SetLineColor(TColor::GetColor("#ff47ca"));
    moduleV->SetTransparency(20);
 
    return moduleV;
}
 
TGeoVolume* CreateFrameForModule_Station66x41(TString frame_name, Double_t dx, Double_t dy, Double_t dz)
{
 
    // frame shapes
    TGeoShape* vertical_frameS =
        new TGeoBBox("vertical_frameS", dx * 0.5, YModuleSize_Station66x41 * 0.5 + dy, dz * 0.5);
    TGeoShape* horizontal_frameS =
        new TGeoBBox("horizontal_frameS", XModuleSize_Station66x41 * 0.5, dy * 0.5, dz * 0.5);
 
    // frame volumes
    TGeoVolume* composite_frameV = new TGeoVolumeAssembly(frame_name);
    TGeoVolume* vertical_frameV = new TGeoVolume(frame_name + "_vertical_frameV", vertical_frameS);
    TGeoVolume* horizontal_frameV = new TGeoVolume(frame_name + "_horizontal_frameV", horizontal_frameS);
 
    // media
    TGeoMedium* frame_medium = pMedCarbon;   // set medium
    if (pMedCarbon) {
        vertical_frameV->SetMedium(frame_medium);
        horizontal_frameV->SetMedium(frame_medium);
    } else
        Fatal("Main", "Invalid medium for frames!");
    if (pMedAir) {
        composite_frameV->SetMedium(pMedAir);
    } else
        Fatal("Main", "Invalid medium for frames!");
 
    // visual parameters
    vertical_frameV->SetLineColor(TColor::GetColor("#9999ff"));
    horizontal_frameV->SetLineColor(TColor::GetColor("#9999ff"));
 
    // frame part positions in the frame
    TGeoCombiTrans* left_vertical_frame_position = new TGeoCombiTrans();
    left_vertical_frame_position->SetTranslation(XModuleSize_Station66x41 * 0.5 + dXFrame_Station66x41 * 0.5, 0.0, 0.0);
 
    TGeoCombiTrans* right_vertical_frame_position = new TGeoCombiTrans();
    right_vertical_frame_position->SetTranslation(-XModuleSize_Station66x41 * 0.5 - dXFrame_Station66x41 * 0.5, 0.0,
                                                  0.0);
 
    TGeoCombiTrans* bottom_horizontal_frame_position = new TGeoCombiTrans();
    bottom_horizontal_frame_position->SetTranslation(0.0, YModuleSize_Station66x41 * 0.5 + dYFrame_Station66x41 * 0.5,
                                                     0.0);
 
    TGeoCombiTrans* lower_horizontal_frame_position = new TGeoCombiTrans();
    lower_horizontal_frame_position->SetTranslation(0.0, -YModuleSize_Station66x41 * 0.5 - dYFrame_Station66x41 * 0.5,
                                                    0.0);
 
    composite_frameV->AddNode(vertical_frameV, 0, left_vertical_frame_position);
    composite_frameV->AddNode(vertical_frameV, 1, right_vertical_frame_position);
    composite_frameV->AddNode(horizontal_frameV, 0, bottom_horizontal_frame_position);
    composite_frameV->AddNode(horizontal_frameV, 1, lower_horizontal_frame_position);
 
    return composite_frameV;
}
 
TGeoVolume* CreateModule_Station163x45(TString module_name,
                                       Double_t xsize,
                                       Double_t ysize,
                                       Double_t zsize,
                                       Double_t hole_radius)
{
 
    gRandom->SetSeed(0);
    Int_t uniqueId = gRandom->Integer(100000);
 
    // module shapes
    TGeoShape* module_partS =
        new TGeoBBox(TString("module_partS") + TString::Itoa(uniqueId, 10), xsize * 0.5, ysize * 0.5, zsize * 0.5);
    TGeoShape* holeS =
        new TGeoTube(TString("holeS") + TString::Itoa(uniqueId, 10), 0.0, hole_radius, zsize * 0.5 + 0.01);
 
    TGeoTranslation* hole_module_trans = new TGeoTranslation(
        TString("hole_module_trans_") + TString::Itoa(uniqueId, 10), -xsize * 0.5, -ysize * 0.5, 0.0);
 
    hole_module_trans->RegisterYourself();
 
    TString bool_expression_module = TString(module_partS->GetName()) + TString(" - ") + TString(holeS->GetName())
                                     + TString(":") + TString(hole_module_trans->GetName());
    TGeoCompositeShape* moduleS = new TGeoCompositeShape(bool_expression_module);
 
    // module volume
    TGeoVolume* moduleV = new TGeoVolume(module_name, moduleS);
 
    // medium
    if (pMedArC4H10_80_20)
        moduleV->SetMedium(pMedArC4H10_80_20);
    else
        Fatal("Main", "Invalid medium for sensitive modules!");
 
    // visual parameters
    moduleV->SetLineColor(TColor::GetColor("#ff47ca"));
    moduleV->SetTransparency(20);
 
    return moduleV;
}
//------------------------------------------------------------------------------
 
TGeoVolume* CreateLayersForModule_Station163x45(TString layer_name)
{
 
    // frame volumes
    TGeoVolume* layers = new TGeoVolumeAssembly(layer_name);
    layers->SetMedium(pMedAir);
 
    // Common parameters for all layers
    Double_t layer_XSize = XModuleSize_Station163x45;   // cm
    Double_t layer_YSize = YModuleSize_Station163x45;   // cm
 
    Double_t layerFront_XSize = 80.85;     // cm
    Double_t layerFront_YSize = 44.2;      // cm
    Double_t layerFront_XShift = -0.375;   // cm
    Double_t layerFront_YShift = -0.2;     // cm
 
    Double_t layerBack_XSize = 80.71;     // cm
    Double_t layerBack_YSize = 44.2;      // cm
    Double_t layerBack_XShift = -0.445;   // cm
    Double_t layerBack_YShift = -0.2;     // cm
 
    // copper layer -------------------------------------------------------------
    Double_t copperLayer_ZSize = 0.00655;    // cm;
    Double_t copperLayerHole_Radius = 6.0;   // cm
 
    // layer shape
    TGeoShape* copperLayerFrontBlankS =
        new TGeoBBox(TString("copperLayerFrontBlankS") += TString("_") + layers->GetName(), layerFront_XSize * 0.5,
                     layerFront_YSize * 0.5, copperLayer_ZSize * 0.5);
    TGeoShape* copperLayerBackBlankS =
        new TGeoBBox(TString("copperLayerBackBlankS") += TString("_") + layers->GetName(), layerBack_XSize * 0.5,
                     layerBack_YSize * 0.5, copperLayer_ZSize * 0.5);
 
    TGeoShape* copperLayerFrontHoleS =
        new TGeoTube(TString("copperLayerFrontHoleS") += TString("_") + layers->GetName(), 0.0, copperLayerHole_Radius,
                     copperLayer_ZSize * 0.5 + 0.01);
    TGeoShape* copperLayerBackHoleS = new TGeoTube(TString("copperLayerBackHoleS") += TString("_") + layers->GetName(),
                                                   0.0, copperLayerHole_Radius, copperLayer_ZSize * 0.5 + 0.01);
 
    TGeoTranslation* copperLayerFrontHole_pos = new TGeoTranslation();
    copperLayerFrontHole_pos->SetName(TString("copperLayerFrontHole_pos") += TString("_") + layers->GetName());
    copperLayerFrontHole_pos->SetDx(-XModuleSize_Station163x45 * 0.5 - layerFront_XShift);
    copperLayerFrontHole_pos->SetDy(-YModuleSize_Station163x45 * 0.5 - layerFront_YShift);
    copperLayerFrontHole_pos->SetDz(0.0);
    copperLayerFrontHole_pos->RegisterYourself();
 
    TGeoTranslation* copperLayerBackHole_pos = new TGeoTranslation();
    copperLayerBackHole_pos->SetName(TString("copperLayerBackHole_pos") += TString("_") + layers->GetName());
    copperLayerBackHole_pos->SetDx(-XModuleSize_Station163x45 * 0.5 - layerBack_XShift);
    copperLayerBackHole_pos->SetDy(-YModuleSize_Station163x45 * 0.5 - layerBack_YShift);
    copperLayerBackHole_pos->SetDz(0.0);
    copperLayerBackHole_pos->RegisterYourself();
 
    TGeoCompositeShape* copperLayerFrontS = new TGeoCompositeShape();
    copperLayerFrontS->SetName(TString("copperLayerFrontS") += TString("_") + layers->GetName());
    {
        TString expression = "copperLayerFrontBlankS";
        expression += TString("_") + layers->GetName();
        expression += "-copperLayerFrontHoleS";
        expression += TString("_") + layers->GetName();
        expression += ":copperLayerFrontHole_pos";
        expression += TString("_") + layers->GetName();
 
        copperLayerFrontS->MakeNode(expression);
        copperLayerFrontS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoCompositeShape* copperLayerBackS = new TGeoCompositeShape();
    copperLayerBackS->SetName(TString("copperLayerBackS") += TString("_") + layers->GetName());
    {
        TString expression = "copperLayerBackBlankS";
        expression += TString("_") + layers->GetName();
        expression += "-copperLayerBackHoleS";
        expression += TString("_") + layers->GetName();
        expression += ":copperLayerBackHole_pos";
        expression += TString("_") + layers->GetName();
 
        copperLayerBackS->MakeNode(expression);
        copperLayerBackS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* copperLayerFrontV =
        new TGeoVolume(TString("copperLayerFrontV") += TString("_") + layers->GetName(), copperLayerFrontS);
    TGeoVolume* copperLayerBackV =
        new TGeoVolume(TString("copperLayerBackV") += TString("_") + layers->GetName(), copperLayerBackS);
 
    // volume medium
    TGeoMedium* copperLayerV_medium = pMedCopper;
    if (copperLayerV_medium) {
        copperLayerFrontV->SetMedium(copperLayerV_medium);
        copperLayerBackV->SetMedium(copperLayerV_medium);
    } else
        Fatal("Main", "Invalid medium for copperLayerV!");
 
    // volume visual property (transparency)
    copperLayerFrontV->SetLineColor(TColor::GetColor("#ff7538"));
    copperLayerFrontV->SetTransparency(0);
    copperLayerBackV->SetLineColor(TColor::GetColor("#ff7538"));
    copperLayerBackV->SetTransparency(0);
 
    TGeoCombiTrans* copperLayerFront_transf[2];
 
    copperLayerFront_transf[0] = new TGeoCombiTrans();
    copperLayerFront_transf[0]->SetDx(+layerFront_XShift);
    copperLayerFront_transf[0]->SetDy(+layerFront_YShift);
    copperLayerFront_transf[0]->SetDz(-ZGEMSize_Station163x45 * 0.5 + copperLayer_ZSize * 0.5);
 
    TGeoCombiTrans* copperLayerBack_transf[1];
 
    copperLayerBack_transf[0] = new TGeoCombiTrans();
    copperLayerBack_transf[0]->SetDx(+layerBack_XShift);
    copperLayerBack_transf[0]->SetDy(+layerBack_YShift);
    copperLayerBack_transf[0]->SetDz(+ZGEMSize_Station163x45 * 0.5 - copperLayer_ZSize * 0.5);
 
    layers->AddNode(copperLayerFrontV, 0, copperLayerFront_transf[0]);
    layers->AddNode(copperLayerBackV, 0, copperLayerBack_transf[0]);
    //--------------------------------------------------------------------------
 
    // glue layer ---------------------------------------------------------------
    Double_t glueLayer_ZSize = 0.01;       // cm; half-thickness of all glue layers
    Double_t glueLayerHole_Radius = 6.0;   // cm
 
    // layer shape
    TGeoShape* glueLayerFrontBlankS =
        new TGeoBBox(TString("glueLayerFrontBlankS") += TString("_") + layers->GetName(), layerFront_XSize * 0.5,
                     layerFront_YSize * 0.5, glueLayer_ZSize * 0.5);
    TGeoShape* glueLayerBackBlankS = new TGeoBBox(TString("glueLayerBackBlankS") += TString("_") + layers->GetName(),
                                                  layerBack_XSize * 0.5, layerBack_YSize * 0.5, glueLayer_ZSize * 0.5);
 
    TGeoShape* glueLayerFrontHoleS = new TGeoTube(TString("glueLayerFrontHoleS") += TString("_") + layers->GetName(),
                                                  0.0, glueLayerHole_Radius, glueLayer_ZSize * 0.5 + 0.01);
    TGeoShape* glueLayerBackHoleS = new TGeoTube(TString("glueLayerBackHoleS") += TString("_") + layers->GetName(), 0.0,
                                                 glueLayerHole_Radius, glueLayer_ZSize * 0.5 + 0.01);
 
    TGeoTranslation* glueLayerFrontHole_pos = new TGeoTranslation();
    glueLayerFrontHole_pos->SetName(TString("glueLayerFrontHole_pos") += TString("_") + layers->GetName());
    glueLayerFrontHole_pos->SetDx(-XModuleSize_Station163x45 * 0.5 - layerFront_XShift);
    glueLayerFrontHole_pos->SetDy(-YModuleSize_Station163x45 * 0.5 - layerFront_YShift);
    glueLayerFrontHole_pos->SetDz(0.0);
    glueLayerFrontHole_pos->RegisterYourself();
 
    TGeoTranslation* glueLayerBackHole_pos = new TGeoTranslation();
    glueLayerBackHole_pos->SetName(TString("glueLayerBackHole_pos") += TString("_") + layers->GetName());
    glueLayerBackHole_pos->SetDx(-XModuleSize_Station163x45 * 0.5 - layerBack_XShift);
    glueLayerBackHole_pos->SetDy(-YModuleSize_Station163x45 * 0.5 - layerBack_YShift);
    glueLayerBackHole_pos->SetDz(0.0);
    glueLayerBackHole_pos->RegisterYourself();
 
    TGeoCompositeShape* glueLayerFrontS = new TGeoCompositeShape();
    glueLayerFrontS->SetName(TString("glueLayerFrontS") += TString("_") + layers->GetName());
    {
        TString expression = "glueLayerFrontBlankS";
        expression += TString("_") + layers->GetName();
        expression += "-glueLayerFrontHoleS";
        expression += TString("_") + layers->GetName();
        expression += ":glueLayerFrontHole_pos";
        expression += TString("_") + layers->GetName();
 
        glueLayerFrontS->MakeNode(expression);
        glueLayerFrontS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoCompositeShape* glueLayerBackS = new TGeoCompositeShape();
    glueLayerBackS->SetName(TString("glueLayerBackS") += TString("_") + layers->GetName());
    {
        TString expression = "glueLayerBackBlankS";
        expression += TString("_") + layers->GetName();
        expression += "-glueLayerBackHoleS";
        expression += TString("_") + layers->GetName();
        expression += ":glueLayerBackHole_pos";
        expression += TString("_") + layers->GetName();
 
        glueLayerBackS->MakeNode(expression);
        glueLayerBackS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* glueLayerFrontV =
        new TGeoVolume(TString("glueLayerFrontV") += TString("_") + layers->GetName(), glueLayerFrontS);
    TGeoVolume* glueLayerBackV =
        new TGeoVolume(TString("glueLayerBackV") += TString("_") + layers->GetName(), glueLayerBackS);
 
    // volume medium
    TGeoMedium* glueLayerV_medium = pMedGlue;
    if (glueLayerV_medium) {
        glueLayerFrontV->SetMedium(glueLayerV_medium);
        glueLayerBackV->SetMedium(glueLayerV_medium);
    } else
        Fatal("Main", "Invalid medium for glueLayerV!");
 
    // volume visual property (transparency)
    glueLayerFrontV->SetLineColor(TColor::GetColor("#1919ff"));
    glueLayerFrontV->SetTransparency(0);
    glueLayerBackV->SetLineColor(TColor::GetColor("#1919ff"));
    glueLayerBackV->SetTransparency(0);
 
    TGeoCombiTrans* glueLayerFront_transf[1];
 
    glueLayerFront_transf[0] = new TGeoCombiTrans();
    glueLayerFront_transf[0]->SetDx(+layerFront_XShift);
    glueLayerFront_transf[0]->SetDy(+layerFront_YShift);
    glueLayerFront_transf[0]->SetDz(-ZGEMSize_Station163x45 * 0.5 + copperLayer_ZSize + glueLayer_ZSize * 0.5);
 
    TGeoCombiTrans* glueLayerBack_transf[1];
 
    glueLayerBack_transf[0] = new TGeoCombiTrans();
    glueLayerBack_transf[0]->SetDx(+layerBack_XShift);
    glueLayerBack_transf[0]->SetDy(+layerBack_YShift);
    glueLayerBack_transf[0]->SetDz(+ZGEMSize_Station163x45 * 0.5 - copperLayer_ZSize - glueLayer_ZSize * 0.5);
 
    layers->AddNode(glueLayerFrontV, 0, glueLayerFront_transf[0]);
    layers->AddNode(glueLayerBackV, 0, glueLayerBack_transf[0]);
    //--------------------------------------------------------------------------
 
    // epoxide layer ------------------------------------------------------------
    Double_t epoxideLayer_ZSize = 0.105;      // cm; half-thickness of all epoxide layers
    Double_t epoxideLayerHole_Radius = 6.0;   // cm
 
    // layer shape
    TGeoShape* epoxideLayerFrontBlankS =
        new TGeoBBox(TString("epoxideLayerFrontBlankS") += TString("_") + layers->GetName(), layerFront_XSize * 0.5,
                     layerFront_YSize * 0.5, epoxideLayer_ZSize * 0.5);
    TGeoShape* epoxideLayerBackBlankS =
        new TGeoBBox(TString("epoxideLayerBackBlankS") += TString("_") + layers->GetName(), layerBack_XSize * 0.5,
                     layerBack_YSize * 0.5, epoxideLayer_ZSize * 0.5);
 
    TGeoShape* epoxideLayerFrontHoleS =
        new TGeoTube(TString("epoxideLayerFrontHoleS") += TString("_") + layers->GetName(), 0.0,
                     epoxideLayerHole_Radius, epoxideLayer_ZSize * 0.5 + 0.01);
    TGeoShape* epoxideLayerBackHoleS =
        new TGeoTube(TString("epoxideLayerBackHoleS") += TString("_") + layers->GetName(), 0.0, epoxideLayerHole_Radius,
                     epoxideLayer_ZSize * 0.5 + 0.01);
 
    TGeoTranslation* epoxideLayerFrontHole_pos = new TGeoTranslation();
    epoxideLayerFrontHole_pos->SetName(TString("epoxideLayerFrontHole_pos") += TString("_") + layers->GetName());
    epoxideLayerFrontHole_pos->SetDx(-XModuleSize_Station163x45 * 0.5 - layerFront_XShift);
    epoxideLayerFrontHole_pos->SetDy(-YModuleSize_Station163x45 * 0.5 - layerFront_YShift);
    epoxideLayerFrontHole_pos->SetDz(0.0);
    epoxideLayerFrontHole_pos->RegisterYourself();
 
    TGeoTranslation* epoxideLayerBackHole_pos = new TGeoTranslation();
    epoxideLayerBackHole_pos->SetName(TString("epoxideLayerBackHole_pos") += TString("_") + layers->GetName());
    epoxideLayerBackHole_pos->SetDx(-XModuleSize_Station163x45 * 0.5 - layerBack_XShift);
    epoxideLayerBackHole_pos->SetDy(-YModuleSize_Station163x45 * 0.5 - layerBack_YShift);
    epoxideLayerBackHole_pos->SetDz(0.0);
    epoxideLayerBackHole_pos->RegisterYourself();
 
    TGeoCompositeShape* epoxideLayerFrontS = new TGeoCompositeShape();
    epoxideLayerFrontS->SetName(TString("epoxideLayerFrontS") += TString("_") + layers->GetName());
    {
        TString expression = "epoxideLayerFrontBlankS";
        expression += TString("_") + layers->GetName();
        expression += "-epoxideLayerFrontHoleS";
        expression += TString("_") + layers->GetName();
        expression += ":epoxideLayerFrontHole_pos";
        expression += TString("_") + layers->GetName();
 
        epoxideLayerFrontS->MakeNode(expression);
        epoxideLayerFrontS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoCompositeShape* epoxideLayerBackS = new TGeoCompositeShape();
    epoxideLayerBackS->SetName(TString("epoxideLayerBackS") += TString("_") + layers->GetName());
    {
        TString expression = "epoxideLayerBackBlankS";
        expression += TString("_") + layers->GetName();
        expression += "-epoxideLayerBackHoleS";
        expression += TString("_") + layers->GetName();
        expression += ":epoxideLayerBackHole_pos";
        expression += TString("_") + layers->GetName();
 
        epoxideLayerBackS->MakeNode(expression);
        epoxideLayerBackS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* epoxideLayerFrontV =
        new TGeoVolume(TString("epoxideLayerFrontV") += TString("_") + layers->GetName(), epoxideLayerFrontS);
    TGeoVolume* epoxideLayerBackV =
        new TGeoVolume(TString("epoxideLayerBackV") += TString("_") + layers->GetName(), epoxideLayerBackS);
 
    // volume medium
    TGeoMedium* epoxideLayerV_medium = pMedEpoxideCompound;
    if (epoxideLayerV_medium) {
        epoxideLayerFrontV->SetMedium(epoxideLayerV_medium);
        epoxideLayerBackV->SetMedium(epoxideLayerV_medium);
    } else
        Fatal("Main", "Invalid medium for epoxideLayerV!");
 
    // volume visual property (transparency)
    epoxideLayerFrontV->SetLineColor(TColor::GetColor("#aa33ff"));
    epoxideLayerFrontV->SetTransparency(0);
    epoxideLayerBackV->SetLineColor(TColor::GetColor("#aa33ff"));
    epoxideLayerBackV->SetTransparency(0);
 
    TGeoCombiTrans* epoxideLayerFront_transf[1];
 
    epoxideLayerFront_transf[0] = new TGeoCombiTrans();
    epoxideLayerFront_transf[0]->SetDx(+layerFront_XShift);
    epoxideLayerFront_transf[0]->SetDy(+layerFront_YShift);
    epoxideLayerFront_transf[0]->SetDz(-ZGEMSize_Station163x45 * 0.5 + copperLayer_ZSize + glueLayer_ZSize
                                       + epoxideLayer_ZSize * 0.5);
 
    TGeoCombiTrans* epoxideLayerBack_transf[1];
 
    epoxideLayerBack_transf[0] = new TGeoCombiTrans();
    epoxideLayerBack_transf[0]->SetDx(+layerBack_XShift);
    epoxideLayerBack_transf[0]->SetDy(+layerBack_YShift);
    epoxideLayerBack_transf[0]->SetDz(+ZGEMSize_Station163x45 * 0.5 - copperLayer_ZSize - glueLayer_ZSize
                                      - epoxideLayer_ZSize * 0.5);
 
    layers->AddNode(epoxideLayerFrontV, 0, epoxideLayerFront_transf[0]);
    layers->AddNode(epoxideLayerBackV, 0, epoxideLayerBack_transf[0]);
    //--------------------------------------------------------------------------
 
    // honeycomb layer ----------------------------------------------------------
    // Double_t honeycombLayer_ZSize = 1.5-(0.05+0.00055+0.005+0.0125); //cm; half-thickness of all honeycomb layers
    Double_t honeycombLayer_ZSize = 1.5 - (0.06805);   // cm; half-thickness of all honeycomb layers
    Double_t honeycombLayerHole_Radius = 6.0;          // cm
 
    // layer shape
    TGeoShape* honeycombLayerFrontBlankS =
        new TGeoBBox(TString("honeycombLayerFrontBlankS") += TString("_") + layers->GetName(), layerFront_XSize * 0.5,
                     layerFront_YSize * 0.5, honeycombLayer_ZSize * 0.5);
    TGeoShape* honeycombLayerBackBlankS =
        new TGeoBBox(TString("honeycombLayerBackBlankS") += TString("_") + layers->GetName(), layerBack_XSize * 0.5,
                     layerBack_YSize * 0.5, honeycombLayer_ZSize * 0.5);
 
    TGeoShape* honeycombLayerFrontHoleS =
        new TGeoTube(TString("honeycombLayerFrontHoleS") += TString("_") + layers->GetName(), 0.0,
                     honeycombLayerHole_Radius, honeycombLayer_ZSize * 0.5 + 0.01);
    TGeoShape* honeycombLayerBackHoleS =
        new TGeoTube(TString("honeycombLayerBackHoleS") += TString("_") + layers->GetName(), 0.0,
                     honeycombLayerHole_Radius, honeycombLayer_ZSize * 0.5 + 0.01);
 
    TGeoTranslation* honeycombLayerFrontHole_pos = new TGeoTranslation();
    honeycombLayerFrontHole_pos->SetName(TString("honeycombLayerFrontHole_pos") += TString("_") + layers->GetName());
    honeycombLayerFrontHole_pos->SetDx(-XModuleSize_Station163x45 * 0.5 - layerFront_XShift);
    honeycombLayerFrontHole_pos->SetDy(-YModuleSize_Station163x45 * 0.5 - layerFront_YShift);
    honeycombLayerFrontHole_pos->SetDz(0.0);
    honeycombLayerFrontHole_pos->RegisterYourself();
 
    TGeoTranslation* honeycombLayerBackHole_pos = new TGeoTranslation();
    honeycombLayerBackHole_pos->SetName(TString("honeycombLayerBackHole_pos") += TString("_") + layers->GetName());
    honeycombLayerBackHole_pos->SetDx(-XModuleSize_Station163x45 * 0.5 - layerBack_XShift);
    honeycombLayerBackHole_pos->SetDy(-YModuleSize_Station163x45 * 0.5 - layerBack_YShift);
    honeycombLayerBackHole_pos->SetDz(0.0);
    honeycombLayerBackHole_pos->RegisterYourself();
 
    TGeoCompositeShape* honeycombLayerFrontS = new TGeoCompositeShape();
    honeycombLayerFrontS->SetName(TString("honeycombLayerFrontS") += TString("_") + layers->GetName());
    {
        TString expression = "honeycombLayerFrontBlankS";
        expression += TString("_") + layers->GetName();
        expression += "-honeycombLayerFrontHoleS";
        expression += TString("_") + layers->GetName();
        expression += ":honeycombLayerFrontHole_pos";
        expression += TString("_") + layers->GetName();
 
        honeycombLayerFrontS->MakeNode(expression);
        honeycombLayerFrontS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoCompositeShape* honeycombLayerBackS = new TGeoCompositeShape();
    honeycombLayerBackS->SetName(TString("honeycombLayerBackS") += TString("_") + layers->GetName());
    {
        TString expression = "honeycombLayerBackBlankS";
        expression += TString("_") + layers->GetName();
        expression += "-honeycombLayerBackHoleS";
        expression += TString("_") + layers->GetName();
        expression += ":honeycombLayerBackHole_pos";
        expression += TString("_") + layers->GetName();
 
        honeycombLayerBackS->MakeNode(expression);
        honeycombLayerBackS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* honeycombLayerFrontV =
        new TGeoVolume(TString("honeycombLayerFrontV") += TString("_") + layers->GetName(), honeycombLayerFrontS);
    TGeoVolume* honeycombLayerBackV =
        new TGeoVolume(TString("honeycombLayerBackV") += TString("_") + layers->GetName(), honeycombLayerBackS);
 
    // volume medium
    TGeoMedium* honeycombLayerV_medium = pMedHoneyComb;
    if (honeycombLayerV_medium) {
        honeycombLayerFrontV->SetMedium(honeycombLayerV_medium);
        honeycombLayerBackV->SetMedium(honeycombLayerV_medium);
    } else
        Fatal("Main", "Invalid medium for honeycombLayerV!");
 
    // volume visual property (transparency)
    honeycombLayerFrontV->SetLineColor(TColor::GetColor("#ffdf84"));
    honeycombLayerFrontV->SetTransparency(0);
    honeycombLayerBackV->SetLineColor(TColor::GetColor("#ffdf84"));
    honeycombLayerBackV->SetTransparency(0);
 
    TGeoCombiTrans* honeycombLayerFront_transf[1];
 
    honeycombLayerFront_transf[0] = new TGeoCombiTrans();
    honeycombLayerFront_transf[0]->SetDx(+layerFront_XShift);
    honeycombLayerFront_transf[0]->SetDy(+layerFront_YShift);
    honeycombLayerFront_transf[0]->SetDz(-ZGEMSize_Station163x45 * 0.5 + copperLayer_ZSize + glueLayer_ZSize
                                         + epoxideLayer_ZSize + honeycombLayer_ZSize * 0.5);
 
    TGeoCombiTrans* honeycombLayerBack_transf[1];
 
    honeycombLayerBack_transf[0] = new TGeoCombiTrans();
    honeycombLayerBack_transf[0]->SetDx(+layerBack_XShift);
    honeycombLayerBack_transf[0]->SetDy(+layerBack_YShift);
    honeycombLayerBack_transf[0]->SetDz(+ZGEMSize_Station163x45 * 0.5 - copperLayer_ZSize - glueLayer_ZSize
                                        - epoxideLayer_ZSize - honeycombLayer_ZSize * 0.5);
 
    layers->AddNode(honeycombLayerFrontV, 0, honeycombLayerFront_transf[0]);
    layers->AddNode(honeycombLayerBackV, 0, honeycombLayerBack_transf[0]);
    //--------------------------------------------------------------------------
 
    // polyamide layer ----------------------------------------------------------
    Double_t polyamideLayer_ZSize = 0.0125;     // cm; half-thickness of all polyamide layers
    Double_t polyamideLayerHole_Radius = 6.0;   // cm
 
    // layer shape
    TGeoShape* polyamideLayerFrontBlankS =
        new TGeoBBox(TString("polyamideLayerFrontBlankS") += TString("_") + layers->GetName(), layerFront_XSize * 0.5,
                     layerFront_YSize * 0.5, polyamideLayer_ZSize * 0.5);
    TGeoShape* polyamideLayerBackBlankS =
        new TGeoBBox(TString("polyamideLayerBackBlankS") += TString("_") + layers->GetName(), layerBack_XSize * 0.5,
                     layerBack_YSize * 0.5, polyamideLayer_ZSize * 0.5);
 
    TGeoShape* polyamideLayerFrontHoleS =
        new TGeoTube(TString("polyamideLayerFrontHoleS") += TString("_") + layers->GetName(), 0.0,
                     polyamideLayerHole_Radius, polyamideLayer_ZSize * 0.5 + 0.01);
    TGeoShape* polyamideLayerBackHoleS =
        new TGeoTube(TString("polyamideLayerBackHoleS") += TString("_") + layers->GetName(), 0.0,
                     polyamideLayerHole_Radius, polyamideLayer_ZSize * 0.5 + 0.01);
 
    TGeoTranslation* polyamideLayerFrontHole_pos = new TGeoTranslation();
    polyamideLayerFrontHole_pos->SetName(TString("polyamideLayerFrontHole_pos") += TString("_") + layers->GetName());
    polyamideLayerFrontHole_pos->SetDx(-XModuleSize_Station163x45 * 0.5 - layerFront_XShift);
    polyamideLayerFrontHole_pos->SetDy(-YModuleSize_Station163x45 * 0.5 - layerFront_YShift);
    polyamideLayerFrontHole_pos->SetDz(0.0);
    polyamideLayerFrontHole_pos->RegisterYourself();
 
    TGeoTranslation* polyamideLayerBackHole_pos = new TGeoTranslation();
    polyamideLayerBackHole_pos->SetName(TString("polyamideLayerBackHole_pos") += TString("_") + layers->GetName());
    polyamideLayerBackHole_pos->SetDx(-XModuleSize_Station163x45 * 0.5 - layerBack_XShift);
    polyamideLayerBackHole_pos->SetDy(-YModuleSize_Station163x45 * 0.5 - layerBack_YShift);
    polyamideLayerBackHole_pos->SetDz(0.0);
    polyamideLayerBackHole_pos->RegisterYourself();
 
    TGeoCompositeShape* polyamideLayerFrontS = new TGeoCompositeShape();
    polyamideLayerFrontS->SetName(TString("polyamideLayerFrontS") += TString("_") + layers->GetName());
    {
        TString expression = "polyamideLayerFrontBlankS";
        expression += TString("_") + layers->GetName();
        expression += "-polyamideLayerFrontHoleS";
        expression += TString("_") + layers->GetName();
        expression += ":polyamideLayerFrontHole_pos";
        expression += TString("_") + layers->GetName();
 
        polyamideLayerFrontS->MakeNode(expression);
        polyamideLayerFrontS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoCompositeShape* polyamideLayerBackS = new TGeoCompositeShape();
    polyamideLayerBackS->SetName(TString("polyamideLayerBackS") += TString("_") + layers->GetName());
    {
        TString expression = "polyamideLayerBackBlankS";
        expression += TString("_") + layers->GetName();
        expression += "-polyamideLayerBackHoleS";
        expression += TString("_") + layers->GetName();
        expression += ":polyamideLayerBackHole_pos";
        expression += TString("_") + layers->GetName();
 
        polyamideLayerBackS->MakeNode(expression);
        polyamideLayerBackS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* polyamideLayerFrontV =
        new TGeoVolume(TString("polyamideLayerFrontV") += TString("_") + layers->GetName(), polyamideLayerFrontS);
    TGeoVolume* polyamideLayerBackV =
        new TGeoVolume(TString("polyamideLayerBackV") += TString("_") + layers->GetName(), polyamideLayerBackS);
 
    // volume medium
    TGeoMedium* polyamideLayerV_medium = pMedPolyamide;
    if (polyamideLayerV_medium) {
        polyamideLayerFrontV->SetMedium(polyamideLayerV_medium);
        polyamideLayerBackV->SetMedium(polyamideLayerV_medium);
    } else
        Fatal("Main", "Invalid medium for polyamideLayerV!");
 
    // volume visual property (transparency)
    polyamideLayerFrontV->SetLineColor(TColor::GetColor("#fefe22"));
    polyamideLayerFrontV->SetTransparency(0);
    polyamideLayerBackV->SetLineColor(TColor::GetColor("#fefe22"));
    polyamideLayerBackV->SetTransparency(0);
 
    TGeoCombiTrans* polyamideLayerFront_transf[1];
 
    polyamideLayerFront_transf[0] = new TGeoCombiTrans();
    polyamideLayerFront_transf[0]->SetDx(+layerFront_XShift);
    polyamideLayerFront_transf[0]->SetDy(+layerFront_YShift);
    polyamideLayerFront_transf[0]->SetDz(-ZGEMSize_Station163x45 * 0.5 + copperLayer_ZSize + glueLayer_ZSize
                                         + epoxideLayer_ZSize + honeycombLayer_ZSize + polyamideLayer_ZSize * 0.5);
 
    TGeoCombiTrans* polyamideLayerBack_transf[1];
 
    polyamideLayerBack_transf[0] = new TGeoCombiTrans();
    polyamideLayerBack_transf[0]->SetDx(+layerBack_XShift);
    polyamideLayerBack_transf[0]->SetDy(+layerBack_YShift);
    polyamideLayerBack_transf[0]->SetDz(+ZGEMSize_Station163x45 * 0.5 - copperLayer_ZSize - glueLayer_ZSize
                                        - epoxideLayer_ZSize - honeycombLayer_ZSize - polyamideLayer_ZSize * 0.5);
 
    layers->AddNode(polyamideLayerFrontV, 0, polyamideLayerFront_transf[0]);
    layers->AddNode(polyamideLayerBackV, 0, polyamideLayerBack_transf[0]);
    //--------------------------------------------------------------------------
 
    return layers;
}
 
//------------------------------------------------------------------------------
TGeoVolume* CreateFrameForModule_Station163x45(TString frame_name,
                                               Double_t dx,
                                               Double_t dy,
                                               Double_t dz,
                                               Double_t hole_radius)
{
 
    // frame shapes
    TGeoShape* vertical_frameS =
        new TGeoBBox("vertical_frameS", dx * 0.5, YModuleSize_Station163x45 * 0.5 + dy, dz * 0.5);
    TGeoShape* horizontal_upper_frameS =
        new TGeoBBox("horizontal_upper_frameS", XModuleSize_Station163x45 * 0.5, dy * 0.5, dz * 0.5);
    TGeoShape* horizontal_bottom_frameS = new TGeoBBox(
        "horizontal_upper_frameS", XModuleSize_Station163x45 * 0.5 - hole_radius * 0.5, dy * 0.5, dz * 0.5);
 
    // frame volumes
    TGeoVolume* composite_frameV = new TGeoVolumeAssembly(frame_name);
    TGeoVolume* vertical_frameV = new TGeoVolume(frame_name + "_vertical_frameV", vertical_frameS);
    TGeoVolume* horizontal_upper_frameV =
        new TGeoVolume(frame_name + "_horizontal_upper_frameV", horizontal_upper_frameS);
    TGeoVolume* horizontal_bottom_frameV =
        new TGeoVolume(frame_name + "_horizontal_bottom_frameV", horizontal_bottom_frameS);
 
    // media
    TGeoMedium* frame_medium = pMedCarbon;   // set medium
    if (pMedCarbon) {
        vertical_frameV->SetMedium(frame_medium);
        horizontal_upper_frameV->SetMedium(frame_medium);
        horizontal_bottom_frameV->SetMedium(frame_medium);
    } else
        Fatal("Main", "Invalid medium for frames!");
    if (pMedAir) {
        composite_frameV->SetMedium(pMedAir);
    } else
        Fatal("Main", "Invalid medium for frames!");
 
    // visual parameters
    vertical_frameV->SetLineColor(TColor::GetColor("#9999ff"));
    horizontal_upper_frameV->SetLineColor(TColor::GetColor("#9999ff"));
    horizontal_bottom_frameV->SetLineColor(TColor::GetColor("#9999ff"));
 
    // frame part positions in the frame
    TGeoCombiTrans* left_vertical_frame_position = new TGeoCombiTrans();
    left_vertical_frame_position->SetTranslation(XModuleSize_Station163x45 * 0.5 + dXFrame_Station163x45 * 0.5, 0.0,
                                                 0.0);
 
    TGeoCombiTrans* upper_horizontal_frame_position = new TGeoCombiTrans();
    upper_horizontal_frame_position->SetTranslation(0.0, YModuleSize_Station163x45 * 0.5 + dYFrame_Station163x45 * 0.5,
                                                    0.0);
 
    TGeoCombiTrans* bottom_horizontal_frame_position = new TGeoCombiTrans();
    bottom_horizontal_frame_position->SetTranslation(
        hole_radius * 0.5, -YModuleSize_Station163x45 * 0.5 - dYFrame_Station163x45 * 0.5, 0.0);
 
    composite_frameV->AddNode(vertical_frameV, 0, left_vertical_frame_position);
    composite_frameV->AddNode(horizontal_upper_frameV, 0, upper_horizontal_frame_position);
    composite_frameV->AddNode(horizontal_bottom_frameV, 0, bottom_horizontal_frame_position);
 
    return composite_frameV;
}
 
//------------------------------------------------------------------------------
TGeoVolume* CreateFrameForHalfPlane_Station163x45(TString frame_name, Int_t istation)
{
 
    // frame volumes
    TGeoVolume* frames = new TGeoVolumeAssembly(frame_name);
    frames->SetMedium(pMedAir);
 
    // back lower epoxide frame -------------------------------------------------
    Double_t backLowerEpoxideFrame_XSize = 161.42;   // cm
    Double_t backLowerEpoxideFrame_YSize = 2.2;      // cm
    Double_t backLowerEpoxideFrame_ZSize = 1.5;      // cm
 
    Double_t backLowerEpoxideFrameArch_Rmin = 4.0;
    Double_t backLowerEpoxideFrameArch_Rmax = 6.0;
    Double_t backLowerEpoxideFrameArch_ZSize = backLowerEpoxideFrame_ZSize;
 
    TGeoShape* backLowerEpoxideFrameHalfS =
        new TGeoBBox(TString("backLowerEpoxideFrameHalfS") += TString("_") + frames->GetName(),
                     (backLowerEpoxideFrame_XSize * 0.5 - backLowerEpoxideFrameArch_Rmin) * 0.5,
                     backLowerEpoxideFrame_YSize * 0.5, backLowerEpoxideFrame_ZSize * 0.5);
    TGeoShape* backLowerEpoxideFrameArchS = new TGeoTubeSeg(
        TString("backLowerEpoxideFrameArchS") += TString("_") + frames->GetName(), backLowerEpoxideFrameArch_Rmin,
        backLowerEpoxideFrameArch_Rmax, backLowerEpoxideFrameArch_ZSize * 0.5, 0.0, 180.0);
 
    TGeoTranslation* backLowerEpoxideFrameHalf_pos1 = new TGeoTranslation();
    backLowerEpoxideFrameHalf_pos1->SetName("backLowerEpoxideFrameHalf1_pos1");
    backLowerEpoxideFrameHalf_pos1->SetDx(
        +((backLowerEpoxideFrame_XSize * 0.5 - backLowerEpoxideFrameArch_Rmin) * 0.5 + backLowerEpoxideFrameArch_Rmin));
    backLowerEpoxideFrameHalf_pos1->SetDy(-(backLowerEpoxideFrame_YSize * 0.5) + 0.2 /*shift*/);
    backLowerEpoxideFrameHalf_pos1->SetDz(0.0);
    backLowerEpoxideFrameHalf_pos1->RegisterYourself();
 
    TGeoTranslation* backLowerEpoxideFrameHalf_pos2 = new TGeoTranslation();
    backLowerEpoxideFrameHalf_pos2->SetName("backLowerEpoxideFrameHalf1_pos2");
    backLowerEpoxideFrameHalf_pos2->SetDx(
        -((backLowerEpoxideFrame_XSize * 0.5 - backLowerEpoxideFrameArch_Rmin) * 0.5 + backLowerEpoxideFrameArch_Rmin));
    backLowerEpoxideFrameHalf_pos2->SetDy(-(backLowerEpoxideFrame_YSize * 0.5) + 0.2 /*shift*/);
    backLowerEpoxideFrameHalf_pos2->SetDz(0.0);
    backLowerEpoxideFrameHalf_pos2->RegisterYourself();
 
    TGeoTranslation* backLowerEpoxideFrameArch_pos = new TGeoTranslation();
    backLowerEpoxideFrameArch_pos->SetName("backLowerEpoxideFrameArch_pos");
    backLowerEpoxideFrameArch_pos->SetDx(0.0);
    backLowerEpoxideFrameArch_pos->SetDy(0.0);
    backLowerEpoxideFrameArch_pos->SetDz(0.0);
    backLowerEpoxideFrameArch_pos->RegisterYourself();
 
    TGeoCompositeShape* backLowerEpoxideFrameS = new TGeoCompositeShape();
    backLowerEpoxideFrameS->SetName(TString("backLowerEpoxideFrameS") += TString("_") + frames->GetName());
    {
        TString expression = "backLowerEpoxideFrameHalfS";
        expression += TString("_") + frames->GetName();
        expression += ":backLowerEpoxideFrameHalf1_pos1";
        expression += "+backLowerEpoxideFrameHalfS";
        expression += TString("_") + frames->GetName();
        expression += ":backLowerEpoxideFrameHalf1_pos2";
        expression += "+backLowerEpoxideFrameArchS";
        expression += TString("_") + frames->GetName();
        expression += ":backLowerEpoxideFrameArch_pos";
 
        backLowerEpoxideFrameS->MakeNode(expression);
        backLowerEpoxideFrameS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* backLowerEpoxideFrameV =
        new TGeoVolume(TString("backLowerEpoxideFrameV") += TString("_") + frames->GetName(), backLowerEpoxideFrameS);
 
    // volume medium
    TGeoMedium* backLowerEpoxideFrameV_medium = pMedEpoxideCompound;
    if (backLowerEpoxideFrameV_medium) {
        backLowerEpoxideFrameV->SetMedium(backLowerEpoxideFrameV_medium);
    } else
        Fatal("Main", "Invalid medium for backLowerEpoxideFrameV!");
 
    // volume visual property (transparency)
    backLowerEpoxideFrameV->SetLineColor(TColor::GetColor("#ccccff"));
    backLowerEpoxideFrameV->SetTransparency(0);
 
    TGeoCombiTrans* backLowerEpoxideFrame_transf[1];
 
    backLowerEpoxideFrame_transf[0] = new TGeoCombiTrans();
    backLowerEpoxideFrame_transf[0]->SetDx(0.0);
    backLowerEpoxideFrame_transf[0]->SetDy(-(YModuleSize_Station163x45 * 0.5));
    backLowerEpoxideFrame_transf[0]->SetDz(+(backLowerEpoxideFrame_ZSize * 0.5 + ZModuleSize_Station163x45 * 0.5));
 
    // alignment corrections
    backLowerEpoxideFrame_transf[0]->SetDx(backLowerEpoxideFrame_transf[0]->GetTranslation()[0] + 0.3);
    backLowerEpoxideFrame_transf[0]->SetDy(backLowerEpoxideFrame_transf[0]->GetTranslation()[1] - 0.5);
    backLowerEpoxideFrame_transf[0]->SetDz(backLowerEpoxideFrame_transf[0]->GetTranslation()[2] - 0.0);
 
    if (istation == 2 || istation == 4) {
        backLowerEpoxideFrame_transf[0]->SetDx(backLowerEpoxideFrame_transf[0]->GetTranslation()[0] + 0.0);
        backLowerEpoxideFrame_transf[0]->SetDy(backLowerEpoxideFrame_transf[0]->GetTranslation()[1] - 0.0);
        backLowerEpoxideFrame_transf[0]->SetDz(backLowerEpoxideFrame_transf[0]->GetTranslation()[2] + 0.12);
    }
 
    frames->AddNode(backLowerEpoxideFrameV, 0, backLowerEpoxideFrame_transf[0]);
    //--------------------------------------------------------------------------
 
    // front lower epoxide frame -------------------------------------------------
    Double_t frontLowerEpoxideFrame_XSize = 161.7;   // cm
    Double_t frontLowerEpoxideFrame_YSize = 2.2;     // cm
    Double_t frontLowerEpoxideFrame_ZSize = 1.6;     // cm
 
    Double_t frontLowerEpoxideFrameArch_Rmin = 4.0;
    Double_t frontLowerEpoxideFrameArch_Rmax = 6.0;
    Double_t frontLowerEpoxideFrameArch_ZSize = frontLowerEpoxideFrame_ZSize;
 
    TGeoShape* frontLowerEpoxideFrameHalfS =
        new TGeoBBox(TString("frontLowerEpoxideFrameHalfS") += TString("_") + frames->GetName(),
                     (frontLowerEpoxideFrame_XSize * 0.5 - frontLowerEpoxideFrameArch_Rmin) * 0.5,
                     frontLowerEpoxideFrame_YSize * 0.5, frontLowerEpoxideFrame_ZSize * 0.5);
    TGeoShape* frontLowerEpoxideFrameArchS = new TGeoTubeSeg(
        TString("frontLowerEpoxideFrameArchS") += TString("_") + frames->GetName(), frontLowerEpoxideFrameArch_Rmin,
        frontLowerEpoxideFrameArch_Rmax, frontLowerEpoxideFrameArch_ZSize * 0.5, 0.0, 180.0);
 
    TGeoTranslation* frontLowerEpoxideFrameHalf_pos1 = new TGeoTranslation();
    frontLowerEpoxideFrameHalf_pos1->SetName("frontLowerEpoxideFrameHalf1_pos1");
    frontLowerEpoxideFrameHalf_pos1->SetDx(
        +((frontLowerEpoxideFrame_XSize * 0.5 - frontLowerEpoxideFrameArch_Rmin) * 0.5
          + frontLowerEpoxideFrameArch_Rmin));
    frontLowerEpoxideFrameHalf_pos1->SetDy(-(frontLowerEpoxideFrame_YSize * 0.5) + 0.2 /*shift*/);
    frontLowerEpoxideFrameHalf_pos1->SetDz(0.0);
    frontLowerEpoxideFrameHalf_pos1->RegisterYourself();
 
    TGeoTranslation* frontLowerEpoxideFrameHalf_pos2 = new TGeoTranslation();
    frontLowerEpoxideFrameHalf_pos2->SetName("frontLowerEpoxideFrameHalf1_pos2");
    frontLowerEpoxideFrameHalf_pos2->SetDx(
        -((frontLowerEpoxideFrame_XSize * 0.5 - frontLowerEpoxideFrameArch_Rmin) * 0.5
          + frontLowerEpoxideFrameArch_Rmin));
    frontLowerEpoxideFrameHalf_pos2->SetDy(-(frontLowerEpoxideFrame_YSize * 0.5) + 0.2 /*shift*/);
    frontLowerEpoxideFrameHalf_pos2->SetDz(0.0);
    frontLowerEpoxideFrameHalf_pos2->RegisterYourself();
 
    TGeoTranslation* frontLowerEpoxideFrameArch_pos = new TGeoTranslation();
    frontLowerEpoxideFrameArch_pos->SetName("frontLowerEpoxideFrameArch_pos");
    frontLowerEpoxideFrameArch_pos->SetDx(0.0);
    frontLowerEpoxideFrameArch_pos->SetDy(0.0);
    frontLowerEpoxideFrameArch_pos->SetDz(0.0);
    frontLowerEpoxideFrameArch_pos->RegisterYourself();
 
    TGeoCompositeShape* frontLowerEpoxideFrameS = new TGeoCompositeShape();
    frontLowerEpoxideFrameS->SetName(TString("frontLowerEpoxideFrameS") += TString("_") + frames->GetName());
    {
        TString expression = "frontLowerEpoxideFrameHalfS";
        expression += TString("_") + frames->GetName();
        expression += ":frontLowerEpoxideFrameHalf1_pos1";
        expression += "+frontLowerEpoxideFrameHalfS";
        expression += TString("_") + frames->GetName();
        expression += ":frontLowerEpoxideFrameHalf1_pos2";
        expression += "+frontLowerEpoxideFrameArchS";
        expression += TString("_") + frames->GetName();
        expression += ":frontLowerEpoxideFrameArch_pos";
 
        frontLowerEpoxideFrameS->MakeNode(expression);
        frontLowerEpoxideFrameS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* frontLowerEpoxideFrameV =
        new TGeoVolume(TString("frontLowerEpoxideFrameV") += TString("_") + frames->GetName(), frontLowerEpoxideFrameS);
 
    // volume medium
    TGeoMedium* frontLowerEpoxideFrameV_medium = pMedEpoxideCompound;
    if (frontLowerEpoxideFrameV_medium) {
        frontLowerEpoxideFrameV->SetMedium(frontLowerEpoxideFrameV_medium);
    } else
        Fatal("Main", "Invalid medium for frontLowerEpoxideFrameV!");
 
    // volume visual property (transparency)
    frontLowerEpoxideFrameV->SetLineColor(TColor::GetColor("#ccccff"));
    frontLowerEpoxideFrameV->SetTransparency(0);
 
    TGeoCombiTrans* frontLowerEpoxideFrame_transf[1];
 
    frontLowerEpoxideFrame_transf[0] = new TGeoCombiTrans();
    frontLowerEpoxideFrame_transf[0]->SetDx(0.0);
    frontLowerEpoxideFrame_transf[0]->SetDy(-(YModuleSize_Station163x45 * 0.5));
    frontLowerEpoxideFrame_transf[0]->SetDz(-(frontLowerEpoxideFrame_ZSize * 0.5 + ZModuleSize_Station163x45 * 0.5));
 
    if (istation == 5) {
        frontLowerEpoxideFrame_transf[0]->SetDx(backLowerEpoxideFrame_transf[0]->GetTranslation()[0] + 0.0);
        frontLowerEpoxideFrame_transf[0]->SetDy(backLowerEpoxideFrame_transf[0]->GetTranslation()[1] - 0.25);
        frontLowerEpoxideFrame_transf[0]->SetDz(backLowerEpoxideFrame_transf[0]->GetTranslation()[2] - 0.0);
    }
 
    if (istation == 1) {
        frontLowerEpoxideFrame_transf[0]->SetDx(0);
        frontLowerEpoxideFrame_transf[0]->SetDy(-(YModuleSize_Station163x39 * 0.5) - 0.0);
        frontLowerEpoxideFrame_transf[0]->SetDz(-(frontLowerEpoxideFrame_ZSize * 0.5 + ZModuleSize_Station163x39 * 0.5)
                                                - 0.007);
    }
 
    if (istation == 0) {
        frontLowerEpoxideFrame_transf[0]->SetDx(0);
        frontLowerEpoxideFrame_transf[0]->SetDy(-(YModuleSize_Station163x39 * 0.5) - 0.0);
        frontLowerEpoxideFrame_transf[0]->SetDz(-(frontLowerEpoxideFrame_ZSize * 0.5 + ZModuleSize_Station163x39 * 0.5)
                                                - 0.008);
    }
    if (istation == 3) {
        frontLowerEpoxideFrame_transf[0]->SetDx(0);
        frontLowerEpoxideFrame_transf[0]->SetDy(-(YModuleSize_Station163x39 * 0.5) - 0.0);
        frontLowerEpoxideFrame_transf[0]->SetDz(-(frontLowerEpoxideFrame_ZSize * 0.5 + ZModuleSize_Station163x39 * 0.5)
                                                - 0.07);
    }
 
    frames->AddNode(frontLowerEpoxideFrameV, 0, frontLowerEpoxideFrame_transf[0]);
    //--------------------------------------------------------------------------
 
    // back side epoxide frame --------------------------------------------------
    Double_t backSideEpoxideFramePart_XSize = 9.2;    // cm
    Double_t backSideEpoxideFramePart_YSize = 52.6;   // cm
    Double_t backSideEpoxideFramePart_ZSize = 1.6;    // cm
 
    Double_t backSideEpoxideFrame_XDistBetweenParts = 162.82;   // cm
 
    TGeoShape* backSideEpoxideFramePartS = new TGeoBBox(
        TString("backSideEpoxideFramePartS") += TString("_") + frames->GetName(), backSideEpoxideFramePart_XSize * 0.5,
        backSideEpoxideFramePart_YSize * 0.5, backSideEpoxideFramePart_ZSize * 0.5);
 
    TGeoTranslation* backSideEpoxideFramePart_pos1 = new TGeoTranslation();
    backSideEpoxideFramePart_pos1->SetName("backSideEpoxideFramePart_pos1");
    backSideEpoxideFramePart_pos1->SetDx(
        +(backSideEpoxideFrame_XDistBetweenParts * 0.5 + backSideEpoxideFramePart_XSize * 0.5));
    backSideEpoxideFramePart_pos1->SetDy(0.0);
    backSideEpoxideFramePart_pos1->SetDz(0.0);
    backSideEpoxideFramePart_pos1->RegisterYourself();
 
    TGeoTranslation* backSideEpoxideFramePart_pos2 = new TGeoTranslation();
    backSideEpoxideFramePart_pos2->SetName("backSideEpoxideFramePart_pos2");
    backSideEpoxideFramePart_pos2->SetDx(
        -(backSideEpoxideFrame_XDistBetweenParts * 0.5 + backSideEpoxideFramePart_XSize * 0.5));
    backSideEpoxideFramePart_pos2->SetDy(0.0);
    backSideEpoxideFramePart_pos2->SetDz(0.0);
    backSideEpoxideFramePart_pos2->RegisterYourself();
 
    TGeoCompositeShape* backSideEpoxideFrameS = new TGeoCompositeShape();
    backSideEpoxideFrameS->SetName(TString("backSideEpoxideFrameS") += TString("_") + frames->GetName());
    {
        TString expression = "backSideEpoxideFramePartS";
        expression += TString("_") + frames->GetName();
        expression += ":backSideEpoxideFramePart_pos1";
        expression += "+backSideEpoxideFramePartS";
        expression += TString("_") + frames->GetName();
        expression += ":backSideEpoxideFramePart_pos2";
 
        backSideEpoxideFrameS->MakeNode(expression);
        backSideEpoxideFrameS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* backSideEpoxideFrameV =
        new TGeoVolume(TString("backSideEpoxideFrameV") += TString("_") + frames->GetName(), backSideEpoxideFrameS);
 
    // volume medium
    TGeoMedium* backSideEpoxideFrameV_medium = pMedEpoxideCompound;
    if (backSideEpoxideFrameV_medium) {
        backSideEpoxideFrameV->SetMedium(backSideEpoxideFrameV_medium);
    } else
        Fatal("Main", "Invalid medium for backSideEpoxideFrameV!");
 
    // volume visual property (transparency)
    backSideEpoxideFrameV->SetLineColor(TColor::GetColor("#ccccff"));
    backSideEpoxideFrameV->SetTransparency(0);
 
    TGeoCombiTrans* backSideEpoxideFrame_transf[1];
 
    backSideEpoxideFrame_transf[0] = new TGeoCombiTrans();
    backSideEpoxideFrame_transf[0]->SetDx(0.0);
    backSideEpoxideFrame_transf[0]->SetDy(+((backSideEpoxideFramePart_YSize - YModuleSize_Station163x45) * 0.5)
                                          - frontLowerEpoxideFrame_YSize + 0.2 /*shift*/);
    backSideEpoxideFrame_transf[0]->SetDz(+(backSideEpoxideFramePart_ZSize * 0.5 + ZModuleSize_Station163x45 * 0.5));
 
    if (istation == 5) {
        backSideEpoxideFrame_transf[0]->SetDx(0.0);
        backSideEpoxideFrame_transf[0]->SetDy(+((backSideEpoxideFramePart_YSize - YModuleSize_Station163x45) * 0.5)
                                              - frontLowerEpoxideFrame_YSize + 0.2 /*shift*/ + 0.0);
        backSideEpoxideFrame_transf[0]->SetDz(+(backSideEpoxideFramePart_ZSize * 0.5 + ZModuleSize_Station163x45 * 0.5)
                                              + 0.1);
    }
 
    if (istation == 2 || istation == 4) {
        backSideEpoxideFrame_transf[0]->SetDx(0.0);
        backSideEpoxideFrame_transf[0]->SetDy(+((backSideEpoxideFramePart_YSize - YModuleSize_Station163x45) * 0.5)
                                              - frontLowerEpoxideFrame_YSize + 0.2 /*shift*/ + 0.0);
        backSideEpoxideFrame_transf[0]->SetDz(+(backSideEpoxideFramePart_ZSize * 0.5 + ZModuleSize_Station163x45 * 0.5)
                                              + 0.06);
    }
 
    frames->AddNode(backSideEpoxideFrameV, 0, backSideEpoxideFrame_transf[0]);
    //--------------------------------------------------------------------------
 
    // back upper epoxide frame -------------------------------------------------
    Double_t backUpperEpoxideFrame_XSize = 161.42;   // cm
    Double_t backUpperEpoxideFrame_YSize = 6.2;      // cm
    Double_t backUpperEpoxideFrame_ZSize = 1.6;      // cm
 
    TGeoShape* backUpperEpoxideFrameS =
        new TGeoBBox("backUpperEpoxideFrameS", backUpperEpoxideFrame_XSize * 0.5, backUpperEpoxideFrame_YSize * 0.5,
                     backUpperEpoxideFrame_ZSize * 0.5);
 
    TGeoVolume* backUpperEpoxideFrameV =
        new TGeoVolume(TString("backUpperEpoxideFrameV") += TString("_") + frames->GetName(), backUpperEpoxideFrameS);
 
    // volume medium
    TGeoMedium* backUpperEpoxideFrameV_medium = pMedEpoxideCompound;
    if (backUpperEpoxideFrameV_medium) {
        backUpperEpoxideFrameV->SetMedium(backUpperEpoxideFrameV_medium);
    } else
        Fatal("Main", "Invalid medium for backUpperEpoxideFrameV!");
 
    // volume visual property (transparency)
    backUpperEpoxideFrameV->SetLineColor(TColor::GetColor("#ccccff"));
    backUpperEpoxideFrameV->SetTransparency(0);
 
    TGeoCombiTrans* backUpperEpoxideFrame_transf[1];
 
    backUpperEpoxideFrame_transf[0] = new TGeoCombiTrans();
    backUpperEpoxideFrame_transf[0]->SetDx(0.0);
    backUpperEpoxideFrame_transf[0]->SetDy(+(backUpperEpoxideFrame_YSize * 0.5 + YModuleSize_Station163x45 * 0.5)
                                           - 0.6 /*shift*/);
    backUpperEpoxideFrame_transf[0]->SetDz(+(backUpperEpoxideFrame_ZSize * 0.5 + ZModuleSize_Station163x45 * 0.5));
 
    if (istation == 5) {
        backUpperEpoxideFrame_transf[0]->SetDx(0.0);
        backUpperEpoxideFrame_transf[0]->SetDy(+(backUpperEpoxideFrame_YSize * 0.5 + YModuleSize_Station163x45 * 0.5)
                                               - 0.6 /*shift*/);
        backUpperEpoxideFrame_transf[0]->SetDz(+(backUpperEpoxideFrame_ZSize * 0.5 + ZModuleSize_Station163x45 * 0.5)
                                               + 0.1);
    }
 
    if (istation == 2 || istation == 4) {
        backUpperEpoxideFrame_transf[0]->SetDx(0.0);
        backUpperEpoxideFrame_transf[0]->SetDy(+(backUpperEpoxideFrame_YSize * 0.5 + YModuleSize_Station163x45 * 0.5)
                                               - 0.6 /*shift*/);
        backUpperEpoxideFrame_transf[0]->SetDz(+(backUpperEpoxideFrame_ZSize * 0.5 + ZModuleSize_Station163x45 * 0.5)
                                               + 0.06);
    }
 
    frames->AddNode(backUpperEpoxideFrameV, 0, backUpperEpoxideFrame_transf[0]);
    //--------------------------------------------------------------------------
 
    // front side epoxide frame -------------------------------------------------
    Double_t frontSideEpoxideFramePart_XSize = 3.5;    // cm
    Double_t frontSideEpoxideFramePart_YSize = 49.9;   // cm
    Double_t frontSideEpoxideFramePart_ZSize = 1.6;    // cm
 
    Double_t frontSideEpoxideFrame_XDistBetweenParts = 162.7;   // cm
 
    TGeoShape* frontSideEpoxideFramePartS =
        new TGeoBBox(TString("frontSideEpoxideFramePartS") += TString("_") + frames->GetName(),
                     frontSideEpoxideFramePart_XSize * 0.5, frontSideEpoxideFramePart_YSize * 0.5,
                     frontSideEpoxideFramePart_ZSize * 0.5);
 
    TGeoTranslation* frontSideEpoxideFramePart_pos1 = new TGeoTranslation();
    frontSideEpoxideFramePart_pos1->SetName("frontSideEpoxideFramePart_pos1");
    frontSideEpoxideFramePart_pos1->SetDx(
        +(frontSideEpoxideFrame_XDistBetweenParts * 0.5 + frontSideEpoxideFramePart_XSize * 0.5));
    frontSideEpoxideFramePart_pos1->SetDy(0.0);
    frontSideEpoxideFramePart_pos1->SetDz(0.0);
    frontSideEpoxideFramePart_pos1->RegisterYourself();
 
    TGeoTranslation* frontSideEpoxideFramePart_pos2 = new TGeoTranslation();
    frontSideEpoxideFramePart_pos2->SetName("frontSideEpoxideFramePart_pos2");
    frontSideEpoxideFramePart_pos2->SetDx(
        -(frontSideEpoxideFrame_XDistBetweenParts * 0.5 + frontSideEpoxideFramePart_XSize * 0.5));
    frontSideEpoxideFramePart_pos2->SetDy(0.0);
    frontSideEpoxideFramePart_pos2->SetDz(0.0);
    frontSideEpoxideFramePart_pos2->RegisterYourself();
 
    TGeoCompositeShape* frontSideEpoxideFrameS = new TGeoCompositeShape();
    frontSideEpoxideFrameS->SetName(TString("frontSideEpoxideFrameS") += TString("_") + frames->GetName());
    {
        TString expression = "frontSideEpoxideFramePartS";
        expression += TString("_") + frames->GetName();
        expression += ":frontSideEpoxideFramePart_pos1";
        expression += "+frontSideEpoxideFramePartS";
        expression += TString("_") + frames->GetName();
        expression += ":frontSideEpoxideFramePart_pos2";
 
        frontSideEpoxideFrameS->MakeNode(expression);
        frontSideEpoxideFrameS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* frontSideEpoxideFrameV =
        new TGeoVolume(TString("frontSideEpoxideFrameV") += TString("_") + frames->GetName(), frontSideEpoxideFrameS);
 
    // volume medium
    TGeoMedium* frontSideEpoxideFrameV_medium = pMedEpoxideCompound;
    if (frontSideEpoxideFrameV_medium) {
        frontSideEpoxideFrameV->SetMedium(frontSideEpoxideFrameV_medium);
    } else
        Fatal("Main", "Invalid medium for frontSideEpoxideFrameV!");
 
    // volume visual property (transparency)
    frontSideEpoxideFrameV->SetLineColor(TColor::GetColor("#ccccff"));
    frontSideEpoxideFrameV->SetTransparency(0);
 
    TGeoCombiTrans* frontSideEpoxideFrame_transf[1];
 
    frontSideEpoxideFrame_transf[0] = new TGeoCombiTrans();
    frontSideEpoxideFrame_transf[0]->SetDx(0.0);
    frontSideEpoxideFrame_transf[0]->SetDy(+((frontSideEpoxideFramePart_YSize - YModuleSize_Station163x45) * 0.5)
                                           - frontLowerEpoxideFrame_YSize + 0.2 /*shift*/);
    frontSideEpoxideFrame_transf[0]->SetDz(-(frontSideEpoxideFramePart_ZSize * 0.5 + ZModuleSize_Station163x45 * 0.5)
                                           - 0.1);
 
    frames->AddNode(frontSideEpoxideFrameV, 0, frontSideEpoxideFrame_transf[0]);
    //--------------------------------------------------------------------------
 
    // front upper epoxide frame -------------------------------------------------
    Double_t frontUpperEpoxideFrame_XSize = 161.7;   // cm
    Double_t frontUpperEpoxideFrame_YSize = 3.5;     // cm
    Double_t frontUpperEpoxideFrame_ZSize = 1.6;     // cm
 
    TGeoShape* frontUpperEpoxideFrameS =
        new TGeoBBox("frontUpperEpoxideFrameS", frontUpperEpoxideFrame_XSize * 0.5, frontUpperEpoxideFrame_YSize * 0.5,
                     frontUpperEpoxideFrame_ZSize * 0.5);
 
    TGeoVolume* frontUpperEpoxideFrameV =
        new TGeoVolume(TString("frontUpperEpoxideFrameV") += TString("_") + frames->GetName(), frontUpperEpoxideFrameS);
 
    // volume medium
    TGeoMedium* frontUpperEpoxideFrameV_medium = pMedEpoxideCompound;
    if (frontUpperEpoxideFrameV_medium) {
        frontUpperEpoxideFrameV->SetMedium(frontUpperEpoxideFrameV_medium);
    } else
        Fatal("Main", "Invalid medium for frontUpperEpoxideFrameV!");
 
    // volume visual property (transparency)
    frontUpperEpoxideFrameV->SetLineColor(TColor::GetColor("#ccccff"));
    frontUpperEpoxideFrameV->SetTransparency(0);
 
    TGeoCombiTrans* frontUpperEpoxideFrame_transf[1];
 
    frontUpperEpoxideFrame_transf[0] = new TGeoCombiTrans();
    frontUpperEpoxideFrame_transf[0]->SetDx(0.0);
    frontUpperEpoxideFrame_transf[0]->SetDy(+(frontUpperEpoxideFrame_YSize * 0.5 + YModuleSize_Station163x45 * 0.5)
                                            - 0.6 /*shift*/);
    frontUpperEpoxideFrame_transf[0]->SetDz(-(frontUpperEpoxideFrame_ZSize * 0.5 + ZModuleSize_Station163x45 * 0.5)
                                            - 0.068);
 
    frames->AddNode(frontUpperEpoxideFrameV, 0, frontUpperEpoxideFrame_transf[0]);
    //--------------------------------------------------------------------------
 
    // alum. cable support ------------------------------------------------------
    Double_t alumCableSupportPart_XSize = 3.0;    // cm
    Double_t alumCableSupportPart_YSize = 41.3;   // cm
    Double_t alumCableSupportPart_ZSize = 0.2;    // cm
 
    Double_t alumCableSupport_XDistBetweenParts = 164.82;   // cm
 
    TGeoShape* alumCableSupportPartS = new TGeoBBox(
        TString("alumCableSupportPartS") += TString("_") + frames->GetName(), alumCableSupportPart_XSize * 0.5,
        alumCableSupportPart_YSize * 0.5, alumCableSupportPart_ZSize * 0.5);
 
    TGeoTranslation* alumCableSupportPart_pos1 = new TGeoTranslation();
    alumCableSupportPart_pos1->SetName("alumCableSupportPart_pos1");
    alumCableSupportPart_pos1->SetDx(+(alumCableSupport_XDistBetweenParts * 0.5 + alumCableSupportPart_XSize * 0.5));
    alumCableSupportPart_pos1->SetDy(0.0);
    alumCableSupportPart_pos1->SetDz(0.0);
    alumCableSupportPart_pos1->RegisterYourself();
 
    TGeoTranslation* alumCableSupportPart_pos2 = new TGeoTranslation();
    alumCableSupportPart_pos2->SetName("alumCableSupportPart_pos2");
    alumCableSupportPart_pos2->SetDx(-(alumCableSupport_XDistBetweenParts * 0.5 + alumCableSupportPart_XSize * 0.5));
    alumCableSupportPart_pos2->SetDy(0.0);
    alumCableSupportPart_pos2->SetDz(0.0);
    alumCableSupportPart_pos2->RegisterYourself();
 
    TGeoCompositeShape* alumCableSupportS = new TGeoCompositeShape();
    alumCableSupportS->SetName(TString("alumCableSupportS") += TString("_") + frames->GetName());
    {
        TString expression = "alumCableSupportPartS";
        expression += TString("_") + frames->GetName();
        expression += ":alumCableSupportPart_pos1";
        expression += "+alumCableSupportPartS";
        expression += TString("_") + frames->GetName();
        expression += ":alumCableSupportPart_pos2";
 
        alumCableSupportS->MakeNode(expression);
        alumCableSupportS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* alumCableSupportV =
        new TGeoVolume(TString("alumCableSupportV") += TString("_") + frames->GetName(), alumCableSupportS);
 
    // volume medium
    TGeoMedium* alumCableSupportV_medium = pMedAluminium;
    if (alumCableSupportV_medium) {
        alumCableSupportV->SetMedium(alumCableSupportV_medium);
    } else
        Fatal("Main", "Invalid medium for alumCableSupportV!");
 
    // volume visual property (transparency)
    alumCableSupportV->SetLineColor(TColor::GetColor("#ffcccc"));
    alumCableSupportV->SetTransparency(0);
 
    TGeoCombiTrans* alumCableSupport_transf[1];
 
    alumCableSupport_transf[0] = new TGeoCombiTrans();
    alumCableSupport_transf[0]->SetDx(0.0);
    alumCableSupport_transf[0]->SetDy(+((alumCableSupportPart_YSize - YModuleSize_Station163x45) * 0.5)
                                      - frontLowerEpoxideFrame_YSize + 0.2 /*shift*/);
    alumCableSupport_transf[0]->SetDz(-(alumCableSupportPart_ZSize * 0.5 + ZModuleSize_Station163x45 * 0.5
                                        + frontSideEpoxideFramePart_ZSize + 1.0 /*shift*/));
 
    frames->AddNode(alumCableSupportV, 0, alumCableSupport_transf[0]);
    //--------------------------------------------------------------------------
 
    // amplifier upper block ----------------------------------------------------
    Double_t amplifierUpperBlockPart1_XSize = 161.44;   // cm
    Double_t amplifierUpperBlockPart1_YSize = 0.6;      // cm
    Double_t amplifierUpperBlockPart1_ZSize = 8.41;     // cm
 
    Double_t amplifierUpperBlockPart2_XSize = 161.16;   // cm
    Double_t amplifierUpperBlockPart2_YSize = 0.6;      // cm
    Double_t amplifierUpperBlockPart2_ZSize = 8.41;     // cm
 
    TGeoShape* amplifierUpperBlockPart1S = new TGeoBBox(
        TString("amplifierUpperBlockPart1S") += TString("_") + frames->GetName(), amplifierUpperBlockPart1_XSize * 0.5,
        amplifierUpperBlockPart1_YSize * 0.5, amplifierUpperBlockPart1_ZSize * 0.5);
    TGeoShape* amplifierUpperBlockPart2S = new TGeoBBox(
        TString("amplifierUpperBlockPart2S") += TString("_") + frames->GetName(), amplifierUpperBlockPart2_XSize * 0.5,
        amplifierUpperBlockPart2_YSize * 0.5, amplifierUpperBlockPart2_ZSize * 0.5);
 
    TGeoTranslation* amplifierUpperBlockPart1_pos = new TGeoTranslation();
    amplifierUpperBlockPart1_pos->SetName("amplifierUpperBlockPart1_pos");
    amplifierUpperBlockPart1_pos->SetDx(0.0);
    amplifierUpperBlockPart1_pos->SetDy(0.0);
    amplifierUpperBlockPart1_pos->SetDz(0.0);
    amplifierUpperBlockPart1_pos->RegisterYourself();
 
    TGeoTranslation* amplifierUpperBlockPart2_pos = new TGeoTranslation();
    amplifierUpperBlockPart2_pos->SetName("amplifierUpperBlockPart2_pos");
    amplifierUpperBlockPart2_pos->SetDx(-0.76);
    amplifierUpperBlockPart2_pos->SetDy(
        -(amplifierUpperBlockPart1_YSize * 0.5 + amplifierUpperBlockPart2_YSize * 0.5 + 0.33));
    amplifierUpperBlockPart2_pos->SetDz(0.0);
    amplifierUpperBlockPart2_pos->RegisterYourself();
 
    TGeoCompositeShape* amplifierUpperBlockS = new TGeoCompositeShape();
    amplifierUpperBlockS->SetName(TString("amplifierUpperBlockS") += TString("_") + frames->GetName());
    {
        TString expression = "amplifierUpperBlockPart1S";
        expression += TString("_") + frames->GetName();
        expression += ":amplifierUpperBlockPart1_pos";
        expression += "+amplifierUpperBlockPart1S";
        expression += TString("_") + frames->GetName();
        expression += ":amplifierUpperBlockPart2_pos";
 
        amplifierUpperBlockS->MakeNode(expression);
        amplifierUpperBlockS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* amplifierUpperBlockV =
        new TGeoVolume(TString("amplifierUpperBlockV") += TString("_") + frames->GetName(), amplifierUpperBlockS);
 
    // volume medium
    TGeoMedium* amplifierUpperBlockV_medium = pMedFiberGlass;
    if (amplifierUpperBlockV_medium) {
        amplifierUpperBlockV->SetMedium(amplifierUpperBlockV_medium);
    } else
        Fatal("Main", "Invalid medium for amplifierUpperBlockV!");
 
    // volume visual property (transparency)
    amplifierUpperBlockV->SetLineColor(TColor::GetColor("#ccffcc"));
    amplifierUpperBlockV->SetTransparency(0);
 
    TGeoCombiTrans* amplifierUpperBlock_transf[1];
 
    amplifierUpperBlock_transf[0] = new TGeoCombiTrans();
    amplifierUpperBlock_transf[0]->SetDx(0.0);
    amplifierUpperBlock_transf[0]->SetDy(+(YModuleSize_Station163x45 * 0.5 + backUpperEpoxideFrame_YSize)
                                         - 1.0 /*shift*/);
    amplifierUpperBlock_transf[0]->SetDz(-(amplifierUpperBlockPart1_ZSize * 0.5) + ZModuleSize_Station163x45 * 0.5);
 
    frames->AddNode(amplifierUpperBlockV, 0, amplifierUpperBlock_transf[0]);
    //--------------------------------------------------------------------------
 
    // amplifier side block -----------------------------------------------------
    Double_t amplifierSideBlockPartBig_XSize = 0.6;    // cm
    Double_t amplifierSideBlockPartBig_YSize = 33.9;   // cm
    Double_t amplifierSideBlockPartBig_ZSize = 8.41;   // cm
 
    Double_t amplifierSideBlockPartSmall_XSize = 0.6;     // cm
    Double_t amplifierSideBlockPartSmall_YSize = 27.05;   // cm
    Double_t amplifierSideBlockPartSmall_ZSize = 8.41;    // cm
 
    TGeoShape* amplifierSideBlockPartBigS =
        new TGeoBBox(TString("amplifierSideBlockPartBigS") += TString("_") + frames->GetName(),
                     amplifierSideBlockPartBig_XSize * 0.5, amplifierSideBlockPartBig_YSize * 0.5,
                     amplifierSideBlockPartBig_ZSize * 0.5);
    TGeoShape* amplifierSideBlockPartSmallS =
        new TGeoBBox(TString("amplifierSideBlockPartSmallS") += TString("_") + frames->GetName(),
                     amplifierSideBlockPartSmall_XSize * 0.5, amplifierSideBlockPartSmall_YSize * 0.5,
                     amplifierSideBlockPartSmall_ZSize * 0.5);
 
    TGeoTranslation* amplifierSideBlockPartBig_pos = new TGeoTranslation();
    amplifierSideBlockPartBig_pos->SetName("amplifierSideBlockPartBig_pos");
    amplifierSideBlockPartBig_pos->SetDx(0.0);
    amplifierSideBlockPartBig_pos->SetDy(0.0);
    amplifierSideBlockPartBig_pos->SetDz(0.0);
    amplifierSideBlockPartBig_pos->RegisterYourself();
 
    TGeoTranslation* amplifierSideBlockPartSmall_pos = new TGeoTranslation();
    amplifierSideBlockPartSmall_pos->SetName("amplifierSideBlockPartSmall_pos");
    amplifierSideBlockPartSmall_pos->SetDx(
        +(amplifierSideBlockPartBig_XSize * 0.5 + amplifierSideBlockPartSmall_XSize * 0.5 + 1.72));
    amplifierSideBlockPartSmall_pos->SetDy(+0.05);
    amplifierSideBlockPartSmall_pos->SetDz(0.0);
    amplifierSideBlockPartSmall_pos->RegisterYourself();
 
    TGeoCompositeShape* amplifierSideBlockS = new TGeoCompositeShape();
    amplifierSideBlockS->SetName(TString("amplifierSideBlockS") += TString("_") + frames->GetName());
    {
        TString expression = "amplifierSideBlockPartBigS";
        expression += TString("_") + frames->GetName();
        expression += ":amplifierSideBlockPartBig_pos";
        expression += "+amplifierSideBlockPartSmallS";
        expression += TString("_") + frames->GetName();
        expression += ":amplifierSideBlockPartSmall_pos";
 
        amplifierSideBlockS->MakeNode(expression);
        amplifierSideBlockS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* amplifierSideBlockV =
        new TGeoVolume(TString("amplifierSideBlockV") += TString("_") + frames->GetName(), amplifierSideBlockS);
 
    // volume medium
    TGeoMedium* amplifierSideBlockV_medium = pMedFiberGlass;
    if (amplifierSideBlockV_medium) {
        amplifierSideBlockV->SetMedium(amplifierSideBlockV_medium);
    } else
        Fatal("Main", "Invalid medium for amplifierSideBlockV!");
 
    // volume visual property (transparency)
    amplifierSideBlockV->SetLineColor(TColor::GetColor("#ccffcc"));
    amplifierSideBlockV->SetTransparency(0);
 
    TGeoCombiTrans* amplifierSideBlock_transf[2];
 
    amplifierSideBlock_transf[0] = new TGeoCombiTrans();
    amplifierSideBlock_transf[0]->SetDx(+86.47 /*shift*/);
    amplifierSideBlock_transf[0]->SetDy(-5.5);
    amplifierSideBlock_transf[0]->SetDz(-(amplifierSideBlockPartBig_ZSize * 0.5) + ZModuleSize_Station163x45 * 0.5
                                        - 0.4 /*shift*/);
 
    amplifierSideBlock_transf[1] = new TGeoCombiTrans();
    amplifierSideBlock_transf[1]->RotateY(180.0);
    amplifierSideBlock_transf[1]->SetDx(-86.47 /*shift*/);
    amplifierSideBlock_transf[1]->SetDy(-5.5);
    amplifierSideBlock_transf[1]->SetDz(-(amplifierSideBlockPartBig_ZSize * 0.5) + ZModuleSize_Station163x45 * 0.5
                                        - 0.4 /*shift*/);
 
    frames->AddNode(amplifierSideBlockV, 0, amplifierSideBlock_transf[0]);
    frames->AddNode(amplifierSideBlockV, 1, amplifierSideBlock_transf[1]);
    //--------------------------------------------------------------------------
 
    // upper PCB (Printed Circuit Board) ----------------------------------------
    Double_t upperPcb_XSize = 177.22;   // cm
    Double_t upperPcb_YSize = 0.3;      // cm
    Double_t upperPcb_ZSize = 4.0;      // cm
 
    TGeoShape* upperPcbS = new TGeoBBox("upperPcbS", upperPcb_XSize * 0.5, upperPcb_YSize * 0.5, upperPcb_ZSize * 0.5);
 
    TGeoVolume* upperPcbV = new TGeoVolume(TString("upperPcbV") += TString("_") + frames->GetName(), upperPcbS);
 
    // volume medium
    TGeoMedium* upperPcbV_medium = pMedFiberGlass;
    if (upperPcbV_medium) {
        upperPcbV->SetMedium(upperPcbV_medium);
    } else
        Fatal("Main", "Invalid medium for upperPcbV!");
 
    // volume visual property (transparency)
    upperPcbV->SetLineColor(TColor::GetColor("#88ff88"));
    upperPcbV->SetTransparency(0);
 
    TGeoCombiTrans* upperPcb_transf[1];
 
    upperPcb_transf[0] = new TGeoCombiTrans();
    upperPcb_transf[0]->SetDx(+0.1 /*shift*/);
    upperPcb_transf[0]->SetDy(+(upperPcb_YSize * 0.5 + 23.5 /*shift*/));
    upperPcb_transf[0]->SetDz(-(upperPcb_ZSize * 0.5) + ZModuleSize_Station163x45 * 0.5 - ZModuleSize_Station163x45
                              - frontSideEpoxideFramePart_ZSize - 0.1);
 
    frames->AddNode(upperPcbV, 0, upperPcb_transf[0]);
    //--------------------------------------------------------------------------
 
    // side PCB -----------------------------------------------------------------
    Double_t sidePcbPart_XSize = 0.3;     // cm
    Double_t sidePcbPart_YSize = 35.05;   // cm
    Double_t sidePcbPart_ZSize = 4.0;     // cm
 
    TGeoShape* sidePcbPartS = new TGeoBBox(TString("sidePcbPartS") += TString("_") + frames->GetName(),
                                           sidePcbPart_XSize * 0.5, sidePcbPart_YSize * 0.5, sidePcbPart_ZSize * 0.5);
 
    TGeoTranslation* sidePcbPart_pos1 = new TGeoTranslation();
    sidePcbPart_pos1->SetName("sidePcbPart_pos1");
    sidePcbPart_pos1->SetDx(0.0);
    sidePcbPart_pos1->SetDy(0.0);
    sidePcbPart_pos1->SetDz(0.0);
    sidePcbPart_pos1->RegisterYourself();
 
    TGeoTranslation* sidePcbPart_pos2 = new TGeoTranslation();
    sidePcbPart_pos2->SetName("sidePcbPart_pos2");
    sidePcbPart_pos2->SetDx(+(sidePcbPart_XSize + 6.1 /*shift*/));
    sidePcbPart_pos2->SetDy(0.0);
    sidePcbPart_pos2->SetDz(0.0);
    sidePcbPart_pos2->RegisterYourself();
 
    TGeoCompositeShape* sidePcbS = new TGeoCompositeShape();
    sidePcbS->SetName(TString("sidePcbS") += TString("_") + frames->GetName());
    {
        TString expression = "sidePcbPartS";
        expression += TString("_") + frames->GetName();
        expression += ":sidePcbPart_pos1";
        expression += "+sidePcbPartS";
        expression += TString("_") + frames->GetName();
        expression += ":sidePcbPart_pos2";
 
        sidePcbS->MakeNode(expression);
        sidePcbS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* sidePcbV = new TGeoVolume(TString("sidePcbV") += TString("_") + frames->GetName(), sidePcbS);
 
    // volume medium
    TGeoMedium* sidePcbV_medium = pMedFiberGlass;
    if (sidePcbV_medium) {
        sidePcbV->SetMedium(sidePcbV_medium);
    } else
        Fatal("Main", "Invalid medium for sidePcbV!");
 
    // volume visual property (transparency)
    sidePcbV->SetLineColor(TColor::GetColor("#88ff88"));
    sidePcbV->SetTransparency(0);
 
    TGeoCombiTrans* sidePcb_transf[2];
 
    sidePcb_transf[0] = new TGeoCombiTrans();
    sidePcb_transf[0]->SetDx(+(sidePcbPart_XSize * 0.5 + 84.1 /*shift*/));
    sidePcb_transf[0]->SetDy(-6.8 /*shift*/);
    sidePcb_transf[0]->SetDz(-(sidePcbPart_ZSize * 0.5 + alumCableSupportPart_ZSize + ZModuleSize_Station163x45 * 0.5
                               + frontSideEpoxideFramePart_ZSize + 1.0 /*shift*/));
 
    sidePcb_transf[1] = new TGeoCombiTrans();
    sidePcb_transf[1]->RotateY(180.0);
    sidePcb_transf[1]->SetDx(-(sidePcbPart_XSize * 0.5 + 84.1 /*shift*/));
    sidePcb_transf[1]->SetDy(-6.8 /*shift*/);
    sidePcb_transf[1]->SetDz(-(sidePcbPart_ZSize * 0.5 + alumCableSupportPart_ZSize + ZModuleSize_Station163x45 * 0.5
                               + frontSideEpoxideFramePart_ZSize + 1.0 /*shift*/));
 
    frames->AddNode(sidePcbV, 0, sidePcb_transf[0]);
    frames->AddNode(sidePcbV, 1, sidePcb_transf[1]);
    //--------------------------------------------------------------------------
 
    // alum. side support -------------------------------------------------------
    Double_t alumSideSupport_XSize = 4.8;    // cm
    Double_t alumSideSupport_YSize = 52.6;   // cm
    Double_t alumSideSupport_ZSize = 0.4;    // cm
 
    TGeoShape* alumSideSupportS = new TGeoBBox("alumSideSupportS", alumSideSupport_XSize * 0.5,
                                               alumSideSupport_YSize * 0.5, alumSideSupport_ZSize * 0.5);
 
    TGeoVolume* alumSideSupportV =
        new TGeoVolume(TString("alumSideSupportV") += TString("_") + frames->GetName(), alumSideSupportS);
 
    // volume medium
    TGeoMedium* alumSideSupportV_medium = pMedAluminium;
    if (alumSideSupportV_medium) {
        alumSideSupportV->SetMedium(alumSideSupportV_medium);
    } else
        Fatal("Main", "Invalid medium for alumSideSupportV!");
 
    // volume visual property (transparency)
    alumSideSupportV->SetLineColor(TColor::GetColor("#ffcccc"));
    alumSideSupportV->SetTransparency(0);
 
    TGeoCombiTrans* alumSideSupport_transf[2];
 
    alumSideSupport_transf[0] = new TGeoCombiTrans();
    alumSideSupport_transf[0]->SetDx(+(alumSideSupport_XSize * 0.5 + 88.31 /*shift*/));
    alumSideSupport_transf[0]->SetDy(+1.8);
    alumSideSupport_transf[0]->SetDz(+(ZModuleSize_Station163x45 * 0.5 - alumSideSupport_ZSize * 0.5));
 
    alumSideSupport_transf[1] = new TGeoCombiTrans();
    alumSideSupport_transf[1]->SetDx(-(alumSideSupport_XSize * 0.5 + 88.31 /*shift*/));
    alumSideSupport_transf[1]->SetDy(+1.8 /*shift*/);
    alumSideSupport_transf[1]->SetDz(+(ZModuleSize_Station163x45 * 0.5 - alumSideSupport_ZSize * 0.5));
 
    frames->AddNode(alumSideSupportV, 0, alumSideSupport_transf[0]);
    frames->AddNode(alumSideSupportV, 1, alumSideSupport_transf[1]);
    //--------------------------------------------------------------------------
 
    // bosch shape --------------------------------------------------------------
    Double_t boschShape_XSize = 2.25;   // cm
    Double_t boschShape_YSize = 52.6;   // cm
    Double_t boschShape_ZSize = 2.25;   // cm
 
    TGeoShape* boschShapeS =
        new TGeoBBox("boschShapeS", boschShape_XSize * 0.5, boschShape_YSize * 0.5, boschShape_ZSize * 0.5);
 
    TGeoVolume* boschShapeV = new TGeoVolume(TString("boschShapeV") += TString("_") + frames->GetName(), boschShapeS);
 
    // volume medium
    TGeoMedium* boschShapeV_medium = pMedAluminium;
    if (boschShapeV_medium) {
        boschShapeV->SetMedium(boschShapeV_medium);
    } else
        Fatal("Main", "Invalid medium for boschShapeV!");
 
    // volume visual property (transparency)
    boschShapeV->SetLineColor(TColor::GetColor("#ffcccc"));
    boschShapeV->SetTransparency(0);
 
    TGeoCombiTrans* boschShape_transf[2];
 
    boschShape_transf[0] = new TGeoCombiTrans();
    boschShape_transf[0]->SetDx(+(boschShape_XSize * 0.5 + 90.86 /*shift*/));
    boschShape_transf[0]->SetDy(+1.8 /*shift*/);
    boschShape_transf[0]->SetDz(+(ZModuleSize_Station163x45 * 0.5 + boschShape_ZSize * 0.5));
 
    boschShape_transf[1] = new TGeoCombiTrans();
    boschShape_transf[1]->SetDx(-(boschShape_XSize * 0.5 + 90.86 /*shift*/));
    boschShape_transf[1]->SetDy(+1.8 /*shift*/);
    boschShape_transf[1]->SetDz(+(ZModuleSize_Station163x45 * 0.5 + boschShape_ZSize * 0.5));
 
    // frames->AddNode(boschShapeV, 0, boschShape_transf[0]);
    // frames->AddNode(boschShapeV, 1, boschShape_transf[1]);
    //--------------------------------------------------------------------------
 
    // polyacetal shape ---------------------------------------------------------
    Double_t insidePolyacetalShapeUpper_XSize = 11.46;                // cm
    Double_t insidePolyacetalShapeUpper_YSize = 52.6;                 // cm
    Double_t insidePolyacetalShapeUpper_ZSize = 3.0;                  // cm
    Double_t insidePolyacetalShapeUpper_XDistBetweenParts = 186.22;   // cm
 
    Double_t outsidePolyacetalShapeUpper_XSize = 23.5;                 // cm
    Double_t outsidePolyacetalShapeUpper_YSize = 5.405;                // cm
    Double_t outsidePolyacetalShapeUpper_ZSize = 3.0;                  // cm
    Double_t outsidePolyacetalShapeUpper_XDistBetweenParts = 209.14;   // cm
 
    Double_t insidePolyacetalShapeLower_XSize = 8.4;                  // cm
    Double_t insidePolyacetalShapeLower_YSize = 52.6;                 // cm
    Double_t insidePolyacetalShapeLower_ZSize = 3.0;                  // cm
    Double_t insidePolyacetalShapeLower_XDistBetweenParts = 186.22;   // cm
 
    Double_t outsidePolyacetalShapeLower_XSize = 26.7;                 // cm
    Double_t outsidePolyacetalShapeLower_YSize = 24.4;                 // cm
    Double_t outsidePolyacetalShapeLower_ZSize = 3.0;                  // cm
    Double_t outsidePolyacetalShapeLower_XDistBetweenParts = 203.02;   // cm
 
    TGeoShape* insidePolyacetalShapeUpperS =
        new TGeoBBox("insidePolyacetalShapeUpperS", insidePolyacetalShapeUpper_XSize * 0.5,
                     insidePolyacetalShapeUpper_YSize * 0.5, insidePolyacetalShapeUpper_ZSize * 0.5);
    TGeoShape* outsidePolyacetalShapeUpperS =
        new TGeoBBox("outsidePolyacetalShapeUpperS", outsidePolyacetalShapeUpper_XSize * 0.5,
                     outsidePolyacetalShapeUpper_YSize * 0.5, outsidePolyacetalShapeUpper_ZSize * 0.5);
    TGeoShape* insidePolyacetalShapeLowerS =
        new TGeoBBox("insidePolyacetalShapeLowerS", insidePolyacetalShapeLower_XSize * 0.5,
                     insidePolyacetalShapeLower_YSize * 0.5, insidePolyacetalShapeLower_ZSize * 0.5);
    TGeoShape* outsidePolyacetalShapeLowerS =
        new TGeoBBox("outsidePolyacetalShapeLowerS", outsidePolyacetalShapeLower_XSize * 0.5,
                     outsidePolyacetalShapeLower_YSize * 0.5, outsidePolyacetalShapeLower_ZSize * 0.5);
 
    TGeoVolume* insidePolyacetalShapeUpperV = new TGeoVolume(
        TString("insidePolyacetalShapeUpperV") += TString("_") + frames->GetName(), insidePolyacetalShapeUpperS);
    TGeoVolume* outsidePolyacetalShapeUpperV = new TGeoVolume(
        TString("outsidePolyacetalShapeUpperV") += TString("_") + frames->GetName(), outsidePolyacetalShapeUpperS);
    TGeoVolume* insidePolyacetalShapeLowerV = new TGeoVolume(
        TString("insidePolyacetalShapeLowerV") += TString("_") + frames->GetName(), insidePolyacetalShapeLowerS);
    TGeoVolume* outsidePolyacetalShapeLowerV = new TGeoVolume(
        TString("outsidePolyacetalShapeLowerV") += TString("_") + frames->GetName(), outsidePolyacetalShapeLowerS);
 
    // volume medium
    TGeoMedium* polyacetalShapeV_medium = pMedPolyacetal;
    if (polyacetalShapeV_medium) {
        insidePolyacetalShapeUpperV->SetMedium(polyacetalShapeV_medium);
        outsidePolyacetalShapeUpperV->SetMedium(polyacetalShapeV_medium);
        insidePolyacetalShapeLowerV->SetMedium(polyacetalShapeV_medium);
        outsidePolyacetalShapeLowerV->SetMedium(polyacetalShapeV_medium);
    } else
        Fatal("Main", "Invalid medium for insidePolyacetalShapeUpperV!");
 
    // volume visual property (transparency)
    insidePolyacetalShapeUpperV->SetLineColor(TColor::GetColor("#ffcccc"));
    insidePolyacetalShapeUpperV->SetTransparency(0);
    outsidePolyacetalShapeUpperV->SetLineColor(TColor::GetColor("#ffcccc"));
    outsidePolyacetalShapeUpperV->SetTransparency(0);
    insidePolyacetalShapeLowerV->SetLineColor(TColor::GetColor("#ffcccc"));
    insidePolyacetalShapeLowerV->SetTransparency(0);
    outsidePolyacetalShapeLowerV->SetLineColor(TColor::GetColor("#ffcccc"));
    outsidePolyacetalShapeLowerV->SetTransparency(0);
 
    TGeoCombiTrans* insidePolyacetalShapeUpper_transf[2];
 
    insidePolyacetalShapeUpper_transf[0] = new TGeoCombiTrans();
    insidePolyacetalShapeUpper_transf[0]->SetDx(
        +(insidePolyacetalShapeUpper_XSize * 0.5 + insidePolyacetalShapeUpper_XDistBetweenParts * 0.5));
    insidePolyacetalShapeUpper_transf[0]->SetDy(+1.8 /*shift*/);
    insidePolyacetalShapeUpper_transf[0]->SetDz(0.0);
 
    insidePolyacetalShapeUpper_transf[1] = new TGeoCombiTrans();
    insidePolyacetalShapeUpper_transf[1]->SetDx(
        -(insidePolyacetalShapeUpper_XSize * 0.5 + insidePolyacetalShapeUpper_XDistBetweenParts * 0.5));
    insidePolyacetalShapeUpper_transf[1]->SetDy(+1.8 /*shift*/);
    insidePolyacetalShapeUpper_transf[1]->SetDz(0.0);
 
    TGeoCombiTrans* outsidePolyacetalShapeUpper_transf[2];
 
    outsidePolyacetalShapeUpper_transf[0] = new TGeoCombiTrans();
    outsidePolyacetalShapeUpper_transf[0]->SetDx(
        +(outsidePolyacetalShapeUpper_XSize * 0.5 + outsidePolyacetalShapeUpper_XDistBetweenParts * 0.5));
    outsidePolyacetalShapeUpper_transf[0]->SetDy(-insidePolyacetalShapeUpper_YSize * 0.5
                                                 + 1.8 /*shift*/ + outsidePolyacetalShapeUpper_YSize * 0.5);
    outsidePolyacetalShapeUpper_transf[0]->SetDz(0.0);
 
    outsidePolyacetalShapeUpper_transf[1] = new TGeoCombiTrans();
    outsidePolyacetalShapeUpper_transf[1]->SetDx(
        -(outsidePolyacetalShapeUpper_XSize * 0.5 + outsidePolyacetalShapeUpper_XDistBetweenParts * 0.5));
    outsidePolyacetalShapeUpper_transf[1]->SetDy(-insidePolyacetalShapeUpper_YSize * 0.5
                                                 + 1.8 /*shift*/ + outsidePolyacetalShapeUpper_YSize * 0.5);
    outsidePolyacetalShapeUpper_transf[1]->SetDz(0.0);
 
    TGeoCombiTrans* insidePolyacetalShapeLower_transf[2];
 
    insidePolyacetalShapeLower_transf[0] = new TGeoCombiTrans();
    insidePolyacetalShapeLower_transf[0]->SetDx(
        +(insidePolyacetalShapeLower_XSize * 0.5 + insidePolyacetalShapeLower_XDistBetweenParts * 0.5));
    insidePolyacetalShapeLower_transf[0]->SetDy(+1.8 /*shift*/);
    insidePolyacetalShapeLower_transf[0]->SetDz(0.0);
 
    insidePolyacetalShapeLower_transf[1] = new TGeoCombiTrans();
    insidePolyacetalShapeLower_transf[1]->SetDx(
        -(insidePolyacetalShapeLower_XSize * 0.5 + insidePolyacetalShapeLower_XDistBetweenParts * 0.5));
    insidePolyacetalShapeLower_transf[1]->SetDy(+1.8 /*shift*/);
    insidePolyacetalShapeLower_transf[1]->SetDz(0.0);
 
    TGeoCombiTrans* outsidePolyacetalShapeLower_transf[2];
 
    outsidePolyacetalShapeLower_transf[0] = new TGeoCombiTrans();
    outsidePolyacetalShapeLower_transf[0]->SetDx(
        +(outsidePolyacetalShapeLower_XSize * 0.5 + outsidePolyacetalShapeLower_XDistBetweenParts * 0.5));
    outsidePolyacetalShapeLower_transf[0]->SetDy(-insidePolyacetalShapeLower_YSize * 0.5
                                                 + 1.8 /*shift*/ + outsidePolyacetalShapeLower_YSize * 0.5);
    outsidePolyacetalShapeLower_transf[0]->SetDz(0.0);
 
    outsidePolyacetalShapeLower_transf[1] = new TGeoCombiTrans();
    outsidePolyacetalShapeLower_transf[1]->SetDx(
        -(outsidePolyacetalShapeLower_XSize * 0.5 + outsidePolyacetalShapeLower_XDistBetweenParts * 0.5));
    outsidePolyacetalShapeLower_transf[1]->SetDy(-insidePolyacetalShapeLower_YSize * 0.5
                                                 + 1.8 /*shift*/ + outsidePolyacetalShapeLower_YSize * 0.5);
    outsidePolyacetalShapeLower_transf[1]->SetDz(0.0);
 
    if (TString(frames->GetName()).Contains("upper")) {
        frames->AddNode(insidePolyacetalShapeUpperV, 0, insidePolyacetalShapeUpper_transf[0]);
        frames->AddNode(insidePolyacetalShapeUpperV, 1, insidePolyacetalShapeUpper_transf[1]);
        frames->AddNode(outsidePolyacetalShapeUpperV, 0, outsidePolyacetalShapeUpper_transf[0]);
        frames->AddNode(outsidePolyacetalShapeUpperV, 1, outsidePolyacetalShapeUpper_transf[1]);
    }
 
    if (TString(frames->GetName()).Contains("lower")) {
        frames->AddNode(insidePolyacetalShapeLowerV, 0, insidePolyacetalShapeLower_transf[0]);
        frames->AddNode(insidePolyacetalShapeLowerV, 1, insidePolyacetalShapeLower_transf[1]);
        frames->AddNode(outsidePolyacetalShapeLowerV, 0, outsidePolyacetalShapeLower_transf[0]);
        frames->AddNode(outsidePolyacetalShapeLowerV, 1, outsidePolyacetalShapeLower_transf[1]);
    }
 
    //--------------------------------------------------------------------------
 
    // lower fiberglass block ---------------------------------------------------
    Double_t lowerFiberGlassBlock_XSize = 165.1;   // cm
    Double_t lowerFiberGlassBlock_YSize = 1.75;    // cm
    Double_t lowerFiberGlassBlock_ZSize = 0.9;     // cm
 
    Double_t lowerFiberGlassBlockArch_Rmin = 4.0;
    Double_t lowerFiberGlassBlockArch_Rmax = 5.75;
    Double_t lowerFiberGlassBlockArch_ZSize = lowerFiberGlassBlock_ZSize;
 
    TGeoShape* lowerFiberGlassBlockHalfS =
        new TGeoBBox(TString("lowerFiberGlassBlockHalfS") += TString("_") + frames->GetName(),
                     (lowerFiberGlassBlock_XSize * 0.5 - lowerFiberGlassBlockArch_Rmin) * 0.5,
                     lowerFiberGlassBlock_YSize * 0.5, lowerFiberGlassBlock_ZSize * 0.5);
    TGeoShape* lowerFiberGlassBlockArchS = new TGeoTubeSeg(
        TString("lowerFiberGlassBlockArchS") += TString("_") + frames->GetName(), lowerFiberGlassBlockArch_Rmin,
        lowerFiberGlassBlockArch_Rmax, lowerFiberGlassBlockArch_ZSize * 0.5, 0.0, 180.0);
 
    TGeoTranslation* lowerFiberGlassBlockHalf_pos1 = new TGeoTranslation();
    lowerFiberGlassBlockHalf_pos1->SetName("lowerFiberGlassBlockHalf_pos1");
    lowerFiberGlassBlockHalf_pos1->SetDx(
        +((lowerFiberGlassBlock_XSize * 0.5 - lowerFiberGlassBlockArch_Rmin) * 0.5 + lowerFiberGlassBlockArch_Rmin));
    lowerFiberGlassBlockHalf_pos1->SetDy(-(lowerFiberGlassBlock_YSize * 0.5));
    lowerFiberGlassBlockHalf_pos1->SetDz(0.0);
    lowerFiberGlassBlockHalf_pos1->RegisterYourself();
 
    TGeoTranslation* lowerFiberGlassBlockHalf_pos2 = new TGeoTranslation();
    lowerFiberGlassBlockHalf_pos2->SetName("lowerFiberGlassBlockHalf_pos2");
    lowerFiberGlassBlockHalf_pos2->SetDx(
        -((lowerFiberGlassBlock_XSize * 0.5 - lowerFiberGlassBlockArch_Rmin) * 0.5 + lowerFiberGlassBlockArch_Rmin));
    lowerFiberGlassBlockHalf_pos2->SetDy(-(lowerFiberGlassBlock_YSize * 0.5));
    lowerFiberGlassBlockHalf_pos2->SetDz(0.0);
    lowerFiberGlassBlockHalf_pos2->RegisterYourself();
 
    TGeoTranslation* lowerFiberGlassBlockArch_pos = new TGeoTranslation();
    lowerFiberGlassBlockArch_pos->SetName("lowerFiberGlassBlockArch_pos");
    lowerFiberGlassBlockArch_pos->SetDx(0.0);
    lowerFiberGlassBlockArch_pos->SetDy(0.0);
    lowerFiberGlassBlockArch_pos->SetDz(0.0);
    lowerFiberGlassBlockArch_pos->RegisterYourself();
 
    TGeoCompositeShape* lowerFiberGlassBlockS = new TGeoCompositeShape();
    lowerFiberGlassBlockS->SetName(TString("lowerFiberGlassBlockS") += TString("_") + frames->GetName());
    {
        TString expression = "lowerFiberGlassBlockHalfS";
        expression += TString("_") + frames->GetName();
        expression += ":lowerFiberGlassBlockHalf_pos1";
        expression += "+lowerFiberGlassBlockHalfS";
        expression += TString("_") + frames->GetName();
        expression += ":lowerFiberGlassBlockHalf_pos2";
        expression += "+lowerFiberGlassBlockArchS";
        expression += TString("_") + frames->GetName();
        expression += ":lowerFiberGlassBlockArch_pos";
 
        lowerFiberGlassBlockS->MakeNode(expression);
        lowerFiberGlassBlockS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* lowerFiberGlassBlockV =
        new TGeoVolume(TString("lowerFiberGlassBlockV") += TString("_") + frames->GetName(), lowerFiberGlassBlockS);
 
    // volume medium
    TGeoMedium* lowerFiberGlassBlockV_medium = pMedFiberGlass;
    if (lowerFiberGlassBlockV_medium) {
        lowerFiberGlassBlockV->SetMedium(lowerFiberGlassBlockV_medium);
    } else
        Fatal("Main", "Invalid medium for lowerFiberGlassBlockV!");
 
    // volume visual property (transparency)
    lowerFiberGlassBlockV->SetLineColor(TColor::GetColor("#88ff88"));
    lowerFiberGlassBlockV->SetTransparency(0);
 
    TGeoCombiTrans* lowerFiberGlassBlock_transf[1];
 
    lowerFiberGlassBlock_transf[0] = new TGeoCombiTrans();
    lowerFiberGlassBlock_transf[0]->SetDx(+0.0);
    lowerFiberGlassBlock_transf[0]->SetDy(-(YModuleSize_Station163x45 * 0.5) - 0.0);
    lowerFiberGlassBlock_transf[0]->SetDz(0.0);
 
    if (istation == 3) {
        lowerFiberGlassBlock_transf[0]->SetDx(lowerFiberGlassBlock_transf[0]->GetTranslation()[0] + 0.1);
        lowerFiberGlassBlock_transf[0]->SetDy(lowerFiberGlassBlock_transf[0]->GetTranslation()[1] - 0.3);
        lowerFiberGlassBlock_transf[0]->SetDz(lowerFiberGlassBlock_transf[0]->GetTranslation()[2] - 0.1);
    }
    if (istation == 1) {
        lowerFiberGlassBlock_transf[0]->SetDx(lowerFiberGlassBlock_transf[0]->GetTranslation()[0] + 0.0);
        lowerFiberGlassBlock_transf[0]->SetDy(lowerFiberGlassBlock_transf[0]->GetTranslation()[1] - 0.0);
        lowerFiberGlassBlock_transf[0]->SetDz(lowerFiberGlassBlock_transf[0]->GetTranslation()[2] - 0.0069);
    }
    if (istation == 5) {
        lowerFiberGlassBlock_transf[0]->SetDx(lowerFiberGlassBlock_transf[0]->GetTranslation()[0] + 0.0);
        lowerFiberGlassBlock_transf[0]->SetDy(lowerFiberGlassBlock_transf[0]->GetTranslation()[1] - 0.0);
        lowerFiberGlassBlock_transf[0]->SetDz(lowerFiberGlassBlock_transf[0]->GetTranslation()[2] + 0.0306);
    }
 
    frames->AddNode(lowerFiberGlassBlockV, 0, lowerFiberGlassBlock_transf[0]);
    //--------------------------------------------------------------------------
 
    // side fiberglass block ----------------------------------------------------
    Double_t sideFiberGlassBlockPart_XSize = 1.45;   // cm
    Double_t sideFiberGlassBlockPart_YSize = 49.4;   // cm
    Double_t sideFiberGlassBlockPart_ZSize = 0.89;   // cm
 
    TGeoShape* sideFiberGlassBlockS =
        new TGeoBBox("sideFiberGlassBlockS", sideFiberGlassBlockPart_XSize * 0.5, sideFiberGlassBlockPart_YSize * 0.5,
                     sideFiberGlassBlockPart_ZSize * 0.5);
 
    TGeoVolume* sideFiberGlassBlockV =
        new TGeoVolume(TString("sideFiberGlassBlockV") += TString("_") + frames->GetName(), sideFiberGlassBlockS);
 
    // volume medium
    TGeoMedium* sideFiberGlassBlockV_medium = pMedFiberGlass;
    if (sideFiberGlassBlockV_medium) {
        sideFiberGlassBlockV->SetMedium(sideFiberGlassBlockV_medium);
    } else
        Fatal("Main", "Invalid medium for sideFiberGlassBlockV!");
 
    // volume visual property (transparency)
    sideFiberGlassBlockV->SetLineColor(TColor::GetColor("#88ff88"));
    sideFiberGlassBlockV->SetTransparency(0);
 
    TGeoCombiTrans* sideFiberGlassBlock_transf[2];
 
    sideFiberGlassBlock_transf[0] = new TGeoCombiTrans();
    sideFiberGlassBlock_transf[0]->SetDx(+(lowerFiberGlassBlock_XSize * 0.5 + sideFiberGlassBlockPart_XSize * 0.5)
                                         + 0.1);
    sideFiberGlassBlock_transf[0]->SetDy(+0.45 /*shift*/);
    sideFiberGlassBlock_transf[0]->SetDz(+0.005 /*shift for fitting-elements as a plane*/);
 
    sideFiberGlassBlock_transf[1] = new TGeoCombiTrans();
    sideFiberGlassBlock_transf[1]->SetDx(-(lowerFiberGlassBlock_XSize * 0.5 + sideFiberGlassBlockPart_XSize * 0.5));
    sideFiberGlassBlock_transf[1]->SetDy(+0.45 /*shift*/);
    sideFiberGlassBlock_transf[1]->SetDz(+0.005 /*shift for fitting-elements as a plane*/);
 
    frames->AddNode(sideFiberGlassBlockV, 0, sideFiberGlassBlock_transf[0]);
    frames->AddNode(sideFiberGlassBlockV, 1, sideFiberGlassBlock_transf[1]);
    //--------------------------------------------------------------------------
 
    // upper fiberglass block ---------------------------------------------------
    Double_t upperFiberGlassBlock_XSize = 165.1;   // cm
    Double_t upperFiberGlassBlock_YSize = 1.45;    // cm
    Double_t upperFiberGlassBlock_ZSize = 0.89;    // cm
 
    TGeoShape* upperFiberGlassBlockS = new TGeoBBox("upperFiberGlassBlockS", upperFiberGlassBlock_XSize * 0.5,
                                                    upperFiberGlassBlock_YSize * 0.5, upperFiberGlassBlock_ZSize * 0.5);
 
    TGeoVolume* upperFiberGlassBlockV =
        new TGeoVolume(TString("upperFiberGlassBlockV") += TString("_") + frames->GetName(), upperFiberGlassBlockS);
 
    // volume medium
    TGeoMedium* upperFiberGlassBlockV_medium = pMedFiberGlass;
    if (upperFiberGlassBlockV_medium) {
        upperFiberGlassBlockV->SetMedium(upperFiberGlassBlockV_medium);
    } else
        Fatal("Main", "Invalid medium for upperFiberGlassBlockV!");
 
    // volume visual property (transparency)
    upperFiberGlassBlockV->SetLineColor(TColor::GetColor("#88ff88"));
    upperFiberGlassBlockV->SetTransparency(0);
 
    TGeoCombiTrans* upperFiberGlassBlock_transf[1];
 
    upperFiberGlassBlock_transf[0] = new TGeoCombiTrans();
    upperFiberGlassBlock_transf[0]->SetDx(0.0);
    upperFiberGlassBlock_transf[0]->SetDy(
        +(sideFiberGlassBlockPart_YSize * 0.5 - upperFiberGlassBlock_YSize * 0.5 + 0.45 /*shift*/));
    upperFiberGlassBlock_transf[0]->SetDz(+0.005 /*shift for fitting-elements as a plane*/);
 
    frames->AddNode(upperFiberGlassBlockV, 0, upperFiberGlassBlock_transf[0]);
    //--------------------------------------------------------------------------
 
    // brass layer --------------------------------------------------------------
    Double_t upperBrasslayer_XSize = 165.1;   // cm
    Double_t upperBrasslayer_YSize = 1.45;    // cm
    Double_t upperBrasslayer_ZSize = 0.005;   // cm
 
    Double_t sideBrasslayer_XSize = 1.45;    // cm
    Double_t sideBrasslayer_YSize = 49.4;    // cm
    Double_t sideBrasslayer_ZSize = 0.005;   // cm
 
    TGeoShape* upperBrasslayerS = new TGeoBBox("upperBrasslayerS", upperBrasslayer_XSize * 0.5,
                                               upperBrasslayer_YSize * 0.5, upperBrasslayer_ZSize * 0.5);
    TGeoShape* sideBrasslayerS = new TGeoBBox("sideBrasslayerS", sideBrasslayer_XSize * 0.5, sideBrasslayer_YSize * 0.5,
                                              sideBrasslayer_ZSize * 0.5);
 
    TGeoVolume* upperBrasslayerV =
        new TGeoVolume(TString("upperBrasslayerV") += TString("_") + frames->GetName(), upperBrasslayerS);
    TGeoVolume* sideBrasslayerV =
        new TGeoVolume(TString("sideBrasslayerV") += TString("_") + frames->GetName(), sideBrasslayerS);
 
    // volume medium
    TGeoMedium* upperBrasslayerV_medium = pMedBrass;
    if (upperBrasslayerV_medium) {
        upperBrasslayerV->SetMedium(upperBrasslayerV_medium);
        sideBrasslayerV->SetMedium(upperBrasslayerV_medium);
    } else
        Fatal("Main", "Invalid medium for upperBrasslayerV_medium!");
 
    // volume visual property (transparency)
    upperBrasslayerV->SetLineColor(TColor::GetColor("#f07422"));
    upperBrasslayerV->SetTransparency(0);
    sideBrasslayerV->SetLineColor(TColor::GetColor("#f07422"));
    sideBrasslayerV->SetTransparency(0);
 
    TGeoCombiTrans* upperBrasslayer_transf[1];
 
    upperBrasslayer_transf[0] = new TGeoCombiTrans();
    upperBrasslayer_transf[0]->SetDx(0.0);
    upperBrasslayer_transf[0]->SetDy(
        +(sideFiberGlassBlockPart_YSize * 0.5 - upperBrasslayer_YSize * 0.5 + 0.45 /*shift*/));
    upperBrasslayer_transf[0]->SetDz(-ZModuleSize_Station163x45 * 0.5 + upperBrasslayer_ZSize * 0.5 + 0.005 /*shift*/);
 
    TGeoCombiTrans* sideBrasslayer_transf[2];
 
    sideBrasslayer_transf[0] = new TGeoCombiTrans();
    sideBrasslayer_transf[0]->SetDx(+(lowerFiberGlassBlock_XSize * 0.5 + sideBrasslayer_XSize * 0.5));
    sideBrasslayer_transf[0]->SetDy(+0.45 /*shift*/);
    sideBrasslayer_transf[0]->SetDz(-ZModuleSize_Station163x45 * 0.5 + sideBrasslayer_ZSize * 0.5 + 0.005 /*shift*/);
 
    sideBrasslayer_transf[1] = new TGeoCombiTrans();
    sideBrasslayer_transf[1]->SetDx(-(lowerFiberGlassBlock_XSize * 0.5 + sideBrasslayer_XSize * 0.5));
    sideBrasslayer_transf[1]->SetDy(+0.45 /*shift*/);
    sideBrasslayer_transf[1]->SetDz(-ZModuleSize_Station163x45 * 0.5 + sideBrasslayer_ZSize * 0.5 + 0.005 /*shift*/);
 
    frames->AddNode(upperBrasslayerV, 0, upperBrasslayer_transf[0]);
 
    frames->AddNode(sideBrasslayerV, 0, sideBrasslayer_transf[0]);
    frames->AddNode(sideBrasslayerV, 1, sideBrasslayer_transf[1]);
    //--------------------------------------------------------------------------
 
    // steel layer --------------------------------------------------------------
    Double_t upperSteellayer_XSize = 165.1;   // cm
    Double_t upperSteellayer_YSize = 1.45;    // cm
    Double_t upperSteellayer_ZSize = 0.005;   // cm
 
    Double_t sideSteellayer_XSize = 1.45;    // cm
    Double_t sideSteellayer_YSize = 49.4;    // cm
    Double_t sideSteellayer_ZSize = 0.005;   // cm
 
    TGeoShape* upperSteellayerS = new TGeoBBox("upperSteellayerS", upperSteellayer_XSize * 0.5,
                                               upperSteellayer_YSize * 0.5, upperSteellayer_ZSize * 0.5);
    TGeoShape* sideSteellayerS = new TGeoBBox("sideSteellayerS", sideSteellayer_XSize * 0.5, sideSteellayer_YSize * 0.5,
                                              sideSteellayer_ZSize * 0.5);
 
    TGeoVolume* upperSteellayerV =
        new TGeoVolume(TString("upperSteellayerV") += TString("_") + frames->GetName(), upperSteellayerS);
    TGeoVolume* sideSteellayerV =
        new TGeoVolume(TString("sideSteellayerV") += TString("_") + frames->GetName(), sideSteellayerS);
 
    // volume medium
    TGeoMedium* upperSteellayerV_medium = pMedSteel;
    if (upperSteellayerV_medium) {
        upperSteellayerV->SetMedium(upperSteellayerV_medium);
        sideSteellayerV->SetMedium(upperSteellayerV_medium);
    } else
        Fatal("Main", "Invalid medium for upperSteellayerV_medium!");
 
    // volume visual property (transparency)
    upperSteellayerV->SetLineColor(TColor::GetColor("#42aaff"));
    upperSteellayerV->SetTransparency(0);
    sideSteellayerV->SetLineColor(TColor::GetColor("#42aaff"));
    sideSteellayerV->SetTransparency(0);
 
    TGeoCombiTrans* upperSteellayer_transf[1];
 
    upperSteellayer_transf[0] = new TGeoCombiTrans();
    upperSteellayer_transf[0]->SetDx(0.0);
    upperSteellayer_transf[0]->SetDy(
        +(sideFiberGlassBlockPart_YSize * 0.5 - upperSteellayer_YSize * 0.5 + 0.45 /*shift*/));
    upperSteellayer_transf[0]->SetDz(-ZModuleSize_Station163x45 * 0.5 + upperSteellayer_ZSize * 0.5 + 0.005 /*shift*/);
 
    TGeoCombiTrans* sideSteellayer_transf[2];
 
    sideSteellayer_transf[0] = new TGeoCombiTrans();
    sideSteellayer_transf[0]->SetDx(+(lowerFiberGlassBlock_XSize * 0.5 + sideSteellayer_XSize * 0.5));
    sideSteellayer_transf[0]->SetDy(+0.45 /*shift*/);
    sideSteellayer_transf[0]->SetDz(-ZModuleSize_Station163x45 * 0.5 + sideSteellayer_ZSize * 0.5);
 
    sideSteellayer_transf[1] = new TGeoCombiTrans();
    sideSteellayer_transf[1]->SetDx(-(lowerFiberGlassBlock_XSize * 0.5 + sideSteellayer_XSize * 0.5));
    sideSteellayer_transf[1]->SetDy(+0.45 /*shift*/);
    sideSteellayer_transf[1]->SetDz(-ZModuleSize_Station163x45 * 0.5 + sideSteellayer_ZSize * 0.5);
 
    frames->AddNode(upperSteellayerV, 0, upperSteellayer_transf[0]);
 
    frames->AddNode(sideSteellayerV, 0, sideSteellayer_transf[0]);
    frames->AddNode(sideSteellayerV, 1, sideSteellayer_transf[1]);
    //--------------------------------------------------------------------------
 
    return frames;
}
//------------------------------------------------------------------------------
 
TGeoVolume* CreateModule_Station163x39(TString module_name,
                                       Double_t xsize,
                                       Double_t ysize,
                                       Double_t zsize,
                                       Double_t hole_radius)
{
 
    return CreateModule_Station163x45(module_name, xsize, ysize, zsize, hole_radius);
}
//------------------------------------------------------------------------------
 
TGeoVolume* CreateLayersForModule_Station163x39(TString layer_name)
{
 
    // frame volumes
    TGeoVolume* layers = new TGeoVolumeAssembly(layer_name);
    layers->SetMedium(pMedAir);
 
    // Common parameters for all layers
    Double_t layer_XSize = XModuleSize_Station163x39;   // cm
    Double_t layer_YSize = YModuleSize_Station163x39;   // cm
 
    Double_t layerFront_XSize = 80.85;     // cm
    Double_t layerFront_YSize = 38.2;      // cm
    Double_t layerFront_XShift = -0.375;   // cm
    Double_t layerFront_YShift = -0.2;     // cm
 
    Double_t layerBack_XSize = 80.71;     // cm
    Double_t layerBack_YSize = 38.2;      // cm
    Double_t layerBack_XShift = -0.445;   // cm
    Double_t layerBack_YShift = -0.2;     // cm
 
    // copper layer -------------------------------------------------------------
    Double_t copperLayer_ZSize = 0.00655;    // cm;
    Double_t copperLayerHole_Radius = 6.0;   // cm
 
    // layer shape
    TGeoShape* copperLayerFrontBlankS =
        new TGeoBBox(TString("copperLayerFrontBlankS") += TString("_") + layers->GetName(), layerFront_XSize * 0.5,
                     layerFront_YSize * 0.5, copperLayer_ZSize * 0.5);
    TGeoShape* copperLayerBackBlankS =
        new TGeoBBox(TString("copperLayerBackBlankS") += TString("_") + layers->GetName(), layerBack_XSize * 0.5,
                     layerBack_YSize * 0.5, copperLayer_ZSize * 0.5);
 
    TGeoShape* copperLayerFrontHoleS =
        new TGeoTube(TString("copperLayerFrontHoleS") += TString("_") + layers->GetName(), 0.0, copperLayerHole_Radius,
                     copperLayer_ZSize * 0.5 + 0.01);
    TGeoShape* copperLayerBackHoleS = new TGeoTube(TString("copperLayerBackHoleS") += TString("_") + layers->GetName(),
                                                   0.0, copperLayerHole_Radius, copperLayer_ZSize * 0.5 + 0.01);
 
    TGeoTranslation* copperLayerFrontHole_pos = new TGeoTranslation();
    copperLayerFrontHole_pos->SetName(TString("copperLayerFrontHole_pos") += TString("_") + layers->GetName());
    copperLayerFrontHole_pos->SetDx(-XModuleSize_Station163x39 * 0.5 - layerFront_XShift);
    copperLayerFrontHole_pos->SetDy(-YModuleSize_Station163x39 * 0.5 - layerFront_YShift);
    copperLayerFrontHole_pos->SetDz(0.0);
    copperLayerFrontHole_pos->RegisterYourself();
 
    TGeoTranslation* copperLayerBackHole_pos = new TGeoTranslation();
    copperLayerBackHole_pos->SetName(TString("copperLayerBackHole_pos") += TString("_") + layers->GetName());
    copperLayerBackHole_pos->SetDx(-XModuleSize_Station163x39 * 0.5 - layerBack_XShift);
    copperLayerBackHole_pos->SetDy(-YModuleSize_Station163x39 * 0.5 - layerBack_YShift);
    copperLayerBackHole_pos->SetDz(0.0);
    copperLayerBackHole_pos->RegisterYourself();
 
    TGeoCompositeShape* copperLayerFrontS = new TGeoCompositeShape();
    copperLayerFrontS->SetName(TString("copperLayerFrontS") += TString("_") + layers->GetName());
    {
        TString expression = "copperLayerFrontBlankS";
        expression += TString("_") + layers->GetName();
        expression += "-copperLayerFrontHoleS";
        expression += TString("_") + layers->GetName();
        expression += ":copperLayerFrontHole_pos";
        expression += TString("_") + layers->GetName();
 
        copperLayerFrontS->MakeNode(expression);
        copperLayerFrontS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoCompositeShape* copperLayerBackS = new TGeoCompositeShape();
    copperLayerBackS->SetName(TString("copperLayerBackS") += TString("_") + layers->GetName());
    {
        TString expression = "copperLayerBackBlankS";
        expression += TString("_") + layers->GetName();
        expression += "-copperLayerBackHoleS";
        expression += TString("_") + layers->GetName();
        expression += ":copperLayerBackHole_pos";
        expression += TString("_") + layers->GetName();
 
        copperLayerBackS->MakeNode(expression);
        copperLayerBackS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* copperLayerFrontV =
        new TGeoVolume(TString("copperLayerFrontV") += TString("_") + layers->GetName(), copperLayerFrontS);
    TGeoVolume* copperLayerBackV =
        new TGeoVolume(TString("copperLayerBackV") += TString("_") + layers->GetName(), copperLayerBackS);
 
    // volume medium
    TGeoMedium* copperLayerV_medium = pMedCopper;
    if (copperLayerV_medium) {
        copperLayerFrontV->SetMedium(copperLayerV_medium);
        copperLayerBackV->SetMedium(copperLayerV_medium);
    } else
        Fatal("Main", "Invalid medium for copperLayerV!");
 
    // volume visual property (transparency)
    copperLayerFrontV->SetLineColor(TColor::GetColor("#ff7538"));
    copperLayerFrontV->SetTransparency(0);
    copperLayerBackV->SetLineColor(TColor::GetColor("#ff7538"));
    copperLayerBackV->SetTransparency(0);
 
    TGeoCombiTrans* copperLayerFront_transf[2];
 
    copperLayerFront_transf[0] = new TGeoCombiTrans();
    copperLayerFront_transf[0]->SetDx(+layerFront_XShift);
    copperLayerFront_transf[0]->SetDy(+layerFront_YShift);
    copperLayerFront_transf[0]->SetDz(-ZGEMSize_Station163x39 * 0.5 + copperLayer_ZSize * 0.5);
 
    TGeoCombiTrans* copperLayerBack_transf[1];
 
    copperLayerBack_transf[0] = new TGeoCombiTrans();
    copperLayerBack_transf[0]->SetDx(+layerBack_XShift);
    copperLayerBack_transf[0]->SetDy(+layerBack_YShift);
    copperLayerBack_transf[0]->SetDz(+ZGEMSize_Station163x39 * 0.5 - copperLayer_ZSize * 0.5);
 
    layers->AddNode(copperLayerFrontV, 0, copperLayerFront_transf[0]);
    layers->AddNode(copperLayerBackV, 0, copperLayerBack_transf[0]);
    //--------------------------------------------------------------------------
 
    // glue layer ---------------------------------------------------------------
    Double_t glueLayer_ZSize = 0.01;       // cm; half-thickness of all glue layers
    Double_t glueLayerHole_Radius = 6.0;   // cm
 
    // layer shape
    TGeoShape* glueLayerFrontBlankS =
        new TGeoBBox(TString("glueLayerFrontBlankS") += TString("_") + layers->GetName(), layerFront_XSize * 0.5,
                     layerFront_YSize * 0.5, glueLayer_ZSize * 0.5);
    TGeoShape* glueLayerBackBlankS = new TGeoBBox(TString("glueLayerBackBlankS") += TString("_") + layers->GetName(),
                                                  layerBack_XSize * 0.5, layerBack_YSize * 0.5, glueLayer_ZSize * 0.5);
 
    TGeoShape* glueLayerFrontHoleS = new TGeoTube(TString("glueLayerFrontHoleS") += TString("_") + layers->GetName(),
                                                  0.0, glueLayerHole_Radius, glueLayer_ZSize * 0.5 + 0.01);
    TGeoShape* glueLayerBackHoleS = new TGeoTube(TString("glueLayerBackHoleS") += TString("_") + layers->GetName(), 0.0,
                                                 glueLayerHole_Radius, glueLayer_ZSize * 0.5 + 0.01);
 
    TGeoTranslation* glueLayerFrontHole_pos = new TGeoTranslation();
    glueLayerFrontHole_pos->SetName(TString("glueLayerFrontHole_pos") += TString("_") + layers->GetName());
    glueLayerFrontHole_pos->SetDx(-XModuleSize_Station163x39 * 0.5 - layerFront_XShift);
    glueLayerFrontHole_pos->SetDy(-YModuleSize_Station163x39 * 0.5 - layerFront_YShift);
    glueLayerFrontHole_pos->SetDz(0.0);
    glueLayerFrontHole_pos->RegisterYourself();
 
    TGeoTranslation* glueLayerBackHole_pos = new TGeoTranslation();
    glueLayerBackHole_pos->SetName(TString("glueLayerBackHole_pos") += TString("_") + layers->GetName());
    glueLayerBackHole_pos->SetDx(-XModuleSize_Station163x39 * 0.5 - layerBack_XShift);
    glueLayerBackHole_pos->SetDy(-YModuleSize_Station163x39 * 0.5 - layerBack_YShift);
    glueLayerBackHole_pos->SetDz(0.0);
    glueLayerBackHole_pos->RegisterYourself();
 
    TGeoCompositeShape* glueLayerFrontS = new TGeoCompositeShape();
    glueLayerFrontS->SetName(TString("glueLayerFrontS") += TString("_") + layers->GetName());
    {
        TString expression = "glueLayerFrontBlankS";
        expression += TString("_") + layers->GetName();
        expression += "-glueLayerFrontHoleS";
        expression += TString("_") + layers->GetName();
        expression += ":glueLayerFrontHole_pos";
        expression += TString("_") + layers->GetName();
 
        glueLayerFrontS->MakeNode(expression);
        glueLayerFrontS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoCompositeShape* glueLayerBackS = new TGeoCompositeShape();
    glueLayerBackS->SetName(TString("glueLayerBackS") += TString("_") + layers->GetName());
    {
        TString expression = "glueLayerBackBlankS";
        expression += TString("_") + layers->GetName();
        expression += "-glueLayerBackHoleS";
        expression += TString("_") + layers->GetName();
        expression += ":glueLayerBackHole_pos";
        expression += TString("_") + layers->GetName();
 
        glueLayerBackS->MakeNode(expression);
        glueLayerBackS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* glueLayerFrontV =
        new TGeoVolume(TString("glueLayerFrontV") += TString("_") + layers->GetName(), glueLayerFrontS);
    TGeoVolume* glueLayerBackV =
        new TGeoVolume(TString("glueLayerBackV") += TString("_") + layers->GetName(), glueLayerBackS);
 
    // volume medium
    TGeoMedium* glueLayerV_medium = pMedGlue;
    if (glueLayerV_medium) {
        glueLayerFrontV->SetMedium(glueLayerV_medium);
        glueLayerBackV->SetMedium(glueLayerV_medium);
    } else
        Fatal("Main", "Invalid medium for glueLayerV!");
 
    // volume visual property (transparency)
    glueLayerFrontV->SetLineColor(TColor::GetColor("#1919ff"));
    glueLayerFrontV->SetTransparency(0);
    glueLayerBackV->SetLineColor(TColor::GetColor("#1919ff"));
    glueLayerBackV->SetTransparency(0);
 
    TGeoCombiTrans* glueLayerFront_transf[1];
 
    glueLayerFront_transf[0] = new TGeoCombiTrans();
    glueLayerFront_transf[0]->SetDx(+layerFront_XShift);
    glueLayerFront_transf[0]->SetDy(+layerFront_YShift);
    glueLayerFront_transf[0]->SetDz(-ZGEMSize_Station163x39 * 0.5 + copperLayer_ZSize + glueLayer_ZSize * 0.5);
 
    TGeoCombiTrans* glueLayerBack_transf[1];
 
    glueLayerBack_transf[0] = new TGeoCombiTrans();
    glueLayerBack_transf[0]->SetDx(+layerBack_XShift);
    glueLayerBack_transf[0]->SetDy(+layerBack_YShift);
    glueLayerBack_transf[0]->SetDz(+ZGEMSize_Station163x39 * 0.5 - copperLayer_ZSize - glueLayer_ZSize * 0.5);
 
    layers->AddNode(glueLayerFrontV, 0, glueLayerFront_transf[0]);
    layers->AddNode(glueLayerBackV, 0, glueLayerBack_transf[0]);
    //--------------------------------------------------------------------------
 
    // epoxide layer ------------------------------------------------------------
    Double_t epoxideLayer_ZSize = 0.105;      // cm; half-thickness of all epoxide layers
    Double_t epoxideLayerHole_Radius = 6.0;   // cm
 
    // layer shape
    TGeoShape* epoxideLayerFrontBlankS =
        new TGeoBBox(TString("epoxideLayerFrontBlankS") += TString("_") + layers->GetName(), layerFront_XSize * 0.5,
                     layerFront_YSize * 0.5, epoxideLayer_ZSize * 0.5);
    TGeoShape* epoxideLayerBackBlankS =
        new TGeoBBox(TString("epoxideLayerBackBlankS") += TString("_") + layers->GetName(), layerBack_XSize * 0.5,
                     layerBack_YSize * 0.5, epoxideLayer_ZSize * 0.5);
 
    TGeoShape* epoxideLayerFrontHoleS =
        new TGeoTube(TString("epoxideLayerFrontHoleS") += TString("_") + layers->GetName(), 0.0,
                     epoxideLayerHole_Radius, epoxideLayer_ZSize * 0.5 + 0.01);
    TGeoShape* epoxideLayerBackHoleS =
        new TGeoTube(TString("epoxideLayerBackHoleS") += TString("_") + layers->GetName(), 0.0, epoxideLayerHole_Radius,
                     epoxideLayer_ZSize * 0.5 + 0.01);
 
    TGeoTranslation* epoxideLayerFrontHole_pos = new TGeoTranslation();
    epoxideLayerFrontHole_pos->SetName(TString("epoxideLayerFrontHole_pos") += TString("_") + layers->GetName());
    epoxideLayerFrontHole_pos->SetDx(-XModuleSize_Station163x39 * 0.5 - layerFront_XShift);
    epoxideLayerFrontHole_pos->SetDy(-YModuleSize_Station163x39 * 0.5 - layerFront_YShift);
    epoxideLayerFrontHole_pos->SetDz(0.0);
    epoxideLayerFrontHole_pos->RegisterYourself();
 
    TGeoTranslation* epoxideLayerBackHole_pos = new TGeoTranslation();
    epoxideLayerBackHole_pos->SetName(TString("epoxideLayerBackHole_pos") += TString("_") + layers->GetName());
    epoxideLayerBackHole_pos->SetDx(-XModuleSize_Station163x39 * 0.5 - layerBack_XShift);
    epoxideLayerBackHole_pos->SetDy(-YModuleSize_Station163x39 * 0.5 - layerBack_YShift);
    epoxideLayerBackHole_pos->SetDz(0.0);
    epoxideLayerBackHole_pos->RegisterYourself();
 
    TGeoCompositeShape* epoxideLayerFrontS = new TGeoCompositeShape();
    epoxideLayerFrontS->SetName(TString("epoxideLayerFrontS") += TString("_") + layers->GetName());
    {
        TString expression = "epoxideLayerFrontBlankS";
        expression += TString("_") + layers->GetName();
        expression += "-epoxideLayerFrontHoleS";
        expression += TString("_") + layers->GetName();
        expression += ":epoxideLayerFrontHole_pos";
        expression += TString("_") + layers->GetName();
 
        epoxideLayerFrontS->MakeNode(expression);
        epoxideLayerFrontS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoCompositeShape* epoxideLayerBackS = new TGeoCompositeShape();
    epoxideLayerBackS->SetName(TString("epoxideLayerBackS") += TString("_") + layers->GetName());
    {
        TString expression = "epoxideLayerBackBlankS";
        expression += TString("_") + layers->GetName();
        expression += "-epoxideLayerBackHoleS";
        expression += TString("_") + layers->GetName();
        expression += ":epoxideLayerBackHole_pos";
        expression += TString("_") + layers->GetName();
 
        epoxideLayerBackS->MakeNode(expression);
        epoxideLayerBackS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* epoxideLayerFrontV =
        new TGeoVolume(TString("epoxideLayerFrontV") += TString("_") + layers->GetName(), epoxideLayerFrontS);
    TGeoVolume* epoxideLayerBackV =
        new TGeoVolume(TString("epoxideLayerBackV") += TString("_") + layers->GetName(), epoxideLayerBackS);
 
    // volume medium
    TGeoMedium* epoxideLayerV_medium = pMedEpoxideCompound;
    if (epoxideLayerV_medium) {
        epoxideLayerFrontV->SetMedium(epoxideLayerV_medium);
        epoxideLayerBackV->SetMedium(epoxideLayerV_medium);
    } else
        Fatal("Main", "Invalid medium for epoxideLayerV!");
 
    // volume visual property (transparency)
    epoxideLayerFrontV->SetLineColor(TColor::GetColor("#aa33ff"));
    epoxideLayerFrontV->SetTransparency(0);
    epoxideLayerBackV->SetLineColor(TColor::GetColor("#aa33ff"));
    epoxideLayerBackV->SetTransparency(0);
 
    TGeoCombiTrans* epoxideLayerFront_transf[1];
 
    epoxideLayerFront_transf[0] = new TGeoCombiTrans();
    epoxideLayerFront_transf[0]->SetDx(+layerFront_XShift);
    epoxideLayerFront_transf[0]->SetDy(+layerFront_YShift);
    epoxideLayerFront_transf[0]->SetDz(-ZGEMSize_Station163x39 * 0.5 + copperLayer_ZSize + glueLayer_ZSize
                                       + epoxideLayer_ZSize * 0.5);
 
    TGeoCombiTrans* epoxideLayerBack_transf[1];
 
    epoxideLayerBack_transf[0] = new TGeoCombiTrans();
    epoxideLayerBack_transf[0]->SetDx(+layerBack_XShift);
    epoxideLayerBack_transf[0]->SetDy(+layerBack_YShift);
    epoxideLayerBack_transf[0]->SetDz(+ZGEMSize_Station163x39 * 0.5 - copperLayer_ZSize - glueLayer_ZSize
                                      - epoxideLayer_ZSize * 0.5);
 
    layers->AddNode(epoxideLayerFrontV, 0, epoxideLayerFront_transf[0]);
    layers->AddNode(epoxideLayerBackV, 0, epoxideLayerBack_transf[0]);
    //--------------------------------------------------------------------------
 
    // honeycomb layer ----------------------------------------------------------
    // Double_t honeycombLayer_ZSize = 1.5-(0.05+0.00055+0.005+0.0125); //cm; half-thickness of all honeycomb layers
    Double_t honeycombLayer_ZSize = 1.5 - (0.06805);   // cm; half-thickness of all honeycomb layers
    Double_t honeycombLayerHole_Radius = 6.0;          // cm
 
    // layer shape
    TGeoShape* honeycombLayerFrontBlankS =
        new TGeoBBox(TString("honeycombLayerFrontBlankS") += TString("_") + layers->GetName(), layerFront_XSize * 0.5,
                     layerFront_YSize * 0.5, honeycombLayer_ZSize * 0.5);
    TGeoShape* honeycombLayerBackBlankS =
        new TGeoBBox(TString("honeycombLayerBackBlankS") += TString("_") + layers->GetName(), layerBack_XSize * 0.5,
                     layerBack_YSize * 0.5, honeycombLayer_ZSize * 0.5);
 
    TGeoShape* honeycombLayerFrontHoleS =
        new TGeoTube(TString("honeycombLayerFrontHoleS") += TString("_") + layers->GetName(), 0.0,
                     honeycombLayerHole_Radius, honeycombLayer_ZSize * 0.5 + 0.01);
    TGeoShape* honeycombLayerBackHoleS =
        new TGeoTube(TString("honeycombLayerBackHoleS") += TString("_") + layers->GetName(), 0.0,
                     honeycombLayerHole_Radius, honeycombLayer_ZSize * 0.5 + 0.01);
 
    TGeoTranslation* honeycombLayerFrontHole_pos = new TGeoTranslation();
    honeycombLayerFrontHole_pos->SetName(TString("honeycombLayerFrontHole_pos") += TString("_") + layers->GetName());
    honeycombLayerFrontHole_pos->SetDx(-XModuleSize_Station163x39 * 0.5 - layerFront_XShift);
    honeycombLayerFrontHole_pos->SetDy(-YModuleSize_Station163x39 * 0.5 - layerFront_YShift);
    honeycombLayerFrontHole_pos->SetDz(0.0);
    honeycombLayerFrontHole_pos->RegisterYourself();
 
    TGeoTranslation* honeycombLayerBackHole_pos = new TGeoTranslation();
    honeycombLayerBackHole_pos->SetName(TString("honeycombLayerBackHole_pos") += TString("_") + layers->GetName());
    honeycombLayerBackHole_pos->SetDx(-XModuleSize_Station163x39 * 0.5 - layerBack_XShift);
    honeycombLayerBackHole_pos->SetDy(-YModuleSize_Station163x39 * 0.5 - layerBack_YShift);
    honeycombLayerBackHole_pos->SetDz(0.0);
    honeycombLayerBackHole_pos->RegisterYourself();
 
    TGeoCompositeShape* honeycombLayerFrontS = new TGeoCompositeShape();
    honeycombLayerFrontS->SetName(TString("honeycombLayerFrontS") += TString("_") + layers->GetName());
    {
        TString expression = "honeycombLayerFrontBlankS";
        expression += TString("_") + layers->GetName();
        expression += "-honeycombLayerFrontHoleS";
        expression += TString("_") + layers->GetName();
        expression += ":honeycombLayerFrontHole_pos";
        expression += TString("_") + layers->GetName();
 
        honeycombLayerFrontS->MakeNode(expression);
        honeycombLayerFrontS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoCompositeShape* honeycombLayerBackS = new TGeoCompositeShape();
    honeycombLayerBackS->SetName(TString("honeycombLayerBackS") += TString("_") + layers->GetName());
    {
        TString expression = "honeycombLayerBackBlankS";
        expression += TString("_") + layers->GetName();
        expression += "-honeycombLayerBackHoleS";
        expression += TString("_") + layers->GetName();
        expression += ":honeycombLayerBackHole_pos";
        expression += TString("_") + layers->GetName();
 
        honeycombLayerBackS->MakeNode(expression);
        honeycombLayerBackS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* honeycombLayerFrontV =
        new TGeoVolume(TString("honeycombLayerFrontV") += TString("_") + layers->GetName(), honeycombLayerFrontS);
    TGeoVolume* honeycombLayerBackV =
        new TGeoVolume(TString("honeycombLayerBackV") += TString("_") + layers->GetName(), honeycombLayerBackS);
 
    // volume medium
    TGeoMedium* honeycombLayerV_medium = pMedHoneyComb;
    if (honeycombLayerV_medium) {
        honeycombLayerFrontV->SetMedium(honeycombLayerV_medium);
        honeycombLayerBackV->SetMedium(honeycombLayerV_medium);
    } else
        Fatal("Main", "Invalid medium for honeycombLayerV!");
 
    // volume visual property (transparency)
    honeycombLayerFrontV->SetLineColor(TColor::GetColor("#ffdf84"));
    honeycombLayerFrontV->SetTransparency(0);
    honeycombLayerBackV->SetLineColor(TColor::GetColor("#ffdf84"));
    honeycombLayerBackV->SetTransparency(0);
 
    TGeoCombiTrans* honeycombLayerFront_transf[1];
 
    honeycombLayerFront_transf[0] = new TGeoCombiTrans();
    honeycombLayerFront_transf[0]->SetDx(+layerFront_XShift);
    honeycombLayerFront_transf[0]->SetDy(+layerFront_YShift);
    honeycombLayerFront_transf[0]->SetDz(-ZGEMSize_Station163x39 * 0.5 + copperLayer_ZSize + glueLayer_ZSize
                                         + epoxideLayer_ZSize + honeycombLayer_ZSize * 0.5);
 
    TGeoCombiTrans* honeycombLayerBack_transf[1];
 
    honeycombLayerBack_transf[0] = new TGeoCombiTrans();
    honeycombLayerBack_transf[0]->SetDx(+layerBack_XShift);
    honeycombLayerBack_transf[0]->SetDy(+layerBack_YShift);
    honeycombLayerBack_transf[0]->SetDz(+ZGEMSize_Station163x39 * 0.5 - copperLayer_ZSize - glueLayer_ZSize
                                        - epoxideLayer_ZSize - honeycombLayer_ZSize * 0.5);
 
    layers->AddNode(honeycombLayerFrontV, 0, honeycombLayerFront_transf[0]);
    layers->AddNode(honeycombLayerBackV, 0, honeycombLayerBack_transf[0]);
    //--------------------------------------------------------------------------
 
    // polyamide layer ----------------------------------------------------------
    Double_t polyamideLayer_ZSize = 0.0125;     // cm; half-thickness of all polyamide layers
    Double_t polyamideLayerHole_Radius = 6.0;   // cm
 
    // layer shape
    TGeoShape* polyamideLayerFrontBlankS =
        new TGeoBBox(TString("polyamideLayerFrontBlankS") += TString("_") + layers->GetName(), layerFront_XSize * 0.5,
                     layerFront_YSize * 0.5, polyamideLayer_ZSize * 0.5);
    TGeoShape* polyamideLayerBackBlankS =
        new TGeoBBox(TString("polyamideLayerBackBlankS") += TString("_") + layers->GetName(), layerBack_XSize * 0.5,
                     layerBack_YSize * 0.5, polyamideLayer_ZSize * 0.5);
 
    TGeoShape* polyamideLayerFrontHoleS =
        new TGeoTube(TString("polyamideLayerFrontHoleS") += TString("_") + layers->GetName(), 0.0,
                     polyamideLayerHole_Radius, polyamideLayer_ZSize * 0.5 + 0.01);
    TGeoShape* polyamideLayerBackHoleS =
        new TGeoTube(TString("polyamideLayerBackHoleS") += TString("_") + layers->GetName(), 0.0,
                     polyamideLayerHole_Radius, polyamideLayer_ZSize * 0.5 + 0.01);
 
    TGeoTranslation* polyamideLayerFrontHole_pos = new TGeoTranslation();
    polyamideLayerFrontHole_pos->SetName(TString("polyamideLayerFrontHole_pos") += TString("_") + layers->GetName());
    polyamideLayerFrontHole_pos->SetDx(-XModuleSize_Station163x39 * 0.5 - layerFront_XShift);
    polyamideLayerFrontHole_pos->SetDy(-YModuleSize_Station163x39 * 0.5 - layerFront_YShift);
    polyamideLayerFrontHole_pos->SetDz(0.0);
    polyamideLayerFrontHole_pos->RegisterYourself();
 
    TGeoTranslation* polyamideLayerBackHole_pos = new TGeoTranslation();
    polyamideLayerBackHole_pos->SetName(TString("polyamideLayerBackHole_pos") += TString("_") + layers->GetName());
    polyamideLayerBackHole_pos->SetDx(-XModuleSize_Station163x39 * 0.5 - layerBack_XShift);
    polyamideLayerBackHole_pos->SetDy(-YModuleSize_Station163x39 * 0.5 - layerBack_YShift);
    polyamideLayerBackHole_pos->SetDz(0.0);
    polyamideLayerBackHole_pos->RegisterYourself();
 
    TGeoCompositeShape* polyamideLayerFrontS = new TGeoCompositeShape();
    polyamideLayerFrontS->SetName(TString("polyamideLayerFrontS") += TString("_") + layers->GetName());
    {
        TString expression = "polyamideLayerFrontBlankS";
        expression += TString("_") + layers->GetName();
        expression += "-polyamideLayerFrontHoleS";
        expression += TString("_") + layers->GetName();
        expression += ":polyamideLayerFrontHole_pos";
        expression += TString("_") + layers->GetName();
 
        polyamideLayerFrontS->MakeNode(expression);
        polyamideLayerFrontS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoCompositeShape* polyamideLayerBackS = new TGeoCompositeShape();
    polyamideLayerBackS->SetName(TString("polyamideLayerBackS") += TString("_") + layers->GetName());
    {
        TString expression = "polyamideLayerBackBlankS";
        expression += TString("_") + layers->GetName();
        expression += "-polyamideLayerBackHoleS";
        expression += TString("_") + layers->GetName();
        expression += ":polyamideLayerBackHole_pos";
        expression += TString("_") + layers->GetName();
 
        polyamideLayerBackS->MakeNode(expression);
        polyamideLayerBackS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* polyamideLayerFrontV =
        new TGeoVolume(TString("polyamideLayerFrontV") += TString("_") + layers->GetName(), polyamideLayerFrontS);
    TGeoVolume* polyamideLayerBackV =
        new TGeoVolume(TString("polyamideLayerBackV") += TString("_") + layers->GetName(), polyamideLayerBackS);
 
    // volume medium
    TGeoMedium* polyamideLayerV_medium = pMedPolyamide;
    if (polyamideLayerV_medium) {
        polyamideLayerFrontV->SetMedium(polyamideLayerV_medium);
        polyamideLayerBackV->SetMedium(polyamideLayerV_medium);
    } else
        Fatal("Main", "Invalid medium for polyamideLayerV!");
 
    // volume visual property (transparency)
    polyamideLayerFrontV->SetLineColor(TColor::GetColor("#fefe22"));
    polyamideLayerFrontV->SetTransparency(0);
    polyamideLayerBackV->SetLineColor(TColor::GetColor("#fefe22"));
    polyamideLayerBackV->SetTransparency(0);
 
    TGeoCombiTrans* polyamideLayerFront_transf[1];
 
    polyamideLayerFront_transf[0] = new TGeoCombiTrans();
    polyamideLayerFront_transf[0]->SetDx(+layerFront_XShift);
    polyamideLayerFront_transf[0]->SetDy(+layerFront_YShift);
    polyamideLayerFront_transf[0]->SetDz(-ZGEMSize_Station163x39 * 0.5 + copperLayer_ZSize + glueLayer_ZSize
                                         + epoxideLayer_ZSize + honeycombLayer_ZSize + polyamideLayer_ZSize * 0.5);
 
    TGeoCombiTrans* polyamideLayerBack_transf[1];
 
    polyamideLayerBack_transf[0] = new TGeoCombiTrans();
    polyamideLayerBack_transf[0]->SetDx(+layerBack_XShift);
    polyamideLayerBack_transf[0]->SetDy(+layerBack_YShift);
    polyamideLayerBack_transf[0]->SetDz(+ZGEMSize_Station163x39 * 0.5 - copperLayer_ZSize - glueLayer_ZSize
                                        - epoxideLayer_ZSize - honeycombLayer_ZSize - polyamideLayer_ZSize * 0.5);
 
    layers->AddNode(polyamideLayerFrontV, 0, polyamideLayerFront_transf[0]);
    layers->AddNode(polyamideLayerBackV, 0, polyamideLayerBack_transf[0]);
    //--------------------------------------------------------------------------
 
    return layers;
}
//------------------------------------------------------------------------------
 
TGeoVolume* CreateFrameForHalfPlane_Station163x39(TString frame_name, Int_t istation)
{
 
    // frame volumes
    TGeoVolume* frames = new TGeoVolumeAssembly(frame_name);
    frames->SetMedium(pMedAir);
 
    // back lower epoxide frame -------------------------------------------------
    Double_t backLowerEpoxideFrame_XSize = 161.42;   // cm
    Double_t backLowerEpoxideFrame_YSize = 2.2;      // cm
    Double_t backLowerEpoxideFrame_ZSize = 1.4;      // cm
 
    Double_t backLowerEpoxideFrameArch_Rmin = 4.0;
    Double_t backLowerEpoxideFrameArch_Rmax = 6.0;
    Double_t backLowerEpoxideFrameArch_ZSize = backLowerEpoxideFrame_ZSize;
 
    TGeoShape* backLowerEpoxideFrameHalfS =
        new TGeoBBox(TString("backLowerEpoxideFrameHalfS") += TString("_") + frames->GetName(),
                     (backLowerEpoxideFrame_XSize * 0.5 - backLowerEpoxideFrameArch_Rmin) * 0.5,
                     backLowerEpoxideFrame_YSize * 0.5, backLowerEpoxideFrame_ZSize * 0.5);
    TGeoShape* backLowerEpoxideFrameArchS = new TGeoTubeSeg(
        TString("backLowerEpoxideFrameArchS") += TString("_") + frames->GetName(), backLowerEpoxideFrameArch_Rmin,
        backLowerEpoxideFrameArch_Rmax, backLowerEpoxideFrameArch_ZSize * 0.5, 0.0, 180.0);
 
    TGeoTranslation* backLowerEpoxideFrameHalf_pos1 = new TGeoTranslation();
    backLowerEpoxideFrameHalf_pos1->SetName("backLowerEpoxideFrameHalf1_pos1");
    backLowerEpoxideFrameHalf_pos1->SetDx(
        +((backLowerEpoxideFrame_XSize * 0.5 - backLowerEpoxideFrameArch_Rmin) * 0.5 + backLowerEpoxideFrameArch_Rmin));
    backLowerEpoxideFrameHalf_pos1->SetDy(-(backLowerEpoxideFrame_YSize * 0.5) + 0.2 /*shift*/);
    backLowerEpoxideFrameHalf_pos1->SetDz(0.0);
    backLowerEpoxideFrameHalf_pos1->RegisterYourself();
 
    TGeoTranslation* backLowerEpoxideFrameHalf_pos2 = new TGeoTranslation();
    backLowerEpoxideFrameHalf_pos2->SetName("backLowerEpoxideFrameHalf1_pos2");
    backLowerEpoxideFrameHalf_pos2->SetDx(
        -((backLowerEpoxideFrame_XSize * 0.5 - backLowerEpoxideFrameArch_Rmin) * 0.5 + backLowerEpoxideFrameArch_Rmin));
    backLowerEpoxideFrameHalf_pos2->SetDy(-(backLowerEpoxideFrame_YSize * 0.5) + 0.2 /*shift*/);
    backLowerEpoxideFrameHalf_pos2->SetDz(0.0);
    backLowerEpoxideFrameHalf_pos2->RegisterYourself();
 
    TGeoTranslation* backLowerEpoxideFrameArch_pos = new TGeoTranslation();
    backLowerEpoxideFrameArch_pos->SetName("backLowerEpoxideFrameArch_pos");
    backLowerEpoxideFrameArch_pos->SetDx(0.0);
    backLowerEpoxideFrameArch_pos->SetDy(0.0);
    backLowerEpoxideFrameArch_pos->SetDz(0.0);
    backLowerEpoxideFrameArch_pos->RegisterYourself();
 
    TGeoCompositeShape* backLowerEpoxideFrameS = new TGeoCompositeShape();
    backLowerEpoxideFrameS->SetName(TString("backLowerEpoxideFrameS") += TString("_") + frames->GetName());
    {
        TString expression = "backLowerEpoxideFrameHalfS";
        expression += TString("_") + frames->GetName();
        expression += ":backLowerEpoxideFrameHalf1_pos1";
        expression += "+backLowerEpoxideFrameHalfS";
        expression += TString("_") + frames->GetName();
        expression += ":backLowerEpoxideFrameHalf1_pos2";
        expression += "+backLowerEpoxideFrameArchS";
        expression += TString("_") + frames->GetName();
        expression += ":backLowerEpoxideFrameArch_pos";
 
        backLowerEpoxideFrameS->MakeNode(expression);
        backLowerEpoxideFrameS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* backLowerEpoxideFrameV =
        new TGeoVolume(TString("backLowerEpoxideFrameV") += TString("_") + frames->GetName(), backLowerEpoxideFrameS);
 
    // volume medium
    TGeoMedium* backLowerEpoxideFrameV_medium = pMedEpoxideCompound;
    if (backLowerEpoxideFrameV_medium) {
        backLowerEpoxideFrameV->SetMedium(backLowerEpoxideFrameV_medium);
    } else
        Fatal("Main", "Invalid medium for backLowerEpoxideFrameV!");
 
    // volume visual property (transparency)
    backLowerEpoxideFrameV->SetLineColor(TColor::GetColor("#ccccff"));
    backLowerEpoxideFrameV->SetTransparency(0);
 
    TGeoCombiTrans* backLowerEpoxideFrame_transf[1];
 
    backLowerEpoxideFrame_transf[0] = new TGeoCombiTrans();
    backLowerEpoxideFrame_transf[0]->SetDx(-0.15);
    backLowerEpoxideFrame_transf[0]->SetDy(-(YModuleSize_Station163x39 * 0.5) - 0.45);
    backLowerEpoxideFrame_transf[0]->SetDz(+(backLowerEpoxideFrame_ZSize * 0.5 + ZModuleSize_Station163x39 * 0.5)
                                           + 0.015);
 
    frames->AddNode(backLowerEpoxideFrameV, 0, backLowerEpoxideFrame_transf[0]);
    //--------------------------------------------------------------------------
 
    // front lower epoxide frame -------------------------------------------------
    Double_t frontLowerEpoxideFrame_XSize = 161.7;   // cm
    Double_t frontLowerEpoxideFrame_YSize = 2.2;     // cm
    Double_t frontLowerEpoxideFrame_ZSize = 1.5;     // cm
 
    Double_t frontLowerEpoxideFrameArch_Rmin = 4.0;
    Double_t frontLowerEpoxideFrameArch_Rmax = 6.0;
    Double_t frontLowerEpoxideFrameArch_ZSize = frontLowerEpoxideFrame_ZSize;
 
    TGeoShape* frontLowerEpoxideFrameHalfS =
        new TGeoBBox(TString("frontLowerEpoxideFrameHalfS") += TString("_") + frames->GetName(),
                     (frontLowerEpoxideFrame_XSize * 0.5 - frontLowerEpoxideFrameArch_Rmin) * 0.5,
                     frontLowerEpoxideFrame_YSize * 0.5, frontLowerEpoxideFrame_ZSize * 0.5);
    TGeoShape* frontLowerEpoxideFrameArchS = new TGeoTubeSeg(
        TString("frontLowerEpoxideFrameArchS") += TString("_") + frames->GetName(), frontLowerEpoxideFrameArch_Rmin,
        frontLowerEpoxideFrameArch_Rmax, frontLowerEpoxideFrameArch_ZSize * 0.5, 0.0, 180.0);
 
    TGeoTranslation* frontLowerEpoxideFrameHalf_pos1 = new TGeoTranslation();
    frontLowerEpoxideFrameHalf_pos1->SetName("frontLowerEpoxideFrameHalf1_pos1");
    frontLowerEpoxideFrameHalf_pos1->SetDx(
        +((frontLowerEpoxideFrame_XSize * 0.5 - frontLowerEpoxideFrameArch_Rmin) * 0.5
          + frontLowerEpoxideFrameArch_Rmin));
    frontLowerEpoxideFrameHalf_pos1->SetDy(-(frontLowerEpoxideFrame_YSize * 0.5) + 0.2 /*shift*/);
    frontLowerEpoxideFrameHalf_pos1->SetDz(0.0);
    frontLowerEpoxideFrameHalf_pos1->RegisterYourself();
 
    TGeoTranslation* frontLowerEpoxideFrameHalf_pos2 = new TGeoTranslation();
    frontLowerEpoxideFrameHalf_pos2->SetName("frontLowerEpoxideFrameHalf1_pos2");
    frontLowerEpoxideFrameHalf_pos2->SetDx(
        -((frontLowerEpoxideFrame_XSize * 0.5 - frontLowerEpoxideFrameArch_Rmin) * 0.5
          + frontLowerEpoxideFrameArch_Rmin));
    frontLowerEpoxideFrameHalf_pos2->SetDy(-(frontLowerEpoxideFrame_YSize * 0.5) + 0.2 /*shift*/);
    frontLowerEpoxideFrameHalf_pos2->SetDz(0.0);
    frontLowerEpoxideFrameHalf_pos2->RegisterYourself();
 
    TGeoTranslation* frontLowerEpoxideFrameArch_pos = new TGeoTranslation();
    frontLowerEpoxideFrameArch_pos->SetName("frontLowerEpoxideFrameArch_pos");
    frontLowerEpoxideFrameArch_pos->SetDx(0.0);
    frontLowerEpoxideFrameArch_pos->SetDy(0.0);
    frontLowerEpoxideFrameArch_pos->SetDz(0.0);
    frontLowerEpoxideFrameArch_pos->RegisterYourself();
 
    TGeoCompositeShape* frontLowerEpoxideFrameS = new TGeoCompositeShape();
    frontLowerEpoxideFrameS->SetName(TString("frontLowerEpoxideFrameS") += TString("_") + frames->GetName());
    {
        TString expression = "frontLowerEpoxideFrameHalfS";
        expression += TString("_") + frames->GetName();
        expression += ":frontLowerEpoxideFrameHalf1_pos1";
        expression += "+frontLowerEpoxideFrameHalfS";
        expression += TString("_") + frames->GetName();
        expression += ":frontLowerEpoxideFrameHalf1_pos2";
        expression += "+frontLowerEpoxideFrameArchS";
        expression += TString("_") + frames->GetName();
        expression += ":frontLowerEpoxideFrameArch_pos";
 
        frontLowerEpoxideFrameS->MakeNode(expression);
        frontLowerEpoxideFrameS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* frontLowerEpoxideFrameV =
        new TGeoVolume(TString("frontLowerEpoxideFrameV") += TString("_") + frames->GetName(), frontLowerEpoxideFrameS);
 
    // volume medium
    TGeoMedium* frontLowerEpoxideFrameV_medium = pMedEpoxideCompound;
    if (frontLowerEpoxideFrameV_medium) {
        frontLowerEpoxideFrameV->SetMedium(frontLowerEpoxideFrameV_medium);
    } else
        Fatal("Main", "Invalid medium for frontLowerEpoxideFrameV!");
 
    // volume visual property (transparency)
    frontLowerEpoxideFrameV->SetLineColor(TColor::GetColor("#ccccff"));
    frontLowerEpoxideFrameV->SetTransparency(0);
 
    TGeoCombiTrans* frontLowerEpoxideFrame_transf[1];
 
    frontLowerEpoxideFrame_transf[0] = new TGeoCombiTrans();
    frontLowerEpoxideFrame_transf[0]->SetDx(-0.2);
    frontLowerEpoxideFrame_transf[0]->SetDy(-(YModuleSize_Station163x39 * 0.5) - 0.2);
    frontLowerEpoxideFrame_transf[0]->SetDz(-(frontLowerEpoxideFrame_ZSize * 0.5 + ZModuleSize_Station163x39 * 0.5)
                                            - 0.1);
 
    frames->AddNode(frontLowerEpoxideFrameV, 0, frontLowerEpoxideFrame_transf[0]);
    //--------------------------------------------------------------------------
 
    // back side epoxide frame --------------------------------------------------
    Double_t backSideEpoxideFramePart_XSize = 9.2;          // cm
    Double_t backSideEpoxideFramePart_YSize = 52.6 - 6.0;   // cm
    Double_t backSideEpoxideFramePart_ZSize = 1.5;          // cm
 
    Double_t backSideEpoxideFrame_XDistBetweenParts = 163.42;   // cm
 
    TGeoShape* backSideEpoxideFramePartS = new TGeoBBox(
        TString("backSideEpoxideFramePartS") += TString("_") + frames->GetName(), backSideEpoxideFramePart_XSize * 0.5,
        backSideEpoxideFramePart_YSize * 0.5, backSideEpoxideFramePart_ZSize * 0.5);
 
    TGeoTranslation* backSideEpoxideFramePart_pos1 = new TGeoTranslation();
    backSideEpoxideFramePart_pos1->SetName("backSideEpoxideFramePart_pos1_low_plane");
    backSideEpoxideFramePart_pos1->SetDx(
        +(backSideEpoxideFrame_XDistBetweenParts * 0.5 + backSideEpoxideFramePart_XSize * 0.5));
    backSideEpoxideFramePart_pos1->SetDy(0.0);
    backSideEpoxideFramePart_pos1->SetDz(0.0);
    backSideEpoxideFramePart_pos1->RegisterYourself();
 
    TGeoTranslation* backSideEpoxideFramePart_pos2 = new TGeoTranslation();
    backSideEpoxideFramePart_pos2->SetName("backSideEpoxideFramePart_pos2_low_plane");
    backSideEpoxideFramePart_pos2->SetDx(
        -(backSideEpoxideFrame_XDistBetweenParts * 0.5 + backSideEpoxideFramePart_XSize * 0.5));
    backSideEpoxideFramePart_pos2->SetDy(0.0);
    backSideEpoxideFramePart_pos2->SetDz(0.0);
    backSideEpoxideFramePart_pos2->RegisterYourself();
 
    TGeoCompositeShape* backSideEpoxideFrameS = new TGeoCompositeShape();
    backSideEpoxideFrameS->SetName(TString("backSideEpoxideFrameS") += TString("_") + frames->GetName());
    {
        TString expression = "backSideEpoxideFramePartS";
        expression += TString("_") + frames->GetName();
        expression += ":backSideEpoxideFramePart_pos1_low_plane";
        expression += "+backSideEpoxideFramePartS";
        expression += TString("_") + frames->GetName();
        expression += ":backSideEpoxideFramePart_pos2_low_plane";
 
        backSideEpoxideFrameS->MakeNode(expression);
        backSideEpoxideFrameS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* backSideEpoxideFrameV =
        new TGeoVolume(TString("backSideEpoxideFrameV") += TString("_") + frames->GetName(), backSideEpoxideFrameS);
 
    // volume medium
    TGeoMedium* backSideEpoxideFrameV_medium = pMedEpoxideCompound;
    if (backSideEpoxideFrameV_medium) {
        backSideEpoxideFrameV->SetMedium(backSideEpoxideFrameV_medium);
    } else
        Fatal("Main", "Invalid medium for backSideEpoxideFrameV!");
 
    // volume visual property (transparency)
    backSideEpoxideFrameV->SetLineColor(TColor::GetColor("#ccccff"));
    backSideEpoxideFrameV->SetTransparency(0);
 
    TGeoCombiTrans* backSideEpoxideFrame_transf[1];
 
    backSideEpoxideFrame_transf[0] = new TGeoCombiTrans();
    backSideEpoxideFrame_transf[0]->SetDx(0.0);
    backSideEpoxideFrame_transf[0]->SetDy(+((backSideEpoxideFramePart_YSize - YModuleSize_Station163x39) * 0.5)
                                          - frontLowerEpoxideFrame_YSize + 0.2 /*shift*/);
    backSideEpoxideFrame_transf[0]->SetDz(+(backSideEpoxideFramePart_ZSize * 0.5 + ZModuleSize_Station163x39 * 0.5));
 
    frames->AddNode(backSideEpoxideFrameV, 0, backSideEpoxideFrame_transf[0]);
    //--------------------------------------------------------------------------
 
    // back upper epoxide frame -------------------------------------------------
    Double_t backUpperEpoxideFrame_XSize = 161.42;   // cm
    Double_t backUpperEpoxideFrame_YSize = 6.2;      // cm
    Double_t backUpperEpoxideFrame_ZSize = 1.6;      // cm
 
    TGeoShape* backUpperEpoxideFrameS =
        new TGeoBBox("backUpperEpoxideFrameS", backUpperEpoxideFrame_XSize * 0.5, backUpperEpoxideFrame_YSize * 0.5,
                     backUpperEpoxideFrame_ZSize * 0.5);
 
    TGeoVolume* backUpperEpoxideFrameV =
        new TGeoVolume(TString("backUpperEpoxideFrameV") += TString("_") + frames->GetName(), backUpperEpoxideFrameS);
 
    // volume medium
    TGeoMedium* backUpperEpoxideFrameV_medium = pMedEpoxideCompound;
    if (backUpperEpoxideFrameV_medium) {
        backUpperEpoxideFrameV->SetMedium(backUpperEpoxideFrameV_medium);
    } else
        Fatal("Main", "Invalid medium for backUpperEpoxideFrameV!");
 
    // volume visual property (transparency)
    backUpperEpoxideFrameV->SetLineColor(TColor::GetColor("#ccccff"));
    backUpperEpoxideFrameV->SetTransparency(0);
 
    TGeoCombiTrans* backUpperEpoxideFrame_transf[1];
 
    backUpperEpoxideFrame_transf[0] = new TGeoCombiTrans();
    backUpperEpoxideFrame_transf[0]->SetDx(0.0);
    backUpperEpoxideFrame_transf[0]->SetDy(+(backUpperEpoxideFrame_YSize * 0.5 + YModuleSize_Station163x39 * 0.5)
                                           - 0.6 /*shift*/ + 0.2);
    backUpperEpoxideFrame_transf[0]->SetDz(+(backUpperEpoxideFrame_ZSize * 0.5 + ZModuleSize_Station163x39 * 0.5)
                                           + 0.005);
 
    frames->AddNode(backUpperEpoxideFrameV, 0, backUpperEpoxideFrame_transf[0]);
    //--------------------------------------------------------------------------
 
    // front side epoxide frame -------------------------------------------------
    Double_t frontSideEpoxideFramePart_XSize = 3.5;          // cm
    Double_t frontSideEpoxideFramePart_YSize = 49.9 - 6.0;   // cm
    Double_t frontSideEpoxideFramePart_ZSize = 1.5;          // cm
 
    Double_t frontSideEpoxideFrame_XDistBetweenParts = 161.7;   // cm
 
    TGeoShape* frontSideEpoxideFramePartS =
        new TGeoBBox(TString("frontSideEpoxideFramePartS") += TString("_") + frames->GetName(),
                     frontSideEpoxideFramePart_XSize * 0.5, frontSideEpoxideFramePart_YSize * 0.5,
                     frontSideEpoxideFramePart_ZSize * 0.5);
 
    TGeoTranslation* frontSideEpoxideFramePart_pos1 = new TGeoTranslation();
    frontSideEpoxideFramePart_pos1->SetName("frontSideEpoxideFramePart_pos1");
    frontSideEpoxideFramePart_pos1->SetDx(
        +(frontSideEpoxideFrame_XDistBetweenParts * 0.5 + frontSideEpoxideFramePart_XSize * 0.5));
    frontSideEpoxideFramePart_pos1->SetDy(0.0);
    frontSideEpoxideFramePart_pos1->SetDz(0.0);
    frontSideEpoxideFramePart_pos1->RegisterYourself();
 
    TGeoTranslation* frontSideEpoxideFramePart_pos2 = new TGeoTranslation();
    frontSideEpoxideFramePart_pos2->SetName("frontSideEpoxideFramePart_pos2");
    frontSideEpoxideFramePart_pos2->SetDx(
        -(frontSideEpoxideFrame_XDistBetweenParts * 0.5 + frontSideEpoxideFramePart_XSize * 0.5));
    frontSideEpoxideFramePart_pos2->SetDy(0.0);
    frontSideEpoxideFramePart_pos2->SetDz(0.0);
    frontSideEpoxideFramePart_pos2->RegisterYourself();
 
    TGeoCompositeShape* frontSideEpoxideFrameS = new TGeoCompositeShape();
    frontSideEpoxideFrameS->SetName(TString("frontSideEpoxideFrameS") += TString("_") + frames->GetName());
    {
        TString expression = "frontSideEpoxideFramePartS";
        expression += TString("_") + frames->GetName();
        expression += ":frontSideEpoxideFramePart_pos1";
        expression += "+frontSideEpoxideFramePartS";
        expression += TString("_") + frames->GetName();
        expression += ":frontSideEpoxideFramePart_pos2";
 
        frontSideEpoxideFrameS->MakeNode(expression);
        frontSideEpoxideFrameS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* frontSideEpoxideFrameV =
        new TGeoVolume(TString("frontSideEpoxideFrameV") += TString("_") + frames->GetName(), frontSideEpoxideFrameS);
 
    // volume medium
    TGeoMedium* frontSideEpoxideFrameV_medium = pMedEpoxideCompound;
    if (frontSideEpoxideFrameV_medium) {
        frontSideEpoxideFrameV->SetMedium(frontSideEpoxideFrameV_medium);
    } else
        Fatal("Main", "Invalid medium for frontSideEpoxideFrameV!");
 
    // volume visual property (transparency)
    frontSideEpoxideFrameV->SetLineColor(TColor::GetColor("#ccccff"));
    frontSideEpoxideFrameV->SetTransparency(0);
 
    TGeoCombiTrans* frontSideEpoxideFrame_transf[1];
 
    frontSideEpoxideFrame_transf[0] = new TGeoCombiTrans();
    frontSideEpoxideFrame_transf[0]->SetDx(0.0);
    frontSideEpoxideFrame_transf[0]->SetDy(+((frontSideEpoxideFramePart_YSize - YModuleSize_Station163x39) * 0.5)
                                           - frontLowerEpoxideFrame_YSize + 0.2 /*shift*/);
    frontSideEpoxideFrame_transf[0]->SetDz(-(frontSideEpoxideFramePart_ZSize * 0.5 + ZModuleSize_Station163x39 * 0.5)
                                           - 0.05);
 
    frames->AddNode(frontSideEpoxideFrameV, 0, frontSideEpoxideFrame_transf[0]);
    //--------------------------------------------------------------------------
 
    // front upper epoxide frame -------------------------------------------------
    Double_t frontUpperEpoxideFrame_XSize = 161.7;   // cm
    Double_t frontUpperEpoxideFrame_YSize = 3.5;     // cm
    Double_t frontUpperEpoxideFrame_ZSize = 1.6;     // cm
 
    TGeoShape* frontUpperEpoxideFrameS =
        new TGeoBBox("frontUpperEpoxideFrameS", frontUpperEpoxideFrame_XSize * 0.5, frontUpperEpoxideFrame_YSize * 0.5,
                     frontUpperEpoxideFrame_ZSize * 0.5);
 
    TGeoVolume* frontUpperEpoxideFrameV =
        new TGeoVolume(TString("frontUpperEpoxideFrameV") += TString("_") + frames->GetName(), frontUpperEpoxideFrameS);
 
    // volume medium
    TGeoMedium* frontUpperEpoxideFrameV_medium = pMedEpoxideCompound;
    if (frontUpperEpoxideFrameV_medium) {
        frontUpperEpoxideFrameV->SetMedium(frontUpperEpoxideFrameV_medium);
    } else
        Fatal("Main", "Invalid medium for frontUpperEpoxideFrameV!");
 
    // volume visual property (transparency)
    frontUpperEpoxideFrameV->SetLineColor(TColor::GetColor("#ccccff"));
    frontUpperEpoxideFrameV->SetTransparency(0);
 
    TGeoCombiTrans* frontUpperEpoxideFrame_transf[1];
 
    frontUpperEpoxideFrame_transf[0] = new TGeoCombiTrans();
    frontUpperEpoxideFrame_transf[0]->SetDx(0.0);
    frontUpperEpoxideFrame_transf[0]->SetDy(+(frontUpperEpoxideFrame_YSize * 0.5 + YModuleSize_Station163x39 * 0.5)
                                            - 0.6 /*shift*/);
    frontUpperEpoxideFrame_transf[0]->SetDz(-(frontUpperEpoxideFrame_ZSize * 0.5 + ZModuleSize_Station163x39 * 0.5)
                                            - 0.05);
 
    frames->AddNode(frontUpperEpoxideFrameV, 0, frontUpperEpoxideFrame_transf[0]);
    //--------------------------------------------------------------------------
 
    // alum. cable support ------------------------------------------------------
    Double_t alumCableSupportPart_XSize = 3.0;          // cm
    Double_t alumCableSupportPart_YSize = 41.3 - 6.0;   // cm
    Double_t alumCableSupportPart_ZSize = 0.2;          // cm
 
    Double_t alumCableSupport_XDistBetweenParts = 164.82;   // cm
 
    TGeoShape* alumCableSupportPartS = new TGeoBBox(
        TString("alumCableSupportPartS") += TString("_") + frames->GetName(), alumCableSupportPart_XSize * 0.5,
        alumCableSupportPart_YSize * 0.5, alumCableSupportPart_ZSize * 0.5);
 
    TGeoTranslation* alumCableSupportPart_pos1 = new TGeoTranslation();
    alumCableSupportPart_pos1->SetName("alumCableSupportPart_pos1");
    alumCableSupportPart_pos1->SetDx(+(alumCableSupport_XDistBetweenParts * 0.5 + alumCableSupportPart_XSize * 0.5));
    alumCableSupportPart_pos1->SetDy(0.0);
    alumCableSupportPart_pos1->SetDz(0.0);
    alumCableSupportPart_pos1->RegisterYourself();
 
    TGeoTranslation* alumCableSupportPart_pos2 = new TGeoTranslation();
    alumCableSupportPart_pos2->SetName("alumCableSupportPart_pos2");
    alumCableSupportPart_pos2->SetDx(-(alumCableSupport_XDistBetweenParts * 0.5 + alumCableSupportPart_XSize * 0.5));
    alumCableSupportPart_pos2->SetDy(0.0);
    alumCableSupportPart_pos2->SetDz(0.0);
    alumCableSupportPart_pos2->RegisterYourself();
 
    TGeoCompositeShape* alumCableSupportS = new TGeoCompositeShape();
    alumCableSupportS->SetName(TString("alumCableSupportS") += TString("_") + frames->GetName());
    {
        TString expression = "alumCableSupportPartS";
        expression += TString("_") + frames->GetName();
        expression += ":alumCableSupportPart_pos1";
        expression += "+alumCableSupportPartS";
        expression += TString("_") + frames->GetName();
        expression += ":alumCableSupportPart_pos2";
 
        alumCableSupportS->MakeNode(expression);
        alumCableSupportS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* alumCableSupportV =
        new TGeoVolume(TString("alumCableSupportV") += TString("_") + frames->GetName(), alumCableSupportS);
 
    // volume medium
    TGeoMedium* alumCableSupportV_medium = pMedAluminium;
    if (alumCableSupportV_medium) {
        alumCableSupportV->SetMedium(alumCableSupportV_medium);
    } else
        Fatal("Main", "Invalid medium for alumCableSupportV!");
 
    // volume visual property (transparency)
    alumCableSupportV->SetLineColor(TColor::GetColor("#ffcccc"));
    alumCableSupportV->SetTransparency(0);
 
    TGeoCombiTrans* alumCableSupport_transf[1];
 
    alumCableSupport_transf[0] = new TGeoCombiTrans();
    alumCableSupport_transf[0]->SetDx(0.0);
    alumCableSupport_transf[0]->SetDy(+((alumCableSupportPart_YSize - YModuleSize_Station163x39) * 0.5)
                                      - frontLowerEpoxideFrame_YSize + 0.2 /*shift*/);
    alumCableSupport_transf[0]->SetDz(-(alumCableSupportPart_ZSize * 0.5 + ZModuleSize_Station163x39 * 0.5
                                        + frontSideEpoxideFramePart_ZSize + 1.0 /*shift*/));
 
    frames->AddNode(alumCableSupportV, 0, alumCableSupport_transf[0]);
    //--------------------------------------------------------------------------
 
    // amplifier upper block ----------------------------------------------------
    Double_t amplifierUpperBlockPart1_XSize = 161.44;   // cm
    Double_t amplifierUpperBlockPart1_YSize = 0.6;      // cm
    Double_t amplifierUpperBlockPart1_ZSize = 8.41;     // cm
 
    Double_t amplifierUpperBlockPart2_XSize = 161.16;   // cm
    Double_t amplifierUpperBlockPart2_YSize = 0.6;      // cm
    Double_t amplifierUpperBlockPart2_ZSize = 8.41;     // cm
 
    TGeoShape* amplifierUpperBlockPart1S = new TGeoBBox(
        TString("amplifierUpperBlockPart1S") += TString("_") + frames->GetName(), amplifierUpperBlockPart1_XSize * 0.5,
        amplifierUpperBlockPart1_YSize * 0.5, amplifierUpperBlockPart1_ZSize * 0.5);
    TGeoShape* amplifierUpperBlockPart2S = new TGeoBBox(
        TString("amplifierUpperBlockPart2S") += TString("_") + frames->GetName(), amplifierUpperBlockPart2_XSize * 0.5,
        amplifierUpperBlockPart2_YSize * 0.5, amplifierUpperBlockPart2_ZSize * 0.5);
 
    TGeoTranslation* amplifierUpperBlockPart1_pos = new TGeoTranslation();
    amplifierUpperBlockPart1_pos->SetName("amplifierUpperBlockPart1_pos");
    amplifierUpperBlockPart1_pos->SetDx(0.0);
    amplifierUpperBlockPart1_pos->SetDy(0.0);
    amplifierUpperBlockPart1_pos->SetDz(0.0);
    amplifierUpperBlockPart1_pos->RegisterYourself();
 
    TGeoTranslation* amplifierUpperBlockPart2_pos = new TGeoTranslation();
    amplifierUpperBlockPart2_pos->SetName("amplifierUpperBlockPart2_pos");
    amplifierUpperBlockPart2_pos->SetDx(-0.76);
    amplifierUpperBlockPart2_pos->SetDy(
        -(amplifierUpperBlockPart1_YSize * 0.5 + amplifierUpperBlockPart2_YSize * 0.5 + 0.33));
    amplifierUpperBlockPart2_pos->SetDz(0.0);
    amplifierUpperBlockPart2_pos->RegisterYourself();
 
    TGeoCompositeShape* amplifierUpperBlockS = new TGeoCompositeShape();
    amplifierUpperBlockS->SetName(TString("amplifierUpperBlockS") += TString("_") + frames->GetName());
    {
        TString expression = "amplifierUpperBlockPart1S";
        expression += TString("_") + frames->GetName();
        expression += ":amplifierUpperBlockPart1_pos";
        expression += "+amplifierUpperBlockPart1S";
        expression += TString("_") + frames->GetName();
        expression += ":amplifierUpperBlockPart2_pos";
 
        amplifierUpperBlockS->MakeNode(expression);
        amplifierUpperBlockS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* amplifierUpperBlockV =
        new TGeoVolume(TString("amplifierUpperBlockV") += TString("_") + frames->GetName(), amplifierUpperBlockS);
 
    // volume medium
    TGeoMedium* amplifierUpperBlockV_medium = pMedFiberGlass;
    if (amplifierUpperBlockV_medium) {
        amplifierUpperBlockV->SetMedium(amplifierUpperBlockV_medium);
    } else
        Fatal("Main", "Invalid medium for amplifierUpperBlockV!");
 
    // volume visual property (transparency)
    amplifierUpperBlockV->SetLineColor(TColor::GetColor("#ccffcc"));
    amplifierUpperBlockV->SetTransparency(0);
 
    TGeoCombiTrans* amplifierUpperBlock_transf[1];
 
    amplifierUpperBlock_transf[0] = new TGeoCombiTrans();
    amplifierUpperBlock_transf[0]->SetDx(0.0);
    amplifierUpperBlock_transf[0]->SetDy(+(YModuleSize_Station163x39 * 0.5 + backUpperEpoxideFrame_YSize)
                                         - 1.0 /*shift*/);
    amplifierUpperBlock_transf[0]->SetDz(-(amplifierUpperBlockPart1_ZSize * 0.5) + ZModuleSize_Station163x39 * 0.5);
 
    frames->AddNode(amplifierUpperBlockV, 0, amplifierUpperBlock_transf[0]);
    //--------------------------------------------------------------------------
 
    // amplifier side block -----------------------------------------------------
    Double_t amplifierSideBlockPartBig_XSize = 0.6;    // cm
    Double_t amplifierSideBlockPartBig_YSize = 33.9;   // cm
    Double_t amplifierSideBlockPartBig_ZSize = 8.41;   // cm
 
    Double_t amplifierSideBlockPartSmall_XSize = 0.6;     // cm
    Double_t amplifierSideBlockPartSmall_YSize = 27.05;   // cm
    Double_t amplifierSideBlockPartSmall_ZSize = 8.41;    // cm
 
    TGeoShape* amplifierSideBlockPartBigS =
        new TGeoBBox(TString("amplifierSideBlockPartBigS") += TString("_") + frames->GetName(),
                     amplifierSideBlockPartBig_XSize * 0.5, amplifierSideBlockPartBig_YSize * 0.5,
                     amplifierSideBlockPartBig_ZSize * 0.5);
    TGeoShape* amplifierSideBlockPartSmallS =
        new TGeoBBox(TString("amplifierSideBlockPartSmallS") += TString("_") + frames->GetName(),
                     amplifierSideBlockPartSmall_XSize * 0.5, amplifierSideBlockPartSmall_YSize * 0.5,
                     amplifierSideBlockPartSmall_ZSize * 0.5);
 
    TGeoTranslation* amplifierSideBlockPartBig_pos = new TGeoTranslation();
    amplifierSideBlockPartBig_pos->SetName("amplifierSideBlockPartBig_pos");
    amplifierSideBlockPartBig_pos->SetDx(0.0);
    amplifierSideBlockPartBig_pos->SetDy(0.0);
    amplifierSideBlockPartBig_pos->SetDz(0.0);
    amplifierSideBlockPartBig_pos->RegisterYourself();
 
    TGeoTranslation* amplifierSideBlockPartSmall_pos = new TGeoTranslation();
    amplifierSideBlockPartSmall_pos->SetName("amplifierSideBlockPartSmall_pos");
    amplifierSideBlockPartSmall_pos->SetDx(
        +(amplifierSideBlockPartBig_XSize * 0.5 + amplifierSideBlockPartSmall_XSize * 0.5 + 1.72));
    amplifierSideBlockPartSmall_pos->SetDy(+0.05);
    amplifierSideBlockPartSmall_pos->SetDz(0.0);
    amplifierSideBlockPartSmall_pos->RegisterYourself();
 
    TGeoCompositeShape* amplifierSideBlockS = new TGeoCompositeShape();
    amplifierSideBlockS->SetName(TString("amplifierSideBlockS") += TString("_") + frames->GetName());
    {
        TString expression = "amplifierSideBlockPartBigS";
        expression += TString("_") + frames->GetName();
        expression += ":amplifierSideBlockPartBig_pos";
        expression += "+amplifierSideBlockPartSmallS";
        expression += TString("_") + frames->GetName();
        expression += ":amplifierSideBlockPartSmall_pos";
 
        amplifierSideBlockS->MakeNode(expression);
        amplifierSideBlockS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* amplifierSideBlockV =
        new TGeoVolume(TString("amplifierSideBlockV") += TString("_") + frames->GetName(), amplifierSideBlockS);
 
    // volume medium
    TGeoMedium* amplifierSideBlockV_medium = pMedFiberGlass;
    if (amplifierSideBlockV_medium) {
        amplifierSideBlockV->SetMedium(amplifierSideBlockV_medium);
    } else
        Fatal("Main", "Invalid medium for amplifierSideBlockV!");
 
    // volume visual property (transparency)
    amplifierSideBlockV->SetLineColor(TColor::GetColor("#ccffcc"));
    amplifierSideBlockV->SetTransparency(0);
 
    TGeoCombiTrans* amplifierSideBlock_transf[2];
 
    amplifierSideBlock_transf[0] = new TGeoCombiTrans();
    amplifierSideBlock_transf[0]->SetDx(+86.47 /*shift*/);
    amplifierSideBlock_transf[0]->SetDy(-5.5 + 6.0 * 0.5);
    amplifierSideBlock_transf[0]->SetDz(-(amplifierSideBlockPartBig_ZSize * 0.5) + ZModuleSize_Station163x39 * 0.5
                                        - 0.4 /*shift*/);
 
    amplifierSideBlock_transf[1] = new TGeoCombiTrans();
    amplifierSideBlock_transf[1]->RotateY(180.0);
    amplifierSideBlock_transf[1]->SetDx(-86.47 /*shift*/);
    amplifierSideBlock_transf[1]->SetDy(-5.5 + 6.0 * 0.5);
    amplifierSideBlock_transf[1]->SetDz(-(amplifierSideBlockPartBig_ZSize * 0.5) + ZModuleSize_Station163x39 * 0.5
                                        - 0.4 /*shift*/);
 
    frames->AddNode(amplifierSideBlockV, 0, amplifierSideBlock_transf[0]);
    frames->AddNode(amplifierSideBlockV, 1, amplifierSideBlock_transf[1]);
    //--------------------------------------------------------------------------
 
    // upper PCB (Printed Circuit Board) ----------------------------------------
    Double_t upperPcb_XSize = 177.22;   // cm
    Double_t upperPcb_YSize = 0.3;      // cm
    Double_t upperPcb_ZSize = 4.0;      // cm
 
    TGeoShape* upperPcbS = new TGeoBBox("upperPcbS", upperPcb_XSize * 0.5, upperPcb_YSize * 0.5, upperPcb_ZSize * 0.5);
 
    TGeoVolume* upperPcbV = new TGeoVolume(TString("upperPcbV") += TString("_") + frames->GetName(), upperPcbS);
 
    // volume medium
    TGeoMedium* upperPcbV_medium = pMedFiberGlass;
    if (upperPcbV_medium) {
        upperPcbV->SetMedium(upperPcbV_medium);
    } else
        Fatal("Main", "Invalid medium for upperPcbV!");
 
    // volume visual property (transparency)
    upperPcbV->SetLineColor(TColor::GetColor("#88ff88"));
    upperPcbV->SetTransparency(0);
 
    TGeoCombiTrans* upperPcb_transf[1];
 
    upperPcb_transf[0] = new TGeoCombiTrans();
    upperPcb_transf[0]->SetDx(+0.1 /*shift*/);
    upperPcb_transf[0]->SetDy(+(upperPcb_YSize * 0.5 + 23.5 /*shift*/) - 6.0 * 0.5);
    upperPcb_transf[0]->SetDz(-(upperPcb_ZSize * 0.5) + ZModuleSize_Station163x39 * 0.5 - ZModuleSize_Station163x39
                              - frontSideEpoxideFramePart_ZSize - 0.2);
 
    frames->AddNode(upperPcbV, 0, upperPcb_transf[0]);
    //--------------------------------------------------------------------------
 
    // side PCB -----------------------------------------------------------------
    Double_t sidePcbPart_XSize = 0.3;     // cm
    Double_t sidePcbPart_YSize = 35.05;   // cm
    Double_t sidePcbPart_ZSize = 4.0;     // cm
 
    TGeoShape* sidePcbPartS = new TGeoBBox(TString("sidePcbPartS") += TString("_") + frames->GetName(),
                                           sidePcbPart_XSize * 0.5, sidePcbPart_YSize * 0.5, sidePcbPart_ZSize * 0.5);
 
    TGeoTranslation* sidePcbPart_pos1 = new TGeoTranslation();
    sidePcbPart_pos1->SetName("sidePcbPart_pos1");
    sidePcbPart_pos1->SetDx(0.0);
    sidePcbPart_pos1->SetDy(0.0);
    sidePcbPart_pos1->SetDz(0.0);
    sidePcbPart_pos1->RegisterYourself();
 
    TGeoTranslation* sidePcbPart_pos2 = new TGeoTranslation();
    sidePcbPart_pos2->SetName("sidePcbPart_pos2");
    sidePcbPart_pos2->SetDx(+(sidePcbPart_XSize + 6.1 /*shift*/));
    sidePcbPart_pos2->SetDy(0.0);
    sidePcbPart_pos2->SetDz(0.0);
    sidePcbPart_pos2->RegisterYourself();
 
    TGeoCompositeShape* sidePcbS = new TGeoCompositeShape();
    sidePcbS->SetName(TString("sidePcbS") += TString("_") + frames->GetName());
    {
        TString expression = "sidePcbPartS";
        expression += TString("_") + frames->GetName();
        expression += ":sidePcbPart_pos1";
        expression += "+sidePcbPartS";
        expression += TString("_") + frames->GetName();
        expression += ":sidePcbPart_pos2";
 
        sidePcbS->MakeNode(expression);
        sidePcbS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* sidePcbV = new TGeoVolume(TString("sidePcbV") += TString("_") + frames->GetName(), sidePcbS);
 
    // volume medium
    TGeoMedium* sidePcbV_medium = pMedFiberGlass;
    if (sidePcbV_medium) {
        sidePcbV->SetMedium(sidePcbV_medium);
    } else
        Fatal("Main", "Invalid medium for sidePcbV!");
 
    // volume visual property (transparency)
    sidePcbV->SetLineColor(TColor::GetColor("#88ff88"));
    sidePcbV->SetTransparency(0);
 
    TGeoCombiTrans* sidePcb_transf[2];
 
    sidePcb_transf[0] = new TGeoCombiTrans();
    sidePcb_transf[0]->SetDx(+(sidePcbPart_XSize * 0.5 + 84.1 /*shift*/));
    sidePcb_transf[0]->SetDy(-6.8 /*shift*/ + 6.0 * 0.5);
    sidePcb_transf[0]->SetDz(-(sidePcbPart_ZSize * 0.5 + alumCableSupportPart_ZSize + ZModuleSize_Station163x39 * 0.5
                               + frontSideEpoxideFramePart_ZSize + 1.0 /*shift*/));
 
    sidePcb_transf[1] = new TGeoCombiTrans();
    sidePcb_transf[1]->RotateY(180.0);
    sidePcb_transf[1]->SetDx(-(sidePcbPart_XSize * 0.5 + 84.1 /*shift*/));
    sidePcb_transf[1]->SetDy(-6.8 /*shift*/ + 6.0 * 0.5);
    sidePcb_transf[1]->SetDz(-(sidePcbPart_ZSize * 0.5 + alumCableSupportPart_ZSize + ZModuleSize_Station163x39 * 0.5
                               + frontSideEpoxideFramePart_ZSize + 1.0 /*shift*/));
 
    frames->AddNode(sidePcbV, 0, sidePcb_transf[0]);
    frames->AddNode(sidePcbV, 1, sidePcb_transf[1]);
    //--------------------------------------------------------------------------
 
    // alum. side support -------------------------------------------------------
    Double_t alumSideSupport_XSize = 4.8;          // cm
    Double_t alumSideSupport_YSize = 52.6 - 6.0;   // cm
    Double_t alumSideSupport_ZSize = 0.4;          // cm
 
    TGeoShape* alumSideSupportS = new TGeoBBox("alumSideSupportS", alumSideSupport_XSize * 0.5,
                                               alumSideSupport_YSize * 0.5, alumSideSupport_ZSize * 0.5);
 
    TGeoVolume* alumSideSupportV =
        new TGeoVolume(TString("alumSideSupportV") += TString("_") + frames->GetName(), alumSideSupportS);
 
    // volume medium
    TGeoMedium* alumSideSupportV_medium = pMedAluminium;
    if (alumSideSupportV_medium) {
        alumSideSupportV->SetMedium(alumSideSupportV_medium);
    } else
        Fatal("Main", "Invalid medium for alumSideSupportV!");
 
    // volume visual property (transparency)
    alumSideSupportV->SetLineColor(TColor::GetColor("#ffcccc"));
    alumSideSupportV->SetTransparency(0);
 
    TGeoCombiTrans* alumSideSupport_transf[2];
 
    alumSideSupport_transf[0] = new TGeoCombiTrans();
    alumSideSupport_transf[0]->SetDx(+(alumSideSupport_XSize * 0.5 + 88.31 /*shift*/));
    alumSideSupport_transf[0]->SetDy(+1.8);
    alumSideSupport_transf[0]->SetDz(+(ZModuleSize_Station163x39 * 0.5 - alumSideSupport_ZSize * 0.5));
 
    alumSideSupport_transf[1] = new TGeoCombiTrans();
    alumSideSupport_transf[1]->SetDx(-(alumSideSupport_XSize * 0.5 + 88.31 /*shift*/));
    alumSideSupport_transf[1]->SetDy(+1.8 /*shift*/);
    alumSideSupport_transf[1]->SetDz(+(ZModuleSize_Station163x39 * 0.5 - alumSideSupport_ZSize * 0.5));
 
    frames->AddNode(alumSideSupportV, 0, alumSideSupport_transf[0]);
    frames->AddNode(alumSideSupportV, 1, alumSideSupport_transf[1]);
    //--------------------------------------------------------------------------
 
    // bosch shape --------------------------------------------------------------
    Double_t boschShape_XSize = 2.25;         // cm
    Double_t boschShape_YSize = 52.6 - 6.0;   // cm
    Double_t boschShape_ZSize = 2.25;         // cm
 
    TGeoShape* boschShapeS =
        new TGeoBBox("boschShapeS", boschShape_XSize * 0.5, boschShape_YSize * 0.5, boschShape_ZSize * 0.5);
 
    TGeoVolume* boschShapeV = new TGeoVolume(TString("boschShapeV") += TString("_") + frames->GetName(), boschShapeS);
 
    // volume medium
    TGeoMedium* boschShapeV_medium = pMedAluminium;
    if (boschShapeV_medium) {
        boschShapeV->SetMedium(boschShapeV_medium);
    } else
        Fatal("Main", "Invalid medium for boschShapeV!");
 
    // volume visual property (transparency)
    boschShapeV->SetLineColor(TColor::GetColor("#ffcccc"));
    boschShapeV->SetTransparency(0);
 
    TGeoCombiTrans* boschShape_transf[2];
 
    boschShape_transf[0] = new TGeoCombiTrans();
    boschShape_transf[0]->SetDx(+(boschShape_XSize * 0.5 + 90.86 /*shift*/));
    boschShape_transf[0]->SetDy(+1.8 /*shift*/);
    boschShape_transf[0]->SetDz(+(ZModuleSize_Station163x39 * 0.5 + boschShape_ZSize * 0.5));
 
    boschShape_transf[1] = new TGeoCombiTrans();
    boschShape_transf[1]->SetDx(-(boschShape_XSize * 0.5 + 90.86 /*shift*/));
    boschShape_transf[1]->SetDy(+1.8 /*shift*/);
    boschShape_transf[1]->SetDz(+(ZModuleSize_Station163x39 * 0.5 + boschShape_ZSize * 0.5));
 
    // frames->AddNode(boschShapeV, 0, boschShape_transf[0]);
    // frames->AddNode(boschShapeV, 1, boschShape_transf[1]);
    //--------------------------------------------------------------------------
 
    // polyacetal shape ---------------------------------------------------------
    Double_t insidePolyacetalShapeUpper_XSize = 11.46;                // cm
    Double_t insidePolyacetalShapeUpper_YSize = 52.6;                 // cm
    Double_t insidePolyacetalShapeUpper_ZSize = 3.0;                  // cm
    Double_t insidePolyacetalShapeUpper_XDistBetweenParts = 186.22;   // cm
 
    Double_t outsidePolyacetalShapeUpper_XSize = 23.5;                 // cm
    Double_t outsidePolyacetalShapeUpper_YSize = 5.405;                // cm
    Double_t outsidePolyacetalShapeUpper_ZSize = 3.0;                  // cm
    Double_t outsidePolyacetalShapeUpper_XDistBetweenParts = 209.14;   // cm
 
    Double_t insidePolyacetalShapeLower_XSize = 8.4;                  // cm
    Double_t insidePolyacetalShapeLower_YSize = 52.6 - 6.0;           // cm
    Double_t insidePolyacetalShapeLower_ZSize = 3.0;                  // cm
    Double_t insidePolyacetalShapeLower_XDistBetweenParts = 186.22;   // cm
 
    Double_t outsidePolyacetalShapeLower_XSize = 26.7;                 // cm
    Double_t outsidePolyacetalShapeLower_YSize = 24.4;                 // cm
    Double_t outsidePolyacetalShapeLower_ZSize = 3.0;                  // cm
    Double_t outsidePolyacetalShapeLower_XDistBetweenParts = 203.02;   // cm
 
    TGeoShape* insidePolyacetalShapeUpperS =
        new TGeoBBox("insidePolyacetalShapeUpperS", insidePolyacetalShapeUpper_XSize * 0.5,
                     insidePolyacetalShapeUpper_YSize * 0.5, insidePolyacetalShapeUpper_ZSize * 0.5);
    TGeoShape* outsidePolyacetalShapeUpperS =
        new TGeoBBox("outsidePolyacetalShapeUpperS", outsidePolyacetalShapeUpper_XSize * 0.5,
                     outsidePolyacetalShapeUpper_YSize * 0.5, outsidePolyacetalShapeUpper_ZSize * 0.5);
    TGeoShape* insidePolyacetalShapeLowerS =
        new TGeoBBox("insidePolyacetalShapeLowerS", insidePolyacetalShapeLower_XSize * 0.5,
                     insidePolyacetalShapeLower_YSize * 0.5, insidePolyacetalShapeLower_ZSize * 0.5);
    TGeoShape* outsidePolyacetalShapeLowerS =
        new TGeoBBox("outsidePolyacetalShapeLowerS", outsidePolyacetalShapeLower_XSize * 0.5,
                     outsidePolyacetalShapeLower_YSize * 0.5, outsidePolyacetalShapeLower_ZSize * 0.5);
 
    TGeoVolume* insidePolyacetalShapeUpperV = new TGeoVolume(
        TString("insidePolyacetalShapeUpperV") += TString("_") + frames->GetName(), insidePolyacetalShapeUpperS);
    TGeoVolume* outsidePolyacetalShapeUpperV = new TGeoVolume(
        TString("outsidePolyacetalShapeUpperV") += TString("_") + frames->GetName(), outsidePolyacetalShapeUpperS);
    TGeoVolume* insidePolyacetalShapeLowerV = new TGeoVolume(
        TString("insidePolyacetalShapeLowerV") += TString("_") + frames->GetName(), insidePolyacetalShapeLowerS);
    TGeoVolume* outsidePolyacetalShapeLowerV = new TGeoVolume(
        TString("outsidePolyacetalShapeLowerV") += TString("_") + frames->GetName(), outsidePolyacetalShapeLowerS);
 
    // volume medium
    TGeoMedium* polyacetalShapeV_medium = pMedPolyacetal;
    if (polyacetalShapeV_medium) {
        insidePolyacetalShapeUpperV->SetMedium(polyacetalShapeV_medium);
        outsidePolyacetalShapeUpperV->SetMedium(polyacetalShapeV_medium);
        insidePolyacetalShapeLowerV->SetMedium(polyacetalShapeV_medium);
        outsidePolyacetalShapeLowerV->SetMedium(polyacetalShapeV_medium);
    } else
        Fatal("Main", "Invalid medium for insidePolyacetalShapeUpperV!");
 
    // volume visual property (transparency)
    insidePolyacetalShapeUpperV->SetLineColor(TColor::GetColor("#ffcccc"));
    insidePolyacetalShapeUpperV->SetTransparency(0);
    outsidePolyacetalShapeUpperV->SetLineColor(TColor::GetColor("#ffcccc"));
    outsidePolyacetalShapeUpperV->SetTransparency(0);
    insidePolyacetalShapeLowerV->SetLineColor(TColor::GetColor("#ffcccc"));
    insidePolyacetalShapeLowerV->SetTransparency(0);
    outsidePolyacetalShapeLowerV->SetLineColor(TColor::GetColor("#ffcccc"));
    outsidePolyacetalShapeLowerV->SetTransparency(0);
 
    TGeoCombiTrans* insidePolyacetalShapeUpper_transf[2];
 
    insidePolyacetalShapeUpper_transf[0] = new TGeoCombiTrans();
    insidePolyacetalShapeUpper_transf[0]->SetDx(
        +(insidePolyacetalShapeUpper_XSize * 0.5 + insidePolyacetalShapeUpper_XDistBetweenParts * 0.5));
    insidePolyacetalShapeUpper_transf[0]->SetDy(+1.8 /*shift*/);
    insidePolyacetalShapeUpper_transf[0]->SetDz(0.0);
 
    insidePolyacetalShapeUpper_transf[1] = new TGeoCombiTrans();
    insidePolyacetalShapeUpper_transf[1]->SetDx(
        -(insidePolyacetalShapeUpper_XSize * 0.5 + insidePolyacetalShapeUpper_XDistBetweenParts * 0.5));
    insidePolyacetalShapeUpper_transf[1]->SetDy(+1.8 /*shift*/);
    insidePolyacetalShapeUpper_transf[1]->SetDz(0.0);
 
    TGeoCombiTrans* outsidePolyacetalShapeUpper_transf[2];
 
    outsidePolyacetalShapeUpper_transf[0] = new TGeoCombiTrans();
    outsidePolyacetalShapeUpper_transf[0]->SetDx(
        +(outsidePolyacetalShapeUpper_XSize * 0.5 + outsidePolyacetalShapeUpper_XDistBetweenParts * 0.5));
    outsidePolyacetalShapeUpper_transf[0]->SetDy(-insidePolyacetalShapeUpper_YSize * 0.5
                                                 + 1.8 /*shift*/ + outsidePolyacetalShapeUpper_YSize * 0.5);
    outsidePolyacetalShapeUpper_transf[0]->SetDz(0.0);
 
    outsidePolyacetalShapeUpper_transf[1] = new TGeoCombiTrans();
    outsidePolyacetalShapeUpper_transf[1]->SetDx(
        -(outsidePolyacetalShapeUpper_XSize * 0.5 + outsidePolyacetalShapeUpper_XDistBetweenParts * 0.5));
    outsidePolyacetalShapeUpper_transf[1]->SetDy(-insidePolyacetalShapeUpper_YSize * 0.5
                                                 + 1.8 /*shift*/ + outsidePolyacetalShapeUpper_YSize * 0.5);
    outsidePolyacetalShapeUpper_transf[1]->SetDz(0.0);
 
    TGeoCombiTrans* insidePolyacetalShapeLower_transf[2];
 
    insidePolyacetalShapeLower_transf[0] = new TGeoCombiTrans();
    insidePolyacetalShapeLower_transf[0]->SetDx(
        +(insidePolyacetalShapeLower_XSize * 0.5 + insidePolyacetalShapeLower_XDistBetweenParts * 0.5));
    insidePolyacetalShapeLower_transf[0]->SetDy(+1.8 /*shift*/);
    insidePolyacetalShapeLower_transf[0]->SetDz(0.0);
 
    insidePolyacetalShapeLower_transf[1] = new TGeoCombiTrans();
    insidePolyacetalShapeLower_transf[1]->SetDx(
        -(insidePolyacetalShapeLower_XSize * 0.5 + insidePolyacetalShapeLower_XDistBetweenParts * 0.5));
    insidePolyacetalShapeLower_transf[1]->SetDy(+1.8 /*shift*/);
    insidePolyacetalShapeLower_transf[1]->SetDz(0.0);
 
    TGeoCombiTrans* outsidePolyacetalShapeLower_transf[2];
 
    outsidePolyacetalShapeLower_transf[0] = new TGeoCombiTrans();
    outsidePolyacetalShapeLower_transf[0]->SetDx(
        +(outsidePolyacetalShapeLower_XSize * 0.5 + outsidePolyacetalShapeLower_XDistBetweenParts * 0.5));
    outsidePolyacetalShapeLower_transf[0]->SetDy(-insidePolyacetalShapeLower_YSize * 0.5
                                                 + 1.8 /*shift*/ + outsidePolyacetalShapeLower_YSize * 0.5);
    outsidePolyacetalShapeLower_transf[0]->SetDz(0.0);
 
    outsidePolyacetalShapeLower_transf[1] = new TGeoCombiTrans();
    outsidePolyacetalShapeLower_transf[1]->SetDx(
        -(outsidePolyacetalShapeLower_XSize * 0.5 + outsidePolyacetalShapeLower_XDistBetweenParts * 0.5));
    outsidePolyacetalShapeLower_transf[1]->SetDy(-insidePolyacetalShapeLower_YSize * 0.5
                                                 + 1.8 /*shift*/ + outsidePolyacetalShapeLower_YSize * 0.5);
    outsidePolyacetalShapeLower_transf[1]->SetDz(0.0);
 
    if (TString(frames->GetName()).Contains("upper")) {
        frames->AddNode(insidePolyacetalShapeUpperV, 0, insidePolyacetalShapeUpper_transf[0]);
        frames->AddNode(insidePolyacetalShapeUpperV, 1, insidePolyacetalShapeUpper_transf[1]);
        frames->AddNode(outsidePolyacetalShapeUpperV, 0, outsidePolyacetalShapeUpper_transf[0]);
        frames->AddNode(outsidePolyacetalShapeUpperV, 1, outsidePolyacetalShapeUpper_transf[1]);
    }
 
    if (TString(frames->GetName()).Contains("lower")) {
        frames->AddNode(insidePolyacetalShapeLowerV, 0, insidePolyacetalShapeLower_transf[0]);
        frames->AddNode(insidePolyacetalShapeLowerV, 1, insidePolyacetalShapeLower_transf[1]);
        frames->AddNode(outsidePolyacetalShapeLowerV, 0, outsidePolyacetalShapeLower_transf[0]);
        frames->AddNode(outsidePolyacetalShapeLowerV, 1, outsidePolyacetalShapeLower_transf[1]);
    }
 
    //--------------------------------------------------------------------------
 
    // lower fiberglass block ---------------------------------------------------
    Double_t lowerFiberGlassBlock_XSize = 165.1;   // cm
    Double_t lowerFiberGlassBlock_YSize = 1.75;    // cm
    Double_t lowerFiberGlassBlock_ZSize = 0.9;     // cm
 
    Double_t lowerFiberGlassBlockArch_Rmin = 4.0;
    Double_t lowerFiberGlassBlockArch_Rmax = 5.75;
    Double_t lowerFiberGlassBlockArch_ZSize = lowerFiberGlassBlock_ZSize;
 
    TGeoShape* lowerFiberGlassBlockHalfS =
        new TGeoBBox(TString("lowerFiberGlassBlockHalfS") += TString("_") + frames->GetName(),
                     (lowerFiberGlassBlock_XSize * 0.5 - lowerFiberGlassBlockArch_Rmin) * 0.5,
                     lowerFiberGlassBlock_YSize * 0.5, lowerFiberGlassBlock_ZSize * 0.5);
    TGeoShape* lowerFiberGlassBlockArchS = new TGeoTubeSeg(
        TString("lowerFiberGlassBlockArchS") += TString("_") + frames->GetName(), lowerFiberGlassBlockArch_Rmin,
        lowerFiberGlassBlockArch_Rmax, lowerFiberGlassBlockArch_ZSize * 0.5, 0.0, 180.0);
 
    TGeoTranslation* lowerFiberGlassBlockHalf_pos1 = new TGeoTranslation();
    lowerFiberGlassBlockHalf_pos1->SetName("lowerFiberGlassBlockHalf_pos1");
    lowerFiberGlassBlockHalf_pos1->SetDx(
        +((lowerFiberGlassBlock_XSize * 0.5 - lowerFiberGlassBlockArch_Rmin) * 0.5 + lowerFiberGlassBlockArch_Rmin));
    lowerFiberGlassBlockHalf_pos1->SetDy(-(lowerFiberGlassBlock_YSize * 0.5));
    lowerFiberGlassBlockHalf_pos1->SetDz(0.0);
    lowerFiberGlassBlockHalf_pos1->RegisterYourself();
 
    TGeoTranslation* lowerFiberGlassBlockHalf_pos2 = new TGeoTranslation();
    lowerFiberGlassBlockHalf_pos2->SetName("lowerFiberGlassBlockHalf_pos2");
    lowerFiberGlassBlockHalf_pos2->SetDx(
        -((lowerFiberGlassBlock_XSize * 0.5 - lowerFiberGlassBlockArch_Rmin) * 0.5 + lowerFiberGlassBlockArch_Rmin));
    lowerFiberGlassBlockHalf_pos2->SetDy(-(lowerFiberGlassBlock_YSize * 0.5));
    lowerFiberGlassBlockHalf_pos2->SetDz(0.0);
    lowerFiberGlassBlockHalf_pos2->RegisterYourself();
 
    TGeoTranslation* lowerFiberGlassBlockArch_pos = new TGeoTranslation();
    lowerFiberGlassBlockArch_pos->SetName("lowerFiberGlassBlockArch_pos");
    lowerFiberGlassBlockArch_pos->SetDx(0.0);
    lowerFiberGlassBlockArch_pos->SetDy(0.0);
    lowerFiberGlassBlockArch_pos->SetDz(0.0);
    lowerFiberGlassBlockArch_pos->RegisterYourself();
 
    TGeoCompositeShape* lowerFiberGlassBlockS = new TGeoCompositeShape();
    lowerFiberGlassBlockS->SetName(TString("lowerFiberGlassBlockS") += TString("_") + frames->GetName());
    {
        TString expression = "lowerFiberGlassBlockHalfS";
        expression += TString("_") + frames->GetName();
        expression += ":lowerFiberGlassBlockHalf_pos1";
        expression += "+lowerFiberGlassBlockHalfS";
        expression += TString("_") + frames->GetName();
        expression += ":lowerFiberGlassBlockHalf_pos2";
        expression += "+lowerFiberGlassBlockArchS";
        expression += TString("_") + frames->GetName();
        expression += ":lowerFiberGlassBlockArch_pos";
 
        lowerFiberGlassBlockS->MakeNode(expression);
        lowerFiberGlassBlockS->ComputeBBox();   // need to compute a bounding box
    }
 
    TGeoVolume* lowerFiberGlassBlockV =
        new TGeoVolume(TString("lowerFiberGlassBlockV") += TString("_") + frames->GetName(), lowerFiberGlassBlockS);
 
    // volume medium
    TGeoMedium* lowerFiberGlassBlockV_medium = pMedFiberGlass;
    if (lowerFiberGlassBlockV_medium) {
        lowerFiberGlassBlockV->SetMedium(lowerFiberGlassBlockV_medium);
    } else
        Fatal("Main", "Invalid medium for lowerFiberGlassBlockV!");
 
    // volume visual property (transparency)
    lowerFiberGlassBlockV->SetLineColor(TColor::GetColor("#88ff88"));
    lowerFiberGlassBlockV->SetTransparency(0);
 
    TGeoCombiTrans* lowerFiberGlassBlock_transf[1];
 
    lowerFiberGlassBlock_transf[0] = new TGeoCombiTrans();
    lowerFiberGlassBlock_transf[0]->SetDx(+0.0);
    lowerFiberGlassBlock_transf[0]->SetDy(-(YModuleSize_Station163x39 * 0.5) - 0.25 - 0.1);
    lowerFiberGlassBlock_transf[0]->SetDz(0.0);
 
    frames->AddNode(lowerFiberGlassBlockV, 0, lowerFiberGlassBlock_transf[0]);
    //--------------------------------------------------------------------------
 
    // side fiberglass block ----------------------------------------------------
    Double_t sideFiberGlassBlockPart_XSize = 1.45;         // cm
    Double_t sideFiberGlassBlockPart_YSize = 49.4 - 6.0;   // cm
    Double_t sideFiberGlassBlockPart_ZSize = 0.89;         // cm
 
    TGeoShape* sideFiberGlassBlockS =
        new TGeoBBox("sideFiberGlassBlockS", sideFiberGlassBlockPart_XSize * 0.5, sideFiberGlassBlockPart_YSize * 0.5,
                     sideFiberGlassBlockPart_ZSize * 0.5);
 
    TGeoVolume* sideFiberGlassBlockV =
        new TGeoVolume(TString("sideFiberGlassBlockV") += TString("_") + frames->GetName(), sideFiberGlassBlockS);
 
    // volume medium
    TGeoMedium* sideFiberGlassBlockV_medium = pMedFiberGlass;
    if (sideFiberGlassBlockV_medium) {
        sideFiberGlassBlockV->SetMedium(sideFiberGlassBlockV_medium);
    } else
        Fatal("Main", "Invalid medium for sideFiberGlassBlockV!");
 
    // volume visual property (transparency)
    sideFiberGlassBlockV->SetLineColor(TColor::GetColor("#88ff88"));
    sideFiberGlassBlockV->SetTransparency(0);
 
    TGeoCombiTrans* sideFiberGlassBlock_transf[2];
 
    sideFiberGlassBlock_transf[0] = new TGeoCombiTrans();
    sideFiberGlassBlock_transf[0]->SetDx(+(lowerFiberGlassBlock_XSize * 0.5 + sideFiberGlassBlockPart_XSize * 0.5));
    sideFiberGlassBlock_transf[0]->SetDy(+0.45 /*shift*/);
    sideFiberGlassBlock_transf[0]->SetDz(+0.005 /*shift for fitting-elements as a plane*/);
 
    sideFiberGlassBlock_transf[1] = new TGeoCombiTrans();
    sideFiberGlassBlock_transf[1]->SetDx(-(lowerFiberGlassBlock_XSize * 0.5 + sideFiberGlassBlockPart_XSize * 0.5));
    sideFiberGlassBlock_transf[1]->SetDy(+0.45 /*shift*/);
    sideFiberGlassBlock_transf[1]->SetDz(+0.005 /*shift for fitting-elements as a plane*/);
 
    frames->AddNode(sideFiberGlassBlockV, 0, sideFiberGlassBlock_transf[0]);
    frames->AddNode(sideFiberGlassBlockV, 1, sideFiberGlassBlock_transf[1]);
    //--------------------------------------------------------------------------
 
    // upper fiberglass block ---------------------------------------------------
    Double_t upperFiberGlassBlock_XSize = 165.1;   // cm
    Double_t upperFiberGlassBlock_YSize = 1.45;    // cm
    Double_t upperFiberGlassBlock_ZSize = 0.89;    // cm
 
    TGeoShape* upperFiberGlassBlockS = new TGeoBBox("upperFiberGlassBlockS", upperFiberGlassBlock_XSize * 0.5,
                                                    upperFiberGlassBlock_YSize * 0.5, upperFiberGlassBlock_ZSize * 0.5);
 
    TGeoVolume* upperFiberGlassBlockV =
        new TGeoVolume(TString("upperFiberGlassBlockV") += TString("_") + frames->GetName(), upperFiberGlassBlockS);
 
    // volume medium
    TGeoMedium* upperFiberGlassBlockV_medium = pMedFiberGlass;
    if (upperFiberGlassBlockV_medium) {
        upperFiberGlassBlockV->SetMedium(upperFiberGlassBlockV_medium);
    } else
        Fatal("Main", "Invalid medium for upperFiberGlassBlockV!");
 
    // volume visual property (transparency)
    upperFiberGlassBlockV->SetLineColor(TColor::GetColor("#88ff88"));
    upperFiberGlassBlockV->SetTransparency(0);
 
    TGeoCombiTrans* upperFiberGlassBlock_transf[1];
 
    upperFiberGlassBlock_transf[0] = new TGeoCombiTrans();
    upperFiberGlassBlock_transf[0]->SetDx(0.0);
    upperFiberGlassBlock_transf[0]->SetDy(
        +(sideFiberGlassBlockPart_YSize * 0.5 - upperFiberGlassBlock_YSize * 0.5 + 0.45 /*shift*/));
    upperFiberGlassBlock_transf[0]->SetDz(+0.005 /*shift for fitting-elements as a plane*/);
 
    frames->AddNode(upperFiberGlassBlockV, 0, upperFiberGlassBlock_transf[0]);
    //--------------------------------------------------------------------------
 
    // brass layer --------------------------------------------------------------
    Double_t upperBrasslayer_XSize = 165.1;   // cm
    Double_t upperBrasslayer_YSize = 1.45;    // cm
    Double_t upperBrasslayer_ZSize = 0.005;   // cm
 
    Double_t sideBrasslayer_XSize = 1.45;         // cm
    Double_t sideBrasslayer_YSize = 49.4 - 6.0;   // cm
    Double_t sideBrasslayer_ZSize = 0.005;        // cm
 
    TGeoShape* upperBrasslayerS = new TGeoBBox("upperBrasslayerS", upperBrasslayer_XSize * 0.5,
                                               upperBrasslayer_YSize * 0.5, upperBrasslayer_ZSize * 0.5);
    TGeoShape* sideBrasslayerS = new TGeoBBox("sideBrasslayerS", sideBrasslayer_XSize * 0.5, sideBrasslayer_YSize * 0.5,
                                              sideBrasslayer_ZSize * 0.5);
 
    TGeoVolume* upperBrasslayerV =
        new TGeoVolume(TString("upperBrasslayerV") += TString("_") + frames->GetName(), upperBrasslayerS);
    TGeoVolume* sideBrasslayerV =
        new TGeoVolume(TString("sideBrasslayerV") += TString("_") + frames->GetName(), sideBrasslayerS);
 
    // volume medium
    TGeoMedium* upperBrasslayerV_medium = pMedBrass;
    if (upperBrasslayerV_medium) {
        upperBrasslayerV->SetMedium(upperBrasslayerV_medium);
        sideBrasslayerV->SetMedium(upperBrasslayerV_medium);
    } else
        Fatal("Main", "Invalid medium for upperBrasslayerV_medium!");
 
    // volume visual property (transparency)
    upperBrasslayerV->SetLineColor(TColor::GetColor("#f07422"));
    upperBrasslayerV->SetTransparency(0);
    sideBrasslayerV->SetLineColor(TColor::GetColor("#f07422"));
    sideBrasslayerV->SetTransparency(0);
 
    TGeoCombiTrans* upperBrasslayer_transf[1];
 
    upperBrasslayer_transf[0] = new TGeoCombiTrans();
    upperBrasslayer_transf[0]->SetDx(0.0);
    upperBrasslayer_transf[0]->SetDy(
        +(sideFiberGlassBlockPart_YSize * 0.5 - upperBrasslayer_YSize * 0.5 + 0.45 /*shift*/));
    upperBrasslayer_transf[0]->SetDz(-ZModuleSize_Station163x39 * 0.5 + upperBrasslayer_ZSize * 0.5 + 0.005 /*shift*/);
 
    TGeoCombiTrans* sideBrasslayer_transf[2];
 
    sideBrasslayer_transf[0] = new TGeoCombiTrans();
    sideBrasslayer_transf[0]->SetDx(+(lowerFiberGlassBlock_XSize * 0.5 + sideBrasslayer_XSize * 0.5));
    sideBrasslayer_transf[0]->SetDy(+0.45 /*shift*/);
    sideBrasslayer_transf[0]->SetDz(-ZModuleSize_Station163x39 * 0.5 + sideBrasslayer_ZSize * 0.5 + 0.005 /*shift*/);
 
    sideBrasslayer_transf[1] = new TGeoCombiTrans();
    sideBrasslayer_transf[1]->SetDx(-(lowerFiberGlassBlock_XSize * 0.5 + sideBrasslayer_XSize * 0.5));
    sideBrasslayer_transf[1]->SetDy(+0.45 /*shift*/);
    sideBrasslayer_transf[1]->SetDz(-ZModuleSize_Station163x39 * 0.5 + sideBrasslayer_ZSize * 0.5 + 0.005 /*shift*/);
 
    frames->AddNode(upperBrasslayerV, 0, upperBrasslayer_transf[0]);
 
    frames->AddNode(sideBrasslayerV, 0, sideBrasslayer_transf[0]);
    frames->AddNode(sideBrasslayerV, 1, sideBrasslayer_transf[1]);
    //--------------------------------------------------------------------------
 
    // steel layer --------------------------------------------------------------
    Double_t upperSteellayer_XSize = 165.1;   // cm
    Double_t upperSteellayer_YSize = 1.45;    // cm
    Double_t upperSteellayer_ZSize = 0.005;   // cm
 
    Double_t sideSteellayer_XSize = 1.45;         // cm
    Double_t sideSteellayer_YSize = 49.4 - 6.0;   // cm
    Double_t sideSteellayer_ZSize = 0.005;        // cm
 
    TGeoShape* upperSteellayerS = new TGeoBBox("upperSteellayerS", upperSteellayer_XSize * 0.5,
                                               upperSteellayer_YSize * 0.5, upperSteellayer_ZSize * 0.5);
    TGeoShape* sideSteellayerS = new TGeoBBox("sideSteellayerS", sideSteellayer_XSize * 0.5, sideSteellayer_YSize * 0.5,
                                              sideSteellayer_ZSize * 0.5);
 
    TGeoVolume* upperSteellayerV =
        new TGeoVolume(TString("upperSteellayerV") += TString("_") + frames->GetName(), upperSteellayerS);
    TGeoVolume* sideSteellayerV =
        new TGeoVolume(TString("sideSteellayerV") += TString("_") + frames->GetName(), sideSteellayerS);
 
    // volume medium
    TGeoMedium* upperSteellayerV_medium = pMedSteel;
    if (upperSteellayerV_medium) {
        upperSteellayerV->SetMedium(upperSteellayerV_medium);
        sideSteellayerV->SetMedium(upperSteellayerV_medium);
    } else
        Fatal("Main", "Invalid medium for upperSteellayerV_medium!");
 
    // volume visual property (transparency)
    upperSteellayerV->SetLineColor(TColor::GetColor("#42aaff"));
    upperSteellayerV->SetTransparency(0);
    sideSteellayerV->SetLineColor(TColor::GetColor("#42aaff"));
    sideSteellayerV->SetTransparency(0);
 
    TGeoCombiTrans* upperSteellayer_transf[1];
 
    upperSteellayer_transf[0] = new TGeoCombiTrans();
    upperSteellayer_transf[0]->SetDx(0.0);
    upperSteellayer_transf[0]->SetDy(
        +(sideFiberGlassBlockPart_YSize * 0.5 - upperSteellayer_YSize * 0.5 + 0.45 /*shift*/));
    upperSteellayer_transf[0]->SetDz(-ZModuleSize_Station163x39 * 0.5 + upperSteellayer_ZSize * 0.5 + 0.005 /*shift*/);
 
    TGeoCombiTrans* sideSteellayer_transf[2];
 
    sideSteellayer_transf[0] = new TGeoCombiTrans();
    sideSteellayer_transf[0]->SetDx(+(lowerFiberGlassBlock_XSize * 0.5 + sideSteellayer_XSize * 0.5));
    sideSteellayer_transf[0]->SetDy(+0.45 /*shift*/);
    sideSteellayer_transf[0]->SetDz(-ZModuleSize_Station163x39 * 0.5 + sideSteellayer_ZSize * 0.5);
 
    sideSteellayer_transf[1] = new TGeoCombiTrans();
    sideSteellayer_transf[1]->SetDx(-(lowerFiberGlassBlock_XSize * 0.5 + sideSteellayer_XSize * 0.5));
    sideSteellayer_transf[1]->SetDy(+0.45 /*shift*/);
    sideSteellayer_transf[1]->SetDz(-ZModuleSize_Station163x39 * 0.5 + sideSteellayer_ZSize * 0.5);
 
    frames->AddNode(upperSteellayerV, 0, upperSteellayer_transf[0]);
 
    frames->AddNode(sideSteellayerV, 0, sideSteellayer_transf[0]);
    frames->AddNode(sideSteellayerV, 1, sideSteellayer_transf[1]);
    //--------------------------------------------------------------------------
 
    return frames;
}
//------------------------------------------------------------------------------
 
TGeoVolume* CreateModule_Station10x10(TString module_name, Double_t xsize, Double_t ysize, Double_t zsize)
{
 
    // module shape
    TGeoShape* moduleS = new TGeoBBox("moduleS", xsize * 0.5, ysize * 0.5, zsize * 0.5);
 
    // module volume
    TGeoVolume* moduleV = new TGeoVolume(module_name, moduleS);
 
    // medium
    if (pMedArC4H10_80_20)
        moduleV->SetMedium(pMedArC4H10_80_20);
    else
        Fatal("Main", "Invalid medium for sensitive modules!");
 
    // visual parameters
    moduleV->SetLineColor(TColor::GetColor("#ff47ca"));
    moduleV->SetTransparency(20);
 
    return moduleV;
}
//------------------------------------------------------------------------------