]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - STRUCT/AliHALL.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | inline | side by side
Missing protection added. (A.Gheata)
[u/mrichter/AliRoot.git] / STRUCT / AliHALL.cxx
diff --git a/STRUCT/AliHALL.cxx b/STRUCT/AliHALL.cxx
index cd6a76cd9cda1fbe76ad2044f01626339aac3f42..c8412a4819c6ea3860a9d007cac2a545b21068f9 100644 (file)
--- a/STRUCT/AliHALL.cxx
+++ b/STRUCT/AliHALL.cxx
@@ -36,12 +36,14 @@
 //                                                                           //
 ///////////////////////////////////////////////////////////////////////////////
 
+#include <TGeoGlobalMagField.h>
 #include <TVirtualMC.h>
 
 #include "AliConst.h"
 #include "AliHALL.h"
 #include "AliMagF.h"
 #include "AliRun.h"
+#include "AliLog.h"
  
 ClassImp(AliHALL)
  
@@ -60,9 +62,9 @@ AliHALL::AliHALL(const char *name, const char *title)
   //
   // Standard constructor for the experimental Hall
   //
-  SetMarkerColor(7);
-  SetMarkerStyle(2);
-  SetMarkerSize(0.4);
+  //PH  SetMarkerColor(7);
+  //PH  SetMarkerStyle(2);
+  //PH  SetMarkerSize(0.4);
 }
  
 //_____________________________________________________________________________
@@ -111,17 +113,17 @@ void AliHALL::CreateGeometry()
   h   = 140.;
   phi = TMath::ACos(h / r);
   xl  = r * TMath::Sin(phi);
-  dr  = 100.;
+  dr  = 1600.;
   dh  = dr * TMath::Cos(phi);
   dl  = dr * TMath::Sin(phi);
   if (gAlice->GetModule("ZDC") == 0) {
     
     //     No ZDC 
-    hullen = 250.;
+    hullen = 370.;
   } else {
     
     //     ZDC is present 
-    hullen = 6400.;
+    hullen = 6520.;
   }
   trdpar[0] = xl + dl;
   trdpar[1] = xl;
@@ -130,7 +132,7 @@ void AliHALL::CreateGeometry()
   AliMatrix(idrotm[1900], 90., 0., 0., 0., 90., 90.);
   AliMatrix(idrotm[1901], 270., 0., 90., 90., 0., 0.);
   gMC->Gsvolu("HUFL", "TRD1", idtmed[1956], trdpar, 4);
-  r2 = hullen + 2020.;
+  r2 = hullen + 1900.;
   gMC->Gspos("HUFL", 1, "ALIC", 70.,-100-trdpar[3] , -r2, idrotm[1900], "ONLY");
   
   //     RB24/26 wall 
@@ -142,22 +144,10 @@ void AliHALL::CreateGeometry()
   tspar[3] = phid - 90.;
   tspar[4] = 270. - phid;
   gMC->Gsvolu("HUWA", "TUBS", idtmed[1956], tspar, 5);
-  gMC->Gspos("HUWA", 1, "ALIC", 70., 40., -2020 - hullen , 0, "ONLY");
+  gMC->Gspos("HUWA", 1, "ALIC", 70., 40., -1900 - hullen , 0, "ONLY");
   
-  //     END WALL 
-  gMC->Gsvolu("HEW1", "BOX ", idtmed[1956], pbox, 0);
-  pbox[0] = 600.;
-  pbox[1] = 418.;
-  pbox[2] = 60.;
-  gMC->Gsposp("HEW1", 1, "ALIC", 0., -pbox[1]-60., -1960, 0, "ONLY", pbox, 3);
-  pbox[1] = 822.;
-  gMC->Gsposp("HEW1", 2, "ALIC", 0.,  pbox[1]+60., -1960, 0, "ONLY", pbox, 3);
-  pbox[0] = 270.;
-  pbox[1] =  60.;
-  gMC->Gsposp("HEW1", 3, "ALIC",  pbox[0]+60.,  0. , -1960, 0, "ONLY", pbox, 3);
-  gMC->Gsposp("HEW1", 4, "ALIC", -pbox[0]-60.,  0. , -1960, 0, "ONLY", pbox, 3);
 
-  //     hall floor 
+  //     Hall floor 
   
   phid      = 16.197;
   trdpar[0] = 700.;
@@ -168,7 +158,7 @@ void AliHALL::CreateGeometry()
   gMC->Gspos("HHF1", 1, "ALIC", 0., -801., 1350., idrotm[1900], "ONLY");
   gMC->Gspos("HHF1", 2, "ALIC", 0., -801.,-1350., idrotm[1900], "ONLY");
   
-  //     hall side walls 
+  //     Hall side walls 
   
   trapar[0] = 550.;
   trapar[1] = 0.;
@@ -281,28 +271,39 @@ void AliHALL::CreateMaterials()
   //
   // Create materials for the experimental hall
   //
-  
-
-  Int_t   isxfld = gAlice->Field()->Integ();
-  Float_t sxmgmx = gAlice->Field()->Max();
+    
+  Int_t   isxfld = 0;
+  Float_t sxmgmx = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Max();
   
   Float_t aconc[10] = { 1.,12.01,15.994,22.99,24.305,26.98,28.086,39.1,40.08,55.85 };
   Float_t zconc[10] = { 1.,6.,8.,11.,12.,13.,14.,19.,20.,26. };
   Float_t wconc[10] = { .01,.001,.529107,.016,.002,.033872,.337021,.013,.044,.014 };
   
+  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;
+
+// Steel  
+  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 };
+
   Float_t epsil, stmin, deemax, tmaxfd, stemax;
 
-  //     FOR CONCRETE 
   
-  AliMaterial(10, "IRON$     ", 55.85, 26., 7.87, 1.76, 17.1);
-  AliMaterial(30, "IRON$     ", 55.85, 26., 7.87, 1.76, 17.1);
-  AliMaterial(50, "IRON$     ", 55.85, 26., 7.87, 1.76, 17.1);
-  AliMaterial(15, "AIR$      ", 14.61, 7.3, .001205, 30423.24, 67500);
-  AliMaterial(35, "AIR$      ", 14.61, 7.3, .001205, 30423.24, 67500);
-  AliMaterial(55, "AIR$      ", 14.61, 7.3, .001205, 30423.24, 67500);
-  AliMixture(17, "CONCRETE$", aconc, zconc, 2.35, 10, wconc);
-  AliMixture(37, "CONCRETE$", aconc, zconc, 2.35, 10, wconc);
-  AliMixture(57, "CONCRETE$", aconc, zconc, 2.35, 10, wconc);
+  AliMixture(15, "AIR0$      ", aAir, zAir, dAir, 4, wAir);
+  AliMixture(35, "AIR1$      ", aAir, zAir, dAir, 4, wAir);
+  AliMixture(55, "AIR2$      ", aAir, zAir, dAir, 4, wAir);
+
+  AliMixture(17, "CONCRETE0$", aconc, zconc, 2.35, 10, wconc);
+  AliMixture(37, "CONCRETE1$", aconc, zconc, 2.35, 10, wconc);
+  AliMixture(57, "CONCRETE2$", aconc, zconc, 2.35, 10, wconc);
+
+  //     Stainless Steel 
+  AliMixture(10, "STAINLESS STEEL1", asteel, zsteel, 7.88, 4, wsteel);
+  AliMixture(30, "STAINLESS STEEL2", asteel, zsteel, 7.88, 4, wsteel);
+  AliMixture(50, "STAINLESS STEEL3", asteel, zsteel, 7.88, 4, wsteel);
   
   // **************** 
   //     Defines tracking media parameters. 
@@ -315,11 +316,10 @@ void AliHALL::CreateMaterials()
   stmin  = -.8;
   // *************** 
   
-  //     IRON 
-  
-  AliMedium(10, "FE_C0             ", 10, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(30, "FE_C1             ", 30, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(50, "FE_C2             ", 50, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  //  Stainless Steel
+  AliMedium(10, "STST_C0             ", 10, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(30, "STST_C1             ", 30, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(50, "STST_C2             ", 50, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
   
   //    Air 
   
@@ -342,7 +342,7 @@ void AliHALL::Init()
   //
   Int_t i;
   //
-  if(fDebug) {
+  if(AliLog::GetGlobalDebugLevel()>0) {
     printf("\n%s: ",ClassName());
     for(i=0;i<35;i++) printf("*");
     printf(" HALL_INIT ");