Shifted platform

Plaform shifted by 40cm to be under the MUCH/RICH. Based upon create_platform_v16.C

macro/passive/create_platform_v22b.C

+/* Copyright (C) 2022 Facility of Antiproton and Ion Research in Europe, Darmstadt
+   SPDX-License-Identifier: GPL-3.0-only
+   Authors: Eoin Clerkin [committer] */
+
+/******************************************************************************
+ ** based on create_platform_v16.C
+ ** 
+ **/
+
+
+#include "TGeoManager.h"
+
+#include <iomanip>
+#include <iostream>
+
+
+using namespace std;
+
+// -------------   Other global variables   -----------------------------------
+// ---> TGeoManager (too lazy to write out 'Manager' all the time
+TGeoManager* gGeoMan = NULL;  // will be set later
+// ----------------------------------------------------------------------------
+
+// ============================================================================
+// ======                         Main function                           =====
+// ============================================================================
+
+void create_platform_v22b()
+{
+
+  // -----   Define platform parts   ----------------------------------------------
+
+  /** For v22b (SIS 100)   **/
+  TString geoTag = "v22b";
+  Double_t sizeX = 725.0;  // symmetric in x
+  Double_t sizeY = 234.0;  // without rails
+  Double_t sizeZ = 297.0;  // short version
+  Double_t endZ  = 410.0;  // ends at z = 450.0
+
+  //  /** For v16b (SIS 300)   **/
+  //  TString geoTag = "v16b";
+  //  Double_t sizeX = 725.0;  // symmetric in x
+  //  Double_t sizeY = 234.0;  // without rails
+  //  Double_t sizeZ = 455.0;  // long version
+  //  Double_t endZ  = 540.0;  // ends at z = 450.0
+
+  Double_t beamY = 570.0;  // nominal beam height
+  Double_t posX  = 0.0;
+  Double_t posY  = -beamY + sizeY / 2.;  // rest on the floor at -beamY
+  Double_t posZ  = endZ - sizeZ / 2.;    // end at MUCH z-end
+
+  // --------------------------------------------------------------------------
+
+
+  // -------   Geometry file name (output)   ----------------------------------
+  TString geoFileName = "platform_";
+  geoFileName         = geoFileName + geoTag + ".geo.root";
+  // --------------------------------------------------------------------------
+
+
+  // -------   Open info file   -----------------------------------------------
+  TString infoFileName = geoFileName;
+  infoFileName.ReplaceAll("root", "info");
+  fstream infoFile;
+  infoFile.open(infoFileName.Data(), fstream::out);
+  infoFile << "Platform geometry created with create_platform_v16.C" << std::endl << std::endl;
+  // --------------------------------------------------------------------------
+
+
+  // -------   Load media from media file   -----------------------------------
+  FairGeoLoader* geoLoad    = new FairGeoLoader("TGeo", "FairGeoLoader");
+  FairGeoInterface* geoFace = geoLoad->getGeoInterface();
+  TString geoPath           = gSystem->Getenv("VMCWORKDIR");
+  TString medFile           = geoPath + "/geometry/media.geo";
+  geoFace->setMediaFile(medFile);
+  geoFace->readMedia();
+  gGeoMan = gGeoManager;
+  // --------------------------------------------------------------------------
+
+
+  // -----------------   Get and create the required media    -----------------
+  FairGeoMedia* geoMedia   = geoFace->getMedia();
+  FairGeoBuilder* geoBuild = geoLoad->getGeoBuilder();
+
+  // ---> aluminium
+  FairGeoMedium* mAluminium = geoMedia->getMedium("aluminium");
+  if (!mAluminium) Fatal("Main", "FairMedium aluminium not found");
+  geoBuild->createMedium(mAluminium);
+  TGeoMedium* aluminium = gGeoMan->GetMedium("aluminium");
+  if (!aluminium) Fatal("Main", "Medium aluminium not found");
+
+  //  // ---> air
+  //  FairGeoMedium* mAir  = geoMedia->getMedium("air");
+  //  if ( ! mAir ) Fatal("Main", "FairMedium air not found");
+  //  geoBuild->createMedium(mAir);
+  //  TGeoMedium* air = gGeoMan->GetMedium("air");
+  //  if ( ! air ) Fatal("Main", "Medium air not found");
+  // --------------------------------------------------------------------------
+
+
+  // --------------   Create geometry and top volume  -------------------------
+  gGeoMan = (TGeoManager*) gROOT->FindObject("FAIRGeom");
+  gGeoMan->SetName("PLATFORMgeom");
+  TGeoVolume* top = new TGeoVolumeAssembly("top");
+  gGeoMan->SetTopVolume(top);
+  TString platformName = "platform_";
+  platformName += geoTag;
+  TGeoVolume* platform = new TGeoVolumeAssembly(platformName.Data());
+  // --------------------------------------------------------------------------
+
+
+  // -----   Create   ---------------------------------------------------------
+  TGeoBBox* platform_base  = new TGeoBBox("", sizeX / 2., sizeY / 2., sizeZ / 2.);
+  TGeoVolume* platform_vol = new TGeoVolume("platform", platform_base, aluminium);
+  platform_vol->SetLineColor(kBlue);
+  platform_vol->SetTransparency(70);
+
+  TGeoTranslation* platform_trans = new TGeoTranslation("", posX, posY, posZ);
+  platform->AddNode(platform_vol, 1, platform_trans);
+
+  infoFile << "sizeX: " << setprecision(2) << sizeX << "sizeY: " << setprecision(2) << sizeY
+           << "sizeZ: " << setprecision(2) << sizeZ << endl;
+  infoFile << "posX : " << setprecision(2) << posX << "posY : " << setprecision(2) << posY
+           << "posZ : " << setprecision(2) << posZ << endl;
+
+  // ---------------   Finish   -----------------------------------------------
+  top->AddNode(platform, 1);
+  cout << endl << endl;
+  gGeoMan->CloseGeometry();
+  gGeoMan->CheckOverlaps(0.001);
+  gGeoMan->PrintOverlaps();
+  gGeoMan->Test();
+
+  platform->Export(geoFileName);  // an alternative way of writing the platform volume
+
+  TFile* geoFile                      = new TFile(geoFileName, "UPDATE");
+  TGeoTranslation* platform_placement = new TGeoTranslation("platform_trans", 0., 0., 0.);
+  platform_placement->Write();
+  geoFile->Close();
+
+  cout << endl;
+  cout << "Geometry " << top->GetName() << " written to " << geoFileName << endl;
+
+  top->Draw("ogl");
+
+  infoFile.close();
+}
+// ============================================================================
+// ======                   End of main function                          =====
+// ============================================================================