New beamshield geometry with increased tolerances and insulation thickness.
[u/mrichter/AliRoot.git] / STRUCT / AliHALL.cxx
index 72c9d61..79d584d 100644 (file)
@@ -1,3 +1,43 @@
+/**************************************************************************
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ *                                                                        *
+ * Author: The ALICE Off-line Project.                                    *
+ * Contributors are mentioned in the code where appropriate.              *
+ *                                                                        *
+ * Permission to use, copy, modify and distribute this software and its   *
+ * documentation strictly for non-commercial purposes is hereby granted   *
+ * without fee, provided that the above copyright notice appears in all   *
+ * copies and that both the copyright notice and this permission notice   *
+ * appear in the supporting documentation. The authors make no claims     *
+ * about the suitability of this software for any purpose. It is          *
+ * provided "as is" without express or implied warranty.                  *
+ **************************************************************************/
+
+/*
+$Log$
+Revision 1.12  2001/06/20 16:08:56  morsch
+Remove some shielding to accomodate compensator magnet.
+
+Revision 1.11  2001/05/16 14:57:22  alibrary
+New files for folders and Stack
+
+Revision 1.10  2000/10/02 21:28:15  fca
+Removal of useless dependecies via forward declarations
+
+Revision 1.9  2000/06/13 14:55:26  morsch
+Unused variables removed
+
+Revision 1.8  2000/06/11 12:35:41  morsch
+Coding rule violations corrected
+
+Revision 1.7  2000/01/12 15:33:28  morsch
+Muon filter removed
+
+Revision 1.6  1999/09/29 09:24:30  fca
+Introduction of the Copyright and cvs Log
+
+*/
+
 ///////////////////////////////////////////////////////////////////////////////
 //                                                                           //
 //  Experimental Hall                                                        //
@@ -21,6 +61,8 @@
 
 #include "AliHALL.h"
 #include "AliRun.h"
+#include "AliMC.h"
+#include "AliMagF.h"
 #include "AliConst.h"
  
 ClassImp(AliHALL)
@@ -77,59 +119,14 @@ void AliHALL::CreateGeometry()
 
   
   Float_t r2;
-  Float_t phid, phim, tpar[3], pbox[3], zfil_out, h, r, tspar[5];
+  Float_t phid, phim, tpar[3], pbox[3], h, r, tspar[5];
   Float_t w1, dh, am, bm, dl,cm, hm, dr, dx, xl;
   Int_t idrotm[1999];
   Float_t trdpar[4], trapar[11], hullen;
-  Float_t dz, phi, par[3], zfil_in;
+  Float_t phi;
   
   Int_t *idtmed = fIdtmed->GetArray()-1899;
   
-  //abs_d   = 90.;  // DEFINES DRIFT LENGTH 
-  //z_nose  = 102.;
-  //z_cone  = 285.;
-  //theta1  = 24.;  // 1. angle defining the front absorber 
-  //theta2  = 5.;   // 2. angle defining the front absorbe 
-  //acc_max = 9.;   // ANGLE POLAIRE MAXIMUM 
-  //acc_min = 2.;   // ANGLE POLAIRE MINIMUM DE DETECTION 
-  //abs_l   = 503.;
-  //d_steel = 1.;   // THICKNESS OF STEEL SUPPORT 
-  //d_poly  = 7.5;
-  //d_pb    = 2.5;
-  //abs_cc  = 315.; // DEFINES LENGTH OF CARBON 
-  //abs_c   = 358.;
-  //abs_s   = 150.; // DEFINES W-SHIELD LENGTH 
-  //abs_n   = 80.;  // START OF NOSE 
-  //r_abs   = 4.;
-  //r_pb    = .1;
-  //epsilon = .01;
-  //theta_r = 3.;
-  //d_rear  = 35.;
-  //theta_open = .75;
-  
-  //z_l3     = 700.;
-  //zmag_in  = 725.;
-  //zmag_out = 1225.;
-  zfil_in  = 1471.;
-  zfil_out = 1591.;
-  //zcon_in  = 1900.;
-  //zcon_out = 2e3;
-  //zcone_e  = 859.0875;
-  //spec_l   = 1800.;
-  //zplug_in = 1780.;
-  //zplug_out= 1900.;
-  
-  //     Chamber position 
-  //      CZ1=515.5 
-  //cz1 = 511.;
-  //cz2 = 686.;
-  //cz3 = 971.;
-  //cz4 = 1245.;
-  //cz5 = 1445.;
-  //cz6 = 1610.;
-  //cz7 = 1710.;
-  
-  
   //     RB24/26 TUNNEL FLOOR 
   
   r   = 220.;
@@ -285,32 +282,23 @@ void AliHALL::CreateGeometry()
   gMC->Gspos("HPIL", 1, "ALIC", 165.,-706+pbox[1] , 1350., 0, "ONLY");
   gMC->Gspos("HPIL", 2, "ALIC",-165.,-706+pbox[1] , 1350., 0, "ONLY");
   
-  //     concrete beam shield 
+  //     simple concrete beam shield 
+  
+  Float_t ppgon[10];
+  ppgon[0] =    45.;
+  ppgon[1] =   360.;
+  ppgon[2] =     4.;
+  ppgon[3] =     2.;
+  ppgon[4] =  -720.;
+  ppgon[5] =   150.;
+  ppgon[6] =   250.;
+  ppgon[7] = -1800.;
+  ppgon[8] =   150.;
+  ppgon[9] =   250.;
+
+  gMC->Gsvolu("HMBS", "PGON", idtmed[1956], ppgon, 10);
+  gMC->Gspos("HMBS", 1, "ALIC", 0., 70., 0., 0, "ONLY");
   
-  pbox[0] = 402.5;
-  pbox[1] = 260.;
-  pbox[2] = 120.;
-  gMC->Gsvolu("HMBS", "BOX ", idtmed[1956], pbox, 3);
-  pbox[0] = 85.;
-  pbox[1] = 100.;
-  gMC->Gsvolu("HBBS", "BOX ", idtmed[1956], pbox, 3);
-  gMC->Gspos("HBBS", 1, "HMBS", -157.5, 0., 0., 0, "ONLY");
-  pbox[0] = 40.;
-  pbox[1] = 130.;
-  gMC->Gsvolu("HPBS", "BOX ", idtmed[1956], pbox, 3);
-  gMC->Gspos("HPBS", 1, "HMBS", 202.5,  30.,    0., 0, "ONLY");
-  gMC->Gspos("HMBS", 1, "ALIC", 157.5, -50., -820., 0, "ONLY");
-  
-  //       MUON FILTER 
-  par[0] = 30.;
-  par[1] = 310.;
-  par[2] = (zfil_out - zfil_in) / 2.;
-  gMC->Gsvolu("HXFI", "TUBE", idtmed[1949], par, 3);
-  dz = (zfil_in + zfil_out) / 2.;
-  par[2] -= 10.;
-  gMC->Gsvolu("HXII", "TUBE", idtmed[1909], par, 3);
-  gMC->Gspos("HXII", 1, "HXFI", 0., 0., 0., 0, "ONLY");
-  gMC->Gspos("HXFI", 1, "ALIC", 0., 0., dz, 0, "ONLY");
 }
 
 //_____________________________________________________________________________
@@ -321,8 +309,8 @@ void AliHALL::CreateMaterials()
   //
   
 
-  Int_t   ISXFLD = gAlice->Field()->Integ();
-  Float_t SXMGMX = gAlice->Field()->Max();
+  Int_t   isxfld = gAlice->Field()->Integ();
+  Float_t sxmgmx = gAlice->Field()->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. };
@@ -355,21 +343,21 @@ void AliHALL::CreateMaterials()
   
   //     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);
+  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);
   
   //    Air 
   
-  AliMedium(15, "AIR_C0           ", 15, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(35, "AIR_C1           ", 35, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(55, "AIR_C2           ", 55, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(15, "AIR_C0           ", 15, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(35, "AIR_C1           ", 35, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(55, "AIR_C2           ", 55, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
   
   //    Concrete 
   
-  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);
+  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);
 }
 
 //_____________________________________________________________________________
@@ -380,15 +368,17 @@ void AliHALL::Init()
   //
   Int_t i;
   //
-  printf("\n");
-  for(i=0;i<35;i++) printf("*");
-  printf(" HALL_INIT ");
-  for(i=0;i<35;i++) printf("*");
-  printf("\n");
-  //
-  // Here the HALL initialisation code (if any!)
-  for(i=0;i<80;i++) printf("*");
-  printf("\n");
+  if(fDebug) {
+    printf("\n%s: ",ClassName());
+    for(i=0;i<35;i++) printf("*");
+    printf(" HALL_INIT ");
+    for(i=0;i<35;i++) printf("*");
+    printf("\n%s: ",ClassName());
+    //
+    // Here the HALL initialisation code (if any!)
+    for(i=0;i<80;i++) printf("*");
+    printf("\n");
+  }
 }
 
 //_____________________________________________________________________________