From ca288bf8ca19fcc3147cd5e892427cac470c185a Mon Sep 17 00:00:00 2001
From: Omveer <osingh@amu.ac.in>
Date: Fri, 16 Apr 2021 18:47:28 +0530
Subject: [PATCH] Clang format for v21 much geometries.
Changes to obey coding rules.
---
.../create_MUCH_geometry_v20b_lmvm_start.C | 404 ++---
.../geometry/create_MUCH_geometry_v21a_jpsi.C | 1522 ++++++++---------
.../geometry/create_MUCH_geometry_v21a_lmvm.C | 1358 ++++++++-------
3 files changed, 1575 insertions(+), 1709 deletions(-)
diff --git a/macro/much/geometry/create_MUCH_geometry_v20b_lmvm_start.C b/macro/much/geometry/create_MUCH_geometry_v20b_lmvm_start.C
index 5279938599..2cc263b0ff 100644
--- a/macro/much/geometry/create_MUCH_geometry_v20b_lmvm_start.C
+++ b/macro/much/geometry/create_MUCH_geometry_v20b_lmvm_start.C
@@ -33,14 +33,12 @@
#include "TGeoVolume.h"
#include "TGeoXtru.h"
#include "TList.h"
+#include "TMath.h"
+#include "TObjArray.h"
#include "TRandom3.h"
#include "TString.h"
#include "TSystem.h"
-#include "TFile.h"
-#include "TMath.h"
-#include "TObjArray.h"
-
#include <cassert>
#include <fstream>
#include <iostream>
@@ -62,13 +60,13 @@ const TString L = "MUCHlead";
const TString I = "MUCHiron";
const TString activemedium = "MUCHargon";
const TString spacermedium = "MUCHnoryl";
-const TString LDcarbon = "MUCHcarbonLD"; //Low density Carbon for Ist Abs
-const TString Concrete = "MUCHconcrete"; //Concrete for Ist Abs
-const TString Al = "MUCHaluminium"; //Al for cooling & support purpose
-const TString copper = "MUCHcopper";
-const TString g10 = "MUCHG10";
-const TString RPCm = "MUCHRPCgas";
-const TString RPCg = "MUCHRPCglass";
+const TString LDcarbon = "MUCHcarbonLD"; //Low density Carbon for Ist Abs
+const TString Concrete = "MUCHconcrete"; //Concrete for Ist Abs
+const TString Al = "MUCHaluminium"; //Al for cooling & support purpose
+const TString copper = "MUCHcopper";
+const TString g10 = "MUCHG10";
+const TString RPCm = "MUCHRPCgas";
+const TString RPCg = "MUCHRPCglass";
// Universal input parameters
@@ -81,10 +79,9 @@ Double_t fAcceptanceTanMax = 0.466; // Acceptance tangent max
// Input parameters for absorbers
//***********************************************************
-const Int_t fAbs = 3;
-const Int_t fNabs = 5; // Number of absorber pieces
-TString AbsMaterial[6] =
- {"LD Graphite", "LD Graphite", "Concrete", "Iron", "Iron", "Iron"};
+const Int_t fAbs = 3;
+const Int_t fNabs = 5; // Number of absorber pieces
+TString AbsMaterial[6] = {"LD Graphite", "LD Graphite", "Concrete", "Iron", "Iron", "Iron"};
// Absorber Zin position [cm] in the cave reference frame
Double_t fAbsorberZ1[6] = {0, 16, 28, 90, 140, 190};
// Absorber thickness [cm]
@@ -100,15 +97,12 @@ const Int_t fNst = 2; // Number of stations
// Number of sectors per layer (should be even for symmetry)
// Needs to be fixed with actual numbers
Int_t fNSectorsPerLayer[4] = {16, 20, 18, 20};
-Double_t fActiveLzSector[4] = {0.3,
- 0.3,
- 0.2,
- 0.2}; // Active volume thickness [cm]
-Double_t fFrameLzSector[4] = {.3, .3, 1., 1.}; // Frame thickness [cm]
+Double_t fActiveLzSector[4] = {0.3, 0.3, 0.2, 0.2}; // Active volume thickness [cm]
+Double_t fFrameLzSector[4] = {.3, .3, 1., 1.}; // Frame thickness [cm]
Double_t fRpcGlassDz[4] = {0.0, 0.0, 0.2, 0.2}; //Rpc Glass thickness [cm]
Double_t fSpacerR = 2.0; // Spacer width in R [cm]
Double_t fSpacerPhi = 2.0; // Spacer width in Phi [cm]
-Double_t fOverlapR = 2.0; // Overlap in R direction [cm]
+Double_t fOverlapR = 2.0; // Overlap in R direction [cm]
// Station Zceneter [cm] in the cave reference frame
@@ -124,14 +118,12 @@ Double_t fG10Dz = 0.3; // 3 mm G10
// spans from 2.9 degree to 5.1 degree
//Inner radius is tan(2.9) + 2 cm, extra 20 mm for clamp connection
-const Int_t fNshs = 4;
-TString ShMaterial[5] = {"Al", "Pb", "Al", "Al", "Al"};
-Double_t fShieldZin[5] = {125, 153.0, 215.0, 265.0, 315.0};
-Double_t fShieldLz[5] = {28, 30, 20, 20, 30};
-Double_t fShield_AcceptanceTanMin[5] =
- {0.043, 0.043, 0.051, 0.051, 0.051}; // Acceptance tangent min for shield
-Double_t fShield_AcceptanceTanMax[5] =
- {0.1, 0.1, 0.1, 0.1, 0.1}; // Acceptance tangent max for shield
+const Int_t fNshs = 4;
+TString ShMaterial[5] = {"Al", "Pb", "Al", "Al", "Al"};
+Double_t fShieldZin[5] = {125, 153.0, 215.0, 265.0, 315.0};
+Double_t fShieldLz[5] = {28, 30, 20, 20, 30};
+Double_t fShield_AcceptanceTanMin[5] = {0.043, 0.043, 0.051, 0.051, 0.051}; // Acceptance tangent min for shield
+Double_t fShield_AcceptanceTanMax[5] = {0.1, 0.1, 0.1, 0.1, 0.1}; // Acceptance tangent max for shield
//***********************************************************
@@ -151,17 +143,15 @@ TGeoVolume* CreateLayers(int istn, int ily);
fstream infoFile;
-void create_MUCH_geometry_v20b_lmvm_start() {
+void create_MUCH_geometry_v20b_lmvm_start()
+{
// ------- Open info file -----------------------------------------------
TString infoFileName = FileNameSim;
infoFileName.ReplaceAll("root", "info");
infoFile.open(infoFileName.Data(), fstream::out);
- infoFile
- << "MUCH geometry created with create_MUCH_geometry_v20b_lmvm_start.C"
- << endl
- << endl;
+ infoFile << "MUCH geometry created with create_MUCH_geometry_v20b_lmvm_start.C" << endl << endl;
infoFile << "Global Variables: " << endl;
infoFile << "MuchCave Zin position = " << fMuchZ1 << " cm " << endl;
infoFile << "Acceptance tangent min = " << fAcceptanceTanMin << endl;
@@ -202,10 +192,8 @@ void create_MUCH_geometry_v20b_lmvm_start() {
infoFile << "Second half is made of Low Density Graphite + Concrete." << endl;
infoFile << "Total No. of Pieces: " << fNabs << endl;
infoFile << endl;
- infoFile << "AbsPieces Position[cm] Thickness[cm] Material"
- << endl;
- infoFile << "--------------------------------------------------------------"
- << endl;
+ infoFile << "AbsPieces Position[cm] Thickness[cm] Material" << endl;
+ infoFile << "--------------------------------------------------------------" << endl;
for (Int_t iabs = 0; iabs < fNabs; iabs++) {
@@ -221,10 +209,8 @@ void create_MUCH_geometry_v20b_lmvm_start() {
infoFile << "No. of Shields: " << fNshs << endl;
infoFile << "Inside the Abs I, Shielding divided into two parts." << endl;
infoFile << endl;
- infoFile << "Shield No. Z_In[cm] Z_Out[cm] R_In[cm] R_Out[cm] Material"
- << endl;
- infoFile << "--------------------------------------------------------------"
- << endl;
+ infoFile << "Shield No. Z_In[cm] Z_Out[cm] R_In[cm] R_Out[cm] Material" << endl;
+ infoFile << "--------------------------------------------------------------" << endl;
for (Int_t ishi = 0; ishi < fNshs; ishi++) {
gModules_shield[ishi] = CreateShields(ishi);
@@ -245,8 +231,7 @@ void create_MUCH_geometry_v20b_lmvm_start() {
"realistic material budget for GEM modules."
<< endl;
infoFile << "#Station #Layers Z[cm] #Sectors ActiveLz[cm]" << endl;
- infoFile << "--------------------------------------------------------------"
- << endl;
+ infoFile << "--------------------------------------------------------------" << endl;
for (Int_t istn = 0; istn < fNst; istn++) { // 4 Stations
@@ -265,9 +250,8 @@ void create_MUCH_geometry_v20b_lmvm_start() {
much->Export(FileNameSim); // an alternative way of writing the much
- TFile* outfile = new TFile(FileNameSim, "UPDATE");
- TGeoTranslation* much_placement =
- new TGeoTranslation("much_trans", 0., 0., 0.);
+ TFile* outfile = new TFile(FileNameSim, "UPDATE");
+ TGeoTranslation* much_placement = new TGeoTranslation("much_trans", 0., 0., 0.);
much_placement->Write();
outfile->Close();
@@ -279,7 +263,8 @@ void create_MUCH_geometry_v20b_lmvm_start() {
infoFile.close();
}
-void create_materials_from_media_file() {
+void create_materials_from_media_file()
+{
// Use the FairRoot geometry interface to load the media which are already defined
FairGeoLoader* geoLoad = new FairGeoLoader("TGeo", "FairGeoLoader");
FairGeoInterface* geoFace = geoLoad->getGeoInterface();
@@ -330,7 +315,8 @@ void create_materials_from_media_file() {
}
-TGeoVolume* CreateShields(int ish) {
+TGeoVolume* CreateShields(int ish)
+{
TGeoMedium* iron = gGeoMan->GetMedium(I);
@@ -353,12 +339,12 @@ TGeoVolume* CreateShields(int ish) {
Double_t rmin2 = globalZ2 * fShield_AcceptanceTanMin[ish] + 2.0;
Double_t rmax2 = globalZ2 * fShield_AcceptanceTanMax[ish] - 0.0001;
- infoFile << " " << ish << "\t\t" << globalZ1 << "\t " << globalZ2 << "\t"
- << rmin1 << "\t " << rmax1 << "\t " << ShMaterial[ish] << endl;
+ infoFile << " " << ish << "\t\t" << globalZ1 << "\t " << globalZ2 << "\t" << rmin1 << "\t " << rmax1 << "\t "
+ << ShMaterial[ish] << endl;
if (ish == 0) {
- TGeoCone* sh = new TGeoCone(conename_sh, dz, rmin1, rmax1, rmin2, rmax2);
+ TGeoCone* sh = new TGeoCone(conename_sh, dz, rmin1, rmax1, rmin2, rmax2);
TGeoVolume* shield = new TGeoVolume("shield", sh, Aluminium);
shield->SetLineColor(7);
@@ -368,7 +354,7 @@ TGeoVolume* CreateShields(int ish) {
}
if (ish == 1) {
- TGeoCone* sh = new TGeoCone(conename_sh, dz, rmin1, rmax1, rmin2, rmax2);
+ TGeoCone* sh = new TGeoCone(conename_sh, dz, rmin1, rmax1, rmin2, rmax2);
TGeoVolume* shield = new TGeoVolume("shield", sh, lead);
shield->SetLineColor(kMagenta);
@@ -380,7 +366,7 @@ TGeoVolume* CreateShields(int ish) {
if (ish > 1) {
- TGeoCone* sh = new TGeoCone(conename_sh, dz, rmin1, rmax1, rmin2, rmax2);
+ TGeoCone* sh = new TGeoCone(conename_sh, dz, rmin1, rmax1, rmin2, rmax2);
TGeoVolume* shield = new TGeoVolume("shield", sh, Aluminium);
shield->SetLineColor(kBlack);
@@ -394,7 +380,8 @@ TGeoVolume* CreateShields(int ish) {
}
-TGeoVolume* CreateAbsorbers(int i) {
+TGeoVolume* CreateAbsorbers(int i)
+{
TGeoMedium* graphite = gGeoMan->GetMedium(LDcarbon);
TGeoMedium* iron = gGeoMan->GetMedium(I);
@@ -420,24 +407,21 @@ TGeoVolume* CreateAbsorbers(int i) {
Double_t rmax1 = globalZ1 * fAcceptanceTanMax + safetyrad[i];
Double_t rmax2 = globalZ2 * fAcceptanceTanMax + safetyrad[i]; //
- infoFile << " " << i + 1 << "\t\t" << globalPos << "\t\t" << 2 * dz
- << "\t\t" << AbsMaterial[i] << endl;
+ infoFile << " " << i + 1 << "\t\t" << globalPos << "\t\t" << 2 * dz << "\t\t" << AbsMaterial[i] << endl;
// 1st part of 1st absorber LD Carbon
//dimensions are hardcoded
if (i == 0) {
printf("absorber %d \n", i);
- TGeoTrd2* trap = new TGeoTrd2(TrapName, 70.0, 70.0, 46.0, 71.0, dz);
- TGeoCone* tube = new TGeoCone(pipename, dz + 0.001, 0., rmin1, 0., rmin2);
- TString expression = TrapName + "-" + pipename;
- TGeoCompositeShape* shSupport =
- new TGeoCompositeShape(supportShapeName, expression);
- TGeoVolume* abs0 = new TGeoVolume("absorber", shSupport, graphite);
+ TGeoTrd2* trap = new TGeoTrd2(TrapName, 70.0, 70.0, 46.0, 71.0, dz);
+ TGeoCone* tube = new TGeoCone(pipename, dz + 0.001, 0., rmin1, 0., rmin2);
+ TString expression = TrapName + "-" + pipename;
+ TGeoCompositeShape* shSupport = new TGeoCompositeShape(supportShapeName, expression);
+ TGeoVolume* abs0 = new TGeoVolume("absorber", shSupport, graphite);
abs0->SetLineColor(kRed);
abs0->SetTransparency(0);
- TGeoTranslation* abs0_trans =
- new TGeoTranslation("", 0., 0., globalZ1 + dz);
+ TGeoTranslation* abs0_trans = new TGeoTranslation("", 0., 0., globalZ1 + dz);
absblock->AddNode(abs0, i, abs0_trans);
}
@@ -445,45 +429,40 @@ TGeoVolume* CreateAbsorbers(int i) {
// 2rd part of 1st absorber box (LD Carbon)
if (i == 1) {
printf("absorber %d \n", i);
- TGeoBBox* box = new TGeoBBox(BoxName, 130.0, 125.0, dz);
- TGeoCone* tube = new TGeoCone(pipename, dz + 0.001, 0., rmin1, 0., rmin2);
- TString expression = BoxName + "-" + pipename;
- TGeoCompositeShape* shSupport =
- new TGeoCompositeShape(supportShapeName, expression);
+ TGeoBBox* box = new TGeoBBox(BoxName, 130.0, 125.0, dz);
+ TGeoCone* tube = new TGeoCone(pipename, dz + 0.001, 0., rmin1, 0., rmin2);
+ TString expression = BoxName + "-" + pipename;
+ TGeoCompositeShape* shSupport = new TGeoCompositeShape(supportShapeName, expression);
TGeoVolume* abs1 = new TGeoVolume("absorber", shSupport, graphite);
abs1->SetLineColor(kRed);
abs1->SetTransparency(0);
- TGeoTranslation* abs1_trans =
- new TGeoTranslation("", 0., 0., globalZ1 + dz);
+ TGeoTranslation* abs1_trans = new TGeoTranslation("", 0., 0., globalZ1 + dz);
absblock->AddNode(abs1, i, abs1_trans);
}
// 3th part of 1st absorber box (Conc)
if (i == 2) {
printf("absorber %d \n", i);
- TGeoBBox* box = new TGeoBBox(BoxName, 130.0, 125.0, dz);
- TGeoCone* tube = new TGeoCone(pipename, dz + 0.001, 0., rmin1, 0., rmin2);
- TString expression = BoxName + "-" + pipename;
- TGeoCompositeShape* shSupport =
- new TGeoCompositeShape(supportShapeName, expression);
+ TGeoBBox* box = new TGeoBBox(BoxName, 130.0, 125.0, dz);
+ TGeoCone* tube = new TGeoCone(pipename, dz + 0.001, 0., rmin1, 0., rmin2);
+ TString expression = BoxName + "-" + pipename;
+ TGeoCompositeShape* shSupport = new TGeoCompositeShape(supportShapeName, expression);
TGeoVolume* abs1 = new TGeoVolume("absorber", shSupport, concrete);
abs1->SetLineColor(kViolet);
abs1->SetTransparency(0);
- TGeoTranslation* abs1_trans =
- new TGeoTranslation("", 0., 0., globalZ1 + dz);
+ TGeoTranslation* abs1_trans = new TGeoTranslation("", 0., 0., globalZ1 + dz);
absblock->AddNode(abs1, i, abs1_trans);
}
//rest of the absorbers
if (i > 2) {
- TGeoBBox* box = new TGeoBBox(BoxName, rmax2, rmax2, dz);
- TGeoCone* tube = new TGeoCone(pipename, dz + 0.001, 0., rmin1, 0., rmin2);
- TString expression = BoxName + "-" + pipename;
- TGeoCompositeShape* shSupport =
- new TGeoCompositeShape(supportShapeName, expression);
+ TGeoBBox* box = new TGeoBBox(BoxName, rmax2, rmax2, dz);
+ TGeoCone* tube = new TGeoCone(pipename, dz + 0.001, 0., rmin1, 0., rmin2);
+ TString expression = BoxName + "-" + pipename;
+ TGeoCompositeShape* shSupport = new TGeoCompositeShape(supportShapeName, expression);
TGeoVolume* abs2 = new TGeoVolume("absorber", shSupport, iron);
@@ -498,12 +477,12 @@ TGeoVolume* CreateAbsorbers(int i) {
return absblock;
}
-TGeoVolume* CreateStations(int ist) {
+TGeoVolume* CreateStations(int ist)
+{
TString stationName = Form("muchstation%02i", ist + 1);
- TGeoVolumeAssembly* station =
- new TGeoVolumeAssembly(stationName); //, shStation, air);
+ TGeoVolumeAssembly* station = new TGeoVolumeAssembly(stationName); //, shStation, air);
TGeoVolume* gLayer[4];
@@ -517,17 +496,15 @@ TGeoVolume* CreateStations(int ist) {
}
-TGeoVolume* CreateLayers(int istn, int ily) {
+TGeoVolume* CreateLayers(int istn, int ily)
+{
- TString layerName = Form("muchstation%02ilayer%i", istn + 1, ily + 1);
+ TString layerName = Form("muchstation%02ilayer%i", istn + 1, ily + 1);
TGeoVolumeAssembly* volayer = new TGeoVolumeAssembly(layerName);
- Double_t stGlobalZ0 =
- fStationZ0[istn] + fMuchZ1; //z position of station center (midplane) [cm]
- Double_t stDz = ((fNlayers[istn] - 1) * fLayersDz[istn] + fSupportLz[istn]
- + 2 * fActiveLzSector[istn])
- / 2.;
+ Double_t stGlobalZ0 = fStationZ0[istn] + fMuchZ1; //z position of station center (midplane) [cm]
+ Double_t stDz = ((fNlayers[istn] - 1) * fLayersDz[istn] + fSupportLz[istn] + 2 * fActiveLzSector[istn]) / 2.;
Double_t stGlobalZ2 = stGlobalZ0 + stDz;
Double_t stGlobalZ1 = stGlobalZ0 - stDz;
@@ -541,29 +518,24 @@ TGeoVolume* CreateLayers(int istn, int ily) {
Double_t RpcGlassDz = fActiveLzSector[istn] / 2. + fRpcGlassDz[istn] / 2.;
Double_t driftDz = RpcGlassDz + fRpcGlassDz[istn] / 2 + fDriftDz / 2.;
Double_t g10Dz = driftDz + fDriftDz / 2. + fG10Dz / 2.;
- Double_t moduleZ = fSupportLz[istn] / 2. + fActiveLzSector[istn] / 2.
- + fRpcGlassDz[istn] + fDriftDz + fG10Dz;
- Double_t RpcGlassZIn = moduleZ - RpcGlassDz;
- Double_t RpcGlassZOut = moduleZ + RpcGlassDz;
- Double_t driftZIn = moduleZ - driftDz;
- Double_t driftZOut = moduleZ + driftDz;
- Double_t g10ZIn = moduleZ - g10Dz;
- Double_t g10ZOut = moduleZ + g10Dz;
-
-
- Double_t phi0 =
- TMath::Pi()
- / fNSectorsPerLayer[istn]; // azimuthal half widh of each module
+ Double_t moduleZ = fSupportLz[istn] / 2. + fActiveLzSector[istn] / 2. + fRpcGlassDz[istn] + fDriftDz + fG10Dz;
+ Double_t RpcGlassZIn = moduleZ - RpcGlassDz;
+ Double_t RpcGlassZOut = moduleZ + RpcGlassDz;
+ Double_t driftZIn = moduleZ - driftDz;
+ Double_t driftZOut = moduleZ + driftDz;
+ Double_t g10ZIn = moduleZ - g10Dz;
+ Double_t g10ZOut = moduleZ + g10Dz;
+
+
+ Double_t phi0 = TMath::Pi() / fNSectorsPerLayer[istn]; // azimuthal half widh of each module
Double_t ymin = rmin + fSpacerR;
Double_t ymax = rmax;
//define the dimensions of the trapezoidal module
- Double_t dy = (ymax - ymin) / 2.; //y (length)
- Double_t dx1 =
- ymin * TMath::Tan(phi0) + fOverlapR / TMath::Cos(phi0); // large x
- Double_t dx2 =
- ymax * TMath::Tan(phi0) + fOverlapR / TMath::Cos(phi0); // small x
- Double_t dz = fActiveLzSector[istn] / 2.; // thickness
+ Double_t dy = (ymax - ymin) / 2.; //y (length)
+ Double_t dx1 = ymin * TMath::Tan(phi0) + fOverlapR / TMath::Cos(phi0); // large x
+ Double_t dx2 = ymax * TMath::Tan(phi0) + fOverlapR / TMath::Cos(phi0); // small x
+ Double_t dz = fActiveLzSector[istn] / 2.; // thickness
//define the spacer dimensions
Double_t tg = (dx2 - dx1) / 2 / dy;
@@ -580,24 +552,21 @@ TGeoVolume* CreateLayers(int istn, int ily) {
Double_t dz_sG = fG10Dz / 2.; // 3mm G10
- infoFile << " " << istn + 1 << "\t " << ily + 1 << "\t\t"
- << layerGlobalZ0 << "\t" << fNSectorsPerLayer[istn] << "\t"
- << fActiveLzSector[istn] << endl;
+ infoFile << " " << istn + 1 << "\t " << ily + 1 << "\t\t" << layerGlobalZ0 << "\t" << fNSectorsPerLayer[istn]
+ << "\t" << fActiveLzSector[istn] << endl;
// Aluminum Plate (Cooling + Support)
TString supportAlName = Form("shStation%02iSupportAl", istn + 1);
TGeoTube* shSupportAl = new TGeoTube(supportAlName, rmin, rmax, SupportDz);
- TString supportName1 =
- Form("muchstation%02ilayer%isupportAl", istn + 1, ily + 1);
- TGeoMedium* coolMat = gGeoMan->GetMedium(Al);
+ TString supportName1 = Form("muchstation%02ilayer%isupportAl", istn + 1, ily + 1);
+ TGeoMedium* coolMat = gGeoMan->GetMedium(Al);
TGeoVolume* voSupport1 = new TGeoVolume(supportName1, shSupportAl, coolMat);
voSupport1->SetLineColor(kCyan);
- TGeoTranslation* support_trans1 =
- new TGeoTranslation("supportName1", 0, 0, layerGlobalZ0);
+ TGeoTranslation* support_trans1 = new TGeoTranslation("supportName1", 0, 0, layerGlobalZ0);
volayer->AddNode(voSupport1, 0, support_trans1);
Double_t pos[9];
@@ -608,10 +577,7 @@ TGeoVolume* CreateLayers(int istn, int ily) {
// Now start adding the GEM modules
for (Int_t iModule = 0; iModule < fNSectorsPerLayer[istn]; iModule++) {
- Double_t phi =
- 2 * phi0
- * (iModule
- + 0.5); // add 0.5 to not overlap with y-axis for left-right layer separation
+ Double_t phi = 2 * phi0 * (iModule + 0.5); // add 0.5 to not overlap with y-axis for left-right layer separation
Bool_t isBack = iModule % 2;
Char_t cside = (isBack == 1) ? 'b' : 'f';
@@ -633,81 +599,44 @@ TGeoVolume* CreateLayers(int istn, int ily) {
if (iModule != 0) iMod = iModule / 2;
- TGeoMedium* argon = gGeoMan->GetMedium(activemedium); // active medium
- TGeoMedium* noryl = gGeoMan->GetMedium(spacermedium); // spacer medium
- TGeoMedium* SpacerRpc = gGeoMan->GetMedium(Al); // spacer medium
- TGeoMedium* copperplate =
- gGeoMan->GetMedium(copper); // copper Drift medium
- TGeoMedium* g10plate = gGeoMan->GetMedium(g10); // G10 medium
- TGeoMedium* RPCglassmat = gGeoMan->GetMedium(RPCg); // RPC Glass
- TGeoMedium* RPCgasmedium = gGeoMan->GetMedium(RPCm); // RPC Gas
+ TGeoMedium* argon = gGeoMan->GetMedium(activemedium); // active medium
+ TGeoMedium* noryl = gGeoMan->GetMedium(spacermedium); // spacer medium
+ TGeoMedium* SpacerRpc = gGeoMan->GetMedium(Al); // spacer medium
+ TGeoMedium* copperplate = gGeoMan->GetMedium(copper); // copper Drift medium
+ TGeoMedium* g10plate = gGeoMan->GetMedium(g10); // G10 medium
+ TGeoMedium* RPCglassmat = gGeoMan->GetMedium(RPCg); // RPC Glass
+ TGeoMedium* RPCgasmedium = gGeoMan->GetMedium(RPCm); // RPC Gas
// Define and place the trapezoidal GEM module in X-Y plane
TGeoVolume *voActive, *voRPCback, *voRPCFront;
- TGeoTrap* shape =
- new TGeoTrap(sdz, 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
- shape->SetName(Form("shStation%02iLayer%i%cModule%03iFrameNoHole",
- istn,
- ily,
- cside,
- iModule));
-
- TGeoTrap* shapeActive =
- new TGeoTrap(dz, 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
- shapeActive->SetName(Form("shStation%02iLayer%i%cModule%03iActiveNoHole",
- istn,
- ily,
- cside,
- iModule));
+ TGeoTrap* shape = new TGeoTrap(sdz, 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
+ shape->SetName(Form("shStation%02iLayer%i%cModule%03iFrameNoHole", istn, ily, cside, iModule));
+
+ TGeoTrap* shapeActive = new TGeoTrap(dz, 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
+ shapeActive->SetName(Form("shStation%02iLayer%i%cModule%03iActiveNoHole", istn, ily, cside, iModule));
if (istn < 2) { //GEM
- TString activeName = Form("muchstation%02ilayer%i%cactive%03igasArgon",
- istn + 1,
- ily + 1,
- cside,
- iMod + 1);
+ TString activeName = Form("muchstation%02ilayer%i%cactive%03igasArgon", istn + 1, ily + 1, cside, iMod + 1);
voActive = new TGeoVolume(activeName, shapeActive, argon);
voActive->SetLineColor(kGreen);
- } else { //RPC Active
- TString activeName = Form("muchstation%02ilayer%i%cactive%03irpcgas",
- istn + 1,
- ily + 1,
- cside,
- iMod + 1);
- voActive = new TGeoVolume(activeName, shapeActive, RPCgasmedium);
+ }
+ else { //RPC Active
+ TString activeName = Form("muchstation%02ilayer%i%cactive%03irpcgas", istn + 1, ily + 1, cside, iMod + 1);
+ voActive = new TGeoVolume(activeName, shapeActive, RPCgasmedium);
voActive->SetLineColor(kRed);
//RPC Glass Front
- TGeoTrap* shapeRPCfront =
- new TGeoTrap(dz, 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
- shapeRPCfront->SetName(
- Form("shStation%02iLayer%i%cModule%03iNoHoleFront",
- istn,
- ily,
- cside,
- iModule));
- TString NameFront = Form("muchstation%02ilayer%i%c%03irpcglassFront",
- istn + 1,
- ily + 1,
- cside,
- iMod + 1);
- voRPCFront = new TGeoVolume(NameFront, shapeRPCfront, RPCglassmat);
+ TGeoTrap* shapeRPCfront = new TGeoTrap(dz, 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
+ shapeRPCfront->SetName(Form("shStation%02iLayer%i%cModule%03iNoHoleFront", istn, ily, cside, iModule));
+ TString NameFront = Form("muchstation%02ilayer%i%c%03irpcglassFront", istn + 1, ily + 1, cside, iMod + 1);
+ voRPCFront = new TGeoVolume(NameFront, shapeRPCfront, RPCglassmat);
voRPCFront->SetLineColor(kRed);
//RPC Glass Back
- TGeoTrap* shapeRPCback =
- new TGeoTrap(dz, 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
- shapeRPCback->SetName(Form("shStation%02iLayer%i%cModule%03iNoHoleback",
- istn,
- ily,
- cside,
- iModule));
- TString Nameback = Form("muchstation%02ilayer%i%c%03irpcglassback",
- istn + 1,
- ily + 1,
- cside,
- iMod + 1);
+ TGeoTrap* shapeRPCback = new TGeoTrap(dz, 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
+ shapeRPCback->SetName(Form("shStation%02iLayer%i%cModule%03iNoHoleback", istn, ily, cside, iModule));
+ TString Nameback = Form("muchstation%02ilayer%i%c%03irpcglassback", istn + 1, ily + 1, cside, iMod + 1);
voRPCback = new TGeoVolume(Nameback, shapeRPCback, RPCglassmat);
voRPCback->SetLineColor(kRed);
}
@@ -721,108 +650,55 @@ TGeoVolume* CreateLayers(int istn, int ily) {
Char_t cpos = (iPos == 0) ? 'i' : 'o';
// Define and place copper drift & G10
if (istn < 2) { //GEM
- Drift[iPos] =
- new TGeoTrap(dz_sD, 0, phi, sdy, sdx1, sdx2, 0, sdy, sdx1, sdx2, 0);
- G10[iPos] =
- new TGeoTrap(dz_sG, 0, phi, sdy, sdx1, sdx2, 0, sdy, sdx1, sdx2, 0);
- } else { //RPC
- Drift[iPos] =
- new TGeoTrap(dz_sD, 0, phi, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
- G10[iPos] =
- new TGeoTrap(dz_sG, 0, phi, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
+ Drift[iPos] = new TGeoTrap(dz_sD, 0, phi, sdy, sdx1, sdx2, 0, sdy, sdx1, sdx2, 0);
+ G10[iPos] = new TGeoTrap(dz_sG, 0, phi, sdy, sdx1, sdx2, 0, sdy, sdx1, sdx2, 0);
+ }
+ else { //RPC
+ Drift[iPos] = new TGeoTrap(dz_sD, 0, phi, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
+ G10[iPos] = new TGeoTrap(dz_sG, 0, phi, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
}
- Drift[iPos]->SetName(Form("shStation%02iLayer%i%cModule%03i%cDrift",
- istn,
- ily,
- cside,
- iModule,
- cpos));
- DriftName[iPos] = Form("muchstation%02ilayer%i%cside%03i%ccopperDrift",
- istn + 1,
- ily + 1,
- cside,
- iMod + 1,
- cpos);
- voDrift[iPos] =
- new TGeoVolume(DriftName[iPos], Drift[iPos], copperplate);
+ Drift[iPos]->SetName(Form("shStation%02iLayer%i%cModule%03i%cDrift", istn, ily, cside, iModule, cpos));
+ DriftName[iPos] =
+ Form("muchstation%02ilayer%i%cside%03i%ccopperDrift", istn + 1, ily + 1, cside, iMod + 1, cpos);
+ voDrift[iPos] = new TGeoVolume(DriftName[iPos], Drift[iPos], copperplate);
voDrift[iPos]->SetLineColor(kRed);
- G10[iPos]->SetName(Form("shStation%02iLayer%i%cModule%03i%cG10",
- istn,
- ily,
- cside,
- iModule,
- cpos));
- G10Name[iPos] = Form("muchstation%02ilayer%i%cside%03i%cG10",
- istn + 1,
- ily + 1,
- cside,
- iMod + 1,
- cpos);
+ G10[iPos]->SetName(Form("shStation%02iLayer%i%cModule%03i%cG10", istn, ily, cside, iModule, cpos));
+ G10Name[iPos] = Form("muchstation%02ilayer%i%cside%03i%cG10", istn + 1, ily + 1, cside, iMod + 1, cpos);
voG10[iPos] = new TGeoVolume(G10Name[iPos], G10[iPos], g10plate);
voG10[iPos]->SetLineColor(28);
}
// Define the trapezoidal spacers
- TGeoTrap* shapeFrame =
- new TGeoTrap(sdz, 0, 0, sdy, sdx1, sdx2, 0, sdy, sdx1, sdx2, 0);
- shapeFrame->SetName(
- Form("shStation%02iLayer%i%cModule%03iFullFrameNoHole",
- istn,
- ily,
- cside,
- iModule));
+ TGeoTrap* shapeFrame = new TGeoTrap(sdz, 0, 0, sdy, sdx1, sdx2, 0, sdy, sdx1, sdx2, 0);
+ shapeFrame->SetName(Form("shStation%02iLayer%i%cModule%03iFullFrameNoHole", istn, ily, cside, iModule));
TString expression;
if (istn < 2) { //GEM
expression = Form("shStation%02iLayer%i%cModule%03iFullFrameNoHole-"
"shStation%02iLayer%i%cModule%03iActiveNoHole",
- istn,
- ily,
- cside,
- iModule,
- istn,
- ily,
- cside,
- iModule);
- } else { //RPC
+ istn, ily, cside, iModule, istn, ily, cside, iModule);
+ }
+ else { //RPC
expression = Form("shStation%02iLayer%i%cModule%03iFullFrameNoHole-"
"shStation%02iLayer%i%cModule%03iFrameNoHole",
- istn,
- ily,
- cside,
- iModule,
- istn,
- ily,
- cside,
- iModule);
+ istn, ily, cside, iModule, istn, ily, cside, iModule);
}
TGeoCompositeShape* shFrame = new TGeoCompositeShape(
- Form("shStation%02iLayer%i%cModule%03iFinalFrameNoHole",
- istn,
- ily,
- cside,
- iModule),
- expression);
- TString frameName = Form("muchstation%02ilayer%i%cframe%03i",
- istn + 1,
- ily + 1,
- cside,
- iMod + 1);
+ Form("shStation%02iLayer%i%cModule%03iFinalFrameNoHole", istn, ily, cside, iModule), expression);
+ TString frameName = Form("muchstation%02ilayer%i%cframe%03i", istn + 1, ily + 1, cside, iMod + 1);
TGeoVolume* voFrame;
if (istn < 2) {
- voFrame =
- new TGeoVolume(frameName, shFrame, noryl); // add a name to the frame
- } else {
- voFrame = new TGeoVolume(
- frameName, shFrame, SpacerRpc); // add a name to the frame
+ voFrame = new TGeoVolume(frameName, shFrame, noryl); // add a name to the frame
+ }
+ else {
+ voFrame = new TGeoVolume(frameName, shFrame, SpacerRpc); // add a name to the frame
}
voFrame->SetLineColor(kMagenta);
// Calculate the phi angle of the sector where it has to be placed
- Double_t angle =
- 180. / TMath::Pi() * phi; // convert angle phi from rad to deg
+ Double_t angle = 180. / TMath::Pi() * phi; // convert angle phi from rad to deg
TGeoRotation* r2 = new TGeoRotation("r2");
diff --git a/macro/much/geometry/create_MUCH_geometry_v21a_jpsi.C b/macro/much/geometry/create_MUCH_geometry_v21a_jpsi.C
index 582ec5191d..2e5b103a32 100644
--- a/macro/much/geometry/create_MUCH_geometry_v21a_jpsi.C
+++ b/macro/much/geometry/create_MUCH_geometry_v21a_jpsi.C
@@ -17,143 +17,135 @@
// in root all sizes are given in cm
-
-#include "TSystem.h"
+#include "TClonesArray.h"
+#include "TDatime.h"
+#include "TFile.h"
+#include "TGeoBBox.h"
+#include "TGeoCompositeShape.h"
+#include "TGeoCone.h"
#include "TGeoManager.h"
-#include "TGeoVolume.h"
#include "TGeoMaterial.h"
+#include "TGeoMatrix.h"
#include "TGeoMedium.h"
#include "TGeoPgon.h"
-#include "TGeoMatrix.h"
-#include "TGeoCompositeShape.h"
-#include "TGeoXtru.h"
-#include "TGeoCone.h"
-#include "TGeoBBox.h"
#include "TGeoTube.h"
-#include "TFile.h"
-#include "TString.h"
+#include "TGeoVolume.h"
+#include "TGeoXtru.h"
#include "TList.h"
-#include "TRandom3.h"
-#include "TDatime.h"
-#include "TClonesArray.h"
-
-#include "TObjArray.h"
-#include "TFile.h"
#include "TMath.h"
+#include "TObjArray.h"
+#include "TRandom3.h"
+#include "TString.h"
+#include "TSystem.h"
-#include <iostream>
-#include <fstream>
#include <cassert>
+#include <fstream>
+#include <iostream>
#include <stdexcept>
-
// Name of output file with geometry
-const TString tagVersion = "_v21a";
-const TString subVersion = "_sis100_1m_jpsi";
-const TString geoVersion = "much";// + tagVersion + subVersion;
-const TString FileNameSim = geoVersion + tagVersion + subVersion+".geo.root";
-const TString FileNameGeo = geoVersion + tagVersion + subVersion+"_geo.root";
+const TString tagVersion = "_v21a";
+const TString subVersion = "_sis100_1m_jpsi";
+const TString geoVersion = "much"; // + tagVersion + subVersion;
+const TString FileNameSim = geoVersion + tagVersion + subVersion + ".geo.root";
+const TString FileNameGeo = geoVersion + tagVersion + subVersion + "_geo.root";
//const TString FileNameInfo = geoVersion + tagVersion + subVersion+".geo.info";
// Names of the different used materials which are used to build the modules
// The materials are defined in the global media.geo file
-const TString KeepingVolumeMedium = "air";
-const TString L= "MUCHlead";
-const TString I= "MUCHiron";
-const TString I_fifth = "MUCHironfifth";
-const TString activemedium="MUCHGEMmixture";
-const TString spacermedium="MUCHnoryl";
-const TString LDcarbon = "MUCHcarbonLD"; //Low density Carbon for Ist Abs
-const TString Concrete = "MUCHconcrete"; //Concrete for Ist Abs
-const TString Al = "MUCHaluminium"; //Al for cooling & support purpose
-const TString copper = "MUCHcopper";
-const TString g10= "MUCHG10";
-const TString RPCm= "MUCHRPCgas";
-const TString RPCg= "MUCHRPCglass";
-const TString Kapton = "MUCHkapton";
+const TString KeepingVolumeMedium = "air";
+const TString L = "MUCHlead";
+const TString I = "MUCHiron";
+const TString I_fifth = "MUCHironfifth";
+const TString activemedium = "MUCHGEMmixture";
+const TString spacermedium = "MUCHnoryl";
+const TString LDcarbon = "MUCHcarbonLD"; //Low density Carbon for Ist Abs
+const TString Concrete = "MUCHconcrete"; //Concrete for Ist Abs
+const TString Al = "MUCHaluminium"; //Al for cooling & support purpose
+const TString copper = "MUCHcopper";
+const TString g10 = "MUCHG10";
+const TString RPCm = "MUCHRPCgas";
+const TString RPCg = "MUCHRPCglass";
+const TString Kapton = "MUCHkapton";
// Universal input parameters
-Double_t fMuchZ1 =125.0; // MuchCave Zin position [cm]
-Double_t fAcceptanceTanMin = 0.1; // Acceptance tangent min
-Double_t fAcceptanceTanMax = 0.466; // Acceptance tangent max
+Double_t fMuchZ1 = 125.0; // MuchCave Zin position [cm]
+Double_t fAcceptanceTanMin = 0.1; // Acceptance tangent min
+Double_t fAcceptanceTanMax = 0.466; // Acceptance tangent max
//************************************************************
- // Input parameters for absorbers
+// Input parameters for absorbers
//***********************************************************
-const Int_t fAbs = 5;
-const Int_t fNabs = 7; // Number of absorber pieces
-TString AbsMaterial[7]={"LD Graphite","LD Graphite","Concrete","Iron","Iron","Iron","Iron_fifth"};
+const Int_t fAbs = 5;
+const Int_t fNabs = 7; // Number of absorber pieces
+TString AbsMaterial[7] = {"LD Graphite", "LD Graphite", "Concrete", "Iron", "Iron", "Iron", "Iron_fifth"};
// Absorber Zin position [cm] in the cave reference frame
-Double_t fAbsorberZ1[7]={0, 16, 28, 90, 140, 190, 252};
+Double_t fAbsorberZ1[7] = {0, 16, 28, 90, 140, 190, 252};
// Absorber thickness [cm]
-Double_t fAbsorberLz[7]= {16, 12, 30, 20, 20, 30, 100};
-Double_t safetyrad[7]={0.0, 0.0, 0.0,30.0,30.0,30.0, 30.0};
+Double_t fAbsorberLz[7] = {16, 12, 30, 20, 20, 30, 100};
+Double_t safetyrad[7] = {0.0, 0.0, 0.0, 30.0, 30.0, 30.0, 30.0};
//..........All the dimensions of the 5th absorber is mentioned here...............
TVector3 posfifthabsorber;
-Int_t noofblocks=8;
-Double_t open_x=110.;
-Double_t open_y=30.0;
-Double_t open_z=50.0;
-Double_t sh_x=30.0;
-Double_t sh_y=30.0;
-Double_t sh_z = 50.;
-Double_t block_x=250.0;
-Double_t block_y=27.5;
-Double_t block_z=50.0;
-Double_t distance_between_blocks=0.3;
-Double_t gap_al_iron= 0.3;
+Int_t noofblocks = 8;
+Double_t open_x = 110.;
+Double_t open_y = 30.0;
+Double_t open_z = 50.0;
+Double_t sh_x = 30.0;
+Double_t sh_y = 30.0;
+Double_t sh_z = 50.;
+Double_t block_x = 250.0;
+Double_t block_y = 27.5;
+Double_t block_z = 50.0;
+Double_t distance_between_blocks = 0.3;
+Double_t gap_al_iron = 0.3;
//........................
// Input parameters for MUCH stations
//********************************************
-const Int_t fNst = 4; // Number of stations
+const Int_t fNst = 4; // Number of stations
// Sector-type module parameters
// Number of sectors per layer (should be even for symmetry)
// Needs to be fixed with actual numbers
Int_t fNSectorsPerLayer[4] = {16, 20, 18, 20};
//Double_t fRpcGlassDz[4] = {0.0,0.0,0.2,0.2}; //Rpc Glass thickness [cm]
-Double_t fSpacerR = 2.0; // Spacer width in R [cm]
-Double_t fSpacerPhi = 2.0; // Spacer width in Phi [cm]
-Double_t fOverlapR = 2.0; // Overlap in R direction [cm]
+Double_t fSpacerR = 2.0; // Spacer width in R [cm]
+Double_t fSpacerPhi = 2.0; // Spacer width in Phi [cm]
+Double_t fOverlapR = 2.0; // Overlap in R direction [cm]
// Station Zceneter [cm] in the cave reference frame
-Double_t fStationZ0[4]={75,125,175,235};
-Int_t fNlayers[4]={3,3,3,3}; // Number of layers
-Int_t fDetType[4]={3,3,4,4}; // Detector type
-Double_t fLayersDz[4]={10,10,10,10};
-
-
-
-
+Double_t fStationZ0[4] = {75, 125, 175, 235};
+Int_t fNlayers[4] = {3, 3, 3, 3}; // Number of layers
+Int_t fDetType[4] = {3, 3, 4, 4}; // Detector type
+Double_t fLayersDz[4] = {10, 10, 10, 10};
// Input parameters for beam pipe shielding
// spans from 2.9 degree to 5.1 degree
//Inner radius is tan(2.9) + 2 cm, extra 20 mm for clamp connection
-const Int_t fNshs=6;
-TString ShMaterial[6]={"Al","Pb","Al","Al","Al", "Al"};
-Double_t fShieldZin[6]={125,153.0, 215.0, 265.0, 315.0, 377.0};
-Double_t fShieldLz[6]={28, 30, 20, 20, 30, 100};
-Double_t fShield_AcceptanceTanMin[6] = {0.043,0.043, 0.051,0.051,0.051, 0.051}; // Acceptance tangent min for shield
-Double_t fShield_AcceptanceTanMax[6] = {0.1,0.1,0.1,0.1,0.1,0.1}; // Acceptance tangent max for shield
+const Int_t fNshs = 6;
+TString ShMaterial[6] = {"Al", "Pb", "Al", "Al", "Al", "Al"};
+Double_t fShieldZin[6] = {125, 153.0, 215.0, 265.0, 315.0, 377.0};
+Double_t fShieldLz[6] = {28, 30, 20, 20, 30, 100};
+Double_t fShield_AcceptanceTanMin[6] = {0.043, 0.043, 0.051, 0.051, 0.051, 0.051}; // Acceptance tangent min for shield
+Double_t fShield_AcceptanceTanMax[6] = {0.1, 0.1, 0.1, 0.1, 0.1, 0.1}; // Acceptance tangent max for shield
//***********************************************************
// some global variables
-TGeoManager* gGeoMan = NULL; // Pointer to TGeoManager instance
-TGeoVolume* gModules[fNabs]; // Global storage for module types
-TGeoVolume* gModules_shield[fNshs]; // Global storage for module types
-TGeoVolume* gModules_station[fNst]; // Global storage for module types
+TGeoManager* gGeoMan = NULL; // Pointer to TGeoManager instance
+TGeoVolume* gModules[fNabs]; // Global storage for module types
+TGeoVolume* gModules_shield[fNshs]; // Global storage for module types
+TGeoVolume* gModules_station[fNst]; // Global storage for module types
// Forward declarations
void create_materials_from_media_file();
@@ -166,128 +158,131 @@ TGeoVolume* CreateRpcLayers(int istn, int ily);
fstream infoFile;
-void create_MUCH_geometry_v21a_jpsi() {
-
-
+void create_MUCH_geometry_v21a_jpsi()
+{
+
+
// ------- Open info file -----------------------------------------------
TString infoFileName = FileNameSim;
infoFileName.ReplaceAll("root", "info");
infoFile.open(infoFileName.Data(), fstream::out);
infoFile << "MUCH geometry created with create_MUCH_geometry_v21a_lmvm.C" << endl << endl;
- infoFile<<"Global Variables: "<<endl;
- infoFile<<"MuchCave Zin position = "<<fMuchZ1<<" cm "<<endl;
- infoFile<<"Acceptance tangent min = "<<fAcceptanceTanMin<<endl;
- infoFile<<"Acceptance tangent max = "<<fAcceptanceTanMax<<endl;
+ infoFile << "Global Variables: " << endl;
+ infoFile << "MuchCave Zin position = " << fMuchZ1 << " cm " << endl;
+ infoFile << "Acceptance tangent min = " << fAcceptanceTanMin << endl;
+ infoFile << "Acceptance tangent max = " << fAcceptanceTanMax << endl;
// Load needed material definition from media.geo file
create_materials_from_media_file();
-
+
// Get the GeoManager for later usage
gGeoMan = (TGeoManager*) gROOT->FindObject("FAIRGeom");
gGeoMan->SetVisLevel(10);
-
+
// Create the top volume
- TGeoBBox* topbox= new TGeoBBox("", 1000., 1000., 2000.);
- TGeoVolume* top = new TGeoVolume("top", topbox, gGeoMan->GetMedium("air"));
+ TGeoBBox* topbox = new TGeoBBox("", 1000., 1000., 2000.);
+ TGeoVolume* top = new TGeoVolume("top", topbox, gGeoMan->GetMedium("air"));
gGeoMan->SetTopVolume(top);
- TString topName=geoVersion+tagVersion+subVersion;
+ TString topName = geoVersion + tagVersion + subVersion;
TGeoVolume* much = new TGeoVolumeAssembly(topName);
top->AddNode(much, 1);
- TGeoVolume *absr = new TGeoVolumeAssembly("absorber");
- much->AddNode(absr,1);
-
- TGeoVolume *shld = new TGeoVolumeAssembly("shield");
- much->AddNode(shld,1);
-
- TGeoVolume *sttn = new TGeoVolumeAssembly("station");
- much->AddNode(sttn,1);
-
- infoFile<<endl;
- infoFile<<" Absorbers "<<endl;
- infoFile<<" -----------"<<endl;
- infoFile<<"Total No. of Absorbers: "<<fAbs<<endl;
- infoFile<<"First abosrber is divided into two halves."<<endl;
- infoFile<<"First half inserted inside the Dipole Magnet."<<endl;
- infoFile<<"Second half is made of Low Density Graphite + Concrete."<<endl;
- infoFile<<"Total No. of Pieces: "<<fNabs<<endl;
- infoFile<<endl;
- infoFile<<"AbsPieces Position[cm] Thickness[cm] Material"<<endl;
- infoFile<<"--------------------------------------------------------------"<<endl;
-
- for (Int_t iabs = 0; iabs <fNabs ; iabs++) {
-
+ TGeoVolume* absr = new TGeoVolumeAssembly("absorber");
+ much->AddNode(absr, 1);
+
+ TGeoVolume* shld = new TGeoVolumeAssembly("shield");
+ much->AddNode(shld, 1);
+
+ TGeoVolume* sttn = new TGeoVolumeAssembly("station");
+ much->AddNode(sttn, 1);
+
+ infoFile << endl;
+ infoFile << " Absorbers " << endl;
+ infoFile << " -----------" << endl;
+ infoFile << "Total No. of Absorbers: " << fAbs << endl;
+ infoFile << "First abosrber is divided into two halves." << endl;
+ infoFile << "First half inserted inside the Dipole Magnet." << endl;
+ infoFile << "Second half is made of Low Density Graphite + Concrete." << endl;
+ infoFile << "Total No. of Pieces: " << fNabs << endl;
+ infoFile << endl;
+ infoFile << "AbsPieces Position[cm] Thickness[cm] Material" << endl;
+ infoFile << "--------------------------------------------------------------" << endl;
+
+ for (Int_t iabs = 0; iabs < fNabs; iabs++) {
+
gModules[iabs] = CreateAbsorbers(iabs);
-
- absr->AddNode(gModules[iabs],iabs);
+
+ absr->AddNode(gModules[iabs], iabs);
}
-
-
- infoFile<<endl;
- infoFile<<" Shielding "<<endl;
- infoFile<<" -----------"<<endl;
- infoFile<<"No. of Shields: "<<fNshs<<endl;
- infoFile<<"Inside the Abs I, Shielding divided into two parts."<<endl;
- infoFile<<endl;
- infoFile<<"Shield No. Z_In[cm] Z_Out[cm] R_In[cm] R_Out[cm] Material"<<endl;
- infoFile<<"--------------------------------------------------------------"<<endl;
- for (Int_t ishi = 0; ishi <fNshs ; ishi++) {
-
- if (ishi<5)
- {
- gModules_shield[ishi] = CreateShields(ishi);
-
- shld->AddNode(gModules_shield[ishi],ishi);
+
+
+ infoFile << endl;
+ infoFile << " Shielding " << endl;
+ infoFile << " -----------" << endl;
+ infoFile << "No. of Shields: " << fNshs << endl;
+ infoFile << "Inside the Abs I, Shielding divided into two parts." << endl;
+ infoFile << endl;
+ infoFile << "Shield No. Z_In[cm] Z_Out[cm] R_In[cm] R_Out[cm] Material" << endl;
+ infoFile << "--------------------------------------------------------------" << endl;
+ for (Int_t ishi = 0; ishi < fNshs; ishi++) {
+
+ if (ishi < 5) {
+ gModules_shield[ishi] = CreateShields(ishi);
+
+ shld->AddNode(gModules_shield[ishi], ishi);
}
-
- else
- {
- gModules_shield[ishi] = CreateShieldfifthabsorber(ishi);
-
- shld->AddNode(gModules_shield[ishi],ishi);
+
+ else {
+ gModules_shield[ishi] = CreateShieldfifthabsorber(ishi);
+
+ shld->AddNode(gModules_shield[ishi], ishi);
}
}
- infoFile<< " 5 377 477 25.00 25.00 Al" <<endl;
- infoFile<<endl;
- infoFile<<" Stations "<<endl;
- infoFile<<" ----------"<<endl;
- infoFile<<"No. of Stations: "<<fNst<<endl;
- infoFile<<"First two stations (1,2) are made up of GEM and last two stations (3,4) are made up of RPC."<<endl;
- infoFile<<"Passive material implemented in GEM modules. Ar:CO2 (70:30) is used as active gas. "<<endl;
- infoFile<<"12 mm thick Al plates are used for support and cooling in the GEM modules."<<endl;
- infoFile<<"2 mm thick Aluminum plates are used for support in the RPC modules behind the active area. 10 mm thick Aluminium at the boundaries as the frame."<<endl;
- infoFile<<"Drift and read-out PCBs (copper coated G10 plates) inserted for realistic material budget for both GEM and RPC modules."<<endl;
- infoFile<<"#Station #Layers Z[cm] #Sectors ActiveLz[cm]"<<endl;
- infoFile<<"--------------------------------------------------------------"<<endl;
- for (Int_t istn = 0; istn < 4; istn++) { // 4 Stations
-
-
+ infoFile << " 5 377 477 25.00 25.00 Al" << endl;
+ infoFile << endl;
+ infoFile << " Stations " << endl;
+ infoFile << " ----------" << endl;
+ infoFile << "No. of Stations: " << fNst << endl;
+ infoFile << "First two stations (1,2) are made up of GEM and last two stations (3,4) are made up of RPC." << endl;
+ infoFile << "Passive material implemented in GEM modules. Ar:CO2 (70:30) is used as active gas. " << endl;
+ infoFile << "12 mm thick Al plates are used for support and cooling in the GEM modules." << endl;
+ infoFile << "2 mm thick Aluminum plates are used for support in the RPC modules behind the active area. 10 mm thick "
+ "Aluminium at the boundaries as the frame."
+ << endl;
+ infoFile << "Drift and read-out PCBs (copper coated G10 plates) inserted for realistic material budget for both GEM "
+ "and RPC modules."
+ << endl;
+ infoFile << "#Station #Layers Z[cm] #Sectors ActiveLz[cm]" << endl;
+ infoFile << "--------------------------------------------------------------" << endl;
+ for (Int_t istn = 0; istn < 4; istn++) { // 4 Stations
+
+
gModules_station[istn] = CreateStations(istn);
-
- sttn->AddNode(gModules_station[istn],istn);
+
+ sttn->AddNode(gModules_station[istn], istn);
}
-
+
gGeoMan->CloseGeometry();
gGeoMan->CheckOverlaps(0.0000001);
gGeoMan->PrintOverlaps();
gGeoMan->CheckOverlaps(0.0001, "s");
gGeoMan->PrintOverlaps();
-
-
- much->Export(FileNameSim); // an alternative way of writing the much
-
- TFile* outfile = new TFile(FileNameSim, "UPDATE");
+
+
+ much->Export(FileNameSim); // an alternative way of writing the much
+
+ TFile* outfile = new TFile(FileNameSim, "UPDATE");
TGeoTranslation* much_placement = new TGeoTranslation("much_trans", 0., 0., 0.);
much_placement->Write();
outfile->Close();
-
- outfile = new TFile(FileNameGeo,"RECREATE");
+
+ outfile = new TFile(FileNameGeo, "RECREATE");
gGeoMan->Write(); // use this if you want GeoManager format in the output
outfile->Close();
-
+
top->Draw("ogl");
infoFile.close();
}
@@ -295,15 +290,15 @@ void create_MUCH_geometry_v21a_jpsi() {
void create_materials_from_media_file()
{
// Use the FairRoot geometry interface to load the media which are already defined
- FairGeoLoader* geoLoad = new FairGeoLoader("TGeo", "FairGeoLoader");
+ FairGeoLoader* geoLoad = new FairGeoLoader("TGeo", "FairGeoLoader");
FairGeoInterface* geoFace = geoLoad->getGeoInterface();
- TString geoPath = gSystem->Getenv("VMCWORKDIR");
- TString geoFile = geoPath + "/geometry/media.geo";
+ TString geoPath = gSystem->Getenv("VMCWORKDIR");
+ TString geoFile = geoPath + "/geometry/media.geo";
geoFace->setMediaFile(geoFile);
geoFace->readMedia();
// Read the required media and create them in the GeoManager
- FairGeoMedia* geoMedia = geoFace->getMedia();
+ FairGeoMedia* geoMedia = geoFace->getMedia();
FairGeoBuilder* geoBuild = geoLoad->getGeoBuilder();
FairGeoMedium* air = geoMedia->getMedium(KeepingVolumeMedium);
@@ -314,7 +309,7 @@ void create_materials_from_media_file()
FairGeoMedium* MUCHironfifth = geoMedia->getMedium(I_fifth);
geoBuild->createMedium(MUCHironfifth);
-
+
FairGeoMedium* MUCHlead = geoMedia->getMedium(L);
geoBuild->createMedium(MUCHlead);
@@ -333,10 +328,10 @@ void create_materials_from_media_file()
FairGeoMedium* MUCHconcrete = geoMedia->getMedium(Concrete);
geoBuild->createMedium(MUCHconcrete);
- FairGeoMedium* copperplate = geoMedia->getMedium(copper);
+ FairGeoMedium* copperplate = geoMedia->getMedium(copper);
geoBuild->createMedium(copperplate);
- FairGeoMedium* g10plate = geoMedia->getMedium(g10); //G10
+ FairGeoMedium* g10plate = geoMedia->getMedium(g10); //G10
geoBuild->createMedium(g10plate);
FairGeoMedium* RPCmedium = geoMedia->getMedium(RPCm);
@@ -350,772 +345,771 @@ void create_materials_from_media_file()
}
+TGeoVolume* CreateShields(int ish)
+{
+
-TGeoVolume* CreateShields(int ish) {
+ TGeoMedium* iron = gGeoMan->GetMedium(I);
+ TGeoMedium* lead = gGeoMan->GetMedium(L);
+ TGeoMedium* Aluminium = gGeoMan->GetMedium(Al);
-
- TGeoMedium* iron = gGeoMan->GetMedium(I);
- TGeoMedium* lead = gGeoMan->GetMedium(L);
- TGeoMedium* Aluminium= gGeoMan->GetMedium(Al);
-
- TString name = Form("shieldblock%d", ish);
+ TString name = Form("shieldblock%d", ish);
TGeoVolumeAssembly* shieldblock = new TGeoVolumeAssembly(name);
-
- TString conename_sh = Form("conesh_%d",ish);
-
-
- Double_t dz = fShieldLz[ish]/2.0 ;
- Double_t globalZ1 = fShieldZin[ish] ;
- Double_t globalZ2 = fShieldZin[ish] + 2 * dz ;
-
-
-
- Double_t rmin1 = globalZ1 * fShield_AcceptanceTanMin[ish]+2.0;
- Double_t rmax1 = globalZ1 * fShield_AcceptanceTanMax[ish]-0.0001;
- Double_t rmin2 = globalZ2 * fShield_AcceptanceTanMin[ish]+2.0;
- Double_t rmax2 = globalZ2 * fShield_AcceptanceTanMax[ish]-0.0001;
-
- infoFile<<" "<<ish<<"\t\t"<<globalZ1<<"\t "<<globalZ2<<"\t"<<rmin1<<"\t "<<rmax1<<"\t "<<ShMaterial[ish]<<endl;
-
- if(ish==0){
-
- TGeoCone * sh =new TGeoCone(conename_sh,dz, rmin1, rmax1, rmin2, rmax2);
+
+ TString conename_sh = Form("conesh_%d", ish);
+
+
+ Double_t dz = fShieldLz[ish] / 2.0;
+ Double_t globalZ1 = fShieldZin[ish];
+ Double_t globalZ2 = fShieldZin[ish] + 2 * dz;
+
+
+ Double_t rmin1 = globalZ1 * fShield_AcceptanceTanMin[ish] + 2.0;
+ Double_t rmax1 = globalZ1 * fShield_AcceptanceTanMax[ish] - 0.0001;
+ Double_t rmin2 = globalZ2 * fShield_AcceptanceTanMin[ish] + 2.0;
+ Double_t rmax2 = globalZ2 * fShield_AcceptanceTanMax[ish] - 0.0001;
+
+ infoFile << " " << ish << "\t\t" << globalZ1 << "\t " << globalZ2 << "\t" << rmin1 << "\t " << rmax1 << "\t "
+ << ShMaterial[ish] << endl;
+
+ if (ish == 0) {
+
+ TGeoCone* sh = new TGeoCone(conename_sh, dz, rmin1, rmax1, rmin2, rmax2);
TGeoVolume* shield = new TGeoVolume("shield", sh, Aluminium);
shield->SetLineColor(7);
shield->SetTransparency(2);
- TGeoTranslation *sh_trans = new TGeoTranslation("", 0., 0., globalZ1+dz);
- shieldblock->AddNode(shield,ish, sh_trans);
-
-
+ TGeoTranslation* sh_trans = new TGeoTranslation("", 0., 0., globalZ1 + dz);
+ shieldblock->AddNode(shield, ish, sh_trans);
}
-
- if(ish==1){
- TGeoCone * sh =new TGeoCone(conename_sh,dz, rmin1, rmax1, rmin2, rmax2);
+
+ if (ish == 1) {
+ TGeoCone* sh = new TGeoCone(conename_sh, dz, rmin1, rmax1, rmin2, rmax2);
TGeoVolume* shield = new TGeoVolume("shield", sh, lead);
-
+
shield->SetLineColor(kMagenta);
shield->SetTransparency(2);
- TGeoTranslation *sh_trans = new TGeoTranslation("", 0., 0., globalZ1+dz);
- shieldblock->AddNode(shield,ish, sh_trans);
+ TGeoTranslation* sh_trans = new TGeoTranslation("", 0., 0., globalZ1 + dz);
+ shieldblock->AddNode(shield, ish, sh_trans);
+ }
+
+ if (ish > 1) {
+
+
+ TGeoCone* sh = new TGeoCone(conename_sh, dz, rmin1, rmax1, rmin2, rmax2);
+ TGeoVolume* shield = new TGeoVolume("shield", sh, Aluminium);
+
+ shield->SetLineColor(kBlack);
+ shield->SetTransparency(2);
+ TGeoTranslation* sh_trans = new TGeoTranslation("", 0., 0., globalZ1 + dz);
+ shieldblock->AddNode(shield, ish, sh_trans);
}
- if(ish>1)
- {
-
-
- TGeoCone * sh =new TGeoCone(conename_sh,dz, rmin1, rmax1, rmin2, rmax2);
- TGeoVolume* shield = new TGeoVolume("shield", sh, Aluminium);
-
- shield->SetLineColor(kBlack);
- shield->SetTransparency(2);
- TGeoTranslation *sh_trans = new TGeoTranslation("", 0., 0., globalZ1+dz);
- shieldblock->AddNode(shield,ish, sh_trans);
- }
-
return shieldblock;
-
}
-TGeoVolume* CreateShieldfifthabsorber(int ish) {
+TGeoVolume* CreateShieldfifthabsorber(int ish)
+{
- TString BoxName1 = Form("Box1_%d",ish);
- TString BoxName3 = Form("Box3_%d",ish);
- TGeoMedium* Aluminium= gGeoMan->GetMedium(Al);
- TString name = Form("shieldblock%d", ish);
+ TString BoxName1 = Form("Box1_%d", ish);
+ TString BoxName3 = Form("Box3_%d", ish);
+ TGeoMedium* Aluminium = gGeoMan->GetMedium(Al);
+ TString name = Form("shieldblock%d", ish);
TGeoVolumeAssembly* shieldblock = new TGeoVolumeAssembly(name);
- TString tubename_sh = Form("tubesh_%d",ish);
- TGeoBBox * box1 = new TGeoBBox(BoxName1,sh_x+0.00001,sh_y+0.00001,sh_z);
- TGeoTube * sh =new TGeoTube(BoxName3,0.,25.0,50.0);
- TString expression1 = BoxName1+"-"+BoxName3;
- TGeoCompositeShape* shSupport = new TGeoCompositeShape("supportShapeName",expression1);
- TGeoVolume* abs2 = new TGeoVolume("shield", shSupport,Aluminium);
+ TString tubename_sh = Form("tubesh_%d", ish);
+ TGeoBBox* box1 = new TGeoBBox(BoxName1, sh_x + 0.00001, sh_y + 0.00001, sh_z);
+ TGeoTube* sh = new TGeoTube(BoxName3, 0., 25.0, 50.0);
+ TString expression1 = BoxName1 + "-" + BoxName3;
+ TGeoCompositeShape* shSupport = new TGeoCompositeShape("supportShapeName", expression1);
+ TGeoVolume* abs2 = new TGeoVolume("shield", shSupport, Aluminium);
- Double_t dz = fShieldLz[ish]/2.0 ;
+ Double_t dz = fShieldLz[ish] / 2.0;
- Double_t globalZ1 = fShieldZin[ish] ;
+ Double_t globalZ1 = fShieldZin[ish];
+
+ Double_t globalZ2 = fShieldZin[ish] + 2 * dz;
- Double_t globalZ2 = fShieldZin[ish] + 2 * dz ;
-
abs2->SetLineColor(kRed);
- TGeoTranslation *sh_trans = new TGeoTranslation("", 0., 0., globalZ1+dz);
- shieldblock->AddNode(abs2,ish, sh_trans);
-
+ TGeoTranslation* sh_trans = new TGeoTranslation("", 0., 0., globalZ1 + dz);
+ shieldblock->AddNode(abs2, ish, sh_trans);
+
return shieldblock;
}
-TGeoVolume* CreateAbsorbers(int i) {
-
- TGeoMedium* graphite = gGeoMan->GetMedium(LDcarbon);
- TGeoMedium* iron = gGeoMan->GetMedium(I);
+TGeoVolume* CreateAbsorbers(int i)
+{
+
+ TGeoMedium* graphite = gGeoMan->GetMedium(LDcarbon);
+ TGeoMedium* iron = gGeoMan->GetMedium(I);
TGeoMedium* iron_fifth = gGeoMan->GetMedium(I_fifth);
- TGeoMedium* concrete = gGeoMan->GetMedium(Concrete);
- TGeoMedium* Aluminium= gGeoMan->GetMedium(Al);
+ TGeoMedium* concrete = gGeoMan->GetMedium(Concrete);
+ TGeoMedium* Aluminium = gGeoMan->GetMedium(Al);
- TString name = Form("absblock%d", i);
+ TString name = Form("absblock%d", i);
TGeoVolumeAssembly* absblock = new TGeoVolumeAssembly(name);
- TString pipename = Form("beampipe_%d",i);
- TString conename = Form("cone_%d",i);
- TString BoxName = Form("Box_%d",i);
- TString BoxName1 = Form("Box1_%d",i);
- TString BoxName2 = Form("Box2_%d",i);
- TString supportShapeName = Form("Support_%d",i);
- TString TrapName = Form("Trap_%d",i);
-
- Double_t dz = fAbsorberLz[i]/2.0 ;
- Double_t globalZ1 = fAbsorberZ1[i] + fMuchZ1;
+ TString pipename = Form("beampipe_%d", i);
+ TString conename = Form("cone_%d", i);
+ TString BoxName = Form("Box_%d", i);
+ TString BoxName1 = Form("Box1_%d", i);
+ TString BoxName2 = Form("Box2_%d", i);
+ TString supportShapeName = Form("Support_%d", i);
+ TString TrapName = Form("Trap_%d", i);
+
+ Double_t dz = fAbsorberLz[i] / 2.0;
+ Double_t globalZ1 = fAbsorberZ1[i] + fMuchZ1;
Double_t globalPos = globalZ1 + dz;
- Double_t globalZ2 = fAbsorberZ1[i] + 2 * dz + fMuchZ1;
+ Double_t globalZ2 = fAbsorberZ1[i] + 2 * dz + fMuchZ1;
Double_t rmin1 = globalZ1 * fAcceptanceTanMin;
Double_t rmin2 = globalZ2 * fAcceptanceTanMin;
Double_t rmax1 = globalZ1 * fAcceptanceTanMax + safetyrad[i];
- Double_t rmax2 = globalZ2 * fAcceptanceTanMax + safetyrad[i];//
+ Double_t rmax2 = globalZ2 * fAcceptanceTanMax + safetyrad[i]; //
- infoFile<<" "<<i+1<<"\t\t"<<globalPos<<"\t\t"<<2*dz<<"\t\t"<<AbsMaterial[i]<<endl;
+ infoFile << " " << i + 1 << "\t\t" << globalPos << "\t\t" << 2 * dz << "\t\t" << AbsMaterial[i] << endl;
// 1st part of 1st absorber LD Carbon
//dimensions are hardcoded
- if(i==0)
- {
- // printf("absorber %d \n",i);
-
- TGeoTrd2 * trap = new TGeoTrd2(TrapName,70.0,70.0,46.0,71.0,dz);
- TGeoCone * tube = new TGeoCone(pipename,dz+0.001,0.,rmin1,0.,rmin2);
- TString expression = TrapName +"-"+pipename;
- TGeoCompositeShape* shSupport = new TGeoCompositeShape(supportShapeName,expression);
- TGeoVolume* abs0 = new TGeoVolume("absorber", shSupport, graphite);
- abs0->SetLineColor(kRed);
- abs0->SetTransparency(0);
- TGeoTranslation *abs0_trans = new TGeoTranslation("", 0., 0., globalZ1+dz);
- absblock->AddNode(abs0, i, abs0_trans);
- }
-
+ if (i == 0) {
+ // printf("absorber %d \n",i);
+
+ TGeoTrd2* trap = new TGeoTrd2(TrapName, 70.0, 70.0, 46.0, 71.0, dz);
+ TGeoCone* tube = new TGeoCone(pipename, dz + 0.001, 0., rmin1, 0., rmin2);
+ TString expression = TrapName + "-" + pipename;
+ TGeoCompositeShape* shSupport = new TGeoCompositeShape(supportShapeName, expression);
+ TGeoVolume* abs0 = new TGeoVolume("absorber", shSupport, graphite);
+ abs0->SetLineColor(kRed);
+ abs0->SetTransparency(0);
+ TGeoTranslation* abs0_trans = new TGeoTranslation("", 0., 0., globalZ1 + dz);
+ absblock->AddNode(abs0, i, abs0_trans);
+ }
// 2rd part of 1st absorber box (LD Carbon)
- if(i==1)
- {
- // printf("absorber %d \n",i);
- TGeoBBox * box = new TGeoBBox(BoxName,130.0,125.0,dz);
- TGeoCone * tube = new TGeoCone(pipename,dz+0.001,0.,rmin1,0.,rmin2);
- TString expression = BoxName +"-"+pipename;
- TGeoCompositeShape* shSupport = new TGeoCompositeShape(supportShapeName,expression);
-
- TGeoVolume* abs1 = new TGeoVolume("absorber", shSupport, graphite);
- abs1->SetLineColor(kRed);
- abs1->SetTransparency(0);
- TGeoTranslation *abs1_trans = new TGeoTranslation("", 0., 0., globalZ1+dz);
- absblock->AddNode(abs1, i, abs1_trans);
-
-
- }
+ if (i == 1) {
+ // printf("absorber %d \n",i);
+ TGeoBBox* box = new TGeoBBox(BoxName, 130.0, 125.0, dz);
+ TGeoCone* tube = new TGeoCone(pipename, dz + 0.001, 0., rmin1, 0., rmin2);
+ TString expression = BoxName + "-" + pipename;
+ TGeoCompositeShape* shSupport = new TGeoCompositeShape(supportShapeName, expression);
+
+ TGeoVolume* abs1 = new TGeoVolume("absorber", shSupport, graphite);
+ abs1->SetLineColor(kRed);
+ abs1->SetTransparency(0);
+ TGeoTranslation* abs1_trans = new TGeoTranslation("", 0., 0., globalZ1 + dz);
+ absblock->AddNode(abs1, i, abs1_trans);
+ }
// 3th part of 1st absorber box (Conc)
- if(i==2)
- {
- // printf("absorber %d \n",i);
- TGeoBBox * box = new TGeoBBox(BoxName,130.0,125.0,dz);
- TGeoCone * tube = new TGeoCone(pipename,dz+0.001,0.,rmin1,0.,rmin2);
- TString expression = BoxName +"-"+pipename;
- TGeoCompositeShape* shSupport = new TGeoCompositeShape(supportShapeName,expression);
-
- TGeoVolume* abs1 = new TGeoVolume("absorber", shSupport, concrete);
- abs1->SetLineColor(kViolet);
- abs1->SetTransparency(0);
- TGeoTranslation *abs1_trans = new TGeoTranslation("", 0., 0., globalZ1+dz);
- absblock->AddNode(abs1, i, abs1_trans);
- }
-
+ if (i == 2) {
+ // printf("absorber %d \n",i);
+ TGeoBBox* box = new TGeoBBox(BoxName, 130.0, 125.0, dz);
+ TGeoCone* tube = new TGeoCone(pipename, dz + 0.001, 0., rmin1, 0., rmin2);
+ TString expression = BoxName + "-" + pipename;
+ TGeoCompositeShape* shSupport = new TGeoCompositeShape(supportShapeName, expression);
+
+ TGeoVolume* abs1 = new TGeoVolume("absorber", shSupport, concrete);
+ abs1->SetLineColor(kViolet);
+ abs1->SetTransparency(0);
+ TGeoTranslation* abs1_trans = new TGeoTranslation("", 0., 0., globalZ1 + dz);
+ absblock->AddNode(abs1, i, abs1_trans);
+ }
//rest of the absorbers
- if (i>2 && i<6)
- {
- TGeoBBox * box = new TGeoBBox(BoxName,rmax2,rmax2,dz);
- TGeoCone * tube = new TGeoCone(pipename,dz+0.001,0.,rmin1,0.,rmin2);
- TString expression = BoxName +"-"+pipename;
- TGeoCompositeShape* shSupport = new TGeoCompositeShape(supportShapeName,expression);
+ if (i > 2 && i < 6) {
+ TGeoBBox* box = new TGeoBBox(BoxName, rmax2, rmax2, dz);
+ TGeoCone* tube = new TGeoCone(pipename, dz + 0.001, 0., rmin1, 0., rmin2);
+ TString expression = BoxName + "-" + pipename;
+ TGeoCompositeShape* shSupport = new TGeoCompositeShape(supportShapeName, expression);
- TGeoVolume* abs2 = new TGeoVolume("absorber", shSupport, iron);
+ TGeoVolume* abs2 = new TGeoVolume("absorber", shSupport, iron);
- abs2->SetLineColor(kBlue);
- abs2->SetTransparency(2);
+ abs2->SetLineColor(kBlue);
+ abs2->SetTransparency(2);
- TGeoTranslation *abs_trans = new TGeoTranslation("", 0., 0., globalZ1+dz);
- absblock->AddNode(abs2,i, abs_trans);
-
+ TGeoTranslation* abs_trans = new TGeoTranslation("", 0., 0., globalZ1 + dz);
+ absblock->AddNode(abs2, i, abs_trans);
+ }
+ if (i == 6) {
+ printf("absorber %d \n", i);
+ TGeoBBox* box = new TGeoBBox(BoxName, open_x, open_y, open_z);
+ TGeoBBox* box1 = new TGeoBBox(BoxName1, 0.0000, 0.0000, 0.0000);
+ TGeoBBox* box2 = new TGeoBBox(BoxName2, block_x, block_y, block_z);
+
+ TGeoVolume* abs1 = new TGeoVolume("absorber", box, iron_fifth);
+
+ TGeoTranslation* abs1_trans = new TGeoTranslation("", 140.3, 0., globalZ1 + dz);
+ TGeoTranslation* abs1_trans1 = new TGeoTranslation("", -140.3, 0., globalZ1 + dz);
+
+ TGeoVolume* abs3 = new TGeoVolume("absorber", box2, iron_fifth);
+
+ abs1->SetLineColor(kBlue);
+ abs3->SetLineColor(kBlue);
+
+ posfifthabsorber[0] = 0.0;
+ posfifthabsorber[1] = 0.0;
+ posfifthabsorber[2] = globalZ1 + dz;
+
+ absblock->AddNode(abs1, i, abs1_trans);
+ absblock->AddNode(abs1, i, abs1_trans1);
+
+ for (int i = 0; i < (noofblocks / 2); i++) {
+ posfifthabsorber[1] = posfifthabsorber[1] + sh_y + block_y + gap_al_iron;
+ TGeoTranslation* trans = new TGeoTranslation("", posfifthabsorber[0], posfifthabsorber[1], posfifthabsorber[2]);
+ TGeoTranslation* trans1 =
+ new TGeoTranslation("", posfifthabsorber[0], -(posfifthabsorber[1]), posfifthabsorber[2]);
+ TGeoHMatrix* pos_fifthabsorber = new TGeoHMatrix("");
+ TGeoHMatrix* pos_fifthabsorber1 = new TGeoHMatrix("");
+ (*pos_fifthabsorber) = (*trans);
+ (*pos_fifthabsorber1) = (*trans1);
+ absblock->AddNode(abs3, i, pos_fifthabsorber);
+ absblock->AddNode(abs3, i, pos_fifthabsorber1);
+ // cout<<"position"<<posfifthabsorber[1]<<endl;
+ posfifthabsorber[1] = posfifthabsorber[1] - sh_y + block_y - gap_al_iron + distance_between_blocks;
}
-
- if(i==6)
- {
- printf("absorber %d \n",i);
- TGeoBBox * box = new TGeoBBox(BoxName,open_x,open_y,open_z);
- TGeoBBox * box1 = new TGeoBBox(BoxName1,0.0000,0.0000,0.0000);
- TGeoBBox * box2 =new TGeoBBox(BoxName2,block_x,block_y,block_z);
-
- TGeoVolume* abs1 = new TGeoVolume("absorber", box, iron_fifth);
-
- TGeoTranslation *abs1_trans = new TGeoTranslation("",140.3,0., globalZ1+dz);
- TGeoTranslation *abs1_trans1 = new TGeoTranslation("",-140.3,0., globalZ1+dz);
-
- TGeoVolume* abs3 = new TGeoVolume("absorber", box2, iron_fifth);
-
- abs1->SetLineColor(kBlue);
- abs3->SetLineColor(kBlue);
-
- posfifthabsorber[0]=0.0;
- posfifthabsorber[1]=0.0;
- posfifthabsorber[2]=globalZ1+dz;
-
- absblock->AddNode(abs1,i,abs1_trans);
- absblock->AddNode(abs1,i,abs1_trans1);
-
- for (int i=0; i<(noofblocks/2); i++)
- {
- posfifthabsorber[1]=posfifthabsorber[1]+sh_y+block_y+gap_al_iron;
- TGeoTranslation*trans=new TGeoTranslation("",posfifthabsorber[0],posfifthabsorber[1],posfifthabsorber[2]);
- TGeoTranslation*trans1=new TGeoTranslation("",posfifthabsorber[0],-(posfifthabsorber[1]),posfifthabsorber[2]);
- TGeoHMatrix *pos_fifthabsorber = new TGeoHMatrix("");
- TGeoHMatrix *pos_fifthabsorber1 = new TGeoHMatrix("");
- (*pos_fifthabsorber) = (*trans);
- (*pos_fifthabsorber1) = (*trans1);
- absblock->AddNode(abs3, i, pos_fifthabsorber);
- absblock->AddNode (abs3,i,pos_fifthabsorber1);
- // cout<<"position"<<posfifthabsorber[1]<<endl;
- posfifthabsorber[1]=posfifthabsorber[1]-sh_y+block_y-gap_al_iron+distance_between_blocks;
-
-
- }
- }
+ }
return absblock;
-
}
-TGeoVolume * CreateStations(int ist){
-
- TString stationName = Form("muchstation%02i",ist+1);
-
- TGeoVolumeAssembly* station = new TGeoVolumeAssembly(stationName);//, shStation, air);
-
- TGeoVolume* gLayer[4];
-
- for (int ii=0;ii<3;ii++){ // 3 Layers
-
- switch(ist){
- case 0:
- case 1: gLayer[ii] = CreateGemLayers(ist, ii); break;
- case 2:
- case 3: gLayer[ii] = CreateRpcLayers(ist, ii); break;}
-
- station->AddNode(gLayer[ii],ii);
-
+TGeoVolume* CreateStations(int ist)
+{
+
+ TString stationName = Form("muchstation%02i", ist + 1);
+
+ TGeoVolumeAssembly* station = new TGeoVolumeAssembly(stationName); //, shStation, air);
+
+ TGeoVolume* gLayer[4];
+
+ for (int ii = 0; ii < 3; ii++) { // 3 Layers
+
+ switch (ist) {
+ case 0:
+ case 1: gLayer[ii] = CreateGemLayers(ist, ii); break;
+ case 2:
+ case 3: gLayer[ii] = CreateRpcLayers(ist, ii); break;
+ }
+
+ station->AddNode(gLayer[ii], ii);
}
-
+
return station;
}
-TGeoVolume * CreateGemLayers(int istn, int ily){
-
- TString layerName = Form("muchstation%02ilayer%i",istn+1,ily+1);
+TGeoVolume* CreateGemLayers(int istn, int ily)
+{
+
+ TString layerName = Form("muchstation%02ilayer%i", istn + 1, ily + 1);
TGeoVolumeAssembly* volayer = new TGeoVolumeAssembly(layerName);
//GEM Parametrs
- Double_t fSupportDz = 1.2; //1.2 cm Al (cooling + support)
- Double_t fCopperDz = 0.0065; // 65 micron copper
- Double_t fCopperSliceDz = 0.0005; //5 micron copper slices
- Double_t fReadoutPlateDz = 0.3; //3mm G10
- Double_t fDriftPlateDz = 0.3; //3 mm G10
- Double_t fPassiveVolumeDz = 0.2; // 2 mm of Argon
- Double_t fActiveVolumeDz = 0.3; //3mm argon
- Double_t fKaptonDz = 0.005; //50 micron Kapton
- Double_t fFrameDz = 0.2; // 2mm frame
-
-
- Double_t stGlobalZ0 = fStationZ0[istn] + fMuchZ1; //z position of station center (midplane) [cm]
-
-
-
- Double_t layerZ0 = (ily - (fNlayers[istn] - 1) / 2.) * fLayersDz[istn];
+ Double_t fSupportDz = 1.2; //1.2 cm Al (cooling + support)
+ Double_t fCopperDz = 0.0065; // 65 micron copper
+ Double_t fCopperSliceDz = 0.0005; //5 micron copper slices
+ Double_t fReadoutPlateDz = 0.3; //3mm G10
+ Double_t fDriftPlateDz = 0.3; //3 mm G10
+ Double_t fPassiveVolumeDz = 0.2; // 2 mm of Argon
+ Double_t fActiveVolumeDz = 0.3; //3mm argon
+ Double_t fKaptonDz = 0.005; //50 micron Kapton
+ Double_t fFrameDz = 0.2; // 2mm frame
+
+
+ Double_t stGlobalZ0 = fStationZ0[istn] + fMuchZ1; //z position of station center (midplane) [cm]
+
+
+ Double_t layerZ0 = (ily - (fNlayers[istn] - 1) / 2.) * fLayersDz[istn];
Double_t layerGlobalZ0 = layerZ0 + stGlobalZ0;
-
+
//Active Gas distance from layer position
- Double_t ModuleZ = fSupportDz/2. + 2*fCopperDz + fReadoutPlateDz + 3*fPassiveVolumeDz + 6*fCopperSliceDz + 3*fKaptonDz + fActiveVolumeDz/2.;
-
+ Double_t ModuleZ = fSupportDz / 2. + 2 * fCopperDz + fReadoutPlateDz + 3 * fPassiveVolumeDz + 6 * fCopperSliceDz
+ + 3 * fKaptonDz + fActiveVolumeDz / 2.;
+
//ReadOut Plate distance from layer position
- Double_t ReadOutZ = fSupportDz/2. + fCopperDz + fReadoutPlateDz/2.;
-
+ Double_t ReadOutZ = fSupportDz / 2. + fCopperDz + fReadoutPlateDz / 2.;
+
//Copper distance from layer position
- Double_t CooperIZ = fSupportDz/2. + fCopperDz/2.;
- Double_t CooperIIZ = ReadOutZ + fReadoutPlateDz/2. + fCopperDz/2.;
- Double_t CooperIIIZ = ModuleZ + fActiveVolumeDz/2. + fCopperDz/2.;
+ Double_t CooperIZ = fSupportDz / 2. + fCopperDz / 2.;
+ Double_t CooperIIZ = ReadOutZ + fReadoutPlateDz / 2. + fCopperDz / 2.;
+ Double_t CooperIIIZ = ModuleZ + fActiveVolumeDz / 2. + fCopperDz / 2.;
Double_t CooperIVZ = CooperIIIZ + fCopperDz + fDriftPlateDz;
-
+
//Drift Plate distance from layer position
- Double_t DriftZ = ModuleZ + fActiveVolumeDz/2. + fCopperDz + fDriftPlateDz/2.;
-
+ Double_t DriftZ = ModuleZ + fActiveVolumeDz / 2. + fCopperDz + fDriftPlateDz / 2.;
+
//Passive Gas distance from layer position
- Double_t PassivePassiveDz = fPassiveVolumeDz + 2*fCopperSliceDz + fKaptonDz;
- Double_t PassiveGasIZ = ReadOutZ + fReadoutPlateDz/2. + fCopperDz + fPassiveVolumeDz/2.;
- Double_t PassiveGasIIZ = PassiveGasIZ + PassivePassiveDz;
- Double_t PassiveGasIIIZ = PassiveGasIIZ + PassivePassiveDz;
+ Double_t PassivePassiveDz = fPassiveVolumeDz + 2 * fCopperSliceDz + fKaptonDz;
+ Double_t PassiveGasIZ = ReadOutZ + fReadoutPlateDz / 2. + fCopperDz + fPassiveVolumeDz / 2.;
+ Double_t PassiveGasIIZ = PassiveGasIZ + PassivePassiveDz;
+ Double_t PassiveGasIIIZ = PassiveGasIIZ + PassivePassiveDz;
//Kapton distance from layer position
- Double_t KaptonPassiveDz = fPassiveVolumeDz/2. + fCopperSliceDz + fKaptonDz/2.;
- Double_t KaptonIZ = PassiveGasIZ + KaptonPassiveDz;
- Double_t KaptonIIZ = PassiveGasIIZ + KaptonPassiveDz;
- Double_t KaptonIIIZ = PassiveGasIIIZ + KaptonPassiveDz;
-
+ Double_t KaptonPassiveDz = fPassiveVolumeDz / 2. + fCopperSliceDz + fKaptonDz / 2.;
+ Double_t KaptonIZ = PassiveGasIZ + KaptonPassiveDz;
+ Double_t KaptonIIZ = PassiveGasIIZ + KaptonPassiveDz;
+ Double_t KaptonIIIZ = PassiveGasIIIZ + KaptonPassiveDz;
+
//Copper Slice distance from layer position
- Double_t PassiveCopperDz = fPassiveVolumeDz/2. + fCopperSliceDz/2.;
- Double_t KaptonCopperDz = fKaptonDz + fCopperSliceDz;
- Double_t CopperSilceIZ = PassiveGasIZ + PassiveCopperDz;
- Double_t CopperSilceIIZ = CopperSilceIZ + KaptonCopperDz;
+ Double_t PassiveCopperDz = fPassiveVolumeDz / 2. + fCopperSliceDz / 2.;
+ Double_t KaptonCopperDz = fKaptonDz + fCopperSliceDz;
+ Double_t CopperSilceIZ = PassiveGasIZ + PassiveCopperDz;
+ Double_t CopperSilceIIZ = CopperSilceIZ + KaptonCopperDz;
Double_t CopperSilceIIIZ = PassiveGasIIZ + PassiveCopperDz;
- Double_t CopperSilceIVZ = CopperSilceIIIZ + KaptonCopperDz;
- Double_t CopperSilceVZ = PassiveGasIIIZ + PassiveCopperDz;
- Double_t CopperSilceVIZ = CopperSilceVZ + KaptonCopperDz;
+ Double_t CopperSilceIVZ = CopperSilceIIIZ + KaptonCopperDz;
+ Double_t CopperSilceVZ = PassiveGasIIIZ + PassiveCopperDz;
+ Double_t CopperSilceVIZ = CopperSilceVZ + KaptonCopperDz;
+
+ //Set Rmin & Rmax
+ Double_t stDz = ((fNlayers[istn] - 1) * fLayersDz[istn]) / 2. + CooperIVZ + fCopperDz / 2.;
- //Set Rmin & Rmax
- Double_t stDz = ((fNlayers[istn] - 1) * fLayersDz[istn])/2. + CooperIVZ + fCopperDz/2.;
-
Double_t stGlobalZ2 = stGlobalZ0 + stDz;
Double_t stGlobalZ1 = stGlobalZ0 - stDz;
-
+
Double_t rmin = stGlobalZ1 * fAcceptanceTanMin;
Double_t rmax = stGlobalZ2 * fAcceptanceTanMax;
// Module dimention calculation
- Double_t phi0 = TMath::Pi()/fNSectorsPerLayer[istn]; // azimuthal half widh of each module
- Double_t ymin = rmin+fSpacerR;
+ Double_t phi0 = TMath::Pi() / fNSectorsPerLayer[istn]; // azimuthal half widh of each module
+ Double_t ymin = rmin + fSpacerR;
Double_t ymax = rmax;
//define the dimensions of the trapezoidal module
- Double_t dy = (ymax-ymin)/2.; //y (length)
- Double_t dx1 = ymin*TMath::Tan(phi0)+fOverlapR/TMath::Cos(phi0); // large x
- Double_t dx2 = ymax*TMath::Tan(phi0)+fOverlapR/TMath::Cos(phi0); // small x
+ Double_t dy = (ymax - ymin) / 2.; //y (length)
+ Double_t dx1 = ymin * TMath::Tan(phi0) + fOverlapR / TMath::Cos(phi0); // large x
+ Double_t dx2 = ymax * TMath::Tan(phi0) + fOverlapR / TMath::Cos(phi0); // small x
//define the spacer dimensions
- Double_t tg = (dx2-dx1)/2/dy;
- Double_t dd1 = fSpacerPhi*tg;
- Double_t dd2 = fSpacerPhi*sqrt(1+tg*tg);
- Double_t sdx1 = dx1+dd2-dd1-0.1; // 0.1 cm to avoid overlaps
- Double_t sdx2 = dx2+dd2+dd1;
- Double_t sdy = dy+fSpacerR;
-
-
-
- infoFile<<" "<<istn+1<<"\t "<<ily+1<<"\t\t"<<layerGlobalZ0<<"\t"<<fNSectorsPerLayer[istn]<<"\t"<<fActiveVolumeDz<<endl;
-
-
+ Double_t tg = (dx2 - dx1) / 2 / dy;
+ Double_t dd1 = fSpacerPhi * tg;
+ Double_t dd2 = fSpacerPhi * sqrt(1 + tg * tg);
+ Double_t sdx1 = dx1 + dd2 - dd1 - 0.1; // 0.1 cm to avoid overlaps
+ Double_t sdx2 = dx2 + dd2 + dd1;
+ Double_t sdy = dy + fSpacerR;
+
+
+ infoFile << " " << istn + 1 << "\t " << ily + 1 << "\t\t" << layerGlobalZ0 << "\t" << fNSectorsPerLayer[istn]
+ << "\t" << fActiveVolumeDz << endl;
+
+
// Aluminum Plate (Cooling + Support)
- TString supportAlName = Form("shStation%02iSupportAl",istn+1);
- TGeoTube* shSupportAl = new TGeoTube(supportAlName,rmin,rmax,fSupportDz/2.);
-
- TString supportName1 = Form("muchstation%02ilayer%isupportAl",istn+1,ily+1);
- TGeoMedium* coolMat = gGeoMan->GetMedium(Al);
-
- TGeoVolume* voSupport1 = new TGeoVolume(supportName1,shSupportAl,coolMat);
+ TString supportAlName = Form("shStation%02iSupportAl", istn + 1);
+ TGeoTube* shSupportAl = new TGeoTube(supportAlName, rmin, rmax, fSupportDz / 2.);
+
+ TString supportName1 = Form("muchstation%02ilayer%isupportAl", istn + 1, ily + 1);
+ TGeoMedium* coolMat = gGeoMan->GetMedium(Al);
+
+ TGeoVolume* voSupport1 = new TGeoVolume(supportName1, shSupportAl, coolMat);
voSupport1->SetLineColor(kCyan);
-
- TGeoTranslation *support_trans1 = new TGeoTranslation("supportName1", 0,0,layerGlobalZ0);
- volayer->AddNode(voSupport1,0,support_trans1);
-
-
-
+
+ TGeoTranslation* support_trans1 = new TGeoTranslation("supportName1", 0, 0, layerGlobalZ0);
+ volayer->AddNode(voSupport1, 0, support_trans1);
+
+
// Now start adding the GEM modules
- for (Int_t iModule=0; iModule<fNSectorsPerLayer[istn]; iModule++){
-
+ for (Int_t iModule = 0; iModule < fNSectorsPerLayer[istn]; iModule++) {
+
Double_t phi = 2 * phi0 * (iModule + 0.5); // add 0.5 to not overlap with y-axis for left-right layer separation
- Bool_t isBack = iModule%2;
- Char_t cside = (isBack==1) ? 'b' : 'f';
- Int_t iMod = iModule/2;
+ Bool_t isBack = iModule % 2;
+ Char_t cside = (isBack == 1) ? 'b' : 'f';
+ Int_t iMod = iModule / 2;
// correct the x, y positions
Double_t pos[21];
- pos[0] = -(ymin+dy)*sin(phi);
- pos[1] = (ymin+dy)*cos(phi);
+ pos[0] = -(ymin + dy) * sin(phi);
+ pos[1] = (ymin + dy) * cos(phi);
// different z positions for odd/even modules
- pos[2] = (isBack ? 1 : -1)*ModuleZ + layerGlobalZ0; //Active volume & Frame
- pos[3] = (isBack ? 1 : -1)*ReadOutZ + layerGlobalZ0;//Readout G10
- pos[4] = (isBack ? 1 : -1)*CooperIZ + layerGlobalZ0; //Copper I
- pos[5] = (isBack ? 1 : -1)*CooperIIZ + layerGlobalZ0; //Copper II
- pos[6] = (isBack ? 1 : -1)*CooperIIIZ + layerGlobalZ0; //Copper III
- pos[7] = (isBack ? 1 : -1)*CooperIVZ + layerGlobalZ0; //Copper IV
- pos[8] = (isBack ? 1 : -1)*DriftZ + layerGlobalZ0; //Drift G10 PassiveGasIZ
- pos[9] = (isBack ? 1 : -1)*PassiveGasIZ + layerGlobalZ0; // Passive Gas I
- pos[10] = (isBack ? 1 : -1)*PassiveGasIIZ + layerGlobalZ0; // Passive Gas II
- pos[11] = (isBack ? 1 : -1)*PassiveGasIIIZ + layerGlobalZ0; // Passive Gas III
- pos[12] = (isBack ? 1 : -1)*KaptonIZ + layerGlobalZ0; // Kapton I
- pos[13] = (isBack ? 1 : -1)*KaptonIIZ + layerGlobalZ0; // Kapton II
- pos[14] = (isBack ? 1 : -1)*KaptonIIIZ + layerGlobalZ0; // Kapton III
- pos[15] = (isBack ? 1 : -1)*CopperSilceIZ + layerGlobalZ0; // Copper Slice I
- pos[16] = (isBack ? 1 : -1)*CopperSilceIIZ + layerGlobalZ0; // Copper Slice II
- pos[17] = (isBack ? 1 : -1)*CopperSilceIIIZ + layerGlobalZ0; // Copper Slice III
- pos[18] = (isBack ? 1 : -1)*CopperSilceIVZ + layerGlobalZ0; // Copper Slice IV
- pos[19] = (isBack ? 1 : -1)*CopperSilceVZ + layerGlobalZ0; // Copper Slice V
- pos[20] = (isBack ? 1 : -1)*CopperSilceVIZ + layerGlobalZ0; // Copper Slice VI
+ pos[2] = (isBack ? 1 : -1) * ModuleZ + layerGlobalZ0; //Active volume & Frame
+ pos[3] = (isBack ? 1 : -1) * ReadOutZ + layerGlobalZ0; //Readout G10
+ pos[4] = (isBack ? 1 : -1) * CooperIZ + layerGlobalZ0; //Copper I
+ pos[5] = (isBack ? 1 : -1) * CooperIIZ + layerGlobalZ0; //Copper II
+ pos[6] = (isBack ? 1 : -1) * CooperIIIZ + layerGlobalZ0; //Copper III
+ pos[7] = (isBack ? 1 : -1) * CooperIVZ + layerGlobalZ0; //Copper IV
+ pos[8] = (isBack ? 1 : -1) * DriftZ + layerGlobalZ0; //Drift G10 PassiveGasIZ
+ pos[9] = (isBack ? 1 : -1) * PassiveGasIZ + layerGlobalZ0; // Passive Gas I
+ pos[10] = (isBack ? 1 : -1) * PassiveGasIIZ + layerGlobalZ0; // Passive Gas II
+ pos[11] = (isBack ? 1 : -1) * PassiveGasIIIZ + layerGlobalZ0; // Passive Gas III
+ pos[12] = (isBack ? 1 : -1) * KaptonIZ + layerGlobalZ0; // Kapton I
+ pos[13] = (isBack ? 1 : -1) * KaptonIIZ + layerGlobalZ0; // Kapton II
+ pos[14] = (isBack ? 1 : -1) * KaptonIIIZ + layerGlobalZ0; // Kapton III
+ pos[15] = (isBack ? 1 : -1) * CopperSilceIZ + layerGlobalZ0; // Copper Slice I
+ pos[16] = (isBack ? 1 : -1) * CopperSilceIIZ + layerGlobalZ0; // Copper Slice II
+ pos[17] = (isBack ? 1 : -1) * CopperSilceIIIZ + layerGlobalZ0; // Copper Slice III
+ pos[18] = (isBack ? 1 : -1) * CopperSilceIVZ + layerGlobalZ0; // Copper Slice IV
+ pos[19] = (isBack ? 1 : -1) * CopperSilceVZ + layerGlobalZ0; // Copper Slice V
+ pos[20] = (isBack ? 1 : -1) * CopperSilceVIZ + layerGlobalZ0; // Copper Slice VI
//define media
- TGeoMedium* argon = gGeoMan->GetMedium(activemedium); // active medium
- TGeoMedium* noryl = gGeoMan->GetMedium(spacermedium); // frame medium
- TGeoMedium* g10plate = gGeoMan->GetMedium(g10); // G10 medium
- TGeoMedium* copperplate = gGeoMan->GetMedium(copper); // copper
- TGeoMedium* kapton = gGeoMan->GetMedium(Kapton); // Kapton
-
-
+ TGeoMedium* argon = gGeoMan->GetMedium(activemedium); // active medium
+ TGeoMedium* noryl = gGeoMan->GetMedium(spacermedium); // frame medium
+ TGeoMedium* g10plate = gGeoMan->GetMedium(g10); // G10 medium
+ TGeoMedium* copperplate = gGeoMan->GetMedium(copper); // copper
+ TGeoMedium* kapton = gGeoMan->GetMedium(Kapton); // Kapton
+
//Active Volume
- TGeoTrap* shapeActive = new TGeoTrap(fActiveVolumeDz/2.,0,0,dy,dx1,dx2,0,dy,dx1,dx2,0);
+ TGeoTrap* shapeActive = new TGeoTrap(fActiveVolumeDz / 2., 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
shapeActive->SetName(Form("shStation%02iLayer%i%cModule%03iActiveNoHole", istn, ily, cside, iModule));
- TString activeName = Form("muchstation%02ilayer%i%cactive%03igasArgon",istn+1,ily+1,cside,iMod+1);
- TGeoVolume* voActive = new TGeoVolume(activeName,shapeActive,argon);
+ TString activeName = Form("muchstation%02ilayer%i%cactive%03igasArgon", istn + 1, ily + 1, cside, iMod + 1);
+ TGeoVolume* voActive = new TGeoVolume(activeName, shapeActive, argon);
voActive->SetLineColor(kGreen);
-
+
//Frame
- TGeoTrap* shapeFrame = new TGeoTrap(fFrameDz/2.,0,0,sdy,sdx1,sdx2,0,sdy,sdx1,sdx2,0);
+ TGeoTrap* shapeFrame = new TGeoTrap(fFrameDz / 2., 0, 0, sdy, sdx1, sdx2, 0, sdy, sdx1, sdx2, 0);
shapeFrame->SetName(Form("shStation%02iLayer%i%cModule%03iFullFrameNoHole", istn, ily, cside, iModule));
-
- TString expression = Form("shStation%02iLayer%i%cModule%03iFullFrameNoHole-shStation%02iLayer%i%cModule%03iActiveNoHole", istn, ily, cside, iModule, istn, ily, cside, iModule);
-
- TGeoCompositeShape* shFrame = new TGeoCompositeShape(Form("shStation%02iLayer%i%cModule%03iFinalFrameNoHole", istn, ily, cside, iModule), expression);
-
- TString frameName = Form("muchstation%02ilayer%i%cframe%03i",istn+1,ily+1,cside,iModule+1);
-
- TGeoVolume* voFrame = new TGeoVolume(frameName,shFrame,noryl); // add a name to the frame
+
+ TString expression =
+ Form("shStation%02iLayer%i%cModule%03iFullFrameNoHole-shStation%02iLayer%i%cModule%03iActiveNoHole", istn, ily,
+ cside, iModule, istn, ily, cside, iModule);
+
+ TGeoCompositeShape* shFrame = new TGeoCompositeShape(
+ Form("shStation%02iLayer%i%cModule%03iFinalFrameNoHole", istn, ily, cside, iModule), expression);
+
+ TString frameName = Form("muchstation%02ilayer%i%cframe%03i", istn + 1, ily + 1, cside, iModule + 1);
+
+ TGeoVolume* voFrame = new TGeoVolume(frameName, shFrame, noryl); // add a name to the frame
voFrame->SetLineColor(kMagenta);
//Readout Plate
- TGeoTrap* shapeReadOut = new TGeoTrap(fReadoutPlateDz/2.0,0,0,sdy,sdx1,sdx2,0,sdy,sdx1,sdx2,0);
+ TGeoTrap* shapeReadOut = new TGeoTrap(fReadoutPlateDz / 2.0, 0, 0, sdy, sdx1, sdx2, 0, sdy, sdx1, sdx2, 0);
shapeReadOut->SetName(Form("shStation%02iLayer%i%cModule%03iReadOut", istn, ily, cside, iModule));
-
- TString ReadOutName = Form("muchstation%02ilayer%i%cReadOut%03i",istn+1,ily+1,cside,iModule+1);
-
- TGeoVolume* voReadOut = new TGeoVolume(ReadOutName,shapeReadOut,g10plate);
+
+ TString ReadOutName = Form("muchstation%02ilayer%i%cReadOut%03i", istn + 1, ily + 1, cside, iModule + 1);
+
+ TGeoVolume* voReadOut = new TGeoVolume(ReadOutName, shapeReadOut, g10plate);
voReadOut->SetLineColor(2);
-
-
+
+
//4 Copper 65 micron
TGeoTrap* shapeCopper[4];
- TString CopperName[4];
+ TString CopperName[4];
TGeoVolume* voCopper[4];
- for(Int_t iCop =0; iCop <4; iCop++){
- shapeCopper[iCop] = new TGeoTrap(fCopperDz/2.0,0,0,sdy,sdx1,sdx2,0,sdy,sdx1,sdx2,0);
+ for (Int_t iCop = 0; iCop < 4; iCop++) {
+ shapeCopper[iCop] = new TGeoTrap(fCopperDz / 2.0, 0, 0, sdy, sdx1, sdx2, 0, sdy, sdx1, sdx2, 0);
shapeCopper[iCop]->SetName(Form("shStation%02iLayer%i%cModule%03iCopper", istn, ily, cside, iModule));
-
- CopperName[iCop] = Form("muchstation%02ilayer%i%cCopper%iModule%03i",istn+1,ily+1,cside,iCop+1,iModule+1);
-
- voCopper[iCop] = new TGeoVolume(CopperName[iCop],shapeCopper[iCop],copperplate);
+
+ CopperName[iCop] =
+ Form("muchstation%02ilayer%i%cCopper%iModule%03i", istn + 1, ily + 1, cside, iCop + 1, iModule + 1);
+
+ voCopper[iCop] = new TGeoVolume(CopperName[iCop], shapeCopper[iCop], copperplate);
voCopper[iCop]->SetLineColor(3);
}
-
+
//Drift Plate
- TGeoTrap* shapeDrift = new TGeoTrap(fDriftPlateDz/2.0,0,0,sdy,sdx1,sdx2,0,sdy,sdx1,sdx2,0);
+ TGeoTrap* shapeDrift = new TGeoTrap(fDriftPlateDz / 2.0, 0, 0, sdy, sdx1, sdx2, 0, sdy, sdx1, sdx2, 0);
shapeDrift->SetName(Form("shStation%02iLayer%i%cModule%03iDrift", istn, ily, cside, iModule));
-
- TString DriftName = Form("muchstation%02ilayer%i%cDrift%03i",istn+1,ily+1,cside,iModule+1);
-
- TGeoVolume* voDrift = new TGeoVolume(DriftName,shapeDrift,g10plate);
+
+ TString DriftName = Form("muchstation%02ilayer%i%cDrift%03i", istn + 1, ily + 1, cside, iModule + 1);
+
+ TGeoVolume* voDrift = new TGeoVolume(DriftName, shapeDrift, g10plate);
voDrift->SetLineColor(2);
-
+
//Passive Gas (2mm Ar)
TGeoTrap* shapePassiveGas[3];
TString PassiveGasName[3];
TGeoVolume* voPassiveGas[3];
- for(Int_t iPGas =0; iPGas<3; iPGas++){
- shapePassiveGas[iPGas] = new TGeoTrap(fPassiveVolumeDz/2.,0,0,dy,dx1,dx2,0,dy,dx1,dx2,0);
- shapePassiveGas[iPGas]->SetName(Form("shStation%02iLayer%i%cModule%03iPassiveGasNoHole", istn, ily, cside, iModule));
-
- PassiveGasName[iPGas] = Form("muchstation%02ilayer%i%cPassiveGas%iModule%03i",istn+1,ily+1,cside,iPGas,iModule+1);
- voPassiveGas[iPGas] = new TGeoVolume(PassiveGasName[iPGas],shapePassiveGas[iPGas],argon);
+ for (Int_t iPGas = 0; iPGas < 3; iPGas++) {
+ shapePassiveGas[iPGas] = new TGeoTrap(fPassiveVolumeDz / 2., 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
+ shapePassiveGas[iPGas]->SetName(
+ Form("shStation%02iLayer%i%cModule%03iPassiveGasNoHole", istn, ily, cside, iModule));
+
+ PassiveGasName[iPGas] =
+ Form("muchstation%02ilayer%i%cPassiveGas%iModule%03i", istn + 1, ily + 1, cside, iPGas, iModule + 1);
+ voPassiveGas[iPGas] = new TGeoVolume(PassiveGasName[iPGas], shapePassiveGas[iPGas], argon);
voPassiveGas[iPGas]->SetLineColor(kGreen);
}
-
+
//kapton (50 micron)
TGeoTrap* shapeKapton[3];
- TString KaptonName[3];
+ TString KaptonName[3];
TGeoVolume* voKapton[3];
- for(Int_t iCop =0; iCop <3; iCop++){
- shapeKapton[iCop] = new TGeoTrap(fKaptonDz/2.0,0,0,dy,dx1,dx2,0,dy,dx1,dx2,0);
+ for (Int_t iCop = 0; iCop < 3; iCop++) {
+ shapeKapton[iCop] = new TGeoTrap(fKaptonDz / 2.0, 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
shapeKapton[iCop]->SetName(Form("shStation%02iLayer%i%cModule%03iKapton", istn, ily, cside, iModule));
-
- KaptonName[iCop] = Form("muchstation%02ilayer%i%cKapton%iModule%03i",istn+1,ily+1,cside,iCop+1,iModule+1);
-
- voKapton[iCop] = new TGeoVolume(KaptonName[iCop],shapeKapton[iCop],kapton);
+
+ KaptonName[iCop] =
+ Form("muchstation%02ilayer%i%cKapton%iModule%03i", istn + 1, ily + 1, cside, iCop + 1, iModule + 1);
+
+ voKapton[iCop] = new TGeoVolume(KaptonName[iCop], shapeKapton[iCop], kapton);
voKapton[iCop]->SetLineColor(4);
}
-
+
//6 Copper Slice (5 micron) fCopperSliceDz
TGeoTrap* shapeCopperSlice[6];
- TString CopperSliceName[6];
+ TString CopperSliceName[6];
TGeoVolume* voCopperSlice[6];
- for(Int_t iCop =0; iCop <6; iCop++){
- shapeCopperSlice[iCop] = new TGeoTrap(fCopperSliceDz/2.0,0,0,dy,dx1,dx2,0,dy,dx1,dx2,0);
+ for (Int_t iCop = 0; iCop < 6; iCop++) {
+ shapeCopperSlice[iCop] = new TGeoTrap(fCopperSliceDz / 2.0, 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
shapeCopperSlice[iCop]->SetName(Form("shStation%02iLayer%i%cModule%03iCopperSlice", istn, ily, cside, iModule));
-
- CopperSliceName[iCop] = Form("muchstation%02ilayer%i%cCopperSlice%iModule%03i",istn+1,ily+1,cside,iCop+1,iModule+1);
-
- voCopperSlice[iCop] = new TGeoVolume(CopperSliceName[iCop],shapeCopperSlice[iCop],copperplate);
+
+ CopperSliceName[iCop] =
+ Form("muchstation%02ilayer%i%cCopperSlice%iModule%03i", istn + 1, ily + 1, cside, iCop + 1, iModule + 1);
+
+ voCopperSlice[iCop] = new TGeoVolume(CopperSliceName[iCop], shapeCopperSlice[iCop], copperplate);
voCopperSlice[iCop]->SetLineColor(3);
}
-
-
-
+
+
// Calculate the phi angle of the sector where it has to be placed
- Double_t angle = 180. / TMath::Pi() * phi; // convert angle phi from rad to deg
-
-
- TGeoRotation *r2 = new TGeoRotation("r2");
- //rotate in the vertical plane (per to z axis) with angle
- r2->RotateZ(angle);
-
-
- TGeoTranslation *trans[20];
- TGeoHMatrix *incline_mod[20];
-
- for(Int_t i=0; i<19; i++){
- trans[i] = new TGeoTranslation("",pos[0],pos[1],pos[i+2]);
- incline_mod[i] = new TGeoHMatrix("");
- (*incline_mod[i]) = (*trans[i]) * (*r2);
- }
- volayer->AddNode(voActive, iModule, incline_mod[0]); // add active volume
+ Double_t angle = 180. / TMath::Pi() * phi; // convert angle phi from rad to deg
+
+
+ TGeoRotation* r2 = new TGeoRotation("r2");
+ //rotate in the vertical plane (per to z axis) with angle
+ r2->RotateZ(angle);
+
+
+ TGeoTranslation* trans[20];
+ TGeoHMatrix* incline_mod[20];
+
+ for (Int_t i = 0; i < 19; i++) {
+ trans[i] = new TGeoTranslation("", pos[0], pos[1], pos[i + 2]);
+ incline_mod[i] = new TGeoHMatrix("");
+ (*incline_mod[i]) = (*trans[i]) * (*r2);
+ }
+ volayer->AddNode(voActive, iModule, incline_mod[0]); // add active volume
// volayer->AddNode(voFrame, iModule, incline_mod[0]); // add frame // Frame removed from Active gas
- volayer->AddNode(voReadOut, iModule, incline_mod[1]); // add Read Out
- for(int iNode =0; iNode<4; iNode++)volayer->AddNode(voCopper[iNode], iModule, incline_mod[iNode+2]); //add Copper
- volayer->AddNode(voDrift, iModule, incline_mod[6]); // add Drift Out
- for(int iNode =0; iNode<3; iNode++)volayer->AddNode(voPassiveGas[iNode], iModule, incline_mod[iNode+7]); //add Passive Gas
- for(int iNode =0; iNode<3; iNode++)volayer->AddNode(voKapton[iNode], iModule, incline_mod[iNode+10]); //add Kapton
- for(int iNode =0; iNode<6; iNode++)volayer->AddNode(voCopperSlice[iNode], iModule, incline_mod[iNode+13]); //add CopperSlices
-
+ volayer->AddNode(voReadOut, iModule, incline_mod[1]); // add Read Out
+ for (int iNode = 0; iNode < 4; iNode++)
+ volayer->AddNode(voCopper[iNode], iModule, incline_mod[iNode + 2]); //add Copper
+ volayer->AddNode(voDrift, iModule, incline_mod[6]); // add Drift Out
+ for (int iNode = 0; iNode < 3; iNode++)
+ volayer->AddNode(voPassiveGas[iNode], iModule, incline_mod[iNode + 7]); //add Passive Gas
+ for (int iNode = 0; iNode < 3; iNode++)
+ volayer->AddNode(voKapton[iNode], iModule, incline_mod[iNode + 10]); //add Kapton
+ for (int iNode = 0; iNode < 6; iNode++)
+ volayer->AddNode(voCopperSlice[iNode], iModule, incline_mod[iNode + 13]); //add CopperSlices
}
-
+
return volayer;
}
+TGeoVolume* CreateRpcLayers(int istn, int ily)
+{
-TGeoVolume * CreateRpcLayers(int istn, int ily){
-
- TString layerName = Form("muchstation%02ilayer%i",istn+1,ily+1);
+ TString layerName = Form("muchstation%02ilayer%i", istn + 1, ily + 1);
TGeoVolumeAssembly* volayer = new TGeoVolumeAssembly(layerName);
-
+
//RPC Parametrs
- Double_t fActiveVolumeDz = 0.2; //2mm RPC
- Double_t fSupportDz = 0.2; //2 mm Al (cooling + support)
- Double_t fRpcGlassDz = 0.2; // 2mm glass
- Double_t fCopperDz = 0.0035; // 35 micron copper
- Double_t fReadoutPlateDz = 0.3; //3mm G10
- Double_t fDriftPlateDz = 0.3; //3 mm G10
- Double_t fFrameDz = 1.0; // 1cm Al frame
-
-
- Double_t stGlobalZ0 = fStationZ0[istn] + fMuchZ1; //z position of station center (midplane) [cm]
-
- Double_t layerZ0 = (ily - (fNlayers[istn] - 1) / 2.) * fLayersDz[istn];
+ Double_t fActiveVolumeDz = 0.2; //2mm RPC
+ Double_t fSupportDz = 0.2; //2 mm Al (cooling + support)
+ Double_t fRpcGlassDz = 0.2; // 2mm glass
+ Double_t fCopperDz = 0.0035; // 35 micron copper
+ Double_t fReadoutPlateDz = 0.3; //3mm G10
+ Double_t fDriftPlateDz = 0.3; //3 mm G10
+ Double_t fFrameDz = 1.0; // 1cm Al frame
+
+
+ Double_t stGlobalZ0 = fStationZ0[istn] + fMuchZ1; //z position of station center (midplane) [cm]
+
+ Double_t layerZ0 = (ily - (fNlayers[istn] - 1) / 2.) * fLayersDz[istn];
Double_t layerGlobalZ0 = layerZ0 + stGlobalZ0;
-
-
+
+
//Active Gas distance from layer position
- Double_t ModuleZ = fSupportDz/2. + fCopperDz + fReadoutPlateDz + fRpcGlassDz + fActiveVolumeDz/2.;
-
+ Double_t ModuleZ = fSupportDz / 2. + fCopperDz + fReadoutPlateDz + fRpcGlassDz + fActiveVolumeDz / 2.;
+
//RPC Glass
- Double_t GlassDz = fRpcGlassDz/2. + fActiveVolumeDz/2.;
- Double_t GlassIZ = ModuleZ - GlassDz;
- Double_t GlassIIZ = ModuleZ + GlassDz;
-
+ Double_t GlassDz = fRpcGlassDz / 2. + fActiveVolumeDz / 2.;
+ Double_t GlassIZ = ModuleZ - GlassDz;
+ Double_t GlassIIZ = ModuleZ + GlassDz;
+
//ReadOut Plate distance from layer position
- Double_t ReadOutZ = fSupportDz/2. + fCopperDz + fReadoutPlateDz/2.;
+ Double_t ReadOutZ = fSupportDz / 2. + fCopperDz + fReadoutPlateDz / 2.;
//Drift Plate distance from layer position
- Double_t DriftZ = GlassIIZ + fRpcGlassDz/2. + fDriftPlateDz/2.;
+ Double_t DriftZ = GlassIIZ + fRpcGlassDz / 2. + fDriftPlateDz / 2.;
//Copper distance from layer position
- Double_t CooperIZ = fSupportDz/2. + fCopperDz/2.;
- Double_t CooperIIZ = DriftZ + fDriftPlateDz/2. + fCopperDz/2.;
+ Double_t CooperIZ = fSupportDz / 2. + fCopperDz / 2.;
+ Double_t CooperIIZ = DriftZ + fDriftPlateDz / 2. + fCopperDz / 2.;
+
-
//Set Rmin & Rmax
- Double_t stDz = ((fNlayers[istn] - 1) * fLayersDz[istn])/2. + CooperIIZ + fCopperDz/2.;
-
+ Double_t stDz = ((fNlayers[istn] - 1) * fLayersDz[istn]) / 2. + CooperIIZ + fCopperDz / 2.;
+
Double_t stGlobalZ2 = stGlobalZ0 + stDz;
Double_t stGlobalZ1 = stGlobalZ0 - stDz;
-
+
Double_t rmin = stGlobalZ1 * fAcceptanceTanMin;
Double_t rmax = stGlobalZ2 * fAcceptanceTanMax;
-
-
+
+
// Module dimention calculation
- Double_t phi0 = TMath::Pi()/fNSectorsPerLayer[istn]; // azimuthal half widh of each module
- Double_t ymin = rmin+fSpacerR;
+ Double_t phi0 = TMath::Pi() / fNSectorsPerLayer[istn]; // azimuthal half widh of each module
+ Double_t ymin = rmin + fSpacerR;
Double_t ymax = rmax;
-
+
//define the dimensions of the trapezoidal module
- Double_t dy = (ymax-ymin)/2.; //y (length)
- Double_t dx1 = ymin*TMath::Tan(phi0)+fOverlapR/TMath::Cos(phi0); // large x
- Double_t dx2 = ymax*TMath::Tan(phi0)+fOverlapR/TMath::Cos(phi0); // small x
-
+ Double_t dy = (ymax - ymin) / 2.; //y (length)
+ Double_t dx1 = ymin * TMath::Tan(phi0) + fOverlapR / TMath::Cos(phi0); // large x
+ Double_t dx2 = ymax * TMath::Tan(phi0) + fOverlapR / TMath::Cos(phi0); // small x
+
//define the spacer dimensions
- Double_t tg = (dx2-dx1)/2/dy;
- Double_t dd1 = fSpacerPhi*tg;
- Double_t dd2 = fSpacerPhi*sqrt(1+tg*tg);
- Double_t sdx1 = dx1+dd2-dd1-0.1; // 0.1 cm to avoid overlaps
- Double_t sdx2 = dx2+dd2+dd1;
- Double_t sdy = dy+fSpacerR;
-
-
-
- infoFile<<" "<<istn+1<<"\t "<<ily+1<<"\t\t"<<layerGlobalZ0<<"\t"<<fNSectorsPerLayer[istn]<<"\t"<<fActiveVolumeDz<<endl;
-
-
+ Double_t tg = (dx2 - dx1) / 2 / dy;
+ Double_t dd1 = fSpacerPhi * tg;
+ Double_t dd2 = fSpacerPhi * sqrt(1 + tg * tg);
+ Double_t sdx1 = dx1 + dd2 - dd1 - 0.1; // 0.1 cm to avoid overlaps
+ Double_t sdx2 = dx2 + dd2 + dd1;
+ Double_t sdy = dy + fSpacerR;
+
+
+ infoFile << " " << istn + 1 << "\t " << ily + 1 << "\t\t" << layerGlobalZ0 << "\t" << fNSectorsPerLayer[istn]
+ << "\t" << fActiveVolumeDz << endl;
+
+
// Aluminum Plate (Cooling + Support)
- TString supportAlName = Form("shStation%02iSupportAl",istn+1);
- TGeoTube* shSupportAl = new TGeoTube(supportAlName,rmin,rmax,fSupportDz/2.);
-
- TString supportName1 = Form("muchstation%02ilayer%isupportAl",istn+1,ily+1);
- TGeoMedium* coolMat = gGeoMan->GetMedium(Al);
-
- TGeoVolume* voSupport1 = new TGeoVolume(supportName1,shSupportAl,coolMat);
+ TString supportAlName = Form("shStation%02iSupportAl", istn + 1);
+ TGeoTube* shSupportAl = new TGeoTube(supportAlName, rmin, rmax, fSupportDz / 2.);
+
+ TString supportName1 = Form("muchstation%02ilayer%isupportAl", istn + 1, ily + 1);
+ TGeoMedium* coolMat = gGeoMan->GetMedium(Al);
+
+ TGeoVolume* voSupport1 = new TGeoVolume(supportName1, shSupportAl, coolMat);
voSupport1->SetLineColor(kCyan);
-
- TGeoTranslation *support_trans1 = new TGeoTranslation("supportName1", 0,0,layerGlobalZ0);
- volayer->AddNode(voSupport1,0,support_trans1);
-
-// Now start adding the GEM modules
- for (Int_t iModule=0; iModule<fNSectorsPerLayer[istn]; iModule++){
-
+
+ TGeoTranslation* support_trans1 = new TGeoTranslation("supportName1", 0, 0, layerGlobalZ0);
+ volayer->AddNode(voSupport1, 0, support_trans1);
+
+ // Now start adding the GEM modules
+ for (Int_t iModule = 0; iModule < fNSectorsPerLayer[istn]; iModule++) {
+
Double_t phi = 2 * phi0 * (iModule + 0.5); // add 0.5 to not overlap with y-axis for left-right layer separation
- Bool_t isBack = iModule%2;
- Char_t cside = (isBack==1) ? 'b' : 'f';
- Int_t iMod = iModule/2;
+ Bool_t isBack = iModule % 2;
+ Char_t cside = (isBack == 1) ? 'b' : 'f';
+ Int_t iMod = iModule / 2;
// correct the x, y positions
Double_t pos[9];
- pos[0] = -(ymin+dy)*sin(phi);
- pos[1] = (ymin+dy)*cos(phi);
+ pos[0] = -(ymin + dy) * sin(phi);
+ pos[1] = (ymin + dy) * cos(phi);
// different z positions for odd/even modules
- pos[2] = (isBack ? 1 : -1)*ModuleZ + layerGlobalZ0; //Active volume & Frame
- pos[3] = (isBack ? 1 : -1)*GlassIZ + layerGlobalZ0; //Glass I
- pos[4] = (isBack ? 1 : -1)*GlassIIZ + layerGlobalZ0; //Glass II
- pos[5] = (isBack ? 1 : -1)*ReadOutZ + layerGlobalZ0; //ReadOut Plate
- pos[6] = (isBack ? 1 : -1)*DriftZ + layerGlobalZ0; //Drift Plate
- pos[7] = (isBack ? 1 : -1)*CooperIZ + layerGlobalZ0; //Copper I Plate
- pos[8] = (isBack ? 1 : -1)*CooperIIZ + layerGlobalZ0; //Copper II Plate
-
-
-
- TGeoMedium* RPCgasmedium = gGeoMan->GetMedium(RPCm); // RPC Gas
- TGeoMedium* RPCglassmat = gGeoMan->GetMedium(RPCg); // RPC Glass
- TGeoMedium* g10plate = gGeoMan->GetMedium(g10); // G10 medium
- TGeoMedium* copperplate = gGeoMan->GetMedium(copper); // copper
- TGeoMedium* AlFrame = gGeoMan->GetMedium(Al); // spacer medium
-
-
+ pos[2] = (isBack ? 1 : -1) * ModuleZ + layerGlobalZ0; //Active volume & Frame
+ pos[3] = (isBack ? 1 : -1) * GlassIZ + layerGlobalZ0; //Glass I
+ pos[4] = (isBack ? 1 : -1) * GlassIIZ + layerGlobalZ0; //Glass II
+ pos[5] = (isBack ? 1 : -1) * ReadOutZ + layerGlobalZ0; //ReadOut Plate
+ pos[6] = (isBack ? 1 : -1) * DriftZ + layerGlobalZ0; //Drift Plate
+ pos[7] = (isBack ? 1 : -1) * CooperIZ + layerGlobalZ0; //Copper I Plate
+ pos[8] = (isBack ? 1 : -1) * CooperIIZ + layerGlobalZ0; //Copper II Plate
+
+
+ TGeoMedium* RPCgasmedium = gGeoMan->GetMedium(RPCm); // RPC Gas
+ TGeoMedium* RPCglassmat = gGeoMan->GetMedium(RPCg); // RPC Glass
+ TGeoMedium* g10plate = gGeoMan->GetMedium(g10); // G10 medium
+ TGeoMedium* copperplate = gGeoMan->GetMedium(copper); // copper
+ TGeoMedium* AlFrame = gGeoMan->GetMedium(Al); // spacer medium
+
+
//Active Volume
- TGeoTrap* shapeActive = new TGeoTrap(fActiveVolumeDz/2.,0,0,dy,dx1,dx2,0,dy,dx1,dx2,0);
+ TGeoTrap* shapeActive = new TGeoTrap(fActiveVolumeDz / 2., 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
shapeActive->SetName(Form("shStation%02iLayer%i%cModule%03iActiveNoHole", istn, ily, cside, iModule));
-
- TString activeName = Form("muchstation%02ilayer%i%cactive%03irpcgas",istn+1,ily+1,cside,iMod+1);
- TGeoVolume* voActive = new TGeoVolume(activeName,shapeActive,RPCgasmedium);
+
+ TString activeName = Form("muchstation%02ilayer%i%cactive%03irpcgas", istn + 1, ily + 1, cside, iMod + 1);
+ TGeoVolume* voActive = new TGeoVolume(activeName, shapeActive, RPCgasmedium);
voActive->SetLineColor(kGreen);
-
+
// RPC Glass
- TGeoTrap* shapeGlass[2];
- TString GlassName[2];
- TGeoVolume* voGlass[2];
- for(Int_t iGlass =0 ; iGlass<2; iGlass++){
- shapeGlass[iGlass] = new TGeoTrap(fRpcGlassDz/2.,0,0,dy,dx1,dx2,0,dy,dx1,dx2,0);
- shapeGlass[iGlass]->SetName(Form("shStation%02iLayer%i%cModule%03iNoHoleGlass", istn, ily, cside, iModule));
- GlassName[iGlass] = Form("muchstation%02ilayer%i%crpcglass%imodule%03i",istn+1,ily+1,cside,iGlass+1,iModule+1);
- voGlass[iGlass] = new TGeoVolume(GlassName[iGlass],shapeGlass[iGlass],RPCglassmat);
- voGlass[iGlass]->SetLineColor(kRed);
- }
-
- //Readout Plate
- TGeoTrap* shapeReadOut = new TGeoTrap(fReadoutPlateDz/2.0,0,0,dy,dx1,dx2,0,dy,dx1,dx2,0);
- shapeReadOut->SetName(Form("shStation%02iLayer%i%cModule%03iReadOut", istn, ily, cside, iModule));
-
- TString ReadOutName = Form("muchstation%02ilayer%i%cReadOut%03i",istn+1,ily+1,cside,iModule+1);
-
- TGeoVolume* voReadOut = new TGeoVolume(ReadOutName,shapeReadOut,g10plate);
- voReadOut->SetLineColor(2);
-
- //Drift Plate
- TGeoTrap* shapeDrift = new TGeoTrap(fDriftPlateDz/2.0,0,0,dy,dx1,dx2,0,dy,dx1,dx2,0);
- shapeDrift->SetName(Form("shStation%02iLayer%i%cModule%03iDrift", istn, ily, cside, iModule));
-
- TString DriftName = Form("muchstation%02ilayer%i%cDrift%03i",istn+1,ily+1,cside,iModule+1);
-
- TGeoVolume* voDrift = new TGeoVolume(DriftName,shapeDrift,g10plate);
- voDrift->SetLineColor(2);
-
-
- //2 Copper 35 micron
- TGeoTrap* shapeCopper[2];
- TString CopperName[2];
- TGeoVolume* voCopper[2];
- for(Int_t iCop =0; iCop <2; iCop++){
- shapeCopper[iCop] = new TGeoTrap(fCopperDz/2.0,0,0,dy,dx1,dx2,0,dy,dx1,dx2,0);
- shapeCopper[iCop]->SetName(Form("shStation%02iLayer%i%cModule%03iCopper", istn, ily, cside, iModule));
-
- CopperName[iCop] = Form("muchstation%02ilayer%i%cCopper%iModule%03i",istn+1,ily+1,cside,iCop+1,iModule+1);
-
- voCopper[iCop] = new TGeoVolume(CopperName[iCop],shapeCopper[iCop],copperplate);
- voCopper[iCop]->SetLineColor(kRed);
- }
-
-
- //Frame (1 cm Al)
-
- TGeoTrap* shape = new TGeoTrap(fFrameDz/2.,0,0,dy,dx1,dx2,0,dy,dx1,dx2,0);
- shape->SetName(Form("shStation%02iLayer%i%cModule%03iFrameNoHole", istn, ily, cside, iModule));
-
- TGeoTrap* shapeFrame = new TGeoTrap(fFrameDz/2.,0,0,sdy,sdx1,sdx2,0,sdy,sdx1,sdx2,0);
- shapeFrame->SetName(Form("shStation%02iLayer%i%cModule%03iFullFrameNoHole", istn, ily, cside, iModule));
-
- TString expression = Form("shStation%02iLayer%i%cModule%03iFullFrameNoHole-shStation%02iLayer%i%cModule%03iFrameNoHole", istn, ily, cside, iModule, istn, ily, cside, iModule);
-
- TGeoCompositeShape* shFrame = new TGeoCompositeShape(Form("shStation%02iLayer%i%cModule%03iFinalFrameNoHole", istn, ily, cside, iModule), expression);
-
- TString frameName = Form("muchstation%02ilayer%i%cframe%03i",istn+1,ily+1,cside,iModule+1);
-
- TGeoVolume* voFrame = new TGeoVolume(frameName,shFrame,AlFrame); // add a name to the frame
- voFrame->SetLineColor(kMagenta);
-
-
-
- // Calculate the phi angle of the sector where it has to be placed
- Double_t angle = 180. / TMath::Pi() * phi; // convert angle phi from rad to deg
-
-
- TGeoRotation *r2 = new TGeoRotation("r2");
- //rotate in the vertical plane (per to z axis) with angle
- r2->RotateZ(angle);
-
-
- TGeoTranslation *trans[8];
- TGeoHMatrix *incline_mod[8];
-
- for(int i=0; i<7; i++){
- trans[i] = new TGeoTranslation("",pos[0],pos[1],pos[i+2]);
-
- incline_mod[i] = new TGeoHMatrix("");
- (*incline_mod[i]) = (*trans[i]) * (*r2);
- }
-
- volayer->AddNode(voActive, iModule, incline_mod[0]); // add active volume
- // volayer->AddNode(voFrame, iModule, incline_mod[0]); // add frame
- for(int iNode =0; iNode<2; iNode++)volayer->AddNode(voGlass[iNode], iModule, incline_mod[iNode+1]); //add Glass
- volayer->AddNode(voReadOut, iModule, incline_mod[3]); // add ReadOut Plate
- volayer->AddNode(voDrift, iModule, incline_mod[4]); // add Drift Plate
- for(int iNode =0; iNode<2; iNode++)volayer->AddNode(voCopper[iNode], iModule, incline_mod[iNode+5]); //add Copper
+ TGeoTrap* shapeGlass[2];
+ TString GlassName[2];
+ TGeoVolume* voGlass[2];
+ for (Int_t iGlass = 0; iGlass < 2; iGlass++) {
+ shapeGlass[iGlass] = new TGeoTrap(fRpcGlassDz / 2., 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
+ shapeGlass[iGlass]->SetName(Form("shStation%02iLayer%i%cModule%03iNoHoleGlass", istn, ily, cside, iModule));
+ GlassName[iGlass] =
+ Form("muchstation%02ilayer%i%crpcglass%imodule%03i", istn + 1, ily + 1, cside, iGlass + 1, iModule + 1);
+ voGlass[iGlass] = new TGeoVolume(GlassName[iGlass], shapeGlass[iGlass], RPCglassmat);
+ voGlass[iGlass]->SetLineColor(kRed);
+ }
+
+ //Readout Plate
+ TGeoTrap* shapeReadOut = new TGeoTrap(fReadoutPlateDz / 2.0, 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
+ shapeReadOut->SetName(Form("shStation%02iLayer%i%cModule%03iReadOut", istn, ily, cside, iModule));
+
+ TString ReadOutName = Form("muchstation%02ilayer%i%cReadOut%03i", istn + 1, ily + 1, cside, iModule + 1);
+
+ TGeoVolume* voReadOut = new TGeoVolume(ReadOutName, shapeReadOut, g10plate);
+ voReadOut->SetLineColor(2);
+
+ //Drift Plate
+ TGeoTrap* shapeDrift = new TGeoTrap(fDriftPlateDz / 2.0, 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
+ shapeDrift->SetName(Form("shStation%02iLayer%i%cModule%03iDrift", istn, ily, cside, iModule));
+
+ TString DriftName = Form("muchstation%02ilayer%i%cDrift%03i", istn + 1, ily + 1, cside, iModule + 1);
+
+ TGeoVolume* voDrift = new TGeoVolume(DriftName, shapeDrift, g10plate);
+ voDrift->SetLineColor(2);
+
+
+ //2 Copper 35 micron
+ TGeoTrap* shapeCopper[2];
+ TString CopperName[2];
+ TGeoVolume* voCopper[2];
+ for (Int_t iCop = 0; iCop < 2; iCop++) {
+ shapeCopper[iCop] = new TGeoTrap(fCopperDz / 2.0, 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
+ shapeCopper[iCop]->SetName(Form("shStation%02iLayer%i%cModule%03iCopper", istn, ily, cside, iModule));
+
+ CopperName[iCop] =
+ Form("muchstation%02ilayer%i%cCopper%iModule%03i", istn + 1, ily + 1, cside, iCop + 1, iModule + 1);
+
+ voCopper[iCop] = new TGeoVolume(CopperName[iCop], shapeCopper[iCop], copperplate);
+ voCopper[iCop]->SetLineColor(kRed);
+ }
+
+
+ //Frame (1 cm Al)
+
+ TGeoTrap* shape = new TGeoTrap(fFrameDz / 2., 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
+ shape->SetName(Form("shStation%02iLayer%i%cModule%03iFrameNoHole", istn, ily, cside, iModule));
+
+ TGeoTrap* shapeFrame = new TGeoTrap(fFrameDz / 2., 0, 0, sdy, sdx1, sdx2, 0, sdy, sdx1, sdx2, 0);
+ shapeFrame->SetName(Form("shStation%02iLayer%i%cModule%03iFullFrameNoHole", istn, ily, cside, iModule));
+
+ TString expression =
+ Form("shStation%02iLayer%i%cModule%03iFullFrameNoHole-shStation%02iLayer%i%cModule%03iFrameNoHole", istn, ily,
+ cside, iModule, istn, ily, cside, iModule);
+
+ TGeoCompositeShape* shFrame = new TGeoCompositeShape(
+ Form("shStation%02iLayer%i%cModule%03iFinalFrameNoHole", istn, ily, cside, iModule), expression);
+
+ TString frameName = Form("muchstation%02ilayer%i%cframe%03i", istn + 1, ily + 1, cside, iModule + 1);
+
+ TGeoVolume* voFrame = new TGeoVolume(frameName, shFrame, AlFrame); // add a name to the frame
+ voFrame->SetLineColor(kMagenta);
+
+
+ // Calculate the phi angle of the sector where it has to be placed
+ Double_t angle = 180. / TMath::Pi() * phi; // convert angle phi from rad to deg
+
+
+ TGeoRotation* r2 = new TGeoRotation("r2");
+ //rotate in the vertical plane (per to z axis) with angle
+ r2->RotateZ(angle);
+
+
+ TGeoTranslation* trans[8];
+ TGeoHMatrix* incline_mod[8];
+
+ for (int i = 0; i < 7; i++) {
+ trans[i] = new TGeoTranslation("", pos[0], pos[1], pos[i + 2]);
+
+ incline_mod[i] = new TGeoHMatrix("");
+ (*incline_mod[i]) = (*trans[i]) * (*r2);
+ }
+
+ volayer->AddNode(voActive, iModule, incline_mod[0]); // add active volume
+ // volayer->AddNode(voFrame, iModule, incline_mod[0]); // add frame
+ for (int iNode = 0; iNode < 2; iNode++)
+ volayer->AddNode(voGlass[iNode], iModule, incline_mod[iNode + 1]); //add Glass
+ volayer->AddNode(voReadOut, iModule, incline_mod[3]); // add ReadOut Plate
+ volayer->AddNode(voDrift, iModule, incline_mod[4]); // add Drift Plate
+ for (int iNode = 0; iNode < 2; iNode++)
+ volayer->AddNode(voCopper[iNode], iModule, incline_mod[iNode + 5]); //add Copper
}
return volayer;
}
diff --git a/macro/much/geometry/create_MUCH_geometry_v21a_lmvm.C b/macro/much/geometry/create_MUCH_geometry_v21a_lmvm.C
index cc1019d213..d9e9f3b63c 100644
--- a/macro/much/geometry/create_MUCH_geometry_v21a_lmvm.C
+++ b/macro/much/geometry/create_MUCH_geometry_v21a_lmvm.C
@@ -17,127 +17,118 @@
// in root all sizes are given in cm
-
-#include "TSystem.h"
+#include "TClonesArray.h"
+#include "TDatime.h"
+#include "TFile.h"
+#include "TGeoBBox.h"
+#include "TGeoCompositeShape.h"
+#include "TGeoCone.h"
#include "TGeoManager.h"
-#include "TGeoVolume.h"
#include "TGeoMaterial.h"
+#include "TGeoMatrix.h"
#include "TGeoMedium.h"
#include "TGeoPgon.h"
-#include "TGeoMatrix.h"
-#include "TGeoCompositeShape.h"
-#include "TGeoXtru.h"
-#include "TGeoCone.h"
-#include "TGeoBBox.h"
#include "TGeoTube.h"
-#include "TFile.h"
-#include "TString.h"
+#include "TGeoVolume.h"
+#include "TGeoXtru.h"
#include "TList.h"
-#include "TRandom3.h"
-#include "TDatime.h"
-#include "TClonesArray.h"
-
-#include "TObjArray.h"
-#include "TFile.h"
#include "TMath.h"
+#include "TObjArray.h"
+#include "TRandom3.h"
+#include "TString.h"
+#include "TSystem.h"
-#include <iostream>
-#include <fstream>
#include <cassert>
+#include <fstream>
+#include <iostream>
#include <stdexcept>
-
// Name of output file with geometry
-const TString tagVersion = "_v21a";
-const TString subVersion = "_sis100_1m_lmvm";
-const TString geoVersion = "much";// + tagVersion + subVersion;
-const TString FileNameSim = geoVersion + tagVersion + subVersion+".geo.root";
-const TString FileNameGeo = geoVersion + tagVersion + subVersion+"_geo.root";
+const TString tagVersion = "_v21a";
+const TString subVersion = "_sis100_1m_lmvm";
+const TString geoVersion = "much"; // + tagVersion + subVersion;
+const TString FileNameSim = geoVersion + tagVersion + subVersion + ".geo.root";
+const TString FileNameGeo = geoVersion + tagVersion + subVersion + "_geo.root";
//const TString FileNameInfo = geoVersion + tagVersion + subVersion+".geo.info";
// Names of the different used materials which are used to build the modules
// The materials are defined in the global media.geo file
-const TString KeepingVolumeMedium = "air";
-const TString L= "MUCHlead";
-const TString I= "MUCHiron";
-const TString activemedium="MUCHGEMmixture";
-const TString spacermedium="MUCHnoryl";
-const TString LDcarbon = "MUCHcarbonLD"; //Low density Carbon for Ist Abs
-const TString Concrete = "MUCHconcrete"; //Concrete for Ist Abs
-const TString Al = "MUCHaluminium"; //Al for cooling & support purpose
-const TString copper = "MUCHcopper";
-const TString g10= "MUCHG10";
-const TString RPCm= "MUCHRPCgas";
-const TString RPCg= "MUCHRPCglass";
-const TString Kapton = "MUCHkapton";
+const TString KeepingVolumeMedium = "air";
+const TString L = "MUCHlead";
+const TString I = "MUCHiron";
+const TString activemedium = "MUCHGEMmixture";
+const TString spacermedium = "MUCHnoryl";
+const TString LDcarbon = "MUCHcarbonLD"; //Low density Carbon for Ist Abs
+const TString Concrete = "MUCHconcrete"; //Concrete for Ist Abs
+const TString Al = "MUCHaluminium"; //Al for cooling & support purpose
+const TString copper = "MUCHcopper";
+const TString g10 = "MUCHG10";
+const TString RPCm = "MUCHRPCgas";
+const TString RPCg = "MUCHRPCglass";
+const TString Kapton = "MUCHkapton";
// Universal input parameters
-Double_t fMuchZ1 =125.0; // MuchCave Zin position [cm]
-Double_t fAcceptanceTanMin = 0.1; // Acceptance tangent min
-Double_t fAcceptanceTanMax = 0.466; // Acceptance tangent max
+Double_t fMuchZ1 = 125.0; // MuchCave Zin position [cm]
+Double_t fAcceptanceTanMin = 0.1; // Acceptance tangent min
+Double_t fAcceptanceTanMax = 0.466; // Acceptance tangent max
//************************************************************
- // Input parameters for absorbers
+// Input parameters for absorbers
//***********************************************************
-const Int_t fAbs = 4;
-const Int_t fNabs = 6; // Number of absorber pieces
-TString AbsMaterial[6]={"LD Graphite","LD Graphite","Concrete","Iron","Iron","Iron"};
+const Int_t fAbs = 4;
+const Int_t fNabs = 6; // Number of absorber pieces
+TString AbsMaterial[6] = {"LD Graphite", "LD Graphite", "Concrete", "Iron", "Iron", "Iron"};
// Absorber Zin position [cm] in the cave reference frame
-Double_t fAbsorberZ1[6]={0, 16, 28, 90, 140, 190};
+Double_t fAbsorberZ1[6] = {0, 16, 28, 90, 140, 190};
// Absorber thickness [cm]
-Double_t fAbsorberLz[6]= {16, 12, 30, 20, 20, 30};
-Double_t safetyrad[6]={0.0, 0.0, 0.0,30.0,30.0,30.0};
-
+Double_t fAbsorberLz[6] = {16, 12, 30, 20, 20, 30};
+Double_t safetyrad[6] = {0.0, 0.0, 0.0, 30.0, 30.0, 30.0};
// Input parameters for MUCH stations
//********************************************
-const Int_t fNst = 4; // Number of stations
+const Int_t fNst = 4; // Number of stations
// Sector-type module parameters
// Number of sectors per layer (should be even for symmetry)
// Needs to be fixed with actual numbers
Int_t fNSectorsPerLayer[4] = {16, 20, 18, 20};
//Double_t fRpcGlassDz[4] = {0.0,0.0,0.2,0.2}; //Rpc Glass thickness [cm]
-Double_t fSpacerR = 2.0; // Spacer width in R [cm]
-Double_t fSpacerPhi = 2.0; // Spacer width in Phi [cm]
-Double_t fOverlapR = 2.0; // Overlap in R direction [cm]
+Double_t fSpacerR = 2.0; // Spacer width in R [cm]
+Double_t fSpacerPhi = 2.0; // Spacer width in Phi [cm]
+Double_t fOverlapR = 2.0; // Overlap in R direction [cm]
// Station Zceneter [cm] in the cave reference frame
-Double_t fStationZ0[4]={75,125,175,235};
-Int_t fNlayers[4]={3,3,3,3}; // Number of layers
-Int_t fDetType[4]={3,3,4,4}; // Detector type
-Double_t fLayersDz[4]={10,10,10,10};
-
-
-
-
+Double_t fStationZ0[4] = {75, 125, 175, 235};
+Int_t fNlayers[4] = {3, 3, 3, 3}; // Number of layers
+Int_t fDetType[4] = {3, 3, 4, 4}; // Detector type
+Double_t fLayersDz[4] = {10, 10, 10, 10};
// Input parameters for beam pipe shielding
// spans from 2.9 degree to 5.1 degree
//Inner radius is tan(2.9) + 2 cm, extra 20 mm for clamp connection
-const Int_t fNshs=5;
-TString ShMaterial[5]={"Al","Pb","Al","Al","Al"};
-Double_t fShieldZin[5]={125,153.0, 215.0, 265.0, 315.0};
-Double_t fShieldLz[5]={28, 30, 20, 20, 30};
-Double_t fShield_AcceptanceTanMin[5] = {0.043,0.043, 0.051,0.051,0.051}; // Acceptance tangent min for shield
-Double_t fShield_AcceptanceTanMax[5] = {0.1,0.1,0.1,0.1,0.1}; // Acceptance tangent max for shield
+const Int_t fNshs = 5;
+TString ShMaterial[5] = {"Al", "Pb", "Al", "Al", "Al"};
+Double_t fShieldZin[5] = {125, 153.0, 215.0, 265.0, 315.0};
+Double_t fShieldLz[5] = {28, 30, 20, 20, 30};
+Double_t fShield_AcceptanceTanMin[5] = {0.043, 0.043, 0.051, 0.051, 0.051}; // Acceptance tangent min for shield
+Double_t fShield_AcceptanceTanMax[5] = {0.1, 0.1, 0.1, 0.1, 0.1}; // Acceptance tangent max for shield
//***********************************************************
// some global variables
-TGeoManager* gGeoMan = NULL; // Pointer to TGeoManager instance
-TGeoVolume* gModules[fNabs]; // Global storage for module types
-TGeoVolume* gModules_shield[fNshs]; // Global storage for module types
-TGeoVolume* gModules_station[fNst]; // Global storage for module types
+TGeoManager* gGeoMan = NULL; // Pointer to TGeoManager instance
+TGeoVolume* gModules[fNabs]; // Global storage for module types
+TGeoVolume* gModules_shield[fNshs]; // Global storage for module types
+TGeoVolume* gModules_station[fNst]; // Global storage for module types
// Forward declarations
void create_materials_from_media_file();
@@ -149,119 +140,123 @@ TGeoVolume* CreateRpcLayers(int istn, int ily);
fstream infoFile;
-void create_MUCH_geometry_v21a_lmvm() {
-
-
+void create_MUCH_geometry_v21a_lmvm()
+{
+
+
// ------- Open info file -----------------------------------------------
TString infoFileName = FileNameSim;
infoFileName.ReplaceAll("root", "info");
infoFile.open(infoFileName.Data(), fstream::out);
infoFile << "MUCH geometry created with create_MUCH_geometry_v21a_lmvm.C" << endl << endl;
- infoFile<<"Global Variables: "<<endl;
- infoFile<<"MuchCave Zin position = "<<fMuchZ1<<" cm "<<endl;
- infoFile<<"Acceptance tangent min = "<<fAcceptanceTanMin<<endl;
- infoFile<<"Acceptance tangent max = "<<fAcceptanceTanMax<<endl;
+ infoFile << "Global Variables: " << endl;
+ infoFile << "MuchCave Zin position = " << fMuchZ1 << " cm " << endl;
+ infoFile << "Acceptance tangent min = " << fAcceptanceTanMin << endl;
+ infoFile << "Acceptance tangent max = " << fAcceptanceTanMax << endl;
// Load needed material definition from media.geo file
create_materials_from_media_file();
-
+
// Get the GeoManager for later usage
gGeoMan = (TGeoManager*) gROOT->FindObject("FAIRGeom");
gGeoMan->SetVisLevel(10);
-
+
// Create the top volume
- TGeoBBox* topbox= new TGeoBBox("", 1000., 1000., 2000.);
- TGeoVolume* top = new TGeoVolume("top", topbox, gGeoMan->GetMedium("air"));
+ TGeoBBox* topbox = new TGeoBBox("", 1000., 1000., 2000.);
+ TGeoVolume* top = new TGeoVolume("top", topbox, gGeoMan->GetMedium("air"));
gGeoMan->SetTopVolume(top);
- TString topName=geoVersion+tagVersion+subVersion;
+ TString topName = geoVersion + tagVersion + subVersion;
TGeoVolume* much = new TGeoVolumeAssembly(topName);
top->AddNode(much, 1);
- TGeoVolume *absr = new TGeoVolumeAssembly("absorber");
- much->AddNode(absr,1);
-
- TGeoVolume *shld = new TGeoVolumeAssembly("shield");
- much->AddNode(shld,1);
-
- TGeoVolume *sttn = new TGeoVolumeAssembly("station");
- much->AddNode(sttn,1);
-
- infoFile<<endl;
- infoFile<<" Absorbers "<<endl;
- infoFile<<" -----------"<<endl;
- infoFile<<"Total No. of Absorbers: "<<fAbs<<endl;
- infoFile<<"First abosrber is divided into two halves."<<endl;
- infoFile<<"First half inserted inside the Dipole Magnet."<<endl;
- infoFile<<"Second half is made of Low Density Graphite + Concrete."<<endl;
- infoFile<<"Total No. of Pieces: "<<fNabs<<endl;
- infoFile<<endl;
- infoFile<<"AbsPieces Position[cm] Thickness[cm] Material"<<endl;
- infoFile<<"--------------------------------------------------------------"<<endl;
-
- for (Int_t iabs = 0; iabs <fNabs ; iabs++) {
-
+ TGeoVolume* absr = new TGeoVolumeAssembly("absorber");
+ much->AddNode(absr, 1);
+
+ TGeoVolume* shld = new TGeoVolumeAssembly("shield");
+ much->AddNode(shld, 1);
+
+ TGeoVolume* sttn = new TGeoVolumeAssembly("station");
+ much->AddNode(sttn, 1);
+
+ infoFile << endl;
+ infoFile << " Absorbers " << endl;
+ infoFile << " -----------" << endl;
+ infoFile << "Total No. of Absorbers: " << fAbs << endl;
+ infoFile << "First abosrber is divided into two halves." << endl;
+ infoFile << "First half inserted inside the Dipole Magnet." << endl;
+ infoFile << "Second half is made of Low Density Graphite + Concrete." << endl;
+ infoFile << "Total No. of Pieces: " << fNabs << endl;
+ infoFile << endl;
+ infoFile << "AbsPieces Position[cm] Thickness[cm] Material" << endl;
+ infoFile << "--------------------------------------------------------------" << endl;
+
+ for (Int_t iabs = 0; iabs < fNabs; iabs++) {
+
gModules[iabs] = CreateAbsorbers(iabs);
-
- absr->AddNode(gModules[iabs],iabs);
+
+ absr->AddNode(gModules[iabs], iabs);
}
-
-
- infoFile<<endl;
- infoFile<<" Shielding "<<endl;
- infoFile<<" -----------"<<endl;
- infoFile<<"No. of Shields: "<<fNshs<<endl;
- infoFile<<"Inside the Abs I, Shielding divided into two parts."<<endl;
- infoFile<<endl;
- infoFile<<"Shield No. Z_In[cm] Z_Out[cm] R_In[cm] R_Out[cm] Material"<<endl;
- infoFile<<"--------------------------------------------------------------"<<endl;
- for (Int_t ishi = 0; ishi <fNshs ; ishi++) {
-
+
+
+ infoFile << endl;
+ infoFile << " Shielding " << endl;
+ infoFile << " -----------" << endl;
+ infoFile << "No. of Shields: " << fNshs << endl;
+ infoFile << "Inside the Abs I, Shielding divided into two parts." << endl;
+ infoFile << endl;
+ infoFile << "Shield No. Z_In[cm] Z_Out[cm] R_In[cm] R_Out[cm] Material" << endl;
+ infoFile << "--------------------------------------------------------------" << endl;
+ for (Int_t ishi = 0; ishi < fNshs; ishi++) {
+
gModules_shield[ishi] = CreateShields(ishi);
-
- shld->AddNode(gModules_shield[ishi],ishi);
-
+
+ shld->AddNode(gModules_shield[ishi], ishi);
}
-
- infoFile<<endl;
- infoFile<<" Stations "<<endl;
- infoFile<<" ----------"<<endl;
- infoFile<<"No. of Stations: "<<fNst<<endl;
- infoFile<<"First two stations (1,2) are made up of GEM and last two stations (3,4) are made up of RPC."<<endl;
- infoFile<<"Passive material implemented in GEM modules. Ar:CO2 (70:30) is used as active gas. "<<endl;
- infoFile<<"12 mm thick Al plates are used for support and cooling in the GEM modules."<<endl;
- infoFile<<"2 mm thick Aluminum plates are used for support in the RPC modules behind the active area. 10 mm thick Aluminium at the boundaries as the frame."<<endl;
- infoFile<<"Drift and read-out PCBs (copper coated G10 plates) inserted for realistic material budget for both GEM and RPC modules."<<endl;
- infoFile<<"#Station #Layers Z[cm] #Sectors ActiveLz[cm]"<<endl;
- infoFile<<"--------------------------------------------------------------"<<endl;
- for (Int_t istn = 0; istn < 4; istn++) { // 4 Stations
-
-
+
+ infoFile << endl;
+ infoFile << " Stations " << endl;
+ infoFile << " ----------" << endl;
+ infoFile << "No. of Stations: " << fNst << endl;
+ infoFile << "First two stations (1,2) are made up of GEM and last two stations (3,4) are made up of RPC." << endl;
+ infoFile << "Passive material implemented in GEM modules. Ar:CO2 (70:30) is used as active gas. " << endl;
+ infoFile << "12 mm thick Al plates are used for support and cooling in the GEM modules." << endl;
+ infoFile << "2 mm thick Aluminum plates are used for support in the RPC modules behind the active area. 10 mm thick "
+ "Aluminium at the boundaries as the frame."
+ << endl;
+ infoFile << "Drift and read-out PCBs (copper coated G10 plates) inserted for realistic material budget for both GEM "
+ "and RPC modules."
+ << endl;
+ infoFile << "#Station #Layers Z[cm] #Sectors ActiveLz[cm]" << endl;
+ infoFile << "--------------------------------------------------------------" << endl;
+ for (Int_t istn = 0; istn < 4; istn++) { // 4 Stations
+
+
gModules_station[istn] = CreateStations(istn);
-
- sttn->AddNode(gModules_station[istn],istn);
+
+ sttn->AddNode(gModules_station[istn], istn);
}
-
+
gGeoMan->CloseGeometry();
gGeoMan->CheckOverlaps(0.0000001);
gGeoMan->PrintOverlaps();
gGeoMan->CheckOverlaps(0.0001, "s");
gGeoMan->PrintOverlaps();
-
-
- much->Export(FileNameSim); // an alternative way of writing the much
-
- TFile* outfile = new TFile(FileNameSim, "UPDATE");
+
+
+ much->Export(FileNameSim); // an alternative way of writing the much
+
+ TFile* outfile = new TFile(FileNameSim, "UPDATE");
TGeoTranslation* much_placement = new TGeoTranslation("much_trans", 0., 0., 0.);
much_placement->Write();
outfile->Close();
-
- outfile = new TFile(FileNameGeo,"RECREATE");
+
+ outfile = new TFile(FileNameGeo, "RECREATE");
gGeoMan->Write(); // use this if you want GeoManager format in the output
outfile->Close();
-
+
top->Draw("ogl");
infoFile.close();
}
@@ -269,15 +264,15 @@ void create_MUCH_geometry_v21a_lmvm() {
void create_materials_from_media_file()
{
// Use the FairRoot geometry interface to load the media which are already defined
- FairGeoLoader* geoLoad = new FairGeoLoader("TGeo", "FairGeoLoader");
+ FairGeoLoader* geoLoad = new FairGeoLoader("TGeo", "FairGeoLoader");
FairGeoInterface* geoFace = geoLoad->getGeoInterface();
- TString geoPath = gSystem->Getenv("VMCWORKDIR");
- TString geoFile = geoPath + "/geometry/media.geo";
+ TString geoPath = gSystem->Getenv("VMCWORKDIR");
+ TString geoFile = geoPath + "/geometry/media.geo";
geoFace->setMediaFile(geoFile);
geoFace->readMedia();
// Read the required media and create them in the GeoManager
- FairGeoMedia* geoMedia = geoFace->getMedia();
+ FairGeoMedia* geoMedia = geoFace->getMedia();
FairGeoBuilder* geoBuild = geoLoad->getGeoBuilder();
FairGeoMedium* air = geoMedia->getMedium(KeepingVolumeMedium);
@@ -304,10 +299,10 @@ void create_materials_from_media_file()
FairGeoMedium* MUCHconcrete = geoMedia->getMedium(Concrete);
geoBuild->createMedium(MUCHconcrete);
- FairGeoMedium* copperplate = geoMedia->getMedium(copper);
+ FairGeoMedium* copperplate = geoMedia->getMedium(copper);
geoBuild->createMedium(copperplate);
- FairGeoMedium* g10plate = geoMedia->getMedium(g10); //G10
+ FairGeoMedium* g10plate = geoMedia->getMedium(g10); //G10
geoBuild->createMedium(g10plate);
FairGeoMedium* RPCmedium = geoMedia->getMedium(RPCm);
@@ -321,701 +316,702 @@ void create_materials_from_media_file()
}
+TGeoVolume* CreateShields(int ish)
+{
+
-TGeoVolume* CreateShields(int ish) {
+ TGeoMedium* iron = gGeoMan->GetMedium(I);
+ TGeoMedium* lead = gGeoMan->GetMedium(L);
+ TGeoMedium* Aluminium = gGeoMan->GetMedium(Al);
-
- TGeoMedium* iron = gGeoMan->GetMedium(I);
- TGeoMedium* lead = gGeoMan->GetMedium(L);
- TGeoMedium* Aluminium= gGeoMan->GetMedium(Al);
-
- TString name = Form("shieldblock%d", ish);
+ TString name = Form("shieldblock%d", ish);
TGeoVolumeAssembly* shieldblock = new TGeoVolumeAssembly(name);
-
- TString conename_sh = Form("conesh_%d",ish);
-
-
- Double_t dz = fShieldLz[ish]/2.0 ;
- Double_t globalZ1 = fShieldZin[ish] ;
- Double_t globalZ2 = fShieldZin[ish] + 2 * dz ;
-
-
-
- Double_t rmin1 = globalZ1 * fShield_AcceptanceTanMin[ish]+2.0;
- Double_t rmax1 = globalZ1 * fShield_AcceptanceTanMax[ish]-0.0001;
- Double_t rmin2 = globalZ2 * fShield_AcceptanceTanMin[ish]+2.0;
- Double_t rmax2 = globalZ2 * fShield_AcceptanceTanMax[ish]-0.0001;
-
- infoFile<<" "<<ish<<"\t\t"<<globalZ1<<"\t "<<globalZ2<<"\t"<<rmin1<<"\t "<<rmax1<<"\t "<<ShMaterial[ish]<<endl;
-
- if(ish==0){
-
- TGeoCone * sh =new TGeoCone(conename_sh,dz, rmin1, rmax1, rmin2, rmax2);
+
+ TString conename_sh = Form("conesh_%d", ish);
+
+
+ Double_t dz = fShieldLz[ish] / 2.0;
+ Double_t globalZ1 = fShieldZin[ish];
+ Double_t globalZ2 = fShieldZin[ish] + 2 * dz;
+
+
+ Double_t rmin1 = globalZ1 * fShield_AcceptanceTanMin[ish] + 2.0;
+ Double_t rmax1 = globalZ1 * fShield_AcceptanceTanMax[ish] - 0.0001;
+ Double_t rmin2 = globalZ2 * fShield_AcceptanceTanMin[ish] + 2.0;
+ Double_t rmax2 = globalZ2 * fShield_AcceptanceTanMax[ish] - 0.0001;
+
+ infoFile << " " << ish << "\t\t" << globalZ1 << "\t " << globalZ2 << "\t" << rmin1 << "\t " << rmax1 << "\t "
+ << ShMaterial[ish] << endl;
+
+ if (ish == 0) {
+
+ TGeoCone* sh = new TGeoCone(conename_sh, dz, rmin1, rmax1, rmin2, rmax2);
TGeoVolume* shield = new TGeoVolume("shield", sh, Aluminium);
shield->SetLineColor(7);
shield->SetTransparency(2);
- TGeoTranslation *sh_trans = new TGeoTranslation("", 0., 0., globalZ1+dz);
- shieldblock->AddNode(shield,ish, sh_trans);
-
-
+ TGeoTranslation* sh_trans = new TGeoTranslation("", 0., 0., globalZ1 + dz);
+ shieldblock->AddNode(shield, ish, sh_trans);
}
-
- if(ish==1){
- TGeoCone * sh =new TGeoCone(conename_sh,dz, rmin1, rmax1, rmin2, rmax2);
+
+ if (ish == 1) {
+ TGeoCone* sh = new TGeoCone(conename_sh, dz, rmin1, rmax1, rmin2, rmax2);
TGeoVolume* shield = new TGeoVolume("shield", sh, lead);
-
+
shield->SetLineColor(kMagenta);
shield->SetTransparency(2);
- TGeoTranslation *sh_trans = new TGeoTranslation("", 0., 0., globalZ1+dz);
- shieldblock->AddNode(shield,ish, sh_trans);
+ TGeoTranslation* sh_trans = new TGeoTranslation("", 0., 0., globalZ1 + dz);
+ shieldblock->AddNode(shield, ish, sh_trans);
+ }
+
+ if (ish > 1) {
+
+
+ TGeoCone* sh = new TGeoCone(conename_sh, dz, rmin1, rmax1, rmin2, rmax2);
+ TGeoVolume* shield = new TGeoVolume("shield", sh, Aluminium);
+
+ shield->SetLineColor(kBlack);
+ shield->SetTransparency(2);
+ TGeoTranslation* sh_trans = new TGeoTranslation("", 0., 0., globalZ1 + dz);
+ shieldblock->AddNode(shield, ish, sh_trans);
}
- if(ish>1)
- {
-
-
- TGeoCone * sh =new TGeoCone(conename_sh,dz, rmin1, rmax1, rmin2, rmax2);
- TGeoVolume* shield = new TGeoVolume("shield", sh, Aluminium);
-
- shield->SetLineColor(kBlack);
- shield->SetTransparency(2);
- TGeoTranslation *sh_trans = new TGeoTranslation("", 0., 0., globalZ1+dz);
- shieldblock->AddNode(shield,ish, sh_trans);
- }
-
return shieldblock;
-
}
-TGeoVolume* CreateAbsorbers(int i) {
-
- TGeoMedium* graphite = gGeoMan->GetMedium(LDcarbon);
- TGeoMedium* iron = gGeoMan->GetMedium(I);
- TGeoMedium* concrete = gGeoMan->GetMedium(Concrete);
- TGeoMedium* Aluminium= gGeoMan->GetMedium(Al);
+TGeoVolume* CreateAbsorbers(int i)
+{
- TString name = Form("absblock%d", i);
+ TGeoMedium* graphite = gGeoMan->GetMedium(LDcarbon);
+ TGeoMedium* iron = gGeoMan->GetMedium(I);
+ TGeoMedium* concrete = gGeoMan->GetMedium(Concrete);
+ TGeoMedium* Aluminium = gGeoMan->GetMedium(Al);
+
+ TString name = Form("absblock%d", i);
TGeoVolumeAssembly* absblock = new TGeoVolumeAssembly(name);
- TString pipename = Form("beampipe_%d",i);
- TString conename = Form("cone_%d",i);
- TString BoxName = Form("Box_%d",i);
- TString supportShapeName = Form("Support_%d",i);
- TString TrapName = Form("Trap_%d",i);
-
- Double_t dz = fAbsorberLz[i]/2.0 ;
- Double_t globalZ1 = fAbsorberZ1[i] + fMuchZ1;
+ TString pipename = Form("beampipe_%d", i);
+ TString conename = Form("cone_%d", i);
+ TString BoxName = Form("Box_%d", i);
+ TString supportShapeName = Form("Support_%d", i);
+ TString TrapName = Form("Trap_%d", i);
+
+ Double_t dz = fAbsorberLz[i] / 2.0;
+ Double_t globalZ1 = fAbsorberZ1[i] + fMuchZ1;
Double_t globalPos = globalZ1 + dz;
- Double_t globalZ2 = fAbsorberZ1[i] + 2 * dz + fMuchZ1;
+ Double_t globalZ2 = fAbsorberZ1[i] + 2 * dz + fMuchZ1;
Double_t rmin1 = globalZ1 * fAcceptanceTanMin;
Double_t rmin2 = globalZ2 * fAcceptanceTanMin;
Double_t rmax1 = globalZ1 * fAcceptanceTanMax + safetyrad[i];
- Double_t rmax2 = globalZ2 * fAcceptanceTanMax + safetyrad[i];//
+ Double_t rmax2 = globalZ2 * fAcceptanceTanMax + safetyrad[i]; //
- infoFile<<" "<<i+1<<"\t\t"<<globalPos<<"\t\t"<<2*dz<<"\t\t"<<AbsMaterial[i]<<endl;
+ infoFile << " " << i + 1 << "\t\t" << globalPos << "\t\t" << 2 * dz << "\t\t" << AbsMaterial[i] << endl;
// 1st part of 1st absorber LD Carbon
//dimensions are hardcoded
- if(i==0)
- {
- // printf("absorber %d \n",i);
-
- TGeoTrd2 * trap = new TGeoTrd2(TrapName,70.0,70.0,46.0,71.0,dz);
- TGeoCone * tube = new TGeoCone(pipename,dz+0.001,0.,rmin1,0.,rmin2);
- TString expression = TrapName +"-"+pipename;
- TGeoCompositeShape* shSupport = new TGeoCompositeShape(supportShapeName,expression);
- TGeoVolume* abs0 = new TGeoVolume("absorber", shSupport, graphite);
- abs0->SetLineColor(kRed);
- abs0->SetTransparency(0);
- TGeoTranslation *abs0_trans = new TGeoTranslation("", 0., 0., globalZ1+dz);
- absblock->AddNode(abs0, i, abs0_trans);
- }
-
+ if (i == 0) {
+ // printf("absorber %d \n",i);
+
+ TGeoTrd2* trap = new TGeoTrd2(TrapName, 70.0, 70.0, 46.0, 71.0, dz);
+ TGeoCone* tube = new TGeoCone(pipename, dz + 0.001, 0., rmin1, 0., rmin2);
+ TString expression = TrapName + "-" + pipename;
+ TGeoCompositeShape* shSupport = new TGeoCompositeShape(supportShapeName, expression);
+ TGeoVolume* abs0 = new TGeoVolume("absorber", shSupport, graphite);
+ abs0->SetLineColor(kRed);
+ abs0->SetTransparency(0);
+ TGeoTranslation* abs0_trans = new TGeoTranslation("", 0., 0., globalZ1 + dz);
+ absblock->AddNode(abs0, i, abs0_trans);
+ }
// 2rd part of 1st absorber box (LD Carbon)
- if(i==1)
- {
- // printf("absorber %d \n",i);
- TGeoBBox * box = new TGeoBBox(BoxName,130.0,125.0,dz);
- TGeoCone * tube = new TGeoCone(pipename,dz+0.001,0.,rmin1,0.,rmin2);
- TString expression = BoxName +"-"+pipename;
- TGeoCompositeShape* shSupport = new TGeoCompositeShape(supportShapeName,expression);
-
- TGeoVolume* abs1 = new TGeoVolume("absorber", shSupport, graphite);
- abs1->SetLineColor(kRed);
- abs1->SetTransparency(0);
- TGeoTranslation *abs1_trans = new TGeoTranslation("", 0., 0., globalZ1+dz);
- absblock->AddNode(abs1, i, abs1_trans);
+ if (i == 1) {
+ // printf("absorber %d \n",i);
+ TGeoBBox* box = new TGeoBBox(BoxName, 130.0, 125.0, dz);
+ TGeoCone* tube = new TGeoCone(pipename, dz + 0.001, 0., rmin1, 0., rmin2);
+ TString expression = BoxName + "-" + pipename;
+ TGeoCompositeShape* shSupport = new TGeoCompositeShape(supportShapeName, expression);
+
+ TGeoVolume* abs1 = new TGeoVolume("absorber", shSupport, graphite);
+ abs1->SetLineColor(kRed);
+ abs1->SetTransparency(0);
+ TGeoTranslation* abs1_trans = new TGeoTranslation("", 0., 0., globalZ1 + dz);
+ absblock->AddNode(abs1, i, abs1_trans);
+ }
+ // 3th part of 1st absorber box (Conc)
+ if (i == 2) {
+ // printf("absorber %d \n",i);
+ TGeoBBox* box = new TGeoBBox(BoxName, 130.0, 125.0, dz);
+ TGeoCone* tube = new TGeoCone(pipename, dz + 0.001, 0., rmin1, 0., rmin2);
+ TString expression = BoxName + "-" + pipename;
+ TGeoCompositeShape* shSupport = new TGeoCompositeShape(supportShapeName, expression);
+
+ TGeoVolume* abs1 = new TGeoVolume("absorber", shSupport, concrete);
+ abs1->SetLineColor(kViolet);
+ abs1->SetTransparency(0);
+ TGeoTranslation* abs1_trans = new TGeoTranslation("", 0., 0., globalZ1 + dz);
+ absblock->AddNode(abs1, i, abs1_trans);
+ }
- }
- // 3th part of 1st absorber box (Conc)
- if(i==2)
- {
- // printf("absorber %d \n",i);
- TGeoBBox * box = new TGeoBBox(BoxName,130.0,125.0,dz);
- TGeoCone * tube = new TGeoCone(pipename,dz+0.001,0.,rmin1,0.,rmin2);
- TString expression = BoxName +"-"+pipename;
- TGeoCompositeShape* shSupport = new TGeoCompositeShape(supportShapeName,expression);
-
- TGeoVolume* abs1 = new TGeoVolume("absorber", shSupport, concrete);
- abs1->SetLineColor(kViolet);
- abs1->SetTransparency(0);
- TGeoTranslation *abs1_trans = new TGeoTranslation("", 0., 0., globalZ1+dz);
- absblock->AddNode(abs1, i, abs1_trans);
- }
+ //rest of the absorbers
+ if (i > 2) {
+ TGeoBBox* box = new TGeoBBox(BoxName, rmax2, rmax2, dz);
+ TGeoCone* tube = new TGeoCone(pipename, dz + 0.001, 0., rmin1, 0., rmin2);
+ TString expression = BoxName + "-" + pipename;
+ TGeoCompositeShape* shSupport = new TGeoCompositeShape(supportShapeName, expression);
+ TGeoVolume* abs2 = new TGeoVolume("absorber", shSupport, iron);
+ abs2->SetLineColor(kBlue);
+ abs2->SetTransparency(2);
- //rest of the absorbers
- if (i>2)
- {
- TGeoBBox * box = new TGeoBBox(BoxName,rmax2,rmax2,dz);
- TGeoCone * tube = new TGeoCone(pipename,dz+0.001,0.,rmin1,0.,rmin2);
- TString expression = BoxName +"-"+pipename;
- TGeoCompositeShape* shSupport = new TGeoCompositeShape(supportShapeName,expression);
- TGeoVolume* abs2 = new TGeoVolume("absorber", shSupport, iron);
+ TGeoTranslation* abs_trans = new TGeoTranslation("", 0., 0., globalZ1 + dz);
+ absblock->AddNode(abs2, i, abs_trans);
+ }
- abs2->SetLineColor(kBlue);
- abs2->SetTransparency(2);
+ return absblock;
+}
+TGeoVolume* CreateStations(int ist)
+{
- TGeoTranslation *abs_trans = new TGeoTranslation("", 0., 0., globalZ1+dz);
- absblock->AddNode(abs2,i, abs_trans);
+ TString stationName = Form("muchstation%02i", ist + 1);
+ TGeoVolumeAssembly* station = new TGeoVolumeAssembly(stationName); //, shStation, air);
- }
+ TGeoVolume* gLayer[4];
- return absblock;
+ for (int ii = 0; ii < 3; ii++) { // 3 Layers
-}
+ switch (ist) {
+ case 0:
+ case 1: gLayer[ii] = CreateGemLayers(ist, ii); break;
+ case 2:
+ case 3: gLayer[ii] = CreateRpcLayers(ist, ii); break;
+ }
-TGeoVolume * CreateStations(int ist){
-
- TString stationName = Form("muchstation%02i",ist+1);
-
- TGeoVolumeAssembly* station = new TGeoVolumeAssembly(stationName);//, shStation, air);
-
- TGeoVolume* gLayer[4];
-
- for (int ii=0;ii<3;ii++){ // 3 Layers
-
- switch(ist){
- case 0:
- case 1: gLayer[ii] = CreateGemLayers(ist, ii); break;
- case 2:
- case 3: gLayer[ii] = CreateRpcLayers(ist, ii); break;}
-
- station->AddNode(gLayer[ii],ii);
-
+ station->AddNode(gLayer[ii], ii);
}
-
+
return station;
}
-TGeoVolume * CreateGemLayers(int istn, int ily){
-
- TString layerName = Form("muchstation%02ilayer%i",istn+1,ily+1);
+TGeoVolume* CreateGemLayers(int istn, int ily)
+{
+
+ TString layerName = Form("muchstation%02ilayer%i", istn + 1, ily + 1);
TGeoVolumeAssembly* volayer = new TGeoVolumeAssembly(layerName);
//GEM Parametrs
- Double_t fSupportDz = 1.2; //1.2 cm Al (cooling + support)
- Double_t fCopperDz = 0.0065; // 65 micron copper
- Double_t fCopperSliceDz = 0.0005; //5 micron copper slices
- Double_t fReadoutPlateDz = 0.3; //3mm G10
- Double_t fDriftPlateDz = 0.3; //3 mm G10
- Double_t fPassiveVolumeDz = 0.2; // 2 mm of Argon
- Double_t fActiveVolumeDz = 0.3; //3mm argon
- Double_t fKaptonDz = 0.005; //50 micron Kapton
- Double_t fFrameDz = 0.2; // 2mm frame
-
-
- Double_t stGlobalZ0 = fStationZ0[istn] + fMuchZ1; //z position of station center (midplane) [cm]
-
-
-
- Double_t layerZ0 = (ily - (fNlayers[istn] - 1) / 2.) * fLayersDz[istn];
+ Double_t fSupportDz = 1.2; //1.2 cm Al (cooling + support)
+ Double_t fCopperDz = 0.0065; // 65 micron copper
+ Double_t fCopperSliceDz = 0.0005; //5 micron copper slices
+ Double_t fReadoutPlateDz = 0.3; //3mm G10
+ Double_t fDriftPlateDz = 0.3; //3 mm G10
+ Double_t fPassiveVolumeDz = 0.2; // 2 mm of Argon
+ Double_t fActiveVolumeDz = 0.3; //3mm argon
+ Double_t fKaptonDz = 0.005; //50 micron Kapton
+ Double_t fFrameDz = 0.2; // 2mm frame
+
+
+ Double_t stGlobalZ0 = fStationZ0[istn] + fMuchZ1; //z position of station center (midplane) [cm]
+
+
+ Double_t layerZ0 = (ily - (fNlayers[istn] - 1) / 2.) * fLayersDz[istn];
Double_t layerGlobalZ0 = layerZ0 + stGlobalZ0;
-
+
//Active Gas distance from layer position
- Double_t ModuleZ = fSupportDz/2. + 2*fCopperDz + fReadoutPlateDz + 3*fPassiveVolumeDz + 6*fCopperSliceDz + 3*fKaptonDz + fActiveVolumeDz/2.;
-
+ Double_t ModuleZ = fSupportDz / 2. + 2 * fCopperDz + fReadoutPlateDz + 3 * fPassiveVolumeDz + 6 * fCopperSliceDz
+ + 3 * fKaptonDz + fActiveVolumeDz / 2.;
+
//ReadOut Plate distance from layer position
- Double_t ReadOutZ = fSupportDz/2. + fCopperDz + fReadoutPlateDz/2.;
-
+ Double_t ReadOutZ = fSupportDz / 2. + fCopperDz + fReadoutPlateDz / 2.;
+
//Copper distance from layer position
- Double_t CooperIZ = fSupportDz/2. + fCopperDz/2.;
- Double_t CooperIIZ = ReadOutZ + fReadoutPlateDz/2. + fCopperDz/2.;
- Double_t CooperIIIZ = ModuleZ + fActiveVolumeDz/2. + fCopperDz/2.;
+ Double_t CooperIZ = fSupportDz / 2. + fCopperDz / 2.;
+ Double_t CooperIIZ = ReadOutZ + fReadoutPlateDz / 2. + fCopperDz / 2.;
+ Double_t CooperIIIZ = ModuleZ + fActiveVolumeDz / 2. + fCopperDz / 2.;
Double_t CooperIVZ = CooperIIIZ + fCopperDz + fDriftPlateDz;
-
+
//Drift Plate distance from layer position
- Double_t DriftZ = ModuleZ + fActiveVolumeDz/2. + fCopperDz + fDriftPlateDz/2.;
-
+ Double_t DriftZ = ModuleZ + fActiveVolumeDz / 2. + fCopperDz + fDriftPlateDz / 2.;
+
//Passive Gas distance from layer position
- Double_t PassivePassiveDz = fPassiveVolumeDz + 2*fCopperSliceDz + fKaptonDz;
- Double_t PassiveGasIZ = ReadOutZ + fReadoutPlateDz/2. + fCopperDz + fPassiveVolumeDz/2.;
- Double_t PassiveGasIIZ = PassiveGasIZ + PassivePassiveDz;
- Double_t PassiveGasIIIZ = PassiveGasIIZ + PassivePassiveDz;
+ Double_t PassivePassiveDz = fPassiveVolumeDz + 2 * fCopperSliceDz + fKaptonDz;
+ Double_t PassiveGasIZ = ReadOutZ + fReadoutPlateDz / 2. + fCopperDz + fPassiveVolumeDz / 2.;
+ Double_t PassiveGasIIZ = PassiveGasIZ + PassivePassiveDz;
+ Double_t PassiveGasIIIZ = PassiveGasIIZ + PassivePassiveDz;
//Kapton distance from layer position
- Double_t KaptonPassiveDz = fPassiveVolumeDz/2. + fCopperSliceDz + fKaptonDz/2.;
- Double_t KaptonIZ = PassiveGasIZ + KaptonPassiveDz;
- Double_t KaptonIIZ = PassiveGasIIZ + KaptonPassiveDz;
- Double_t KaptonIIIZ = PassiveGasIIIZ + KaptonPassiveDz;
-
+ Double_t KaptonPassiveDz = fPassiveVolumeDz / 2. + fCopperSliceDz + fKaptonDz / 2.;
+ Double_t KaptonIZ = PassiveGasIZ + KaptonPassiveDz;
+ Double_t KaptonIIZ = PassiveGasIIZ + KaptonPassiveDz;
+ Double_t KaptonIIIZ = PassiveGasIIIZ + KaptonPassiveDz;
+
//Copper Slice distance from layer position
- Double_t PassiveCopperDz = fPassiveVolumeDz/2. + fCopperSliceDz/2.;
- Double_t KaptonCopperDz = fKaptonDz + fCopperSliceDz;
- Double_t CopperSilceIZ = PassiveGasIZ + PassiveCopperDz;
- Double_t CopperSilceIIZ = CopperSilceIZ + KaptonCopperDz;
+ Double_t PassiveCopperDz = fPassiveVolumeDz / 2. + fCopperSliceDz / 2.;
+ Double_t KaptonCopperDz = fKaptonDz + fCopperSliceDz;
+ Double_t CopperSilceIZ = PassiveGasIZ + PassiveCopperDz;
+ Double_t CopperSilceIIZ = CopperSilceIZ + KaptonCopperDz;
Double_t CopperSilceIIIZ = PassiveGasIIZ + PassiveCopperDz;
- Double_t CopperSilceIVZ = CopperSilceIIIZ + KaptonCopperDz;
- Double_t CopperSilceVZ = PassiveGasIIIZ + PassiveCopperDz;
- Double_t CopperSilceVIZ = CopperSilceVZ + KaptonCopperDz;
+ Double_t CopperSilceIVZ = CopperSilceIIIZ + KaptonCopperDz;
+ Double_t CopperSilceVZ = PassiveGasIIIZ + PassiveCopperDz;
+ Double_t CopperSilceVIZ = CopperSilceVZ + KaptonCopperDz;
+
+ //Set Rmin & Rmax
+ Double_t stDz = ((fNlayers[istn] - 1) * fLayersDz[istn]) / 2. + CooperIVZ + fCopperDz / 2.;
- //Set Rmin & Rmax
- Double_t stDz = ((fNlayers[istn] - 1) * fLayersDz[istn])/2. + CooperIVZ + fCopperDz/2.;
-
Double_t stGlobalZ2 = stGlobalZ0 + stDz;
Double_t stGlobalZ1 = stGlobalZ0 - stDz;
-
+
Double_t rmin = stGlobalZ1 * fAcceptanceTanMin;
Double_t rmax = stGlobalZ2 * fAcceptanceTanMax;
// Module dimention calculation
- Double_t phi0 = TMath::Pi()/fNSectorsPerLayer[istn]; // azimuthal half widh of each module
- Double_t ymin = rmin+fSpacerR;
+ Double_t phi0 = TMath::Pi() / fNSectorsPerLayer[istn]; // azimuthal half widh of each module
+ Double_t ymin = rmin + fSpacerR;
Double_t ymax = rmax;
//define the dimensions of the trapezoidal module
- Double_t dy = (ymax-ymin)/2.; //y (length)
- Double_t dx1 = ymin*TMath::Tan(phi0)+fOverlapR/TMath::Cos(phi0); // large x
- Double_t dx2 = ymax*TMath::Tan(phi0)+fOverlapR/TMath::Cos(phi0); // small x
+ Double_t dy = (ymax - ymin) / 2.; //y (length)
+ Double_t dx1 = ymin * TMath::Tan(phi0) + fOverlapR / TMath::Cos(phi0); // large x
+ Double_t dx2 = ymax * TMath::Tan(phi0) + fOverlapR / TMath::Cos(phi0); // small x
//define the spacer dimensions
- Double_t tg = (dx2-dx1)/2/dy;
- Double_t dd1 = fSpacerPhi*tg;
- Double_t dd2 = fSpacerPhi*sqrt(1+tg*tg);
- Double_t sdx1 = dx1+dd2-dd1-0.1; // 0.1 cm to avoid overlaps
- Double_t sdx2 = dx2+dd2+dd1;
- Double_t sdy = dy+fSpacerR;
-
-
-
- infoFile<<" "<<istn+1<<"\t "<<ily+1<<"\t\t"<<layerGlobalZ0<<"\t"<<fNSectorsPerLayer[istn]<<"\t"<<fActiveVolumeDz<<endl;
-
-
+ Double_t tg = (dx2 - dx1) / 2 / dy;
+ Double_t dd1 = fSpacerPhi * tg;
+ Double_t dd2 = fSpacerPhi * sqrt(1 + tg * tg);
+ Double_t sdx1 = dx1 + dd2 - dd1 - 0.1; // 0.1 cm to avoid overlaps
+ Double_t sdx2 = dx2 + dd2 + dd1;
+ Double_t sdy = dy + fSpacerR;
+
+
+ infoFile << " " << istn + 1 << "\t " << ily + 1 << "\t\t" << layerGlobalZ0 << "\t" << fNSectorsPerLayer[istn]
+ << "\t" << fActiveVolumeDz << endl;
+
+
// Aluminum Plate (Cooling + Support)
- TString supportAlName = Form("shStation%02iSupportAl",istn+1);
- TGeoTube* shSupportAl = new TGeoTube(supportAlName,rmin,rmax,fSupportDz/2.);
-
- TString supportName1 = Form("muchstation%02ilayer%isupportAl",istn+1,ily+1);
- TGeoMedium* coolMat = gGeoMan->GetMedium(Al);
-
- TGeoVolume* voSupport1 = new TGeoVolume(supportName1,shSupportAl,coolMat);
+ TString supportAlName = Form("shStation%02iSupportAl", istn + 1);
+ TGeoTube* shSupportAl = new TGeoTube(supportAlName, rmin, rmax, fSupportDz / 2.);
+
+ TString supportName1 = Form("muchstation%02ilayer%isupportAl", istn + 1, ily + 1);
+ TGeoMedium* coolMat = gGeoMan->GetMedium(Al);
+
+ TGeoVolume* voSupport1 = new TGeoVolume(supportName1, shSupportAl, coolMat);
voSupport1->SetLineColor(kCyan);
-
- TGeoTranslation *support_trans1 = new TGeoTranslation("supportName1", 0,0,layerGlobalZ0);
- volayer->AddNode(voSupport1,0,support_trans1);
-
-
-
+
+ TGeoTranslation* support_trans1 = new TGeoTranslation("supportName1", 0, 0, layerGlobalZ0);
+ volayer->AddNode(voSupport1, 0, support_trans1);
+
+
// Now start adding the GEM modules
- for (Int_t iModule=0; iModule<fNSectorsPerLayer[istn]; iModule++){
-
+ for (Int_t iModule = 0; iModule < fNSectorsPerLayer[istn]; iModule++) {
+
Double_t phi = 2 * phi0 * (iModule + 0.5); // add 0.5 to not overlap with y-axis for left-right layer separation
- Bool_t isBack = iModule%2;
- Char_t cside = (isBack==1) ? 'b' : 'f';
- Int_t iMod = iModule/2;
+ Bool_t isBack = iModule % 2;
+ Char_t cside = (isBack == 1) ? 'b' : 'f';
+ Int_t iMod = iModule / 2;
// correct the x, y positions
Double_t pos[21];
- pos[0] = -(ymin+dy)*sin(phi);
- pos[1] = (ymin+dy)*cos(phi);
+ pos[0] = -(ymin + dy) * sin(phi);
+ pos[1] = (ymin + dy) * cos(phi);
// different z positions for odd/even modules
- pos[2] = (isBack ? 1 : -1)*ModuleZ + layerGlobalZ0; //Active volume & Frame
- pos[3] = (isBack ? 1 : -1)*ReadOutZ + layerGlobalZ0;//Readout G10
- pos[4] = (isBack ? 1 : -1)*CooperIZ + layerGlobalZ0; //Copper I
- pos[5] = (isBack ? 1 : -1)*CooperIIZ + layerGlobalZ0; //Copper II
- pos[6] = (isBack ? 1 : -1)*CooperIIIZ + layerGlobalZ0; //Copper III
- pos[7] = (isBack ? 1 : -1)*CooperIVZ + layerGlobalZ0; //Copper IV
- pos[8] = (isBack ? 1 : -1)*DriftZ + layerGlobalZ0; //Drift G10 PassiveGasIZ
- pos[9] = (isBack ? 1 : -1)*PassiveGasIZ + layerGlobalZ0; // Passive Gas I
- pos[10] = (isBack ? 1 : -1)*PassiveGasIIZ + layerGlobalZ0; // Passive Gas II
- pos[11] = (isBack ? 1 : -1)*PassiveGasIIIZ + layerGlobalZ0; // Passive Gas III
- pos[12] = (isBack ? 1 : -1)*KaptonIZ + layerGlobalZ0; // Kapton I
- pos[13] = (isBack ? 1 : -1)*KaptonIIZ + layerGlobalZ0; // Kapton II
- pos[14] = (isBack ? 1 : -1)*KaptonIIIZ + layerGlobalZ0; // Kapton III
- pos[15] = (isBack ? 1 : -1)*CopperSilceIZ + layerGlobalZ0; // Copper Slice I
- pos[16] = (isBack ? 1 : -1)*CopperSilceIIZ + layerGlobalZ0; // Copper Slice II
- pos[17] = (isBack ? 1 : -1)*CopperSilceIIIZ + layerGlobalZ0; // Copper Slice III
- pos[18] = (isBack ? 1 : -1)*CopperSilceIVZ + layerGlobalZ0; // Copper Slice IV
- pos[19] = (isBack ? 1 : -1)*CopperSilceVZ + layerGlobalZ0; // Copper Slice V
- pos[20] = (isBack ? 1 : -1)*CopperSilceVIZ + layerGlobalZ0; // Copper Slice VI
+ pos[2] = (isBack ? 1 : -1) * ModuleZ + layerGlobalZ0; //Active volume & Frame
+ pos[3] = (isBack ? 1 : -1) * ReadOutZ + layerGlobalZ0; //Readout G10
+ pos[4] = (isBack ? 1 : -1) * CooperIZ + layerGlobalZ0; //Copper I
+ pos[5] = (isBack ? 1 : -1) * CooperIIZ + layerGlobalZ0; //Copper II
+ pos[6] = (isBack ? 1 : -1) * CooperIIIZ + layerGlobalZ0; //Copper III
+ pos[7] = (isBack ? 1 : -1) * CooperIVZ + layerGlobalZ0; //Copper IV
+ pos[8] = (isBack ? 1 : -1) * DriftZ + layerGlobalZ0; //Drift G10 PassiveGasIZ
+ pos[9] = (isBack ? 1 : -1) * PassiveGasIZ + layerGlobalZ0; // Passive Gas I
+ pos[10] = (isBack ? 1 : -1) * PassiveGasIIZ + layerGlobalZ0; // Passive Gas II
+ pos[11] = (isBack ? 1 : -1) * PassiveGasIIIZ + layerGlobalZ0; // Passive Gas III
+ pos[12] = (isBack ? 1 : -1) * KaptonIZ + layerGlobalZ0; // Kapton I
+ pos[13] = (isBack ? 1 : -1) * KaptonIIZ + layerGlobalZ0; // Kapton II
+ pos[14] = (isBack ? 1 : -1) * KaptonIIIZ + layerGlobalZ0; // Kapton III
+ pos[15] = (isBack ? 1 : -1) * CopperSilceIZ + layerGlobalZ0; // Copper Slice I
+ pos[16] = (isBack ? 1 : -1) * CopperSilceIIZ + layerGlobalZ0; // Copper Slice II
+ pos[17] = (isBack ? 1 : -1) * CopperSilceIIIZ + layerGlobalZ0; // Copper Slice III
+ pos[18] = (isBack ? 1 : -1) * CopperSilceIVZ + layerGlobalZ0; // Copper Slice IV
+ pos[19] = (isBack ? 1 : -1) * CopperSilceVZ + layerGlobalZ0; // Copper Slice V
+ pos[20] = (isBack ? 1 : -1) * CopperSilceVIZ + layerGlobalZ0; // Copper Slice VI
//define media
- TGeoMedium* argon = gGeoMan->GetMedium(activemedium); // active medium
- TGeoMedium* noryl = gGeoMan->GetMedium(spacermedium); // frame medium
- TGeoMedium* g10plate = gGeoMan->GetMedium(g10); // G10 medium
- TGeoMedium* copperplate = gGeoMan->GetMedium(copper); // copper
- TGeoMedium* kapton = gGeoMan->GetMedium(Kapton); // Kapton
-
-
+ TGeoMedium* argon = gGeoMan->GetMedium(activemedium); // active medium
+ TGeoMedium* noryl = gGeoMan->GetMedium(spacermedium); // frame medium
+ TGeoMedium* g10plate = gGeoMan->GetMedium(g10); // G10 medium
+ TGeoMedium* copperplate = gGeoMan->GetMedium(copper); // copper
+ TGeoMedium* kapton = gGeoMan->GetMedium(Kapton); // Kapton
+
//Active Volume
- TGeoTrap* shapeActive = new TGeoTrap(fActiveVolumeDz/2.,0,0,dy,dx1,dx2,0,dy,dx1,dx2,0);
+ TGeoTrap* shapeActive = new TGeoTrap(fActiveVolumeDz / 2., 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
shapeActive->SetName(Form("shStation%02iLayer%i%cModule%03iActiveNoHole", istn, ily, cside, iModule));
- TString activeName = Form("muchstation%02ilayer%i%cactive%03igasArgon",istn+1,ily+1,cside,iMod+1);
- TGeoVolume* voActive = new TGeoVolume(activeName,shapeActive,argon);
+ TString activeName = Form("muchstation%02ilayer%i%cactive%03igasArgon", istn + 1, ily + 1, cside, iMod + 1);
+ TGeoVolume* voActive = new TGeoVolume(activeName, shapeActive, argon);
voActive->SetLineColor(kGreen);
-
+
//Frame
- TGeoTrap* shapeFrame = new TGeoTrap(fFrameDz/2.,0,0,sdy,sdx1,sdx2,0,sdy,sdx1,sdx2,0);
+ TGeoTrap* shapeFrame = new TGeoTrap(fFrameDz / 2., 0, 0, sdy, sdx1, sdx2, 0, sdy, sdx1, sdx2, 0);
shapeFrame->SetName(Form("shStation%02iLayer%i%cModule%03iFullFrameNoHole", istn, ily, cside, iModule));
-
- TString expression = Form("shStation%02iLayer%i%cModule%03iFullFrameNoHole-shStation%02iLayer%i%cModule%03iActiveNoHole", istn, ily, cside, iModule, istn, ily, cside, iModule);
-
- TGeoCompositeShape* shFrame = new TGeoCompositeShape(Form("shStation%02iLayer%i%cModule%03iFinalFrameNoHole", istn, ily, cside, iModule), expression);
-
- TString frameName = Form("muchstation%02ilayer%i%cframe%03i",istn+1,ily+1,cside,iModule+1);
-
- TGeoVolume* voFrame = new TGeoVolume(frameName,shFrame,noryl); // add a name to the frame
+
+ TString expression =
+ Form("shStation%02iLayer%i%cModule%03iFullFrameNoHole-shStation%02iLayer%i%cModule%03iActiveNoHole", istn, ily,
+ cside, iModule, istn, ily, cside, iModule);
+
+ TGeoCompositeShape* shFrame = new TGeoCompositeShape(
+ Form("shStation%02iLayer%i%cModule%03iFinalFrameNoHole", istn, ily, cside, iModule), expression);
+
+ TString frameName = Form("muchstation%02ilayer%i%cframe%03i", istn + 1, ily + 1, cside, iModule + 1);
+
+ TGeoVolume* voFrame = new TGeoVolume(frameName, shFrame, noryl); // add a name to the frame
voFrame->SetLineColor(kMagenta);
//Readout Plate
- TGeoTrap* shapeReadOut = new TGeoTrap(fReadoutPlateDz/2.0,0,0,sdy,sdx1,sdx2,0,sdy,sdx1,sdx2,0);
+ TGeoTrap* shapeReadOut = new TGeoTrap(fReadoutPlateDz / 2.0, 0, 0, sdy, sdx1, sdx2, 0, sdy, sdx1, sdx2, 0);
shapeReadOut->SetName(Form("shStation%02iLayer%i%cModule%03iReadOut", istn, ily, cside, iModule));
-
- TString ReadOutName = Form("muchstation%02ilayer%i%cReadOut%03i",istn+1,ily+1,cside,iModule+1);
-
- TGeoVolume* voReadOut = new TGeoVolume(ReadOutName,shapeReadOut,g10plate);
+
+ TString ReadOutName = Form("muchstation%02ilayer%i%cReadOut%03i", istn + 1, ily + 1, cside, iModule + 1);
+
+ TGeoVolume* voReadOut = new TGeoVolume(ReadOutName, shapeReadOut, g10plate);
voReadOut->SetLineColor(2);
-
-
+
+
//4 Copper 65 micron
TGeoTrap* shapeCopper[4];
- TString CopperName[4];
+ TString CopperName[4];
TGeoVolume* voCopper[4];
- for(Int_t iCop =0; iCop <4; iCop++){
- shapeCopper[iCop] = new TGeoTrap(fCopperDz/2.0,0,0,sdy,sdx1,sdx2,0,sdy,sdx1,sdx2,0);
+ for (Int_t iCop = 0; iCop < 4; iCop++) {
+ shapeCopper[iCop] = new TGeoTrap(fCopperDz / 2.0, 0, 0, sdy, sdx1, sdx2, 0, sdy, sdx1, sdx2, 0);
shapeCopper[iCop]->SetName(Form("shStation%02iLayer%i%cModule%03iCopper", istn, ily, cside, iModule));
-
- CopperName[iCop] = Form("muchstation%02ilayer%i%cCopper%iModule%03i",istn+1,ily+1,cside,iCop+1,iModule+1);
-
- voCopper[iCop] = new TGeoVolume(CopperName[iCop],shapeCopper[iCop],copperplate);
+
+ CopperName[iCop] =
+ Form("muchstation%02ilayer%i%cCopper%iModule%03i", istn + 1, ily + 1, cside, iCop + 1, iModule + 1);
+
+ voCopper[iCop] = new TGeoVolume(CopperName[iCop], shapeCopper[iCop], copperplate);
voCopper[iCop]->SetLineColor(3);
}
-
+
//Drift Plate
- TGeoTrap* shapeDrift = new TGeoTrap(fDriftPlateDz/2.0,0,0,sdy,sdx1,sdx2,0,sdy,sdx1,sdx2,0);
+ TGeoTrap* shapeDrift = new TGeoTrap(fDriftPlateDz / 2.0, 0, 0, sdy, sdx1, sdx2, 0, sdy, sdx1, sdx2, 0);
shapeDrift->SetName(Form("shStation%02iLayer%i%cModule%03iDrift", istn, ily, cside, iModule));
-
- TString DriftName = Form("muchstation%02ilayer%i%cDrift%03i",istn+1,ily+1,cside,iModule+1);
-
- TGeoVolume* voDrift = new TGeoVolume(DriftName,shapeDrift,g10plate);
+
+ TString DriftName = Form("muchstation%02ilayer%i%cDrift%03i", istn + 1, ily + 1, cside, iModule + 1);
+
+ TGeoVolume* voDrift = new TGeoVolume(DriftName, shapeDrift, g10plate);
voDrift->SetLineColor(2);
-
+
//Passive Gas (2mm Ar)
TGeoTrap* shapePassiveGas[3];
TString PassiveGasName[3];
TGeoVolume* voPassiveGas[3];
- for(Int_t iPGas =0; iPGas<3; iPGas++){
- shapePassiveGas[iPGas] = new TGeoTrap(fPassiveVolumeDz/2.,0,0,dy,dx1,dx2,0,dy,dx1,dx2,0);
- shapePassiveGas[iPGas]->SetName(Form("shStation%02iLayer%i%cModule%03iPassiveGasNoHole", istn, ily, cside, iModule));
-
- PassiveGasName[iPGas] = Form("muchstation%02ilayer%i%cPassiveGas%iModule%03i",istn+1,ily+1,cside,iPGas,iModule+1);
- voPassiveGas[iPGas] = new TGeoVolume(PassiveGasName[iPGas],shapePassiveGas[iPGas],argon);
+ for (Int_t iPGas = 0; iPGas < 3; iPGas++) {
+ shapePassiveGas[iPGas] = new TGeoTrap(fPassiveVolumeDz / 2., 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
+ shapePassiveGas[iPGas]->SetName(
+ Form("shStation%02iLayer%i%cModule%03iPassiveGasNoHole", istn, ily, cside, iModule));
+
+ PassiveGasName[iPGas] =
+ Form("muchstation%02ilayer%i%cPassiveGas%iModule%03i", istn + 1, ily + 1, cside, iPGas, iModule + 1);
+ voPassiveGas[iPGas] = new TGeoVolume(PassiveGasName[iPGas], shapePassiveGas[iPGas], argon);
voPassiveGas[iPGas]->SetLineColor(kGreen);
}
-
+
//kapton (50 micron)
TGeoTrap* shapeKapton[3];
- TString KaptonName[3];
+ TString KaptonName[3];
TGeoVolume* voKapton[3];
- for(Int_t iCop =0; iCop <3; iCop++){
- shapeKapton[iCop] = new TGeoTrap(fKaptonDz/2.0,0,0,dy,dx1,dx2,0,dy,dx1,dx2,0);
+ for (Int_t iCop = 0; iCop < 3; iCop++) {
+ shapeKapton[iCop] = new TGeoTrap(fKaptonDz / 2.0, 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
shapeKapton[iCop]->SetName(Form("shStation%02iLayer%i%cModule%03iKapton", istn, ily, cside, iModule));
-
- KaptonName[iCop] = Form("muchstation%02ilayer%i%cKapton%iModule%03i",istn+1,ily+1,cside,iCop+1,iModule+1);
-
- voKapton[iCop] = new TGeoVolume(KaptonName[iCop],shapeKapton[iCop],kapton);
+
+ KaptonName[iCop] =
+ Form("muchstation%02ilayer%i%cKapton%iModule%03i", istn + 1, ily + 1, cside, iCop + 1, iModule + 1);
+
+ voKapton[iCop] = new TGeoVolume(KaptonName[iCop], shapeKapton[iCop], kapton);
voKapton[iCop]->SetLineColor(4);
}
-
+
//6 Copper Slice (5 micron) fCopperSliceDz
TGeoTrap* shapeCopperSlice[6];
- TString CopperSliceName[6];
+ TString CopperSliceName[6];
TGeoVolume* voCopperSlice[6];
- for(Int_t iCop =0; iCop <6; iCop++){
- shapeCopperSlice[iCop] = new TGeoTrap(fCopperSliceDz/2.0,0,0,dy,dx1,dx2,0,dy,dx1,dx2,0);
+ for (Int_t iCop = 0; iCop < 6; iCop++) {
+ shapeCopperSlice[iCop] = new TGeoTrap(fCopperSliceDz / 2.0, 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
shapeCopperSlice[iCop]->SetName(Form("shStation%02iLayer%i%cModule%03iCopperSlice", istn, ily, cside, iModule));
-
- CopperSliceName[iCop] = Form("muchstation%02ilayer%i%cCopperSlice%iModule%03i",istn+1,ily+1,cside,iCop+1,iModule+1);
-
- voCopperSlice[iCop] = new TGeoVolume(CopperSliceName[iCop],shapeCopperSlice[iCop],copperplate);
+
+ CopperSliceName[iCop] =
+ Form("muchstation%02ilayer%i%cCopperSlice%iModule%03i", istn + 1, ily + 1, cside, iCop + 1, iModule + 1);
+
+ voCopperSlice[iCop] = new TGeoVolume(CopperSliceName[iCop], shapeCopperSlice[iCop], copperplate);
voCopperSlice[iCop]->SetLineColor(3);
}
-
-
-
+
+
// Calculate the phi angle of the sector where it has to be placed
- Double_t angle = 180. / TMath::Pi() * phi; // convert angle phi from rad to deg
-
-
- TGeoRotation *r2 = new TGeoRotation("r2");
- //rotate in the vertical plane (per to z axis) with angle
- r2->RotateZ(angle);
-
-
- TGeoTranslation *trans[20];
- TGeoHMatrix *incline_mod[20];
-
- for(Int_t i=0; i<19; i++){
- trans[i] = new TGeoTranslation("",pos[0],pos[1],pos[i+2]);
- incline_mod[i] = new TGeoHMatrix("");
- (*incline_mod[i]) = (*trans[i]) * (*r2);
- }
- volayer->AddNode(voActive, iModule, incline_mod[0]); // add active volume
+ Double_t angle = 180. / TMath::Pi() * phi; // convert angle phi from rad to deg
+
+
+ TGeoRotation* r2 = new TGeoRotation("r2");
+ //rotate in the vertical plane (per to z axis) with angle
+ r2->RotateZ(angle);
+
+
+ TGeoTranslation* trans[20];
+ TGeoHMatrix* incline_mod[20];
+
+ for (Int_t i = 0; i < 19; i++) {
+ trans[i] = new TGeoTranslation("", pos[0], pos[1], pos[i + 2]);
+ incline_mod[i] = new TGeoHMatrix("");
+ (*incline_mod[i]) = (*trans[i]) * (*r2);
+ }
+ volayer->AddNode(voActive, iModule, incline_mod[0]); // add active volume
// volayer->AddNode(voFrame, iModule, incline_mod[0]); // add frame // Frame removed from Active gas
- volayer->AddNode(voReadOut, iModule, incline_mod[1]); // add Read Out
- for(int iNode =0; iNode<4; iNode++)volayer->AddNode(voCopper[iNode], iModule, incline_mod[iNode+2]); //add Copper
- volayer->AddNode(voDrift, iModule, incline_mod[6]); // add Drift Out
- for(int iNode =0; iNode<3; iNode++)volayer->AddNode(voPassiveGas[iNode], iModule, incline_mod[iNode+7]); //add Passive Gas
- for(int iNode =0; iNode<3; iNode++)volayer->AddNode(voKapton[iNode], iModule, incline_mod[iNode+10]); //add Kapton
- for(int iNode =0; iNode<6; iNode++)volayer->AddNode(voCopperSlice[iNode], iModule, incline_mod[iNode+13]); //add CopperSlices
-
+ volayer->AddNode(voReadOut, iModule, incline_mod[1]); // add Read Out
+ for (int iNode = 0; iNode < 4; iNode++)
+ volayer->AddNode(voCopper[iNode], iModule, incline_mod[iNode + 2]); //add Copper
+ volayer->AddNode(voDrift, iModule, incline_mod[6]); // add Drift Out
+ for (int iNode = 0; iNode < 3; iNode++)
+ volayer->AddNode(voPassiveGas[iNode], iModule, incline_mod[iNode + 7]); //add Passive Gas
+ for (int iNode = 0; iNode < 3; iNode++)
+ volayer->AddNode(voKapton[iNode], iModule, incline_mod[iNode + 10]); //add Kapton
+ for (int iNode = 0; iNode < 6; iNode++)
+ volayer->AddNode(voCopperSlice[iNode], iModule, incline_mod[iNode + 13]); //add CopperSlices
}
-
+
return volayer;
}
+TGeoVolume* CreateRpcLayers(int istn, int ily)
+{
-TGeoVolume * CreateRpcLayers(int istn, int ily){
-
- TString layerName = Form("muchstation%02ilayer%i",istn+1,ily+1);
+ TString layerName = Form("muchstation%02ilayer%i", istn + 1, ily + 1);
TGeoVolumeAssembly* volayer = new TGeoVolumeAssembly(layerName);
-
+
//RPC Parametrs
- Double_t fActiveVolumeDz = 0.2; //2mm RPC
- Double_t fSupportDz = 0.2; //2 mm Al (cooling + support)
- Double_t fRpcGlassDz = 0.2; // 2mm glass
- Double_t fCopperDz = 0.0035; // 35 micron copper
- Double_t fReadoutPlateDz = 0.3; //3mm G10
- Double_t fDriftPlateDz = 0.3; //3 mm G10
- Double_t fFrameDz = 1.0; // 1cm Al frame
-
-
- Double_t stGlobalZ0 = fStationZ0[istn] + fMuchZ1; //z position of station center (midplane) [cm]
-
- Double_t layerZ0 = (ily - (fNlayers[istn] - 1) / 2.) * fLayersDz[istn];
+ Double_t fActiveVolumeDz = 0.2; //2mm RPC
+ Double_t fSupportDz = 0.2; //2 mm Al (cooling + support)
+ Double_t fRpcGlassDz = 0.2; // 2mm glass
+ Double_t fCopperDz = 0.0035; // 35 micron copper
+ Double_t fReadoutPlateDz = 0.3; //3mm G10
+ Double_t fDriftPlateDz = 0.3; //3 mm G10
+ Double_t fFrameDz = 1.0; // 1cm Al frame
+
+
+ Double_t stGlobalZ0 = fStationZ0[istn] + fMuchZ1; //z position of station center (midplane) [cm]
+
+ Double_t layerZ0 = (ily - (fNlayers[istn] - 1) / 2.) * fLayersDz[istn];
Double_t layerGlobalZ0 = layerZ0 + stGlobalZ0;
-
-
+
+
//Active Gas distance from layer position
- Double_t ModuleZ = fSupportDz/2. + fCopperDz + fReadoutPlateDz + fRpcGlassDz + fActiveVolumeDz/2.;
-
+ Double_t ModuleZ = fSupportDz / 2. + fCopperDz + fReadoutPlateDz + fRpcGlassDz + fActiveVolumeDz / 2.;
+
//RPC Glass
- Double_t GlassDz = fRpcGlassDz/2. + fActiveVolumeDz/2.;
- Double_t GlassIZ = ModuleZ - GlassDz;
- Double_t GlassIIZ = ModuleZ + GlassDz;
-
+ Double_t GlassDz = fRpcGlassDz / 2. + fActiveVolumeDz / 2.;
+ Double_t GlassIZ = ModuleZ - GlassDz;
+ Double_t GlassIIZ = ModuleZ + GlassDz;
+
//ReadOut Plate distance from layer position
- Double_t ReadOutZ = fSupportDz/2. + fCopperDz + fReadoutPlateDz/2.;
+ Double_t ReadOutZ = fSupportDz / 2. + fCopperDz + fReadoutPlateDz / 2.;
//Drift Plate distance from layer position
- Double_t DriftZ = GlassIIZ + fRpcGlassDz/2. + fDriftPlateDz/2.;
+ Double_t DriftZ = GlassIIZ + fRpcGlassDz / 2. + fDriftPlateDz / 2.;
//Copper distance from layer position
- Double_t CooperIZ = fSupportDz/2. + fCopperDz/2.;
- Double_t CooperIIZ = DriftZ + fDriftPlateDz/2. + fCopperDz/2.;
+ Double_t CooperIZ = fSupportDz / 2. + fCopperDz / 2.;
+ Double_t CooperIIZ = DriftZ + fDriftPlateDz / 2. + fCopperDz / 2.;
+
-
//Set Rmin & Rmax
- Double_t stDz = ((fNlayers[istn] - 1) * fLayersDz[istn])/2. + CooperIIZ + fCopperDz/2.;
-
+ Double_t stDz = ((fNlayers[istn] - 1) * fLayersDz[istn]) / 2. + CooperIIZ + fCopperDz / 2.;
+
Double_t stGlobalZ2 = stGlobalZ0 + stDz;
Double_t stGlobalZ1 = stGlobalZ0 - stDz;
-
+
Double_t rmin = stGlobalZ1 * fAcceptanceTanMin;
Double_t rmax = stGlobalZ2 * fAcceptanceTanMax;
-
-
+
+
// Module dimention calculation
- Double_t phi0 = TMath::Pi()/fNSectorsPerLayer[istn]; // azimuthal half widh of each module
- Double_t ymin = rmin+fSpacerR;
+ Double_t phi0 = TMath::Pi() / fNSectorsPerLayer[istn]; // azimuthal half widh of each module
+ Double_t ymin = rmin + fSpacerR;
Double_t ymax = rmax;
-
+
//define the dimensions of the trapezoidal module
- Double_t dy = (ymax-ymin)/2.; //y (length)
- Double_t dx1 = ymin*TMath::Tan(phi0)+fOverlapR/TMath::Cos(phi0); // large x
- Double_t dx2 = ymax*TMath::Tan(phi0)+fOverlapR/TMath::Cos(phi0); // small x
-
+ Double_t dy = (ymax - ymin) / 2.; //y (length)
+ Double_t dx1 = ymin * TMath::Tan(phi0) + fOverlapR / TMath::Cos(phi0); // large x
+ Double_t dx2 = ymax * TMath::Tan(phi0) + fOverlapR / TMath::Cos(phi0); // small x
+
//define the spacer dimensions
- Double_t tg = (dx2-dx1)/2/dy;
- Double_t dd1 = fSpacerPhi*tg;
- Double_t dd2 = fSpacerPhi*sqrt(1+tg*tg);
- Double_t sdx1 = dx1+dd2-dd1-0.1; // 0.1 cm to avoid overlaps
- Double_t sdx2 = dx2+dd2+dd1;
- Double_t sdy = dy+fSpacerR;
-
-
-
- infoFile<<" "<<istn+1<<"\t "<<ily+1<<"\t\t"<<layerGlobalZ0<<"\t"<<fNSectorsPerLayer[istn]<<"\t"<<fActiveVolumeDz<<endl;
-
-
+ Double_t tg = (dx2 - dx1) / 2 / dy;
+ Double_t dd1 = fSpacerPhi * tg;
+ Double_t dd2 = fSpacerPhi * sqrt(1 + tg * tg);
+ Double_t sdx1 = dx1 + dd2 - dd1 - 0.1; // 0.1 cm to avoid overlaps
+ Double_t sdx2 = dx2 + dd2 + dd1;
+ Double_t sdy = dy + fSpacerR;
+
+
+ infoFile << " " << istn + 1 << "\t " << ily + 1 << "\t\t" << layerGlobalZ0 << "\t" << fNSectorsPerLayer[istn]
+ << "\t" << fActiveVolumeDz << endl;
+
+
// Aluminum Plate (Cooling + Support)
- TString supportAlName = Form("shStation%02iSupportAl",istn+1);
- TGeoTube* shSupportAl = new TGeoTube(supportAlName,rmin,rmax,fSupportDz/2.);
-
- TString supportName1 = Form("muchstation%02ilayer%isupportAl",istn+1,ily+1);
- TGeoMedium* coolMat = gGeoMan->GetMedium(Al);
-
- TGeoVolume* voSupport1 = new TGeoVolume(supportName1,shSupportAl,coolMat);
+ TString supportAlName = Form("shStation%02iSupportAl", istn + 1);
+ TGeoTube* shSupportAl = new TGeoTube(supportAlName, rmin, rmax, fSupportDz / 2.);
+
+ TString supportName1 = Form("muchstation%02ilayer%isupportAl", istn + 1, ily + 1);
+ TGeoMedium* coolMat = gGeoMan->GetMedium(Al);
+
+ TGeoVolume* voSupport1 = new TGeoVolume(supportName1, shSupportAl, coolMat);
voSupport1->SetLineColor(kCyan);
-
- TGeoTranslation *support_trans1 = new TGeoTranslation("supportName1", 0,0,layerGlobalZ0);
- volayer->AddNode(voSupport1,0,support_trans1);
-
-// Now start adding the GEM modules
- for (Int_t iModule=0; iModule<fNSectorsPerLayer[istn]; iModule++){
-
+
+ TGeoTranslation* support_trans1 = new TGeoTranslation("supportName1", 0, 0, layerGlobalZ0);
+ volayer->AddNode(voSupport1, 0, support_trans1);
+
+ // Now start adding the GEM modules
+ for (Int_t iModule = 0; iModule < fNSectorsPerLayer[istn]; iModule++) {
+
Double_t phi = 2 * phi0 * (iModule + 0.5); // add 0.5 to not overlap with y-axis for left-right layer separation
- Bool_t isBack = iModule%2;
- Char_t cside = (isBack==1) ? 'b' : 'f';
- Int_t iMod = iModule/2;
+ Bool_t isBack = iModule % 2;
+ Char_t cside = (isBack == 1) ? 'b' : 'f';
+ Int_t iMod = iModule / 2;
// correct the x, y positions
Double_t pos[9];
- pos[0] = -(ymin+dy)*sin(phi);
- pos[1] = (ymin+dy)*cos(phi);
+ pos[0] = -(ymin + dy) * sin(phi);
+ pos[1] = (ymin + dy) * cos(phi);
// different z positions for odd/even modules
- pos[2] = (isBack ? 1 : -1)*ModuleZ + layerGlobalZ0; //Active volume & Frame
- pos[3] = (isBack ? 1 : -1)*GlassIZ + layerGlobalZ0; //Glass I
- pos[4] = (isBack ? 1 : -1)*GlassIIZ + layerGlobalZ0; //Glass II
- pos[5] = (isBack ? 1 : -1)*ReadOutZ + layerGlobalZ0; //ReadOut Plate
- pos[6] = (isBack ? 1 : -1)*DriftZ + layerGlobalZ0; //Drift Plate
- pos[7] = (isBack ? 1 : -1)*CooperIZ + layerGlobalZ0; //Copper I Plate
- pos[8] = (isBack ? 1 : -1)*CooperIIZ + layerGlobalZ0; //Copper II Plate
-
-
-
- TGeoMedium* RPCgasmedium = gGeoMan->GetMedium(RPCm); // RPC Gas
- TGeoMedium* RPCglassmat = gGeoMan->GetMedium(RPCg); // RPC Glass
- TGeoMedium* g10plate = gGeoMan->GetMedium(g10); // G10 medium
- TGeoMedium* copperplate = gGeoMan->GetMedium(copper); // copper
- TGeoMedium* AlFrame = gGeoMan->GetMedium(Al); // spacer medium
-
-
+ pos[2] = (isBack ? 1 : -1) * ModuleZ + layerGlobalZ0; //Active volume & Frame
+ pos[3] = (isBack ? 1 : -1) * GlassIZ + layerGlobalZ0; //Glass I
+ pos[4] = (isBack ? 1 : -1) * GlassIIZ + layerGlobalZ0; //Glass II
+ pos[5] = (isBack ? 1 : -1) * ReadOutZ + layerGlobalZ0; //ReadOut Plate
+ pos[6] = (isBack ? 1 : -1) * DriftZ + layerGlobalZ0; //Drift Plate
+ pos[7] = (isBack ? 1 : -1) * CooperIZ + layerGlobalZ0; //Copper I Plate
+ pos[8] = (isBack ? 1 : -1) * CooperIIZ + layerGlobalZ0; //Copper II Plate
+
+
+ TGeoMedium* RPCgasmedium = gGeoMan->GetMedium(RPCm); // RPC Gas
+ TGeoMedium* RPCglassmat = gGeoMan->GetMedium(RPCg); // RPC Glass
+ TGeoMedium* g10plate = gGeoMan->GetMedium(g10); // G10 medium
+ TGeoMedium* copperplate = gGeoMan->GetMedium(copper); // copper
+ TGeoMedium* AlFrame = gGeoMan->GetMedium(Al); // spacer medium
+
+
//Active Volume
- TGeoTrap* shapeActive = new TGeoTrap(fActiveVolumeDz/2.,0,0,dy,dx1,dx2,0,dy,dx1,dx2,0);
+ TGeoTrap* shapeActive = new TGeoTrap(fActiveVolumeDz / 2., 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
shapeActive->SetName(Form("shStation%02iLayer%i%cModule%03iActiveNoHole", istn, ily, cside, iModule));
-
- TString activeName = Form("muchstation%02ilayer%i%cactive%03irpcgas",istn+1,ily+1,cside,iMod+1);
- TGeoVolume* voActive = new TGeoVolume(activeName,shapeActive,RPCgasmedium);
+
+ TString activeName = Form("muchstation%02ilayer%i%cactive%03irpcgas", istn + 1, ily + 1, cside, iMod + 1);
+ TGeoVolume* voActive = new TGeoVolume(activeName, shapeActive, RPCgasmedium);
voActive->SetLineColor(kGreen);
-
+
// RPC Glass
- TGeoTrap* shapeGlass[2];
- TString GlassName[2];
- TGeoVolume* voGlass[2];
- for(Int_t iGlass =0 ; iGlass<2; iGlass++){
- shapeGlass[iGlass] = new TGeoTrap(fRpcGlassDz/2.,0,0,dy,dx1,dx2,0,dy,dx1,dx2,0);
- shapeGlass[iGlass]->SetName(Form("shStation%02iLayer%i%cModule%03iNoHoleGlass", istn, ily, cside, iModule));
- GlassName[iGlass] = Form("muchstation%02ilayer%i%crpcglass%imodule%03i",istn+1,ily+1,cside,iGlass+1,iModule+1);
- voGlass[iGlass] = new TGeoVolume(GlassName[iGlass],shapeGlass[iGlass],RPCglassmat);
- voGlass[iGlass]->SetLineColor(kRed);
- }
-
- //Readout Plate
- TGeoTrap* shapeReadOut = new TGeoTrap(fReadoutPlateDz/2.0,0,0,dy,dx1,dx2,0,dy,dx1,dx2,0);
- shapeReadOut->SetName(Form("shStation%02iLayer%i%cModule%03iReadOut", istn, ily, cside, iModule));
-
- TString ReadOutName = Form("muchstation%02ilayer%i%cReadOut%03i",istn+1,ily+1,cside,iModule+1);
-
- TGeoVolume* voReadOut = new TGeoVolume(ReadOutName,shapeReadOut,g10plate);
- voReadOut->SetLineColor(2);
-
- //Drift Plate
- TGeoTrap* shapeDrift = new TGeoTrap(fDriftPlateDz/2.0,0,0,dy,dx1,dx2,0,dy,dx1,dx2,0);
- shapeDrift->SetName(Form("shStation%02iLayer%i%cModule%03iDrift", istn, ily, cside, iModule));
-
- TString DriftName = Form("muchstation%02ilayer%i%cDrift%03i",istn+1,ily+1,cside,iModule+1);
-
- TGeoVolume* voDrift = new TGeoVolume(DriftName,shapeDrift,g10plate);
- voDrift->SetLineColor(2);
-
-
- //2 Copper 35 micron
- TGeoTrap* shapeCopper[2];
- TString CopperName[2];
- TGeoVolume* voCopper[2];
- for(Int_t iCop =0; iCop <2; iCop++){
- shapeCopper[iCop] = new TGeoTrap(fCopperDz/2.0,0,0,dy,dx1,dx2,0,dy,dx1,dx2,0);
- shapeCopper[iCop]->SetName(Form("shStation%02iLayer%i%cModule%03iCopper", istn, ily, cside, iModule));
-
- CopperName[iCop] = Form("muchstation%02ilayer%i%cCopper%iModule%03i",istn+1,ily+1,cside,iCop+1,iModule+1);
-
- voCopper[iCop] = new TGeoVolume(CopperName[iCop],shapeCopper[iCop],copperplate);
- voCopper[iCop]->SetLineColor(kRed);
- }
-
-
- //Frame (1 cm Al)
-
- TGeoTrap* shape = new TGeoTrap(fFrameDz/2.,0,0,dy,dx1,dx2,0,dy,dx1,dx2,0);
- shape->SetName(Form("shStation%02iLayer%i%cModule%03iFrameNoHole", istn, ily, cside, iModule));
-
- TGeoTrap* shapeFrame = new TGeoTrap(fFrameDz/2.,0,0,sdy,sdx1,sdx2,0,sdy,sdx1,sdx2,0);
- shapeFrame->SetName(Form("shStation%02iLayer%i%cModule%03iFullFrameNoHole", istn, ily, cside, iModule));
-
- TString expression = Form("shStation%02iLayer%i%cModule%03iFullFrameNoHole-shStation%02iLayer%i%cModule%03iFrameNoHole", istn, ily, cside, iModule, istn, ily, cside, iModule);
-
- TGeoCompositeShape* shFrame = new TGeoCompositeShape(Form("shStation%02iLayer%i%cModule%03iFinalFrameNoHole", istn, ily, cside, iModule), expression);
-
- TString frameName = Form("muchstation%02ilayer%i%cframe%03i",istn+1,ily+1,cside,iModule+1);
-
- TGeoVolume* voFrame = new TGeoVolume(frameName,shFrame,AlFrame); // add a name to the frame
- voFrame->SetLineColor(kMagenta);
-
-
-
- // Calculate the phi angle of the sector where it has to be placed
- Double_t angle = 180. / TMath::Pi() * phi; // convert angle phi from rad to deg
-
-
- TGeoRotation *r2 = new TGeoRotation("r2");
- //rotate in the vertical plane (per to z axis) with angle
- r2->RotateZ(angle);
-
-
- TGeoTranslation *trans[8];
- TGeoHMatrix *incline_mod[8];
-
- for(int i=0; i<7; i++){
- trans[i] = new TGeoTranslation("",pos[0],pos[1],pos[i+2]);
-
- incline_mod[i] = new TGeoHMatrix("");
- (*incline_mod[i]) = (*trans[i]) * (*r2);
- }
-
- volayer->AddNode(voActive, iModule, incline_mod[0]); // add active volume
- // volayer->AddNode(voFrame, iModule, incline_mod[0]); // add frame
- for(int iNode =0; iNode<2; iNode++)volayer->AddNode(voGlass[iNode], iModule, incline_mod[iNode+1]); //add Glass
- volayer->AddNode(voReadOut, iModule, incline_mod[3]); // add ReadOut Plate
- volayer->AddNode(voDrift, iModule, incline_mod[4]); // add Drift Plate
- for(int iNode =0; iNode<2; iNode++)volayer->AddNode(voCopper[iNode], iModule, incline_mod[iNode+5]); //add Copper
+ TGeoTrap* shapeGlass[2];
+ TString GlassName[2];
+ TGeoVolume* voGlass[2];
+ for (Int_t iGlass = 0; iGlass < 2; iGlass++) {
+ shapeGlass[iGlass] = new TGeoTrap(fRpcGlassDz / 2., 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
+ shapeGlass[iGlass]->SetName(Form("shStation%02iLayer%i%cModule%03iNoHoleGlass", istn, ily, cside, iModule));
+ GlassName[iGlass] =
+ Form("muchstation%02ilayer%i%crpcglass%imodule%03i", istn + 1, ily + 1, cside, iGlass + 1, iModule + 1);
+ voGlass[iGlass] = new TGeoVolume(GlassName[iGlass], shapeGlass[iGlass], RPCglassmat);
+ voGlass[iGlass]->SetLineColor(kRed);
+ }
+
+ //Readout Plate
+ TGeoTrap* shapeReadOut = new TGeoTrap(fReadoutPlateDz / 2.0, 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
+ shapeReadOut->SetName(Form("shStation%02iLayer%i%cModule%03iReadOut", istn, ily, cside, iModule));
+
+ TString ReadOutName = Form("muchstation%02ilayer%i%cReadOut%03i", istn + 1, ily + 1, cside, iModule + 1);
+
+ TGeoVolume* voReadOut = new TGeoVolume(ReadOutName, shapeReadOut, g10plate);
+ voReadOut->SetLineColor(2);
+
+ //Drift Plate
+ TGeoTrap* shapeDrift = new TGeoTrap(fDriftPlateDz / 2.0, 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
+ shapeDrift->SetName(Form("shStation%02iLayer%i%cModule%03iDrift", istn, ily, cside, iModule));
+
+ TString DriftName = Form("muchstation%02ilayer%i%cDrift%03i", istn + 1, ily + 1, cside, iModule + 1);
+
+ TGeoVolume* voDrift = new TGeoVolume(DriftName, shapeDrift, g10plate);
+ voDrift->SetLineColor(2);
+
+
+ //2 Copper 35 micron
+ TGeoTrap* shapeCopper[2];
+ TString CopperName[2];
+ TGeoVolume* voCopper[2];
+ for (Int_t iCop = 0; iCop < 2; iCop++) {
+ shapeCopper[iCop] = new TGeoTrap(fCopperDz / 2.0, 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
+ shapeCopper[iCop]->SetName(Form("shStation%02iLayer%i%cModule%03iCopper", istn, ily, cside, iModule));
+
+ CopperName[iCop] =
+ Form("muchstation%02ilayer%i%cCopper%iModule%03i", istn + 1, ily + 1, cside, iCop + 1, iModule + 1);
+
+ voCopper[iCop] = new TGeoVolume(CopperName[iCop], shapeCopper[iCop], copperplate);
+ voCopper[iCop]->SetLineColor(kRed);
+ }
+
+
+ //Frame (1 cm Al)
+
+ TGeoTrap* shape = new TGeoTrap(fFrameDz / 2., 0, 0, dy, dx1, dx2, 0, dy, dx1, dx2, 0);
+ shape->SetName(Form("shStation%02iLayer%i%cModule%03iFrameNoHole", istn, ily, cside, iModule));
+
+ TGeoTrap* shapeFrame = new TGeoTrap(fFrameDz / 2., 0, 0, sdy, sdx1, sdx2, 0, sdy, sdx1, sdx2, 0);
+ shapeFrame->SetName(Form("shStation%02iLayer%i%cModule%03iFullFrameNoHole", istn, ily, cside, iModule));
+
+ TString expression =
+ Form("shStation%02iLayer%i%cModule%03iFullFrameNoHole-shStation%02iLayer%i%cModule%03iFrameNoHole", istn, ily,
+ cside, iModule, istn, ily, cside, iModule);
+
+ TGeoCompositeShape* shFrame = new TGeoCompositeShape(
+ Form("shStation%02iLayer%i%cModule%03iFinalFrameNoHole", istn, ily, cside, iModule), expression);
+
+ TString frameName = Form("muchstation%02ilayer%i%cframe%03i", istn + 1, ily + 1, cside, iModule + 1);
+
+ TGeoVolume* voFrame = new TGeoVolume(frameName, shFrame, AlFrame); // add a name to the frame
+ voFrame->SetLineColor(kMagenta);
+
+
+ // Calculate the phi angle of the sector where it has to be placed
+ Double_t angle = 180. / TMath::Pi() * phi; // convert angle phi from rad to deg
+
+
+ TGeoRotation* r2 = new TGeoRotation("r2");
+ //rotate in the vertical plane (per to z axis) with angle
+ r2->RotateZ(angle);
+
+
+ TGeoTranslation* trans[8];
+ TGeoHMatrix* incline_mod[8];
+
+ for (int i = 0; i < 7; i++) {
+ trans[i] = new TGeoTranslation("", pos[0], pos[1], pos[i + 2]);
+
+ incline_mod[i] = new TGeoHMatrix("");
+ (*incline_mod[i]) = (*trans[i]) * (*r2);
+ }
+
+ volayer->AddNode(voActive, iModule, incline_mod[0]); // add active volume
+ // volayer->AddNode(voFrame, iModule, incline_mod[0]); // add frame
+ for (int iNode = 0; iNode < 2; iNode++)
+ volayer->AddNode(voGlass[iNode], iModule, incline_mod[iNode + 1]); //add Glass
+ volayer->AddNode(voReadOut, iModule, incline_mod[3]); // add ReadOut Plate
+ volayer->AddNode(voDrift, iModule, incline_mod[4]); // add Drift Plate
+ for (int iNode = 0; iNode < 2; iNode++)
+ volayer->AddNode(voCopper[iNode], iModule, incline_mod[iNode + 5]); //add Copper
}
return volayer;
}
--
GitLab