#include <TSystem.h>
#include <TVirtualMC.h>
-
+#include <TString.h>
+#include <TGeoManager.h>
#include "AliFRAMEv2.h"
#include "AliMagF.h"
#include "AliRun.h"
#include "AliConst.h"
#include "AliMC.h"
+#include "AliLog.h"
ClassImp(AliFRAMEv2)
Float_t dx, dy, dz;
Int_t i, j, jmod;
+ jmod = 0;
//
// Constants
const Float_t kEps = 0.01;
ppgon[9] = ppgon[6];
gMC->Gsvolu("BREF", "PGON", kAir, ppgon, 10);
gMC->Gspos("BREF", 1, "B077", 0., 0., 0., 0, "ONLY");
+ gGeoManager->GetVolume("BREF")->SetVisibility(kFALSE);
//
// The outer Frame
//
// The internal frame
//
- char* module[3] = {"B071\0", "B074\0", "B075\0"};
//
//
// Mother Volumes
//
ptrd1[0] = 50.18;
ptrd1[1] = 70.35;
-// ptrd1[0] = 50.10;
-// ptrd1[1] = 70.20;
-
ptrd1[2] = 375.5;
ptrd1[3] = 57.2;
- for (jmod = 0; jmod < 3; jmod++)
- {
- gMC->Gsvolu(module[jmod], "TRD1", kAir, ptrd1, 4);
- }
-// Position of Holes for PHOS (P) and RICH (R) starting at 6h
-// P P P - - R R R - - - - - - - - P P
- Int_t mod[18] = {1, 1, 1, 0, 0, 2, 2, 2, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1};
-
Float_t r = 341.8;
Float_t rout1 = 410.564;
Float_t rout2 = 415.2;
- Int_t modcount[3] = {0, 0, 0};
+ TString module[18];
+// Position of Holes for PHOS (P) and RICH (R) starting at 6h
+// P P P - - R R R - - - - - - - - P P
+// Int_t mod[18] = {0, 0, 0, 0, 0, 2, 2, 2, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1};
- for (i=0; i<18; i++) {
- Float_t phi = i*20.;
- Float_t phi2 = 270+phi;
- if (phi2 >= 360.) phi2-=360.;
+ for (i = 0; i < 18; i++) {
+ // Create volume i
+ char name[16];
+ sprintf(name, "BSEGMO%d", i);
+ gMC->Gsvolu(name, "TRD1", kAir, ptrd1, 4);
+ gGeoManager->GetVolume(name)->SetVisibility(kFALSE);
+
+ module[i] = name;
+ // Place volume i
+ Float_t phi = i * 20.;
+ Float_t phi2 = 270 + phi;
+ if (phi2 >= 360.) phi2 -= 360.;
dx = TMath::Sin(phi*kdeg2rad)*r;
dy = -TMath::Cos(phi*kdeg2rad)*r;
- modcount[mod[i]]++;
AliMatrix(idrotm[2034+i], 90.0, phi, 0., 0., 90., phi2);
- gMC->Gspos(module[mod[i]], modcount[mod[i]], "B077", dx, dy, 0., idrotm[2034+i], "ONLY");
+ gMC->Gspos(name, 1, "B077", dx, dy, 0., idrotm[2034+i], "ONLY");
//
// Position elements of outer Frame
//
dz = -iFrH/2.+ringH/2.+kEps;
- for (i = 0; i< 3; i++)
+ for (jmod = 0; jmod< 18; jmod++)
{
// ring bars
- for (jmod = 0; jmod<3; jmod++) {
- gMC->Gspos("B072", 6*i+jmod+1, module[jmod], 0, dymodL[i], dz, 0, "ONLY");
- gMC->Gspos("B072", 6*i+jmod+4, module[jmod], 0, -dymodL[i], dz, idrotm[2070], "ONLY");
+ for (i = 0; i < 3; i++) {
+ gMC->Gspos("B072", 6*jmod+i+1, module[jmod], 0, dymodL[i], dz, 0, "ONLY");
+ gMC->Gspos("B072", 6*jmod+i+4, module[jmod], 0, -dymodL[i], dz, idrotm[2070], "ONLY");
// 0-deg web
- if (i==2) {
- gMC->Gspos("B263", 12*i+jmod+1, module[jmod], 60.0732, dymodU[i], 4.6669, idrotm[2072], "ONLY");
- gMC->Gspos("B263", 12*i+jmod+4, module[jmod], 60.0732, -dymodU[i], 4.6669, idrotm[2071], "ONLY");
- gMC->Gspos("B263", 12*i+jmod+7, module[jmod], -60.0732, dymodU[i], 4.6669, idrotm[2074], "ONLY");
- gMC->Gspos("B263", 12*i+jmod+10, module[jmod], -60.0732, -dymodU[i], 4.6669, idrotm[2073], "ONLY");
+ if (i == 2) {
+ gMC->Gspos("B263", 4*jmod+1, module[jmod], 60.0732, dymodU[2], 4.6669, idrotm[2072], "ONLY");
+ gMC->Gspos("B263", 4*jmod+2, module[jmod], 60.0732, -dymodU[2], 4.6669, idrotm[2071], "ONLY");
+ gMC->Gspos("B263", 4*jmod+3, module[jmod], -60.0732, dymodU[2], 4.6669, idrotm[2074], "ONLY");
+ gMC->Gspos("B263", 4*jmod+4, module[jmod], -60.0732, -dymodU[2], 4.6669, idrotm[2073], "ONLY");
}
}
}
// outer diagonal web
dy = (dymodU[0]+dymodL[0])/2.;
- for (jmod = 0; jmod<3; jmod++) {
+ for (jmod = 0; jmod < 18; jmod++) {
gMC->Gspos("B063", 4*jmod+1, module[jmod], 60.0732, dy, 4.6669, idrotm[2086], "ONLY");
gMC->Gspos("B063", 4*jmod+2, module[jmod], 60.0732, -dy, 4.6669, idrotm[2087], "ONLY");
gMC->Gspos("B063", 4*jmod+3, module[jmod], -60.0732, dy, 4.6669, idrotm[2088], "ONLY");
}
dy = (dymodU[1]+dymodL[1])/2.;
- for (jmod = 0; jmod<3; jmod++) {
+ for (jmod = 0; jmod < 18; jmod++) {
gMC->Gspos("B163", 4*jmod+1, module[jmod], 60.0732, dy, 4.6669, idrotm[2080], "ONLY");
gMC->Gspos("B163", 4*jmod+2, module[jmod], 60.0732, -dy, 4.6669, idrotm[2079], "ONLY");
gMC->Gspos("B163", 4*jmod+3, module[jmod], -60.0732, dy, 4.6669, idrotm[2082], "ONLY");
dy = (dymodL[2]+dymodU[1])/2.-dwl/2.;
- for (jmod = 0; jmod<3; jmod++) {
+ for (jmod = 0; jmod < 18; jmod++) {
gMC->Gspos("B065", 4*jmod+1, module[jmod], 60.0732, dy, 4.6669, idrotm[2076], "ONLY");
gMC->Gspos("B065", 4*jmod+2, module[jmod], 60.0732, -dy, 4.6669, idrotm[2075], "ONLY");
gMC->Gspos("B065", 4*jmod+3, module[jmod], -60.0732, dy, 4.6669, idrotm[2078], "ONLY");
dy = (dymodL[1]+dymodU[0])/2.;
- for (jmod = 0; jmod<3; jmod++) {
+ for (jmod = 0; jmod < 18; jmod++) {
gMC->Gspos("B067", 4*jmod+1, module[jmod], 60.0732, dy, 4.6669, idrotm[2076], "ONLY");
gMC->Gspos("B067", 4*jmod+2, module[jmod], 60.0732, -dy, 4.6669, idrotm[2075], "ONLY");
gMC->Gspos("B067", 4*jmod+3, module[jmod], -60.0732, dy, 4.6669, idrotm[2078], "ONLY");
dz = -iFrH/2.+ringH+longH/2.;
- for (jmod = 0; jmod<3; jmod++) {
+ for (jmod = 0; jmod < 18; jmod++) {
gMC->Gspos("BA59", 2*jmod+1, module[jmod], 49.6476, 0.0, dz, idrotm[2084], "ONLY");
gMC->Gspos("BA59", 2*jmod+2, module[jmod], -49.6476, 0.0, dz, idrotm[2083], "ONLY");
}
//
//
// Spacer between TRD and TOF modules
+/*
ptrap[0] = 0.99;
ptrap[1] = 10.00;
ptrap[2] = 0.00;
ptrap[9] = ptrap[4];
ptrap[10] = 0;
gMC->Gsvolu("B056", "TRAP", kSteel, ptrap, 11);
- for (jmod = 0; jmod<3; jmod++) {
+ for (jmod = 0; jmod < 18; jmod++) {
gMC->Gspos("B056", 2*jmod+1, module[jmod], 61.9607, 0.0, 27.2, 0 , "ONLY");
gMC->Gspos("B056", 2*jmod+2, module[jmod], -61.9607, 0.0, 27.2, idrotm[2085], "ONLY");
}
-//
-// Mother volumes for TRD and TOF
-//
- if (!fHoles) {
- printf("\n FRAME Version without Holes !");
- ptrd1[0] = 49.8065;
- ptrd1[1] = 62.8535;
- ptrd1[2] = 375.5;
- ptrd1[3] = 37.;
- gMC->Gsvolu("BTR1", "TRD1", kAir, ptrd1, 4);
- gMC->Gsvolu("BTR2", "TRD1", kAir, ptrd1, 4);
- gMC->Gsvolu("BTR3", "TRD1", kAir, ptrd1, 4);
-
- ptrd1[0] = 63.2061;
- ptrd1[1] = 68.3192;
- ptrd1[2] = 375.5;
- ptrd1[3] = 14.5;
- gMC->Gsvolu("BTO1", "TRD1", kAir, ptrd1, 4);
- gMC->Gsvolu("BTO2", "TRD1", kAir, ptrd1, 4);
- gMC->Gsvolu("BTO3", "TRD1", kAir, ptrd1, 4);
-
- gMC->Gspos("BTR1", 1, "B071", 0., 0., -10.8, 0, "ONLY");
- gMC->Gspos("BTR2", 1, "B074", 0., 0., -10.8, 0, "ONLY");
- gMC->Gspos("BTR3", 1, "B075", 0., 0., -10.8, 0, "ONLY");
-
- gMC->Gspos("BTO1", 1, "B071", 0., 0., 42.69, 0, "ONLY");
- gMC->Gspos("BTO2", 1, "B074", 0., 0., 42.69, 0, "ONLY");
- gMC->Gspos("BTO3", 1, "B075", 0., 0., 42.69, 0, "ONLY");
- } else {
- printf("\n FRAME Version with Holes !");
- ptrd1[0] = 49.8065;
- ptrd1[1] = 62.8535;
- ptrd1[2] = 375.5;
- ptrd1[3] = 37;
- gMC->Gsvolu("BTR1", "TRD1", kAir, ptrd1, 4);
- ptrd1[2] = 156.75;
- gMC->Gsvolu("BTR2", "TRD1", kAir, ptrd1, 4);
- ptrd1[2] = 79.75;
- gMC->Gsvolu("BTR3", "TRD1", kAir, ptrd1, 4);
-
- ptrd1[0] = 63.2061;
- ptrd1[1] = 68.3192;
- ptrd1[2] = 375.5;
- ptrd1[3] = 14.5;
- gMC->Gsvolu("BTO1", "TRD1", kAir, ptrd1, 4);
- ptrd1[2] = 156.75;
- gMC->Gsvolu("BTO2", "TRD1", kAir, ptrd1, 4);
- ptrd1[2] = 79.75;
- gMC->Gsvolu("BTO3", "TRD1", kAir, ptrd1, 4);
-
+*/
+ //
+ // TRD mother volumes
+ //
- gMC->Gspos("BTR1", 1, "B071", 0., 0.00, -10.8, 0, "ONLY");
- gMC->Gspos("BTR2", 1, "B074", 0., -218.75, -10.8, idrotm[2070], "ONLY");
- gMC->Gspos("BTR2", 2, "B074", 0., 218.75, -10.8, 0, "ONLY");
- gMC->Gspos("BTR3", 1, "B075", 0., -295.75, -10.8, idrotm[2070], "ONLY");
- gMC->Gspos("BTR3", 2, "B075", 0., 295.75, -10.8, 0, "ONLY");
-
- gMC->Gspos("BTO1", 1, "B071", 0., 0.00, 42.69, 0, "ONLY");
- gMC->Gspos("BTO2", 1, "B074", 0., -218.75, 42.69, idrotm[2070], "ONLY");
- gMC->Gspos("BTO2", 2, "B074", 0., 218.75, 42.69, 0, "ONLY");
- gMC->Gspos("BTO3", 1, "B075", 0., -295.75, 42.69, idrotm[2070], "ONLY");
- gMC->Gspos("BTO3", 2, "B075", 0., 295.75, 42.69, 0, "ONLY");
+ ptrd1[0] = 49.8065;
+ ptrd1[1] = 63.0060;
+ ptrd1[2] = 375.5;
+ ptrd1[3] = 37.43;
+
+// ptrd1[0] = 49.8065;
+// ptrd1[1] = 62.8535;
+// ptrd1[2] = 375.5;
+// ptrd1[3] = 37.;
+
+ for (i = 0; i < 18; i++) {
+ char nameCh[16];
+ sprintf(nameCh, "BTRD%d",i);
+ char nameMo[16];
+ sprintf(nameMo, "BSEGMO%d",i);
+ gMC->Gsvolu(nameCh, "TRD1", kAir, ptrd1, 4);
+ gGeoManager->GetVolume(nameCh)->SetVisibility(kFALSE);
+
+ gMC->Gspos(nameCh, 1, nameMo, 0., 0., -10.37, 0, "ONLY");
}
-
+
+//
+// TOF mother volumes
+//
+ ptrd1[0] = 63.2061;
+ ptrd1[1] = 68.3192;
+ ptrd1[2] = 375.5;
+ ptrd1[3] = 14.5;
+ for (i = 0; i < 18; i++) {
+ char nameCh[16];
+ sprintf(nameCh, "BTOF%d",i);
+ char nameMo[16];
+ sprintf(nameMo, "BSEGMO%d",i);
+ gMC->Gsvolu(nameCh, "TRD1", kAir, ptrd1, 4);
+ gGeoManager->GetVolume(nameCh)->SetVisibility(kFALSE);
+ gMC->Gspos(nameCh, 1, nameMo, 0., 0., 42.69, 0, "ONLY");
+ }
+
//
// Geometry of Rails starts here
//
gMC->Gspos("BRS2", 1, "BRS1", 0., -27.5+3.75, 0., 0, "ONLY");
gMC->Gspos("BRS2", 2, "BRS1", 0., 27.5-3.75, 0., 0, "ONLY");
gMC->Gspos("BRS3", 1, "BRS1", 0., 0., 0., 0, "ONLY");
- gMC->Gspos("BRS1", 1, "ALIC", -430.-3., -180., 0., 0, "ONLY");
- gMC->Gspos("BRS1", 2, "ALIC", 430.+3., -180., 0., 0, "ONLY");
+ gMC->Gspos("BRS1", 1, "ALIC", -430.-3., -190., 0., 0, "ONLY");
+ gMC->Gspos("BRS1", 2, "ALIC", 430.+3., -190., 0., 0, "ONLY");
rbox[0] = 3.0;
- rbox[1] = 145./2.;
+ rbox[1] = 145./4.;
rbox[2] = 25.0;
gMC->Gsvolu("BRS4", "BOX", kSteel, rbox, 3);
- gMC->Gspos("BRS4", 1, "ALIC", 430.+3., -180.+55./2.+rbox[1], 224., 0, "ONLY");
- gMC->Gspos("BRS4", 2, "ALIC", 430.+3., -180.+55./2.+rbox[1], -224., 0, "ONLY");
- gMC->Gspos("BRS4", 3, "ALIC", -430.-3., -180.+55./2.+rbox[1], 224., 0, "ONLY");
- gMC->Gspos("BRS4", 4, "ALIC", -430.-3., -180.+55./2.+rbox[1], -224., 0, "ONLY");
+ gMC->Gspos("BRS4", 1, "ALIC", 430.+3., -190.+55./2.+rbox[1], 224., 0, "ONLY");
+ gMC->Gspos("BRS4", 2, "ALIC", 430.+3., -190.+55./2.+rbox[1], -224., 0, "ONLY");
+// gMC->Gspos("BRS4", 3, "ALIC", -430.+3, -180.+55./2.+rbox[1], 224., 0, "ONLY");
+// gMC->Gspos("BRS4", 4, "ALIC", -430.+3, -180.+55./2.+rbox[1], -224., 0, "ONLY");
stmin = -.8;
Int_t isxfld = gAlice->Field()->Integ();
Float_t sxmgmx = gAlice->Field()->Max();
-
+
+
Float_t asteel[4] = { 55.847,51.9961,58.6934,28.0855 };
Float_t zsteel[4] = { 26.,24.,28.,14. };
Float_t wsteel[4] = { .715,.18,.1,.005 };
+
+ //Air
+
+ Float_t aAir[4]={12.0107,14.0067,15.9994,39.948};
+ Float_t zAir[4]={6.,7.,8.,18.};
+ Float_t wAir[4]={0.000124,0.755267,0.231781,0.012827};
+ Float_t dAir = 1.20479E-3;
AliMixture(65, "STAINLESS STEEL$", asteel, zsteel, 7.88, 4, wsteel);
- AliMaterial(5, "AIR$ ", 14.61, 7.3, .001205, 30423.24, 67500.);
+ AliMixture(5, "AIR$ ", aAir, zAir, dAir,4, wAir);
+
AliMedium(65, "Stainless Steel", 65, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
- AliMedium( 5, "Air ", 5, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+ AliMedium( 5, "Air", 5, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+
}
//_____________________________________________________________________________
//
// Initialise the module after the geometry has been defined
//
- if(fDebug) {
+ if(AliLog::GetGlobalDebugLevel()>0) {
printf("%s: **************************************"
" FRAME "
"**************************************\n",ClassName());
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff