]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - STRUCT/AliHALL.cxx
added slewing correction by data
[u/mrichter/AliRoot.git] / STRUCT / AliHALL.cxx
index bc1cdbb06ade19b0fb2e6a3ffc1a4f78aca51616..3b7a229479ef5b41fe1ffe30c768c2f83291ce27 100644 (file)
@@ -36,7 +36,9 @@
 //                                                                           //
 ///////////////////////////////////////////////////////////////////////////////
 
+#include <TGeoGlobalMagField.h>
 #include <TVirtualMC.h>
+#include <TArrayI.h>
 
 #include "AliConst.h"
 #include "AliHALL.h"
@@ -61,9 +63,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);
 }
  
 //_____________________________________________________________________________
@@ -145,20 +147,8 @@ void AliHALL::CreateGeometry()
   gMC->Gsvolu("HUWA", "TUBS", idtmed[1956], tspar, 5);
   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.;
@@ -169,7 +159,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.;
@@ -284,7 +274,7 @@ void AliHALL::CreateMaterials()
   //
     
   Int_t   isxfld = 0;
-  Float_t sxmgmx = gAlice->Field()->Max();
+  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. };
@@ -295,21 +285,30 @@ void AliHALL::CreateMaterials()
   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); 
+  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(15, "AIR$      ", aAir, zAir, dAir, 4, wAir);
-  AliMixture(35, "AIR$      ", aAir, zAir, dAir, 4, wAir);
-  AliMixture(55, "AIR$      ", 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);
 
-  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);
+  //     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);
+  // Iron
+  //
+  AliMaterial(52, "IRON", 55.85, 26., 7.87, 1.76, 17.1);
+  
   
   // **************** 
   //     Defines tracking media parameters. 
@@ -322,11 +321,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 
   
@@ -339,6 +337,9 @@ void AliHALL::CreateMaterials()
   AliMedium(17, "CC_C0            ", 17, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
   AliMedium(37, "CC_C1            ", 37, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
   AliMedium(57, "CC_C2            ", 57, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+
+  // Iron
+  AliMedium(52, "FE_C2            ", 52, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 }
 
 //_____________________________________________________________________________
@@ -362,58 +363,3 @@ void AliHALL::Init()
   }
 }
 
-//_____________________________________________________________________________
-void AliHALL::DrawModule() const
-{
-  //
-  // Draw a shaded view of Experimental Hall
-  //
-
-  // Set everything unseen
-  gMC->Gsatt("*", "seen", -1);
-  // 
-  // Set ALIC mother transparent
-  gMC->Gsatt("ALIC","SEEN",0);
-  //
-  // Set the volumes visible
-  gMC->Gsatt("HUFL","seen",1);
-  gMC->Gsatt("HUWA","seen",1);
-  gMC->Gsatt("HUP2","seen",1);
-  gMC->Gsatt("HEW1","seen",1);
-  gMC->Gsatt("HHF1","seen",1);
-  gMC->Gsatt("HHW1","seen",1);
-  gMC->Gsatt("HBW1","seen",1);
-  gMC->Gsatt("HHW2","seen",1);
-  gMC->Gsatt("HHW3","seen",1);
-  gMC->Gsatt("HHC1","seen",1);
-  gMC->Gsatt("HHF2","seen",1);
-  gMC->Gsatt("HPIL","seen",1);
-  gMC->Gsatt("HMBS","seen",1);
-  gMC->Gsatt("HBBS","seen",1);
-  gMC->Gsatt("HPBS","seen",1);
-  gMC->Gsatt("HXFI","seen",1);
-  gMC->Gsatt("HXII","seen",1);
-  //
-  gMC->Gdopt("hide", "on");
-  gMC->Gdopt("shad", "on");
-  gMC->Gsatt("*", "fill", 7);
-  gMC->SetClipBox(".");
-  if (gAlice->GetModule("ZDC") == 0) {
-    //
-    // ZDC is not present
-    //
-    gMC->SetClipBox("*", 0, 3000, -3000, 3000, -6000, 6000);
-    gMC->DefaultRange();
-    gMC->Gdraw("alic", 40, 30, 0, 12, 7.5, .005, .005);
-  } else {
-    //
-    // ZDC is present
-    //
-    gMC->SetClipBox("*", 0, 2000, -2000, 2000, -2000, 16000);
-    gMC->DefaultRange();
-    gMC->Gdraw("alic", 40, 30, 0, 17.5, 10, .0019, .0019);
-  }
-  gMC->Gdhead(1111, "Experimental Hall");
-  gMC->Gdman(18, 2, "MAN");
-}