]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - STRUCT/AliDIPOv2.cxx
Details of RB24/2 and RB24/3 geometry for correct BCM simulation.
[u/mrichter/AliRoot.git] / STRUCT / AliDIPOv2.cxx
index 46c99f35c18a82aeab1cec2bd91707e4aed265b5..90ca5be99aa1f56b9498ff71332e6d4da41eb0d4 100644 (file)
@@ -57,9 +57,6 @@ AliDIPOv2::AliDIPOv2(const char *name, const char *title)
 {
   //
   // Standard constructor for the magnetic dipole version 2    
-   SetMarkerColor(7);
-   SetMarkerStyle(2);
-   SetMarkerSize(0.4);
 }
 
 void AliDIPOv2::CreateGeometry()
@@ -98,33 +95,28 @@ void AliDIPOv2::CreateSpectrometerDipole()
   //       DIPOLE MAGNET 
   const Float_t kZDipole = 975; 
 
-  tpar[0] = 0.; 
-  tpar[1] = 360.;
-  tpar[2] = 5.; 
+  tpar[ 0] = 0.; 
+  tpar[ 1] = 360.;
+  tpar[ 2] = 4.; 
   //
-  tpar[3] = -250.55+kZDipole;
-  tpar[4] = 30.5;
-  tpar[5] = 527.34; 
+  tpar[ 3] = -250.55 + kZDipole;
+  tpar[ 4] =  30.1;
+  tpar[ 5] = 527.34; 
   //
-  tpar[6] = -160.7+kZDipole;
-  tpar[7] = 30.5;
-  tpar[8] = 527.34; 
+  tpar[ 6] = 37. + kZDipole;
+  tpar[ 7] =  30.1;
+  tpar[ 8] = 527.34;
   //
-  tpar[9] = 30.+kZDipole;
-  tpar[10] = 30.5;
+  tpar[ 9] = 37. + kZDipole;
+  tpar[10] = tpar[ 9] * TMath::Tan(2. * TMath::Pi() / 180.);
   tpar[11] = 527.34;
-
-  tpar[12] = 150.8+kZDipole;
-  tpar[13] = 182.4;
+  //
+  tpar[12] = 260.55 + kZDipole;
+  tpar[13] = tpar[12] * TMath::Tan(2. * TMath::Pi() / 180.);
   tpar[14] = 527.34;
+  gMC->Gsvolu("DDIP", "PCON", idtmed[1874], tpar, 15);
   //
-  tpar[15] = 260.55+kZDipole;
-  tpar[16] = 182.4;
-  tpar[17] = 527.34;
-
-
-  gMC->Gsvolu("DDIP", "PCON", idtmed[1814], tpar, 18);  
-  //       COILS 
+  //       Coils 
   // air - m.f. 
   cpar[0] = 207.;
   cpar[1] = 274.;
@@ -621,12 +613,7 @@ void AliDIPOv2::CreateSpectrometerDipole()
   
   AliMatrix(idrotm[1811], 270., 0., 90., 90., 180., 0.);
   gMC->Gspos("DDIP", 1, "ALIC", 0., 0., 0., idrotm[1811], "ONLY");
-
   gMC->Gsatt("DDIP", "SEEN", 0);
-//  gMC->Gsatt("DC21", "SEEN", 0);
-//  gMC->Gsatt("DC22", "SEEN", 0);
-//  gMC->Gsatt("DC3 ", "SEEN", 0);
-//  gMC->Gsatt("DC4 ", "SEEN", 0);
 }
 
 
@@ -634,70 +621,89 @@ void AliDIPOv2::CreateCompensatorDipole()
 {
     //
     //  Geometry of the compensator Dipole MBWMD (was MCB @ SPS)
-    // 
+    //  LAB I/EA Note 74.10
+    //  6/5/1974
+    //
+    const Float_t kHCoil       =  22.;  // Coil Height
+    const Float_t kWCoil       =  12.;  // Coil Width
+    const Float_t kLCoilH      = 250.;  // Hor. Coil Length
+    const Float_t kRCoilC      =  31.;  // Circ Coil Radius
+    const Float_t kWBase       = 125.;  // Base Width
+    const Float_t kHBase       =  30.;  // Base Height
+    
+    const Float_t kWUYoke      =  16.;
+    const Float_t kHUYoke      =  31.;
+
+    const Float_t kWLYoke      =  50.0;
+    const Float_t kHLYoke      =  61.0;
+    const Float_t kLLYoke      =  kLCoilH + kRCoilC;
+
+    const Float_t kWApperture  =  12.;
+    const Float_t kDCoil       =  kHUYoke + kHLYoke - 6. - 2. * kRCoilC;
+    
+    const Float_t kH           =  kHBase +  kHUYoke +  kHLYoke;
+    
     Int_t *idtmed = fIdtmed->GetArray()-1799;
     Int_t idrotm[1899];
 //
-    Float_t pbox[3] = {63., 63., 170.};
-    
+    Float_t pbox[3];
+    Float_t ptube[3];
+    ptube[0] = 4.;
+    ptube[1] = 65. * TMath::Sqrt(2.);
+    ptube[2] = (kLCoilH + 2. * kRCoilC) / 2.;
 //  Mother volumes
-    gMC->Gsvolu("DCM0", "BOX", idtmed[1814], pbox, 3);
-
+    gMC->Gsvolu("DCM0", "TUBE", idtmed[1814 + 40], ptube, 3);
+    gMC->Gsatt ("DCM0", "SEEN", 0);
+    
 //
 //  Mother volume containing lower coil
-    pbox[0] = 58.5/2.;
-    pbox[1] = 30.0;
-    pbox[2] = 152.5;
+    pbox[0] = kWLYoke / 2.;
+    pbox[1] = kHLYoke / 2.;
+    pbox[2] = kLLYoke / 2.;
     
-    gMC->Gsvolu("DCML", "BOX", idtmed[1809], pbox, 3);
+    gMC->Gsvolu("DCML", "BOX", idtmed[1809 + 40], pbox, 3);
 //
 // Base
-    pbox[0] = 62.5;
-    pbox[1] = 15.0;
-    gMC->Gsvolu("DCBA", "BOX", idtmed[1809], pbox, 3);
+    pbox[0] = kWBase / 2.;
+    pbox[1] = kHBase / 2.;
+    gMC->Gsvolu("DCBA", "BOX", idtmed[1809 + 40], pbox, 3);
 //
 // Coil: straight sections, horizontal
-    pbox[0] =   6.;
-    pbox[1] =  11.;
-    pbox[2] = 135.;
-    gMC->Gsvolu("DCH1", "BOX", idtmed[1816], pbox, 3);
+    pbox[0] = kWCoil  / 2.;
+    pbox[1] = kHCoil  / 2.;
+    pbox[2] = kLCoilH / 2.;
+    gMC->Gsvolu("DCH1", "BOX", idtmed[1816 + 40], pbox, 3);
 //
 // Coil: straight sections, horizontal
-    pbox[0] =   6.;
-    pbox[1] =  11.;
-    pbox[2] = 135.;
-    gMC->Gsvolu("DCH2", "BOX", idtmed[1816], pbox, 3);
+    pbox[0] = kWCoil  / 2.;
+    pbox[1] = kHCoil  / 2.;
+    pbox[2] = kLCoilH / 2.;
+    gMC->Gsvolu("DCH2", "BOX", idtmed[1816 + 40], pbox, 3);
 
 //
 // Mother volume containing upper coil
-    pbox[0] =    8.0;
-    pbox[1] =   17.5;
-    pbox[2] =  135.0;
-    gMC->Gsvolu("DCMU", "BOX", idtmed[1809], pbox, 3);
+    pbox[0] =  kWUYoke / 2.;
+    pbox[1] =  kHUYoke / 2.;
+    pbox[2] =  kLCoilH / 2.;
+    gMC->Gsvolu("DCMU", "BOX", idtmed[1809 + 40], pbox, 3);
 
 //
 // Coil: straight sections, vertical
-    pbox[0] =  6.0;
-    pbox[1] =  9.5;
-    pbox[2] = 11.0;
+    pbox[0] = kWCoil / 2.;
+    pbox[1] = kDCoil / 2.;
+    pbox[2] = kHCoil / 2.;
     
-    gMC->Gsvolu("DCCV", "BOX", idtmed[1816], pbox, 3);
+    gMC->Gsvolu("DCCV", "BOX", idtmed[1816 + 40], pbox, 3);
 //
 // Coil: circular section 
 
     Float_t ptubs[5];
-    ptubs[0] =  0.;
-    ptubs[1] = 35.;
-    ptubs[2] =  8.;
+    ptubs[0] = kRCoilC - kHCoil;
+    ptubs[1] = kRCoilC;
+    ptubs[2] = kWCoil / 2.;
     ptubs[3] =  0.;
     ptubs[4] = 90.;
-//    gMC->Gsvolu("DCC1", "TUBS", idtmed[1809], ptubs, 5);
-    ptubs[0] = 13.;
-    ptubs[1] = 35.;
-    ptubs[2] =  6.;
-    ptubs[3] =  0.;
-    ptubs[4] = 90.;
-    gMC->Gsvolu("DCC1", "TUBS", idtmed[1816], ptubs, 5);
+    gMC->Gsvolu("DCC1", "TUBS", idtmed[1816 + 40], ptubs, 5);
 //
 // Clamps
     Float_t ppgon[10];
@@ -711,7 +717,7 @@ void AliDIPOv2::CreateCompensatorDipole()
     ppgon[7] =  1.;
     ppgon[8] =  0.;
     ppgon[9] = 24.75;
-    gMC->Gsvolu("DCLA", "PGON", idtmed[1809], ppgon, 10);
+    gMC->Gsvolu("DCLA", "PGON", idtmed[1809 + 40], ppgon, 10);
 //
 // Assemble all
 //
@@ -720,48 +726,83 @@ void AliDIPOv2::CreateCompensatorDipole()
     AliMatrix(idrotm[1813], 180., 0., 90., 90.,  90., 0.);
     AliMatrix(idrotm[1814],   0., 180., 90., 270.,  90., 0.);
     AliMatrix(idrotm[1815], 180., 180., 90., 270.,  90., 0.);  
-       
-    gMC->Gspos("DCH1", 1, "DCML", 23.25, -13., -17.5, 0, "ONLY");
-    gMC->Gspos("DCCV", 1, "DCM0",  12., 19., -159., 0, "ONLY");
-    gMC->Gspos("DCCV", 2, "DCM0", -12., 19., -159., 0, "ONLY");
-    gMC->Gspos("DCCV", 3, "DCML", 23.25, 20.5, 141.5, 0, "ONLY");
-
-    gMC->Gspos("DCML", 1, "DCM0", -33.25, -2.5, 17.5, 0, "ONLY");
-    gMC->Gspos("DCML", 2, "DCM0",  33.25, -2.5, 17.5, idrotm[1811], "ONLY");
-
-
-    gMC->Gspos("DCH2", 1, "DCMU", 2., 6.5, 0., 0, "ONLY");
-    gMC->Gspos("DCMU", 1, "DCM0", -12., 45., 0., 0, "ONLY");
-    gMC->Gspos("DCMU", 2, "DCM0",  12., 45., 0., idrotm[1811], "ONLY");
-
-//    gMC->Gspos("DCC2", 1, "DCC1", 0., 0., 0., 0, "ONLY");
     
-    gMC->Gspos("DCC1", 1, "DCM0", -12., 27.5,  135., idrotm[1812], "ONLY");
-    gMC->Gspos("DCC1", 2, "DCM0",  12., 27.5,  135., idrotm[1812], "ONLY");
-    gMC->Gspos("DCC1", 3, "DCM0", -12., 27.5, -135., idrotm[1813], "ONLY");
-    gMC->Gspos("DCC1", 4, "DCM0",  12., 27.5, -135., idrotm[1813], "ONLY");
+    Float_t dx, dy, dz;
+    Float_t dy0 = 0.;
 
-    gMC->Gspos("DCC1", 5, "DCM0",  12., 27.5-32.+13., -135., idrotm[1815], "ONLY");
-    gMC->Gspos("DCC1", 6, "DCM0", -12., 27.5-32.+13., -135., idrotm[1815], "ONLY");
+    dy0 = -kH / 2. + kHBase/2.;
+    gMC->Gspos("DCBA", 1, "DCM0",  0., dy0, 15.0, 0, "ONLY");
+    
+    // Lower coil
+    dx = ( kWLYoke - kWCoil) / 2.;
+    dy = (-kHLYoke + kHCoil) / 2. + 6.;
+    gMC->Gspos("DCH1", 1, "DCML",  dx, dy,  -kRCoilC / 2., 0, "ONLY");
+    // Lower mother volume
+    dx = (kWLYoke + kWApperture) / 2.;
+    dy0 += (kHBase +  kHLYoke) / 2.;
+    gMC->Gspos("DCML", 1, "DCM0", -dx, dy0, kRCoilC / 2., 0, "ONLY");
+    gMC->Gspos("DCML", 2, "DCM0", +dx, dy0, kRCoilC / 2., idrotm[1811], "ONLY");
+    
+    dx = (kWUYoke - kWCoil) / 2.;
+    dy = (kHUYoke - kHCoil) / 2;
+    // Upper coil
+    gMC->Gspos("DCH2", 1, "DCMU",   dx,  dy, 0., 0, "ONLY");
+    // Upper mother volume
+    dx = (kWUYoke + kWApperture) / 2.;
+    dy0 +=  (kHLYoke + kHUYoke) / 2.;
+    gMC->Gspos("DCMU", 1, "DCM0", -dx, dy0, 0., 0, "ONLY");
+    gMC->Gspos("DCMU", 2, "DCM0", +dx, dy0, 0., idrotm[1811], "ONLY");
+
+    // Vertical coils
+    dx =  (kWCoil +  kWApperture) / 2.;
+    dy =  kH / 2. - kDCoil / 2. - kRCoilC;
+    dz =  (kLCoilH - kHCoil) / 2. + kRCoilC;
+    gMC->Gspos("DCCV", 1, "DCM0",  dx,  dy, -dz, 0, "ONLY");
+    gMC->Gspos("DCCV", 2, "DCM0", -dx,  dy, -dz, 0, "ONLY");
+
+    dx = (kWLYoke - kWCoil) / 2.;
+    dy = -kHLYoke / 2. + kDCoil / 2. + 6. + kRCoilC;
+    dz =  kLLYoke / 2. - kHCoil / 2.;
+    
+    gMC->Gspos("DCCV", 3, "DCML", dx, dy,  dz, 0, "ONLY");
+
+
+
+    // Circular coil
+    dx =  (kWCoil +  kWApperture) / 2.;
+    dy = dy0 + kHUYoke / 2. - kRCoilC;
+    dz =  kLCoilH / 2.;
+    gMC->Gspos("DCC1", 1, "DCM0", -dx, dy,  dz, idrotm[1812], "ONLY");
+    gMC->Gspos("DCC1", 2, "DCM0", +dx, dy,  dz, idrotm[1812], "ONLY");
+    gMC->Gspos("DCC1", 3, "DCM0", +dx, dy, -dz, idrotm[1813], "ONLY");
+    gMC->Gspos("DCC1", 4, "DCM0", -dx, dy, -dz, idrotm[1813], "ONLY");
+    dy = -kH / 2. + kHBase + 6. + kRCoilC;
+    gMC->Gspos("DCC1", 5, "DCM0", +dx, dy, -dz, idrotm[1815], "ONLY");
+    gMC->Gspos("DCC1", 6, "DCM0", -dx, dy, -dz, idrotm[1815], "ONLY");
+
+    dx = ( kWLYoke - kWCoil) / 2.;
+    dy =  -kHLYoke / 2. + 6. + kRCoilC;
+    dz =  kLLYoke / 2. - kRCoilC;
+    gMC->Gspos("DCC1", 7, "DCML", dx, dy, dz, idrotm[1814], "ONLY");
+
+//  Clamps
+    dx = kWApperture / 2. + kWUYoke;
+    dy = -kH / 2. + kHLYoke + kHBase;
     
-    gMC->Gspos("DCC1", 7, "DCML", 23.25, -13+13.+11., 117.5, idrotm[1814], "ONLY");
 
-    gMC->Gspos("DCLA", 1, "DCM0",  20., 27.5, -134., 0, "ONLY");
-    gMC->Gspos("DCLA", 2, "DCM0",  20., 27.5,  -44., 0, "ONLY");
-    gMC->Gspos("DCLA", 3, "DCM0",  20., 27.5,   46., 0, "ONLY");
-    gMC->Gspos("DCLA", 4, "DCM0",  20., 27.5,  134., 0, "ONLY");
+    gMC->Gspos("DCLA", 1, "DCM0",  dx, dy, -119., 0, "ONLY");
+    gMC->Gspos("DCLA", 2, "DCM0",  dx, dy,  -44., 0, "ONLY");
+    gMC->Gspos("DCLA", 3, "DCM0",  dx, dy,   46., 0, "ONLY");
+    gMC->Gspos("DCLA", 4, "DCM0",  dx, dy,  119., 0, "ONLY");
 
-    gMC->Gspos("DCLA", 5, "DCM0",  -20., 27.5, -134., idrotm[1811], "ONLY");
-    gMC->Gspos("DCLA", 6, "DCM0",  -20., 27.5,  -44., idrotm[1811], "ONLY");
-    gMC->Gspos("DCLA", 7, "DCM0",  -20., 27.5,   46., idrotm[1811], "ONLY");
-    gMC->Gspos("DCLA", 8, "DCM0",  -20., 27.5,  134., idrotm[1811], "ONLY");
+    gMC->Gspos("DCLA", 5, "DCM0",  -dx, dy, -119., idrotm[1811], "ONLY");
+    gMC->Gspos("DCLA", 6, "DCM0",  -dx, dy,  -44., idrotm[1811], "ONLY");
+    gMC->Gspos("DCLA", 7, "DCM0",  -dx, dy,   46., idrotm[1811], "ONLY");
+    gMC->Gspos("DCLA", 8, "DCM0",  -dx, dy,  119., idrotm[1811], "ONLY");
 
 
-    gMC->Gspos("DCBA", 1, "DCM0",  0., -47.5 , 17.5, 0, "ONLY");
     AliMatrix(idrotm[1816], 270., 0., 90., 90.,  180., 0.);  
-    gMC->Gspos("DCM0", 1, "ALIC",  0., -6.75,  975., idrotm[1816], "ONLY");
-
-
+    gMC->Gspos("DCM0", 1, "ALIC",  0., -12.,  1075., idrotm[1816], "ONLY");
 }
 
 //_____________________________________________________________________________
@@ -819,13 +860,9 @@ void AliDIPOv2::CreateMaterials()
   // Create Materials for Magnetic Dipole version 2
   //
   
-  Int_t isxfld   = gAlice->Field()->Integ();
-  Float_t sxmgmx = gAlice->Field()->Max();
-  
-  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;
+  Int_t isxfld1   = gAlice->Field()->Integ();
+  Int_t isxfld2   = gAlice->Field()->PrecInteg();
+  Float_t sxmgmx  = gAlice->Field()->Max();
   
   
   Float_t asteel[4] = { 55.847,51.9961,58.6934,28.0855 };
@@ -844,61 +881,61 @@ void AliDIPOv2::CreateMaterials()
   Float_t zG10[5] = { 1.,6.,8.,14.,35. };
   Float_t wG10[5] = { .02089,.22338,.28493,.41342,.05738 };
 
-  Float_t aAlCon[2] = { 14.61, 26.98};
-  Float_t zAlCon[2] = { 7.3, 13.};
-  Float_t wAlCon[2] = { .0004,.9996};
+    // 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;
+  Float_t dAir1 = 1.20479E-10;
 
   
   Float_t epsil, stmin, deemax, tmaxfd, stemax;
   
   // --- Define the various materials for GEANT --- 
   //     Aluminum 
-  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);
+  AliMaterial( 9, "ALUMINIUM0$", 26.98, 13., 2.7, 8.9, 37.2);
+  AliMaterial(29, "ALUMINIUM1$", 26.98, 13., 2.7, 8.9, 37.2);
+  AliMaterial(49, "ALUMINIUM2$", 26.98, 13., 2.7, 8.9, 37.2);
   
   //     Iron 
-  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(10, "IRON0$     ", 55.85, 26., 7.87, 1.76, 17.1);
+  AliMaterial(30, "IRON1$     ", 55.85, 26., 7.87, 1.76, 17.1);
+  AliMaterial(50, "IRON2$     ", 55.85, 26., 7.87, 1.76, 17.1);
   //     Copper
-  AliMaterial(17, "COPPER$   ", 63.55, 29., 8.96, 1.43, 15.1);
-  AliMaterial(37, "COPPER$   ", 63.55, 29., 8.96, 1.43, 15.1);
-  AliMaterial(57, "COPPER$   ", 63.55, 29., 8.96, 1.43, 15.1);
+  AliMaterial(17, "COPPER0$   ", 63.55, 29., 8.96, 1.43, 15.1);
+  AliMaterial(37, "COPPER1$   ", 63.55, 29., 8.96, 1.43, 15.1);
+  AliMaterial(57, "COPPER2$   ", 63.55, 29., 8.96, 1.43, 15.1);
   //     Air 
-  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(15, "AIR0$      ", aAir, zAir, dAir, 4, wAir);
+  AliMixture(35, "AIR1$      ", aAir, zAir, dAir, 4, wAir);
+  AliMixture(55, "AIR2$      ", aAir, zAir, dAir, 4, wAir);
+  AliMixture(75, "AIR_MUON   ", aAir, zAir, dAir, 4, wAir);
   //     Vacuum 
-  AliMaterial(16, "VACUUM$ ", 1e-16, 1e-16, 1e-16, 1e16, 1e16);
-  AliMaterial(36, "VACUUM$ ", 1e-16, 1e-16, 1e-16, 1e16, 1e16);
-  AliMaterial(56, "VACUUM$ ", 1e-16, 1e-16, 1e-16, 1e16, 1e16);
+  AliMixture(16, "VACUUM0$ ", aAir, zAir, dAir1, 4, wAir);
+  AliMixture(36, "VACUUM1$ ", aAir, zAir, dAir1, 4, wAir);
+  AliMixture(56, "VACUUM2$ ", aAir, zAir, dAir1, 4, wAir);
   
   //     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);
+  AliMixture(19, "STAINLESS STEEL0$", asteel, zsteel, 7.88, 4, wsteel);
+  AliMixture(39, "STAINLESS STEEL1$", asteel, zsteel, 7.88, 4, wsteel);
+  AliMixture(59, "STAINLESS STEEL2$", asteel, zsteel, 7.88, 4, wsteel);
   
   //     Coil 
-  AliMixture(14, "Al$", acoil, zcoil, 2.122, 3, wcoil);
-  AliMixture(34, "Al$", acoil, zcoil, 2.122, 3, wcoil);
-  AliMixture(54, "Al$", acoil, zcoil, 2.122, 3, wcoil);
+  AliMixture(14, "Al0$", acoil, zcoil, 2.122, 3, wcoil);
+  AliMixture(34, "Al1$", acoil, zcoil, 2.122, 3, wcoil);
+  AliMixture(54, "Al2$", acoil, zcoil, 2.122, 3, wcoil);
 
   //RESIN
-  AliMixture(13, "RESIN$", aresi, zresi, 1.05, 3, wresi);
-  AliMixture(33, "RESIN$", aresi, zresi, 1.05, 3, wresi);
-  AliMixture(53, "RESIN$", aresi, zresi, 1.05, 3, wresi);  
+  AliMixture(13, "RESIN0$", aresi, zresi, 1.05, 3, wresi);
+  AliMixture(33, "RESIN1$", aresi, zresi, 1.05, 3, wresi);
+  AliMixture(53, "RESIN2$", aresi, zresi, 1.05, 3, wresi);  
 
   //G10
-  AliMixture(11, "G10$", aG10, zG10, 1.7, 5, wG10);
-  AliMixture(31, "G10$", aG10, zG10, 1.7, 5, wG10);
-  AliMixture(51, "G10$", aG10, zG10, 1.7, 5, wG10); 
+  AliMixture(11, "G100$", aG10, zG10, 1.7, 5, wG10);
+  AliMixture(31, "G101$", aG10, zG10, 1.7, 5, wG10);
+  AliMixture(51, "G102$", aG10, zG10, 1.7, 5, wG10); 
  
-  //Aluminium Conductor
-  AliMixture(12, "AlCond$", aAlCon, zAlCon, 1.3506, 2, wAlCon);
-  AliMixture(32, "AlCond$", aAlCon, zAlCon, 1.3506, 2, wAlCon);
-  AliMixture(52, "AlCond$", aAlCon, zAlCon, 1.3506, 2, wAlCon);  
-
   // **************** 
   //     Defines tracking media parameters. 
   //     Les valeurs sont commentees pour laisser le defaut 
@@ -911,55 +948,50 @@ void AliDIPOv2::CreateMaterials()
   // *************** 
   
   //    Aluminum 
-  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);
+  AliMedium( 9, "ALU_C0          ",  9, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(29, "ALU_C1          ", 29, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(49, "ALU_C2          ", 49, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
   
   //    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, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(30, "FE_C1           ", 30, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(50, "FE_C2           ", 50, 0, isxfld1, 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, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(35, "AIR_C1          ", 35, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(55, "AIR_C2          ", 55, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(75, "AIR_MUON        ", 75, 0, isxfld2, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
   
   //    Vacuum 
-  AliMedium(16, "VA_C0           ", 16, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(36, "VA_C1           ", 36, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(56, "VA_C2           ", 56, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(16, "VA_C0           ", 16, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(36, "VA_C1           ", 36, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(56, "VA_C2           ", 56, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
   
   //    Steel 
-  AliMedium(19, "ST_C0           ", 19, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(39, "ST_C1           ", 39, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(59, "ST_C3           ", 59, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(19, "ST_C0           ", 19, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(39, "ST_C1           ", 39, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(59, "ST_C3           ", 59, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
   
   //    Coil 
-  AliMedium(14, "Coil_C1         ", 14, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(34, "Coil_C2         ", 34, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(54, "Coil_C3         ", 54, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(14, "Coil_C1         ", 14, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(34, "Coil_C2         ", 34, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(54, "Coil_C3         ", 54, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 
   //    Resin 
-  AliMedium(13, "RESIN_C0         ", 13, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(33, "RESIN_C1         ", 33, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(53, "RESIN_C2         ", 53, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(13, "RESIN_C0        ", 13, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(33, "RESIN_C1        ", 33, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(53, "RESIN_C2        ", 53, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 
   //    G10 
-  AliMedium(11, "G10_C0         ", 11, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(31, "G10_C1         ", 31, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(51, "G10_C2         ", 51, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
-
-  //Aluminium Contactor
-  AliMedium(12, "AlCond_C0         ", 12, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(32, "AlCond_C1         ", 32, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(52, "AlCond_C2         ", 52, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
-
+  AliMedium(11, "G10_C0          ", 11, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(31, "G10_C1          ", 31, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(51, "G10_C2          ", 51, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
   //
   //    Copper
-  AliMedium(17, "Cu_C0            ", 17, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(37, "Cu_C1            ", 37, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
-  AliMedium(57, "Cu_C2            ", 57, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(17, "Cu_C0           ", 17, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(37, "Cu_C1           ", 37, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
+  AliMedium(57, "Cu_C2           ", 57, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 
 }