Base class only
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 12 Jan 2000 15:37:57 +0000 (15:37 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 12 Jan 2000 15:37:57 +0000 (15:37 +0000)
STRUCT/AliABSO.cxx
STRUCT/AliSHIL.cxx

index 4de512d0fc3ea085695aa3192edf257ebdeab36c..5c4c337b5861c356b033c4ae07eb320bffac4f21 100644 (file)
@@ -15,6 +15,9 @@
 
 /*
 $Log$
+Revision 1.6  1999/09/29 09:24:29  fca
+Introduction of the Copyright and cvs Log
+
 */
 
 ///////////////////////////////////////////////////////////////////////////////
@@ -70,340 +73,6 @@ void AliABSO::CreateGeometry()
   //
   // Creation of the geometry of the muon absorber
   //
-  //Begin_Html
-  /*
-    <img src="picts/AliABSOTree.gif">
-  */
-  //End_Html
-  //Begin_Html
-  /*
-    <img src="picts/AliABSO.gif">
-  */
-  //End_Html
-
-  Int_t *idtmed = fIdtmed->GetArray()-1599;
-  
-  Float_t d_pb, cpar[5], dpar[12], tpar[3], zpos,
-    cpar1[5], cpar2[5], cpar3[5], cpar4[5], cpar5[12], 
-    cpar7[5], cpar8[5], cpar9[5], abs_c, abs_d,
-    abs_l, cpar10[5], r_abs;
-  Float_t theta1, theta2, abs_cc, d_rear, dz, zr,
-    z_cone, d_poly, z_nose, theta_open;
-  Float_t acc_min, acc_max, par[50], d_steel, z_w,
-    theta_r, epsilon;
-  //
-  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.;
-  //
-  // --- Outer shape of front absorber 
-  par[0] = 0.;
-  par[1] = 360.;
-  par[2] = 4.;
-  par[3] = abs_d;
-  par[4] = 0.;
-  par[5] = abs_d * TMath::Tan(theta1 * kDegrad);
-  par[6] = z_nose;
-  par[7] = 0.;
-  par[8] = par[6] * TMath::Tan(theta1 * kDegrad);
-  par[9] = z_cone;
-  par[10] = 0.;
-  par[11] = par[8] + (par[9] - par[6]) * TMath::Tan(theta2 * kDegrad);
-  par[12] = abs_l;
-  par[13] = 0.;
-  par[14] = par[11] + (par[12] - par[9]) * TMath::Tan(acc_max * kDegrad);
-  gMC->Gsvolu("ABSM", "PCON", idtmed[1605], par, 15);
-  //
-  // --- Now define all elements of the absorber 
-  //
-  //       TUNGSTEN NOSE SEGMENT BETWEEN Z=90 AND 112 CM 
-  //       SHAPED ALONG A 24 DEG LINE 
-  //
-  cpar1[0] = (z_nose - abs_d) / 2.;
-  cpar1[1] = abs_d * TMath::Tan(acc_max * kDegrad) + d_steel;
-  cpar1[2] = abs_d * TMath::Tan(theta1 * kDegrad);
-  cpar1[3] = z_nose * TMath::Tan(acc_max * kDegrad) + d_steel;
-  cpar1[4] = z_nose * TMath::Tan(theta1 * kDegrad);
-  gMC->Gsvolu("ANOS", "CONE", idtmed[1611], cpar1, 5);
-  //
-  dz = cpar1[0] + abs_d;
-  gMC->Gspos("ANOS", 1, "ABSM", 0., 0., dz, 0, "ONLY");
-  //
-  //       IRON  SUPPORT STRUCTURE 
-  //
-  cpar2[0] = (abs_l - abs_d) / 2.;
-  cpar2[1] = abs_d * TMath::Tan(acc_max * kDegrad);
-  cpar2[2] = cpar2[1] + d_steel;
-  cpar2[3] = abs_l * TMath::Tan(acc_max * kDegrad);
-  cpar2[4] = cpar2[3] + d_steel;
-  gMC->Gsvolu("ASST", "CONE", idtmed[1658], cpar2, 5);
-  dz = cpar2[0] + abs_d;
-  gMC->Gspos("ASST", 1, "ABSM", 0., 0., dz, 0, "ONLY");
-  //
-  //       PB FRONT SHIELD INNER SEGMENT, ALSO POLYETHYLENE WAS 
-  //       CONSIDERED FOR THIS REGION 
-  //
-  cpar3[0] = (z_cone - z_nose) / 2.;
-  cpar3[1] = cpar1[3];
-  cpar3[2] = cpar1[3] + d_poly;
-  cpar3[3] = z_cone * TMath::Tan(acc_max * kDegrad) + d_steel;
-  cpar3[4] = cpar3[3] + d_poly;
-  gMC->Gsvolu("AWFS", "CONE", idtmed[1652], cpar3, 5);
-  dz = cpar3[0] + z_nose;
-  gMC->Gspos("AWFS", 1, "ABSM", 0., 0., dz, 0, "ONLY");
-  //
-  //       PB OUTER SURFACE 
-  //
-  cpar5[0] = 0.;
-  cpar5[1] = 360.;
-  cpar5[2] = 3.;
-  cpar5[3] = z_nose;
-  cpar5[4] = z_nose * TMath::Tan(acc_max * kDegrad) + d_steel + d_poly;
-  cpar5[5] = z_nose * TMath::Tan(theta1 * kDegrad);
-  cpar5[6] = z_cone;
-  cpar5[7] = z_cone * TMath::Tan(acc_max * kDegrad) + d_steel + d_poly;
-  cpar5[8] = cpar5[7] + d_pb;
-  cpar5[9] = abs_l;
-  cpar5[10] = abs_l * TMath::Tan(acc_max * kDegrad) + d_steel + d_poly;
-  cpar5[11] = cpar5[10] + d_pb;
-  gMC->Gsvolu("APBS", "PCON", idtmed[1612], cpar5, 12);
-  dz = 0.;
-  gMC->Gspos("APBS", 1, "ABSM", 0., 0., dz, 0, "ONLY");
-  //
-  //     POLYETHYLEN LAYER 
-  //
-  cpar4[0] = (abs_l - z_cone) / 2.;
-  cpar4[1] = z_cone * TMath::Tan(acc_max * kDegrad) + d_steel;
-  cpar4[2] = cpar4[1] + d_poly;
-  cpar4[3] = abs_l * TMath::Tan(acc_max * kDegrad) + d_steel;
-  cpar4[4] = cpar4[3] + d_poly;
-  gMC->Gsvolu("APOL", "CONE", idtmed[1657], cpar4, 5);
-  dz = cpar4[0] + z_cone;
-  gMC->Gspos("APOL", 1, "ABSM", 0., 0., dz, 0, "ONLY");
-  //
-  //     LEAD INNER SHIELD (inner radius const up to z=abs_c) 
-  //
-  z_w      = r_abs / TMath::Tan(acc_min * kDegrad);
-  cpar8[0] = (abs_c - z_w) / 2.;
-  cpar8[1] = r_abs;
-  cpar8[2] = r_abs + epsilon;
-  cpar8[3] = r_abs;
-  cpar8[4] = abs_c * TMath::Tan(acc_min * kDegrad);
-  gMC->Gsvolu("AWI1", "CONE", idtmed[1652], cpar8, 5);
-  dz = cpar8[0] + z_w;
-  gMC->Gspos("AWI1", 1, "ABSM", 0., 0., dz, 0, "ONLY");
-  //
-  //     TUNGSTEN OPENING CONE UP TO THE END 
-  //
-  cpar8[0] = (abs_l - abs_c) / 2.;
-  cpar8[1] = r_abs;
-  cpar8[2] = abs_c * TMath::Tan(acc_min * kDegrad);
-  cpar8[3] = cpar8[1] + cpar8[0] * 2. * TMath::Tan(theta_open * kDegrad);
-  cpar8[4] = abs_l * TMath::Tan(acc_min * kDegrad);
-  gMC->Gsvolu("AWI2", "CONE", idtmed[1651], cpar8, 5);
-  dz = cpar8[0] + abs_c;
-  gMC->Gspos("AWI2", 1, "ABSM", 0., 0., dz, 0, "ONLY");
-  //
-  //     CONCRETE CONE 
-  //
-  cpar7[0] = (abs_l - d_rear - abs_cc) / 2.;
-  cpar7[1] = abs_cc * TMath::Tan(acc_min * kDegrad);
-  cpar7[2] = abs_cc * TMath::Tan(acc_max * kDegrad);
-  cpar7[3] = (abs_l - d_rear) * TMath::Tan(acc_min * kDegrad);
-  cpar7[4] = (abs_l - d_rear) * TMath::Tan(acc_max * kDegrad);
-  gMC->Gsvolu("ACON", "CONE", idtmed[1656], cpar7, 5);
-  dz = cpar7[0] + abs_cc;
-  gMC->Gspos("ACON", 1, "ABSM", 0., 0., dz, 0, "ONLY");
-  //
-  //     REAR SHIELD 
-  //
-  zr = abs_l - d_rear;
-  cpar9[0] = 2.5;
-  cpar9[1] = zr * TMath::Tan(theta_r * kDegrad);
-  cpar9[2] = zr * TMath::Tan(acc_max * kDegrad);
-  cpar9[3] = cpar9[1] + TMath::Tan(theta_r * kDegrad) * 5.;
-  cpar9[4] = cpar9[2] + TMath::Tan(acc_max * kDegrad) * 5.;
-  gMC->Gsvolu("ARE1", "CONE", idtmed[1652], cpar9, 5);
-  dz  = cpar9[0] + zr;
-  zr += 5.;
-  gMC->Gspos("ARE1", 1, "ABSM", 0., 0., dz, 0, "ONLY");
-  //
-  cpar9[1] = zr * TMath::Tan(theta_r * kDegrad);
-  cpar9[2] = zr * TMath::Tan(acc_max * kDegrad);
-  cpar9[3] = cpar9[1] + TMath::Tan(theta_r * kDegrad) * 5.;
-  cpar9[4] = cpar9[2] + TMath::Tan(acc_max * kDegrad) * 5.;
-  gMC->Gsvolu("ARE2", "CONE", idtmed[1657], cpar9, 5);
-  dz  = cpar9[0] + zr;
-  zr += 5.;
-  gMC->Gspos("ARE2", 1, "ABSM", 0., 0., dz, 0, "ONLY");
-  //
-  cpar9[1] = zr * TMath::Tan(theta_r * kDegrad);
-  cpar9[2] = zr * TMath::Tan(acc_max * kDegrad);
-  cpar9[3] = cpar9[1] + TMath::Tan(theta_r * kDegrad) * 5.;
-  cpar9[4] = cpar9[2] + TMath::Tan(acc_max * kDegrad) * 5.;
-  gMC->Gsvolu("ARE3", "CONE", idtmed[1652], cpar9, 5);
-  dz  = cpar9[0] + zr;
-  zr += 5.;
-  gMC->Gspos("ARE3", 1, "ABSM", 0., 0., dz, 0, "ONLY");
-  //
-  cpar9[1] = zr * TMath::Tan(theta_r * kDegrad);
-  cpar9[2] = zr * TMath::Tan(acc_max * kDegrad);
-  cpar9[3] = cpar9[1] + TMath::Tan(theta_r * kDegrad) * 5.;
-  cpar9[4] = cpar9[2] + TMath::Tan(acc_max * kDegrad) * 5.;
-  gMC->Gsvolu("ARE4", "CONE", idtmed[1657], cpar9, 5);
-  dz  = cpar9[0] + zr;
-  zr += 5.;
-  gMC->Gspos("ARE4", 1, "ABSM", 0., 0., dz, 0, "ONLY");
-  //
-  cpar9[1] = zr * TMath::Tan(theta_r * kDegrad);
-  cpar9[2] = zr * TMath::Tan(acc_max * kDegrad);
-  cpar9[3] = cpar9[1] + TMath::Tan(theta_r * kDegrad) * 5.;
-  cpar9[4] = cpar9[2] + TMath::Tan(acc_max * kDegrad) * 5.;
-  gMC->Gsvolu("ARE5", "CONE", idtmed[1652], cpar9, 5);
-  dz  = cpar9[0] + zr;
-  zr += 5.;
-  gMC->Gspos("ARE5", 1, "ABSM", 0., 0., dz, 0, "ONLY");
-  //
-  cpar9[1] = zr * TMath::Tan(theta_r * kDegrad);
-  cpar9[2] = zr * TMath::Tan(acc_max * kDegrad);
-  cpar9[3] = cpar9[1] + TMath::Tan(theta_r * kDegrad) * 5.;
-  cpar9[4] = cpar9[2] + TMath::Tan(acc_max * kDegrad) * 5.;
-  gMC->Gsvolu("ARE6", "CONE", idtmed[1657], cpar9, 5);
-  dz  = cpar9[0] + zr;
-  zr += 5.;
-  gMC->Gspos("ARE6", 1, "ABSM", 0., 0., dz, 0, "ONLY");
-  //
-  cpar9[1] = zr * TMath::Tan(theta_r * kDegrad);
-  cpar9[2] = zr * TMath::Tan(acc_max * kDegrad);
-  cpar9[3] = cpar9[1] + TMath::Tan(theta_r * kDegrad) * 5.;
-  cpar9[4] = cpar9[2] + TMath::Tan(acc_max * kDegrad) * 5.;
-  gMC->Gsvolu("ARE7", "CONE", idtmed[1612], cpar9, 5);
-  dz  = cpar9[0] + zr;
-  zr += 5.;
-  gMC->Gspos("ARE7", 1, "ABSM", 0., 0., dz, 0, "ONLY");
-  //
-  //     TUNGSTEN REAR SHIELD INNER PART 
-  //
-  zr = abs_l - d_rear;
-  cpar10[0] = d_rear / 2.;
-  cpar10[1] = zr * TMath::Tan(acc_min * kDegrad);
-  cpar10[2] = zr * TMath::Tan(theta_r * kDegrad);
-  cpar10[3] = cpar10[1] + d_rear * TMath::Tan(acc_min * kDegrad);
-  cpar10[4] = cpar10[2] + d_rear * TMath::Tan(theta_r * kDegrad);
-  gMC->Gsvolu("ARIN", "CONE", idtmed[1611], cpar10, 5);
-  dz = cpar10[0] + zr;
-  gMC->Gspos("ARIN", 1, "ABSM", 0., 0., dz, 0, "ONLY");
-  //
-  //     ELEMENTS OF THE BEAM PIPE TO BE POSITIONED INTO THE ABSORBER 
-  //
-  //     MOTHER VOLUME 1. SEGMENT 
-  //
-  tpar[0] = 0.;
-  tpar[1] = r_abs;
-  tpar[2] = (abs_c - abs_d) / 2.;
-  gMC->Gsvolu("AATU", "TUBE", idtmed[1655], tpar, 3);
-  //
-  tpar[1] = r_abs - .8;
-  tpar[0] = tpar[1] - .2;
-  tpar[2] = (abs_c - abs_d) / 2.;
-  gMC->Gsvolu("ATUB", "TUBE", idtmed[1649], tpar, 3);
-  dz = 0.;
-  gMC->Gspos("ATUB", 1, "AATU", 0., 0., dz, 0, "ONLY");
-  //
-  dz = (abs_c - abs_d) / 2. + abs_d;
-  gMC->Gspos("AATU", 1, "ABSM", 0., 0., dz, 0, "ONLY");
-  //
-  //     MOTHER VOLUME 2. SEGMENT 
-  //
-  cpar[0] = (abs_l - abs_c) / 2.;
-  cpar[1] = 0.;
-  cpar[2] = r_abs;
-  cpar[3] = 0.;
-  cpar[4] = cpar[2] + cpar[0] * 2. * TMath::Tan(theta_open * kDegrad);
-  gMC->Gsvolu("AAT1", "CONE", idtmed[1655], cpar, 5);
-  //
-  cpar[0]  = (abs_l - abs_c) / 2.;
-  cpar[2] += -.8;
-  cpar[1]  = cpar[2] - .2;
-  cpar[4] += -.8;
-  cpar[3]  = cpar[4] - .2;
-  gMC->Gsvolu("ATU1", "CONE", idtmed[1649], cpar, 5);
-  dz = 0.;
-  gMC->Gspos("ATU1", 1, "AAT1", 0., 0., dz, 0, "ONLY");
-  //
-  dz = (abs_l - abs_c) / 2. + abs_c;
-  gMC->Gspos("AAT1", 1, "ABSM", 0., 0., dz, 0, "ONLY");
-  //
-  gMC->Gspos("ABSM", 1, "ALIC", 0., 0., 0., 0, "ONLY");
-  //
-  //       absorber support structure 
-  //
-  //       attention this element is positioned into ALIC 
-  //
-  dpar[0] = 0.;
-  dpar[1] = 360.;
-  dpar[2] = 3.;
-  dpar[3] = abs_l;
-  dpar[4] = abs_l * TMath::Tan(acc_max * kDegrad);
-  dpar[5] = dpar[4] + 4. / TMath::Cos(acc_max * kDegrad);
-  dpar[6] = 600.;
-  dpar[7] = TMath::Tan(acc_max * kDegrad) * 600;
-  dpar[8] = dpar[7] + 4. / TMath::Cos(acc_max * kDegrad);
-  dpar[9] = 670.;
-  dpar[10] = 159.;
-  dpar[11] = 163.5;
-  gMC->Gsvolu("ASUP", "PCON", idtmed[1618], dpar, 12);
-  dz = 0.;
-  gMC->Gspos("ASUP", 1, "ALIC", 0., 0., dz, 0, "ONLY");
-  //
-  //     Flange at the entrance of the absorber 
-  //
-  tpar[0] = 3.;
-  tpar[1] = 5.7;
-  tpar[2] = 2.;
-  gMC->Gsvolu("AF63", "TUBE", idtmed[1618], tpar, 3);
-  zpos = abs_d + tpar[2];
-  gMC->Gspos("AF63", 1, "ABSM", 0., 0., zpos, 0, "ONLY");
 }
 
 //_____________________________________________________________________________
@@ -412,50 +81,9 @@ void AliABSO::DrawModule()
   //
   // Draw a shaded view of the muon absorber
   //
-
-  // Set everything unseen
-  gMC->Gsatt("*", "seen", -1);
-  // 
-  // Set ALIC mother transparent
-  gMC->Gsatt("ALIC","SEEN",0);
-  //
-  // Set the volumes visible
-  gMC->Gsatt("ABSM","seen",1);
-  gMC->Gsatt("ANOS","seen",1);
-  gMC->Gsatt("ASST","seen",1);
-  gMC->Gsatt("AWFS","seen",1);
-  gMC->Gsatt("APBS","seen",1);
-  gMC->Gsatt("APOL","seen",1);
-  gMC->Gsatt("AWI1","seen",1);
-  gMC->Gsatt("AWI2","seen",1);
-  gMC->Gsatt("ACON","seen",1);
-  gMC->Gsatt("ARE1","seen",1);
-  gMC->Gsatt("ARE2","seen",1);
-  gMC->Gsatt("ARE3","seen",1);
-  gMC->Gsatt("ARE4","seen",1);
-  gMC->Gsatt("ARE5","seen",1);
-  gMC->Gsatt("ARE6","seen",1);
-  gMC->Gsatt("ARE7","seen",1);
-  gMC->Gsatt("ARIN","seen",1);
-  gMC->Gsatt("AATU","seen",1);
-  gMC->Gsatt("ATUB","seen",1);
-  gMC->Gsatt("AAT1","seen",1);
-  gMC->Gsatt("ATU1","seen",1);
-  gMC->Gsatt("ASUP","seen",1);
-  gMC->Gsatt("AF63","seen",1);
-  //
-  gMC->Gdopt("hide", "on");
-  gMC->Gdopt("shad", "on");
-  gMC->Gsatt("*", "fill", 7);
-  gMC->SetClipBox(".");
-  gMC->SetClipBox("*", 0, 3000, -3000, 3000, -6000, 6000);
-  gMC->DefaultRange();
-  gMC->Gdraw("alic", 40, 30, 0, 21.5, 15, .04, .04);
-  gMC->Gdhead(1111, "Muon Absorber");
-  gMC->Gdman(16, 6, "MAN");
-  gMC->Gdopt("hide","off");
 }
 
+
 //_____________________________________________________________________________
 void AliABSO::CreateMaterials()
 {
@@ -476,16 +104,27 @@ void AliABSO::CreateMaterials()
   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 };
-  
+//
+//
+// Ni-Cu-W alloy
+  Float_t aniwcu[3] ={58.6934, 183.84, 63.546};
+  Float_t zniwcu[3] ={28., 74., 29};
+  Float_t wniwcu[3] ={0.015,0.95,0.035};
+//
+// Insulation powder
+  Float_t ains[3] ={28.0855, 15.9994};
+  Float_t zins[3] ={14.,      8.};
+  Float_t wins[3] ={1.,       2.};
+//
   Float_t epsil, stmin, tmaxfd, deemax, stemax;
   //
   //     Carbon 
-  AliMaterial(6, "CARBON$   ", 12.01, 6., 2.265, 18.8, 49.9);
+  AliMaterial(6,  "CARBON$   ", 12.01, 6., 2.265, 18.8, 49.9);
   AliMaterial(26, "CARBON$   ", 12.01, 6., 2.265, 18.8, 49.9);
   AliMaterial(46, "CARBON$   ", 12.01, 6., 2.265, 18.8, 49.9);
   //
   //     Aluminum 
-  AliMaterial(9, "ALUMINIUM$", 26.98, 13., 2.7, 8.9, 37.2);
+  AliMaterial(9,  "ALUMINIUM$", 26.98, 13., 2.7, 8.9, 37.2);
   AliMaterial(29, "ALUMINIUM$", 26.98, 13., 2.7, 8.9, 37.2);
   AliMaterial(49, "ALUMINIUM$", 26.98, 13., 2.7, 8.9, 37.2);
   //
@@ -499,6 +138,11 @@ void AliABSO::CreateMaterials()
   AliMaterial(32, "TUNGSTEN$ ", 183.85, 74., 19.3, .35, 10.3);
   AliMaterial(52, "TUNGSTEN$ ", 183.85, 74., 19.3, .35, 10.3);
   //
+  //     Ni-W-Cu 
+  AliMixture(21, "Ni-W-Cu$", aniwcu, zniwcu, 18.78, 3, wniwcu);
+  AliMixture(41, "Ni-W-Cu$", aniwcu, zniwcu, 18.78, 3, wniwcu);
+  AliMixture(61, "Ni-W-Cu$", aniwcu, zniwcu, 18.78, 3, wniwcu);
+  //
   //     Lead 
   AliMaterial(13, "LEAD$     ", 207.19, 82., 11.35, .56, 18.5);
   AliMaterial(33, "LEAD$     ", 207.19, 82., 11.35, .56, 18.5);
@@ -534,24 +178,30 @@ void AliABSO::CreateMaterials()
   AliMixture(19, "STAINLESS STEEL$", asteel, zsteel, 7.88, 4, wsteel);
   AliMixture(39, "STAINLESS STEEL$", asteel, zsteel, 7.88, 4, wsteel);
   AliMixture(59, "STAINLESS STEEL$", asteel, zsteel, 7.88, 4, wsteel);
+  //
+  //     Insulation powder 
+  AliMixture(14, "INSULATION$", ains, zins, 0.41, -2, wins);
+  AliMixture(34, "INSULATION$", ains, zins, 0.41, 2, wins);
+  AliMixture(54, "INSULATION$", ains, zins, 0.41, 2, wins);
+
   //
   // **************** 
   //     Defines tracking media parameters. 
   //
-  epsil  = .001;  // Tracking precision, 
-  stemax = -1.;   // Maximum displacement for multiple scat 
-  tmaxfd = -20.;  // Maximum angle due to field deflection 
-  deemax = -.3;   // Maximum fractional energy loss, DLS 
+  epsil  = .001;    // Tracking precision, 
+  stemax = -0.01;   // Maximum displacement for multiple scat 
+  tmaxfd = -20.;    // Maximum angle due to field deflection 
+  deemax = -.3;     // Maximum fractional energy loss, DLS 
   stmin  = -.8;
   // *************** 
   //
   //    Carbon 
-  AliMedium(6, "C_C0             ", 6, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(6,  "C_C0             ",  6, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
   AliMedium(26, "C_C1             ", 26, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
   AliMedium(46, "C_C2             ", 46, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
   //
   //    Aluminum 
-  AliMedium(9, "ALU_C0          ", 9, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(9,  "ALU_C0          ",  9, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
   AliMedium(29, "ALU_C1          ", 29, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
   AliMedium(49, "ALU_C2          ", 49, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
   //
@@ -562,14 +212,23 @@ void AliABSO::CreateMaterials()
   //
   //    Tungsten 
   AliMedium(12, "W_C0            ", 12, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(32, "W_C1            ", 32, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(32, "W_C1            ", 32, 0, ISXFLD, SXMGMX, tmaxfd, -stemax, deemax, epsil, stmin);
   AliMedium(52, "W_C2            ", 52, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
+  //    Ni/Tungsten 
+  AliMedium(21, "Ni/W0           ", 21, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(41, "Ni/W1           ", 41, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(61, "Ni/W3           ", 61, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
   //
   //    Lead 
   AliMedium(13, "PB_C0           ", 13, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(33, "PB_C1           ", 33, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(33, "PB_C1           ", 33, 0, ISXFLD, SXMGMX, tmaxfd, -stemax, deemax, epsil, stmin);
   AliMedium(53, "PB_C2           ", 53, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
   //
+  //    Insulation Powder 
+  AliMedium(14, "INS_C0          ", 14, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(34, "INS_C1          ", 34, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(54, "INS_C2          ", 54, 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);
@@ -586,7 +245,7 @@ void AliABSO::CreateMaterials()
   AliMedium(57, "CC_C2           ", 57, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
   //
   //    Polyethilene 
-  AliMedium(18, "CH2_C0 B        ", 18, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(18, "CH2_C0          ", 18, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
   AliMedium(38, "CH2_C1          ", 38, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
   AliMedium(58, "CH2_C2          ", 58, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
   //
index edde3c30d0e0df04a16e9680cf84d978e1ff0ee8..749b5fabac1f6d7da88d0f2c93412052f3e5a49f 100644 (file)
@@ -15,6 +15,9 @@
 
 /*
 $Log$
+Revision 1.6  1999/09/29 09:24:30  fca
+Introduction of the Copyright and cvs Log
+
 */
 
 ///////////////////////////////////////////////////////////////////////////////
@@ -63,340 +66,6 @@ void AliSHIL::CreateGeometry()
   //
   // Build muon shield geometry
   //
-  // Origin N.Van Eijndhoven
-  //
-  //Begin_Html
-  /*
-    <img src="picts/AliSHIL.gif">
-  */
-  //End_Html
-  //Begin_Html
-  /*
-    <img src="picts/AliSHILTree.gif">
-  */
-  //End_Html
-
-  Float_t r_pb, cpar[5], parm[12], tpar[3], zmag_out,
-    abs_c, abs_l, r_abs, pcpar[12];
-  Float_t dr, dz, 
-    zs, theta_open;
-  Float_t cz1, cz2;
-  Float_t dzcoch1, dzcoch2, acc_min, dzs, zcone_e, zmag_in,
-    z_close, z_l3, shl1, shl4;
-    
-  Int_t *idtmed = fIdtmed->GetArray()-1699;
-  
-  //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.;
-  
-  
-  //  the mother of all shields 
-  
-  parm[0] = 0.;
-  parm[1] = 360.;
-  parm[2] = 3.;
-  parm[3] = abs_l;
-  parm[4] = 0.;
-  parm[5] = abs_l * TMath::Tan(acc_min * kDegrad);
-  parm[6] = zcone_e;
-  parm[7] = 0.;
-  parm[8] = 30.;
-  parm[9] = 1900.;
-  parm[10] = 0.;
-  parm[11] = 30.;
-  gMC->Gsvolu("YMOT", "PCON", idtmed[1715], parm, 12);
-  
-  //       beam shield ouside absorber inside l3 
-  
-  shl1    = (z_l3 - abs_l) / 2.;
-  cpar[0] = shl1;
-  cpar[1] = 0.;
-  cpar[2] = abs_l * TMath::Tan(acc_min * kDegrad);
-  cpar[3] = 0.;
-  cpar[4] = z_l3 * TMath::Tan(acc_min * kDegrad);
-  gMC->Gsvolu("YMS1", "CONE", idtmed[1759], cpar, 5);
-  
-  //       OUTSIDE STEEL TUBE 
-  // Pb/W 
-  cpar[1] = cpar[2] - 4.;
-  cpar[3] = cpar[4] - 4.;
-  gMC->Gsvolu("YSH1", "CONE", idtmed[1718], cpar, 5);
-  gMC->Gspos("YSH1", 1, "YMS1", 0., 0., 0., 0, "ONLY");
-  
-  //      COCHES FOR CHAMBERS 
-  dzcoch1 = cz1 + 4. - (abs_l + shl1);
-  dzcoch2 = cz2 + 4. - (abs_l + shl1);
-  
-  //       AIR 
-  cpar[0] = 7.5;
-  dr      = cpar[0] * 2. * TMath::Tan(acc_min * kDegrad);
-  cpar[2] = (cz1 + 4. - cpar[0]) * TMath::Tan(acc_min * kDegrad);
-  cpar[1] = cpar[2] - 3.;
-  cpar[3] = cpar[1] + dr;
-  cpar[4] = cpar[2] + dr;
-  gMC->Gsvolu("YAC1", "CONE", idtmed[1714], cpar, 5);
-  gMC->Gspos("YAC1", 1, "YMS1", 0., 0., dzcoch1, 0, "ONLY");
-  cpar[2] = (cz2 + 4. - cpar[0]) * TMath::Tan(acc_min * kDegrad);
-  cpar[1] = cpar[2] - 3.;
-  cpar[3] = cpar[1] + dr;
-  cpar[4] = cpar[2] + dr;
-  gMC->Gsvolu("YAC2", "CONE", idtmed[1714], cpar, 5);
-  gMC->Gspos("YAC2", 1, "YMS1", 0., 0., dzcoch2, 0, "ONLY");
-  
-  //       STEEL 
-  cpar[0] = (cz1 + 4. + 20. - abs_l) / 2.;
-  dr      = cpar[0] * 2. * TMath::Tan(acc_min * kDegrad);
-  //      CPAR(3)=ABS_L*TAN(kDegrad*ACC_MIN)-3. 
-  cpar[2] = abs_l * TMath::Tan(acc_min * kDegrad) - 4.;
-  cpar[1] = cpar[2] - 3.;
-  cpar[3] = cpar[1] + dr;
-  cpar[4] = cpar[2] + dr;
-  dzcoch1 = abs_l + cpar[0] - (abs_l + shl1);
-  gMC->Gsvolu("YSC1", "CONE", idtmed[1718], cpar, 5);
-  gMC->Gspos("YSC1", 1, "YMS1", 0., 0., dzcoch1, 0, "ONLY");
-  
-  cpar[0] = 16.;
-  dr      = cpar[0] * 2. * TMath::Tan(acc_min * kDegrad);
-  //      CPAR(3)=(CZ2+4.-CPAR(1))*TAN(kDegrad*ACC_MIN)-3. 
-  cpar[2] = (cz2 + 4. - cpar[0]) * TMath::Tan(acc_min * kDegrad) - 4.;
-  cpar[1] = cpar[2] - 3.;
-  cpar[3] = cpar[1] + dr;
-  cpar[4] = cpar[2] + dr;
-  gMC->Gsvolu("YSC2", "CONE", idtmed[1718], cpar, 5);
-  gMC->Gspos("YSC2", 1, "YMS1", -4., 0., dzcoch2, 0, "ONLY");
-  
-  //       ... beam pipe 
-  cpar[0] = shl1;
-  cpar[1] = 0.;
-  cpar[2] = r_abs + (abs_l - abs_c) * TMath::Tan(theta_open * kDegrad);
-  cpar[3] = 0.;
-  cpar[4] = cpar[2] + shl1 * 2. * TMath::Tan(theta_open * kDegrad);
-  gMC->Gsvolu("YMB1", "CONE", idtmed[1755], cpar, 5);
-  
-  cpar[0]  = shl1;
-  cpar[2] += -.8;
-  cpar[1]  = cpar[2] - .2;
-  cpar[4] += -.8;
-  cpar[3]  = cpar[4] - .2;
-  gMC->Gsvolu("YBS1", "CONE", idtmed[1749], cpar, 5);
-  gMC->Gspos("YBS1", 1, "YMB1", 0., 0., 0., 0, "ONLY");
-  
-  
-  gMC->Gspos("YMB1", 1, "YMS1", 0., 0., 0., 0, "ONLY");
-  dz = shl1 + abs_l;
-  gMC->Gspos("YMS1", 1, "YMOT", 0., 0., dz, 0, "ONLY");
-  
-  //       BEAM SHIELD OUTSIDE L3 
-  
-  //     L3->DIPOLE 
-  cpar[0] = (zmag_in - z_l3) / 2.;
-  cpar[1] = 0.;
-  cpar[2] = z_l3 * TMath::Tan(acc_min * kDegrad);
-  cpar[3] = 0.;
-  cpar[4] = zmag_in * TMath::Tan(acc_min * kDegrad);
-  gMC->Gsvolu("YMS2", "CONE", idtmed[1759], cpar, 5);
-  
-  //       OUTSIDE STEEL TUBE 
-  
-  // Pb/W 
-  cpar[1] = cpar[2] - 4.;
-  cpar[3] = cpar[4] - 4.;
-  gMC->Gsvolu("YSH2", "CONE", idtmed[1718], cpar, 5);
-  gMC->Gspos("YSH2", 1, "YMS2", 0., 0., 0., 0, "ONLY");
-  
-  //       ... beam pipe 
-
-  zs      = z_l3 - abs_c;
-  cpar[1] = 0.;
-  cpar[2] = r_abs + zs * TMath::Tan(theta_open * kDegrad);
-  cpar[3] = 0.;
-  cpar[4] = cpar[2] + cpar[0] * 2. * TMath::Tan(theta_open * kDegrad);
-  gMC->Gsvolu("YMB2", "CONE", idtmed[1755], cpar, 5);
-  
-  cpar[2] += -.8;
-  cpar[1]  = cpar[2] - .2;
-  cpar[4] += -.8;
-  cpar[3]  = cpar[4] - .2;
-  gMC->Gsvolu("YBS2", "CONE", idtmed[1749], cpar, 5);
-  gMC->Gspos("YBS2", 1, "YMB2", 0., 0., 0., 0, "ONLY");
-  
-  
-  gMC->Gspos("YMB2", 1, "YMS2", 0., 0., 0., 0, "ONLY");
-  dz = cpar[0] + z_l3;
-  gMC->Gspos("YMS2", 1, "YMOT", 0., 0., dz, 0, "ONLY");
-  pcpar[0] = 0.;
-  pcpar[1] = 360.;
-  pcpar[2] = 3.;
-  pcpar[3] = 0.;
-  pcpar[4] = 0.;
-  pcpar[5] = zmag_in * TMath::Tan(acc_min * kDegrad);
-  pcpar[6] = zcone_e - zmag_in;
-  pcpar[7] = 0.;
-  pcpar[8] = 30.;
-  pcpar[9] = zmag_out - zmag_in;
-  pcpar[10] = 0.;
-  pcpar[11] = 30.;
-  gMC->Gsvolu("YMS3", "PCON", idtmed[1759], pcpar, 12);
-  
-  //       OUTSIDE STEEL TUBE 
-  
-  // Pb/W 
-  pcpar[4]  = pcpar[5] - 4.;
-  pcpar[7]  = pcpar[8] - 4.;
-  pcpar[10] = pcpar[11] - 4.;
-  gMC->Gsvolu("YSH3", "PCON", idtmed[1718], pcpar, 12);
-  gMC->Gspos("YSH3", 1, "YMS3", 0., 0., 0., 0, "MANY");
-  
-  //       ... beam pipe up to closing cone 
-  
-  zs      = zmag_in - abs_c;
-  z_close = 804.;
-  
-  cpar[0] = (z_close - zmag_in) / 2.;
-  cpar[1] = 0.;
-  cpar[2] = r_abs + zs * TMath::Tan(theta_open * kDegrad);
-  cpar[3] = 0.;
-  cpar[4] = cpar[2] + cpar[0] * 2. * TMath::Tan(theta_open * kDegrad);
-  gMC->Gsvolu("YMB3", "CONE", idtmed[1755], cpar, 5);
-  
-  cpar[2] += -.8;
-  cpar[1]  = cpar[2] - .2;
-  cpar[4] += -.8;
-  cpar[3]  = cpar[4] - .2;
-  gMC->Gsvolu("YBS3", "CONE", idtmed[1749], cpar, 5);
-  gMC->Gspos("YBS3", 1, "YMB3", 0., 0., 0.,      0, "ONLY");
-  gMC->Gspos("YMB3", 1, "YMS3", 0., 0., cpar[0], 0, "ONLY");
-  
-  //       .closing cone 
-  
-  dzs = cpar[0] * 2.;
-  zs  = z_close - abs_c;
-  
-  cpar[0] = 25.;
-  cpar[1] = 0.;
-  cpar[2] = r_abs + zs * TMath::Tan(theta_open * kDegrad);
-  cpar[3] = 0.;
-  cpar[4] = r_abs;
-  gMC->Gsvolu("YMB5", "CONE", idtmed[1755], cpar, 5);
-  
-  cpar[2] += -.8;
-  cpar[1]  = cpar[2] - .2;
-  cpar[4] += -.8;
-  cpar[3]  = cpar[4] - .2;
-  gMC->Gsvolu("YBS5", "CONE", idtmed[1749], cpar, 5);
-  gMC->Gspos("YBS5", 1, "YMB5", 0., 0., 0., 0, "ONLY");
-  
-  dzs += cpar[0];
-  gMC->Gspos("YMB5", 1, "YMS3", 0., 0., dzs, 0, "ONLY");
-  dzs += cpar[0];
-  
-  //       OUTSIDE PB-TUBE 
-  
-  tpar[0] = 30. - r_pb - 4.;
-  tpar[1] = 30. - r_pb;
-  tpar[2] = (zmag_out - z_close - 50.) / 2.;
-  dzs    += tpar[2];
-  gMC->Gsvolu("YNW1", "TUBE", idtmed[1752], tpar, 3);
-  gMC->Gspos("YNW1", 1, "YMS3", 0., 0., dzs, 0, "MANY");
-  
-  //       constant beam pipe up to end of magnet 
-  
-  tpar[0] = 0.;
-  tpar[1] = r_abs;
-  tpar[2] = (zmag_out - z_close - 50.) / 2.;
-  gMC->Gsvolu("YMB6", "TUBE", idtmed[1755], tpar, 3);
-  
-  tpar[1] = r_abs - .8;
-  tpar[0] = tpar[1] - .2;
-  gMC->Gsvolu("YBS6", "TUBE", idtmed[1749], tpar, 3);
-  gMC->Gspos("YBS6", 1, "YMB6", 0., 0., 0.,  0, "ONLY");
-  gMC->Gspos("YMB6", 1, "YMS3", 0., 0., dzs, 0, "ONLY");
-  
-  dz = zmag_in;
-  gMC->Gspos("YMS3", 1, "YMOT", 0., 0., dz,  0, "ONLY");
-  
-  //       DIPOLE-> 
-  
-  tpar[0] = 0.;
-  tpar[1] = 30.;
-  tpar[2] = (1900. - zmag_out) / 2.;
-  shl4 = tpar[2];
-  gMC->Gsvolu("YMS4", "TUBE", idtmed[1759], tpar, 3);
-  //      CALL GSVOLU('YMS4','TUBE',IDTMED(1752),TPAR,3,IL3) ! W 
-  
-  //       OUTSIDE STEEL TUBE 
-  
-  // Pb/W 
-  tpar[0] = tpar[1] - 4.;
-  gMC->Gsvolu("YSH4", "TUBE", idtmed[1718], tpar, 3);
-  gMC->Gspos("YSH4", 1, "YMS4", 0., 0., 0., 0, "MANY");
-  
-  //       OUTSIDE PB-TUBE 
-  
-  tpar[0] = 30. - r_pb - 4.;
-  tpar[1] = 30. - r_pb;
-  gMC->Gsvolu("YNW2", "TUBE", idtmed[1752], tpar, 3);
-  gMC->Gspos("YNW2", 1, "YMS4", 0., 0., 0., 0, "MANY");
-  
-  //       ... beam pipe 
-  
-  tpar[0] = 0.;
-  tpar[1] = r_abs;
-  tpar[2] = shl4;
-  gMC->Gsvolu("YMB4", "TUBE", idtmed[1755], tpar, 3);
-  
-  tpar[2] = shl4;
-  tpar[1] = r_abs - .8;
-  tpar[0] = tpar[1] - .2;
-  gMC->Gsvolu("YBS4", "TUBE", idtmed[1749], tpar, 3);
-  gMC->Gspos("YBS4", 1, "YMB4", 0., 0., 0., 0, "ONLY");
-  
-  
-  dz = zmag_out + shl4;
-  gMC->Gspos("YMB4", 1, "YMS4", 0., 0., 0., 0, "ONLY");
-  gMC->Gspos("YMS4", 1, "YMOT", 0., 0., dz, 0, "ONLY");
-  gMC->Gspos("YMOT", 1, "ALIC", 0., 0., 0., 0, "ONLY");
 }
 
 //_____________________________________________________________________________
@@ -405,18 +74,31 @@ void AliSHIL::CreateMaterials()
   //
   // Defines materials for the muon shield
   //
-  // Origin N. Van Eijndhoven
-  //
-  
   Int_t   ISXFLD = gAlice->Field()->Integ();
   Float_t SXMGMX = gAlice->Field()->Max();
-  
+// 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 };
+// PbW
   Float_t apbw[2]   = { 207.2,183.85 };
   Float_t zpbw[2]   = { 82.,74. };
   Float_t wpbw[2]   = { .5,.5 };
+// Concrete
+  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 };
+// Ni-Cu-W alloy
+  Float_t aniwcu[3] ={58.6934, 183.84, 63.546};
+  Float_t zniwcu[3] ={28., 74., 29};
+  Float_t wniwcu[3] ={0.015,0.95,0.035};
+//
+// Insulation powder
+  Float_t ains[3] ={28.0855, 15.9994};
+  Float_t zins[3] ={14.,      8.};
+  Float_t wins[3] ={1.,       2.};
   
   Float_t epsil, stmin, tmaxfd, deemax, stemax;
   
@@ -457,7 +139,7 @@ void AliSHIL::CreateMaterials()
   AliMaterial(36, "VACUUM$ ", 1e-16, 1e-16, 1e-16, 1e16, 1e16);
   AliMaterial(56, "VACUUM$ ", 1e-16, 1e-16, 1e-16, 1e16, 1e16);
   
-  //     stainless Steel 
+  //     Stainless Steel 
   AliMixture(19, "STAINLESS STEEL$", asteel, zsteel, 7.88, 4, wsteel);
   AliMixture(39, "STAINLESS STEEL$", asteel, zsteel, 7.88, 4, wsteel);
   AliMixture(59, "STAINLESS STEEL$", asteel, zsteel, 7.88, 4, wsteel);
@@ -466,6 +148,21 @@ void AliSHIL::CreateMaterials()
   AliMixture(20, "LEAD/TUNGSTEN$", apbw, zpbw, 15.325, 2, wpbw);
   AliMixture(40, "LEAD/TUNGSTEN$", apbw, zpbw, 15.325, 2, wpbw);
   AliMixture(60, "LEAD/TUNGSTEN$", apbw, zpbw, 15.325, 2, wpbw);
+
+  //     Ni-W-Cu 
+  AliMixture(21, "Ni-W-Cu$", aniwcu, zniwcu, 18.78, 3, wniwcu);
+  AliMixture(41, "Ni-W-Cu$", aniwcu, zniwcu, 18.78, 3, wniwcu);
+  AliMixture(61, "Ni-W-Cu$", aniwcu, zniwcu, 18.78, 3, wniwcu);
+
+  //     Concrete 
+  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);
+
+  //     Insulation powder 
+  AliMixture(14, "INSULATION$", ains, zins, 0.41, -2, wins);
+  AliMixture(34, "INSULATION$", ains, zins, 0.41, 2, wins);
+  AliMixture(54, "INSULATION$", ains, zins, 0.41, 2, wins);
   
   // **************** 
   //     Defines tracking media parameters. 
@@ -479,7 +176,7 @@ void AliSHIL::CreateMaterials()
   // *************** 
   
   //    Aluminum 
-  AliMedium(9, "ALU_C0          ", 9, 0,  ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(9,  "ALU_C0          ", 9, 0,  ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
   AliMedium(29, "ALU_C1          ", 29, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
   AliMedium(49, "ALU_C2          ", 49, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
   
@@ -497,6 +194,11 @@ void AliSHIL::CreateMaterials()
   AliMedium(13, "PB_C0           ", 13, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
   AliMedium(33, "PB_C1           ", 33, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
   AliMedium(53, "PB_C2           ", 53, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
+
+  //    Insulation Powder 
+  AliMedium(14, "INS_C0          ", 14, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(34, "INS_C1          ", 34, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(54, "INS_C2          ", 54, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
   
   //    Air 
   AliMedium(15, "AIR_C0          ", 15, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
@@ -517,59 +219,23 @@ void AliSHIL::CreateMaterials()
   AliMedium(20, "PB/W0           ", 20, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
   AliMedium(40, "PB/W1           ", 40, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
   AliMedium(60, "PB/W3           ", 60, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
+
+  //    Ni/Tungsten 
+  AliMedium(21, "Ni/W0           ", 21, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(41, "Ni/W1           ", 41, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(61, "Ni/W3           ", 61, 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);
+
 }
 
 //_____________________________________________________________________________
 void AliSHIL::DrawModule () 
 {
-  //
-  // Draw a shaded view of the muon shield
-  //
-
-  // Set everything unseen
-  gMC->Gsatt("*", "seen", -1);
-  // 
-  // Set ALIC mother transparent
-  gMC->Gsatt("ALIC","SEEN",0);
-  //
-  // Set the volumes visible
-  gMC->Gsatt("YMOT","seen",1);
-  gMC->Gsatt("YMS1","seen",1);
-  gMC->Gsatt("YSH1","seen",1);
-  gMC->Gsatt("YAC1","seen",1);
-  gMC->Gsatt("YAC2","seen",1);
-  gMC->Gsatt("YSC1","seen",1);
-  gMC->Gsatt("YSC2","seen",1);
-  gMC->Gsatt("YMB1","seen",1);
-  gMC->Gsatt("YBS1","seen",1);
-  gMC->Gsatt("YMS2","seen",1);
-  gMC->Gsatt("YSH2","seen",1);
-  gMC->Gsatt("YMB2","seen",1);
-  gMC->Gsatt("YBS2","seen",1);
-  gMC->Gsatt("YMS3","seen",1);
-  gMC->Gsatt("YSH3","seen",1);
-  gMC->Gsatt("YMB3","seen",1);
-  gMC->Gsatt("YBS3","seen",1);
-  gMC->Gsatt("YMB5","seen",1);
-  gMC->Gsatt("YBS5","seen",1);
-  gMC->Gsatt("YNW1","seen",1);
-  gMC->Gsatt("YMB6","seen",1);
-  gMC->Gsatt("YBS6","seen",1);
-  gMC->Gsatt("YMS4","seen",1);
-  gMC->Gsatt("YSH4","seen",1);
-  gMC->Gsatt("YNW2","seen",1);
-  gMC->Gsatt("YMB4","seen",1);
-  gMC->Gsatt("YBS4","seen",1);
-  //
-  gMC->Gdopt("hide", "on");
-  gMC->Gdopt("shad", "on");
-  gMC->Gsatt("*", "fill", 7);
-  gMC->SetClipBox(".");
-  gMC->SetClipBox("*", 0, 3000, -3000, 3000, -6000, 6000);
-  gMC->DefaultRange();
-  gMC->Gdraw("alic", 30, 30, 0, 26.5, 18, .03, .03);
-  gMC->Gdhead(1111, "Muon Shield");
-  gMC->Gdman(16, 6, "MAN");
 }
 
 //_____________________________________________________________________________
@@ -586,7 +252,7 @@ void AliSHIL::Init()
   for(i=0;i<35;i++) printf("*");
   printf("\n");
   //
-  // Here the ABSO initialisation code (if any!)
+  // Here the SHIL initialisation code (if any!)
   for(i=0;i<80;i++) printf("*");
   printf("\n");
 }