// //
///////////////////////////////////////////////////////////////////////////////
+#include <TGeoGlobalMagField.h>
+#include <TArrayI.h>
+
#include "AliABSO.h"
#include "AliRun.h"
#include "AliMagF.h"
//
// Default constructor
//
+ for (Int_t i = 0; i < 2; i++)
+ {
+ fNLayers[i] = 0;
+ for (Int_t j = 0; j < 15; j++)
+ {
+ fZLayers[i][j] = 0.;
+ fMLayers[i][j] = 0;
+ }
+ }
}
//_____________________________________________________________________________
//
// Standard constructor
//
- SetMarkerColor(7);
- SetMarkerStyle(2);
- SetMarkerSize(0.4);
+ for (Int_t i = 0; i < 2; i++)
+ {
+ fNLayers[i] = 0;
+ for (Int_t j = 0; j < 15; j++)
+ {
+ fZLayers[i][j] = 0.;
+ fMLayers[i][j] = 0;
+ }
+ }
}
//_____________________________________________________________________________
//
}
-//_____________________________________________________________________________
-void AliABSO::DrawModule() const
-{
- //
- // Draw a shaded view of the muon absorber
- //
-}
-
-
//_____________________________________________________________________________
void AliABSO::CreateMaterials()
{
//
// Define materials for muon absorber
//
- Int_t isxfld = gAlice->Field()->Integ();
- Float_t sxmgmx = gAlice->Field()->Max();
+ Int_t isxfld = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Integ();
+ Float_t sxmgmx = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Max();
//
// Air
//
//
// 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 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 wconc[10] = { .01,.001,.529107,.016,.002,.033872, .337021,.013,.044,.014 };
//
// Steel
//
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};
+//
// Poly Concrete
// H Li F C Al Si Ca Pb O
Float_t aPolyCc[9] = {1. , 6.941, 18.998, 12.01, 26.98, 28.086, 40.078, 207.2, 15.999};
Float_t epsil, stmin, tmaxfd, deemax, stemax;
//
// Carbon
- AliMaterial( 6, "CARBON0$ ", 12.01, 6., 1.75, 24.4, 49.9);
- AliMaterial(26, "CARBON1$ ", 12.01, 6., 1.75, 24.4, 49.9);
- AliMaterial(46, "CARBON2$ ", 12.01, 6., 1.75, 24.4, 49.9);
+ AliMaterial( 6, "CARBON0$ ", 12.01, 6., 1.75, 24.4, 49.9);
+ AliMaterial(26, "CARBON1$ ", 12.01, 6., 1.75, 24.4, 49.9);
+ AliMaterial(46, "CARBON2$ ", 12.01, 6., 1.75, 24.4, 49.9);
+ //
+ // Magnesium
+ AliMaterial( 7, "MAGNESIUM$ ", 24.31, 12., 1.74, 25.3, 46.0);
//
// Aluminum
AliMaterial(9, "ALUMINIUM0$", 26.98, 13., 2.7, 8.9, 37.2);
AliMaterial(49, "ALUMINIUM2$", 26.98, 13., 2.7, 8.9, 37.2);
//
// Iron
- 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);
+ 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(11, "COPPER0$ ", 63.55, 29., 8.96, 1.43, 15.1);
- AliMaterial(31, "COPPER1$ ", 63.55, 29., 8.96, 1.43, 15.1);
- AliMaterial(51, "COPPER2$ ", 63.55, 29., 8.96, 1.43, 15.1);
+ AliMaterial(11, "COPPER0$ ", 63.55, 29., 8.96, 1.43, 15.1);
+ AliMaterial(31, "COPPER1$ ", 63.55, 29., 8.96, 1.43, 15.1);
+ AliMaterial(51, "COPPER2$ ", 63.55, 29., 8.96, 1.43, 15.1);
//
// Tungsten
- AliMaterial(12, "TUNGSTEN0$", 183.85, 74., 19.3, .35, 10.3);
- AliMaterial(32, "TUNGSTEN1$", 183.85, 74., 19.3, .35, 10.3);
- AliMaterial(52, "TUNGSTEN2$", 183.85, 74., 19.3, .35, 10.3);
+ AliMaterial(12, "TUNGSTEN0$ ", 183.85, 74., 19.3, .35, 10.3);
+ AliMaterial(32, "TUNGSTEN1$ ", 183.85, 74., 19.3, .35, 10.3);
+ AliMaterial(52, "TUNGSTEN2$ ", 183.85, 74., 19.3, .35, 10.3);
//
// Ni-W-Cu
AliMixture(21, "Ni-W-Cu0$", aniwcu, zniwcu, 18.78, 3, wniwcu);
AliMixture(61, "Ni-W-Cu2$", aniwcu, zniwcu, 18.78, 3, wniwcu);
//
// Lead
- AliMaterial(13, "LEAD0$ ", 207.19, 82., 11.35, .56, 18.5);
- AliMaterial(33, "LEAD1$ ", 207.19, 82., 11.35, .56, 18.5);
- AliMaterial(53, "LEAD2$ ", 207.19, 82., 11.35, .56, 18.5);
+ AliMaterial(13, "LEAD0$ ", 207.19, 82., 11.35, .56, 18.5);
+ AliMaterial(33, "LEAD1$ ", 207.19, 82., 11.35, .56, 18.5);
+ AliMaterial(53, "LEAD2$ ", 207.19, 82., 11.35, .56, 18.5);
//
// Air
- 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, "AIR0$ ", aAir, zAir, dAir, 4, wAir);
+ AliMixture(35, "AIR1$ ", aAir, zAir, dAir, 4, wAir);
+ AliMixture(55, "AIR2$ ", aAir, zAir, dAir, 4, wAir);
//
// Vacuum
- AliMixture(16, "VACUUM0$", aAir, zAir, dAir1, 4, wAir);
- AliMixture(36, "VACUUM1$", aAir, zAir, dAir1, 4, wAir);
- AliMixture(56, "VACUUM2$", aAir, zAir, dAir1, 4, wAir);
+ AliMixture(16, "VACUUM0$ ", aAir, zAir, dAir1, 4, wAir);
+ AliMixture(36, "VACUUM1$ ", aAir, zAir, dAir1, 4, wAir);
+ AliMixture(56, "VACUUM2$ ", aAir, zAir, dAir1, 4, wAir);
//
// Concrete
AliMixture(17, "CONCRETE0$", aconc, zconc, 2.35, 10, wconc);
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);
//
+ // Magnesium
+ AliMedium(7, "MG_C0 ", 7, 0, isxfld, 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(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(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);