[u/mrichter/AliRoot.git] / STRUCT / AliMAG.cxx

/**************************************************************************
 * 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.                  *
 **************************************************************************/

/* $Id$ */

///////////////////////////////////////////////////////////////////////////////
//                                                                           //
//  L3 Magnet                                                                //
//                                                                           //
//Begin_Html
/*
<img src="picts/AliMAGClass.gif">
</pre>
<br clear=left>
<font size=+2 color=red>
<p>The responsible person for this module is
<a href="mailto:andreas.morsch@cern.ch">Andreas Morsch</a>.
</font>
<pre>

*/
//End_Html
//                                                                           //
///////////////////////////////////////////////////////////////////////////////
 
#include <TVirtualMC.h>
#include <TGeoMedium.h>
#include <TGeoVolume.h>
#include <TGeoMatrix.h>
#include <TGeoPgon.h>
#include <TGeoXtru.h>
#include <TGeoCompositeShape.h>
#include <TGeoManager.h>

#include "AliMAG.h"
#include "AliMagF.h"
#include "AliRun.h"
 
ClassImp(AliMAG)
 
//_____________________________________________________________________________
AliMAG::AliMAG()
{
  //
  // Default constructor for L3 magnet
  //
}
 
//_____________________________________________________________________________
AliMAG::AliMAG(const char *name, const char *title)
  : AliModule(name,title)
{
  //
  // Standard constructor for L3 magnet
  //
  //Begin_Html
  /*
    <img src="picts/aliMAG.gif">
  */
  //End_Html
  
  SetMarkerColor(7);
  SetMarkerStyle(2);
  SetMarkerSize(0.4);
}

//_____________________________________________________________________________
void AliMAG::CreateGeometry()
{
  //
  // Create geometry for L3 magnet
  //
  //Begin_Html
  /*
    <img src="picts/mag.gif">
  */
  //End_Html
    
  //Begin_Html
  /*
    <img src="picts/tree_mag.gif">
    <br> Dimensions taken from drawing: ALIL3___00010
  //End_Html
  */
// Octagon
    const Int_t   kNSides              =    8;
    const Float_t kStartAngle          =   22.5; // deg
    const Float_t kFullAngle           =  360.0; // deg
//  Mother volume 
    const Float_t kRBMotherInner       = 560.00; // cm
    const Float_t kRBMotherOuter       = 790.50; // cm
    const Float_t kLBMother            = 706.00; // cm
// Yoke     
    const Float_t kRYokeInner          = 703.50; // cm
    const Float_t kRYokeOuter          = 790.50; // cm
    const Float_t kLYoke               = 620.00; // cm
// Coil
    const Float_t kRCoilInner          = 593.00; // cm
    const Float_t kRCoilOuter          = 682.00; // cm
    const Float_t kLCoil               = 588.00; // cm
// Cooling
    const Float_t kRCoolingOuter       = 1.70; // cm
    const Float_t kRCoolingInner       = 1.00; // cm
// Thermal Shield    
    const Float_t kRThermalShieldInner = 566.00; // cm
    const Float_t kRThermalShieldOuter = 571.00; // cm
// Crown    
    const Float_t kRCrownInner         = 560.00; // cm    
    const Float_t kRCrownOuter         = 785.50; // cm
    const Float_t kLCrown1             = 605.00; // cm
    const Float_t kLCrown2             = 620.00; // cm
    const Float_t kLCrown3             = 706.00; // cm
// Door
    const Float_t kRDoorOuter          = 560.00; // cm
    const Float_t kRPlugInner          = 183.50; // cm
    const Float_t kLDoor1              = 615.50; // cm
    const Float_t kLDoor2              = 714.60; // cm
//
    const Float_t kDegRad              = TMath::Pi()/180.;
    
    
  //
  // Top volume 
  TGeoVolume* top = gGeoManager->GetVolume("ALIC");
  // Media 
  TGeoMedium* medAir    = gGeoManager->GetMedium("MAG_AIR_C1");
  TGeoMedium* medAlu    = gGeoManager->GetMedium("MAG_ALU_C1");  
  TGeoMedium* medAluI   = gGeoManager->GetMedium("MAG_ALU_C0");
  TGeoMedium* medSteel  = gGeoManager->GetMedium("MAG_ST_C1");
  TGeoMedium* medWater  = gGeoManager->GetMedium("MAG_WATER");
  //
  // Offset between LHC and LEP axis
  Float_t os = -30.;

  //
  //  Define Barrel Mother 
  //  
  TGeoPgon* shBMother = new TGeoPgon(kStartAngle, kFullAngle, kNSides, 2);
  shBMother->DefineSection(0, -kLBMother, kRBMotherInner, kRBMotherOuter);
  shBMother->DefineSection(1,  kLBMother, kRBMotherInner, kRBMotherOuter);  
  // 
  TGeoVolume* voBMother = new TGeoVolume("L3BM", shBMother, medAir);
  //
  // Define Thermal Shield
  //
  // Only one layer
  // This can be improved: replace by (protection - shield - insulation) !
  //
  TGeoPgon* shThermSh = new TGeoPgon(kStartAngle, kFullAngle, kNSides, 2);
  shThermSh->DefineSection(0, -kLCoil, kRThermalShieldInner, kRThermalShieldOuter);
  shThermSh->DefineSection(1,  kLCoil, kRThermalShieldInner, kRThermalShieldOuter);  
  // 
  TGeoVolume* voThermSh = new TGeoVolume("L3TS", shThermSh, medAluI);
  voBMother->AddNode(voThermSh, 1, new TGeoTranslation(0., 0., 0.));
  //  
  // Define Coils and cooling circuits
  //
  TGeoPgon* shCoilMother = new TGeoPgon(kStartAngle, kFullAngle, kNSides, 2);
  shCoilMother->DefineSection(0, -kLCoil, kRCoilInner - 2. * kRCoolingOuter, kRCoilOuter + 2. * kRCoolingOuter);
  shCoilMother->DefineSection(1,  kLCoil, kRCoilInner - 2. * kRCoolingOuter, kRCoilOuter + 2. * kRCoolingOuter);  
  // 
  // Coils
  TGeoVolume* voCoilMother = new TGeoVolume("L3CM", shCoilMother, medAir);
  voBMother->AddNode(voCoilMother, 1, new TGeoTranslation(0., 0., 0.));
  // Devide into the 168 turns
  TGeoVolume* voCoilTurn   = voCoilMother->Divide("L3CD", 3, 168, 0., 0.);
  TGeoPgon* shCoils = new TGeoPgon(kStartAngle, kFullAngle, kNSides, 2);
  shCoils->DefineSection(0, -3., kRCoilInner, kRCoilOuter);
  shCoils->DefineSection(1,  3., kRCoilInner, kRCoilOuter); 
  // 
  TGeoVolume* voCoils = new TGeoVolume("L3C0", shCoils, medAlu);
  voCoilTurn->AddNode(voCoils, 1, new TGeoTranslation(0., 0., 0.));
  //
  // Hexagonal Cooling circuits
  // 
  const Float_t kRCC = kRCoolingOuter;
  const Float_t kRCW = kRCoolingInner;
  const Float_t kRCL = kRCC * TMath::Tan(30. / 180. * TMath::Pi());
  const Float_t kRWL = kRCW * TMath::Tan(30. / 180. * TMath::Pi());
  // Outer Circuits
  //
  // Pipe
  TGeoPgon* shCoolingPipeO = new TGeoPgon(kStartAngle, kFullAngle, kNSides, 4);
  shCoolingPipeO->DefineSection(0, -kRCC, kRCoilOuter + kRCC, kRCoilOuter + kRCC + 0.01);
  shCoolingPipeO->DefineSection(1, -kRCL, kRCoilOuter, kRCoilOuter + 2. * kRCC);
  shCoolingPipeO->DefineSection(2,  kRCL, kRCoilOuter, kRCoilOuter + 2. * kRCC);
  shCoolingPipeO->DefineSection(3,  kRCC, kRCoilOuter + kRCC, kRCoilOuter + kRCC + 0.01);
  //
  TGeoVolume* voCoolingPipeO = new TGeoVolume("L3CCO", shCoolingPipeO, medAlu);
  voCoilTurn->AddNode(voCoolingPipeO, 1, new TGeoTranslation(0., 0., 0.));
  //
  TGeoPgon* shCoolingWaterO = new TGeoPgon(kStartAngle, kFullAngle, kNSides, 4);
  shCoolingWaterO->DefineSection(0, -kRCW, kRCoilOuter + kRCC, kRCoilOuter + kRCC + 0.01);
  shCoolingWaterO->DefineSection(1, -kRWL, kRCoilOuter + (kRCC - kRCW), kRCoilOuter + kRCC + kRCW);
  shCoolingWaterO->DefineSection(2,  kRWL, kRCoilOuter + (kRCC - kRCW), kRCoilOuter + kRCC + kRCW);
  shCoolingWaterO->DefineSection(3,  kRCW, kRCoilOuter + kRCC, kRCoilOuter + kRCC + 0.01);
  //
  TGeoVolume* voCoolingWaterO = new TGeoVolume("L3CWO", shCoolingWaterO, medWater);
  voCoolingPipeO->AddNode(voCoolingWaterO, 1, new TGeoTranslation(0., 0., 0.));

  // Inner Circuits
  //
  // Pipe
  TGeoPgon* shCoolingPipeI = new TGeoPgon(kStartAngle, kFullAngle, kNSides, 4);
  shCoolingPipeI->DefineSection(0, -kRCC, kRCoilInner - kRCC, kRCoilInner - kRCC + 0.01);
  shCoolingPipeI->DefineSection(1, -kRCL, kRCoilInner - 2. * kRCC, kRCoilInner);
  shCoolingPipeI->DefineSection(2,  kRCL, kRCoilInner - 2. * kRCC, kRCoilInner);
  shCoolingPipeI->DefineSection(3,  kRCC, kRCoilInner - kRCC, kRCoilInner - kRCC + 0.01);
  //
  TGeoVolume* voCoolingPipeI = new TGeoVolume("L3CCI", shCoolingPipeI, medAlu);
  voCoilTurn->AddNode(voCoolingPipeI, 1, new TGeoTranslation(0., 0., 0.));
  //
  TGeoPgon* shCoolingWaterI = new TGeoPgon(kStartAngle, kFullAngle, kNSides, 4);
  shCoolingWaterI->DefineSection(0, -kRCW, kRCoilInner - kRCC, kRCoilInner - kRCC + 0.01);
  shCoolingWaterI->DefineSection(1, -kRWL, kRCoilInner - kRCC - kRCW, kRCoilInner - (kRCC - kRCW));
  shCoolingWaterI->DefineSection(2,  kRWL, kRCoilInner - kRCC - kRCW, kRCoilInner - (kRCC - kRCW));
  shCoolingWaterI->DefineSection(3,  kRCW, kRCoilInner - kRCC, kRCoilInner - kRCC + 0.01);
  //
  TGeoVolume* voCoolingWaterI = new TGeoVolume("L3CWI", shCoolingWaterI, medWater);
  voCoolingPipeI->AddNode(voCoolingWaterI, 1, new TGeoTranslation(0., 0., 0.));

  //
  // Define Yoke
  //
  TGeoPgon* shYoke = new TGeoPgon(kStartAngle, kFullAngle, kNSides, 2);
  shYoke->DefineSection(0, -kLYoke, kRYokeInner, kRYokeOuter);
  shYoke->DefineSection(1, +kLYoke, kRYokeInner, kRYokeOuter);  
  // 
  TGeoVolume* voYoke = new TGeoVolume("L3YO", shYoke, medSteel);
  voBMother->AddNode(voYoke, 1, new TGeoTranslation(0., 0., 0.));

  //
  // Define Crown
  //
  TGeoPgon* shCrown = new TGeoPgon(kStartAngle, kFullAngle, kNSides, 4);
  shCrown->DefineSection(0,  kLCrown1, kRCrownInner, kRYokeInner);
  shCrown->DefineSection(1,  kLCrown2, kRCrownInner, kRYokeInner);  
  shCrown->DefineSection(2,  kLCrown2, kRCrownInner, kRCrownOuter);  
  shCrown->DefineSection(3,  kLCrown3, kRCrownInner, kRCrownOuter);  
  // 
  TGeoVolume* voCrown = new TGeoVolume("L3CR", shCrown, medSteel);

  //
  // Door including "Plug"
  //
  Float_t slo = 2. * kRDoorOuter * TMath::Tan(22.5 * kDegRad);
  Float_t sli = 2. * kRPlugInner * TMath::Tan(22.5 * kDegRad);
  Double_t xpol1[12], xpol2[12], ypol1[12], ypol2[12];
  
  xpol1[ 0] = 2.         ;  ypol1[ 0] = kRDoorOuter;
  xpol1[ 1] = slo/2.     ;  ypol1[ 1] = kRDoorOuter;
  xpol1[ 2] = kRDoorOuter;  ypol1[ 2] = slo/2.;
  xpol1[ 3] = kRDoorOuter;  ypol1[ 3] = -slo/2.;
  xpol1[ 4] = slo/2.     ;  ypol1[ 4] = -kRDoorOuter;
  xpol1[ 5] = 2.         ;  ypol1[ 5] = -kRDoorOuter;
  xpol1[ 6] = 2.         ;  ypol1[ 6] = -kRPlugInner - os;
  xpol1[ 7] = sli/2.     ;  ypol1[ 7] = -kRPlugInner - os;
  xpol1[ 8] = kRPlugInner;  ypol1[ 8] = -sli/2.      - os;
  xpol1[ 9] = kRPlugInner;  ypol1[ 9] =  sli/2.      - os;
  xpol1[10] = sli/2.     ;  ypol1[10] = kRPlugInner  - os;
  xpol1[11] = 2.         ;  ypol1[11] = kRPlugInner  - os;

  TGeoXtru* shL3DoorR = new TGeoXtru(2);
  shL3DoorR->DefinePolygon(12, xpol1, ypol1);
  shL3DoorR->DefineSection(0, kLDoor1);
  shL3DoorR->DefineSection(1, kLDoor2);  
  TGeoVolume* voL3DoorR = new TGeoVolume("L3DoorR", shL3DoorR, medSteel);

  for (Int_t i = 0; i < 12; i++) {
      xpol2[i] = - xpol1[11 - i];
      ypol2[i] =   ypol1[11 - i];
  }
   
  TGeoXtru* shL3DoorL = new TGeoXtru(2);
  shL3DoorL->DefinePolygon(12, xpol2, ypol2);
  shL3DoorL->DefineSection(0, kLDoor1);
  shL3DoorL->DefineSection(1, kLDoor2);  
  TGeoVolume* voL3DoorL = new TGeoVolume("L3DoorL", shL3DoorL, medSteel);
  //
  // Plug support plate
  // 
  Float_t ro = kRPlugInner + 50.;
  slo = 2. * ro * TMath::Tan(22.5 * kDegRad);
  
  xpol1[ 0] = 2.         ;  ypol1[ 0] = ro     - os;
  xpol1[ 1] = slo/2.     ;  ypol1[ 1] = ro     - os;
  xpol1[ 2] = ro         ;  ypol1[ 2] = slo/2. - os;
  xpol1[ 3] = ro         ;  ypol1[ 3] = -slo/2.- os;
  xpol1[ 4] = slo/2.     ;  ypol1[ 4] = -ro    - os;
  xpol1[ 5] = 2.         ;  ypol1[ 5] = -ro    - os;

  for (Int_t i = 0; i < 12; i++) {
      xpol2[i] = - xpol1[11 - i];
      ypol2[i] =   ypol1[11 - i];
  }


  TGeoXtru* shL3PlugSPR = new TGeoXtru(2);
  shL3PlugSPR->DefinePolygon(12, xpol1, ypol1);
  shL3PlugSPR->DefineSection(0, kLDoor1-10.);
  shL3PlugSPR->DefineSection(1, kLDoor1);  
  TGeoVolume* voL3PlugSPR = new TGeoVolume("L3PlugSPR", shL3PlugSPR, medSteel);

  TGeoXtru* shL3PlugSPL = new TGeoXtru(2);
  shL3PlugSPL->DefinePolygon(12, xpol2, ypol2);
  shL3PlugSPL->DefineSection(0, kLDoor1-10.);
  shL3PlugSPL->DefineSection(1, kLDoor1);  
  TGeoVolume* voL3PlugSPL = new TGeoVolume("L3PlugSPL", shL3PlugSPL, medSteel);
  

  // Position crown and door
  TGeoRotation* rotxz = new TGeoRotation("rotxz",  90., 0., 90., 90., 180., 0.);

  TGeoVolumeAssembly *l3 = new TGeoVolumeAssembly("L3MO");
  voBMother->AddNode(voCrown, 1, new TGeoTranslation(0., 0., 0.));  
  voBMother->AddNode(voCrown, 2, new TGeoCombiTrans(0., 0., 0., rotxz));
  l3->AddNode(voBMother, 1, new TGeoTranslation(0.,0.,0.));
  l3->AddNode(voL3DoorR,  1, new TGeoTranslation(0., 0., 0.));  
  l3->AddNode(voL3DoorR,  2, new TGeoCombiTrans(0., 0., 0., rotxz));
  l3->AddNode(voL3DoorL,  1, new TGeoTranslation(0., 0., 0.));  
  l3->AddNode(voL3DoorL,  2, new TGeoCombiTrans(0., 0., 0., rotxz));
  l3->AddNode(voL3PlugSPR,  1, new TGeoTranslation(0., 0., 0.));  
  l3->AddNode(voL3PlugSPR,  2, new TGeoCombiTrans(0., 0., 0., rotxz));
  l3->AddNode(voL3PlugSPL,  1, new TGeoTranslation(0., 0., 0.));  
  l3->AddNode(voL3PlugSPL,  2, new TGeoCombiTrans(0., 0., 0., rotxz));
  top->AddNode(l3, 1, new TGeoTranslation(0., os, 0.));
}

//_____________________________________________________________________________
void AliMAG::CreateMaterials()
{
  //
  // Create materials for L3 magnet
  //
  
  Int_t   isxfld = gAlice->Field()->Integ();
  Float_t sxmgmx = gAlice->Field()->Max();
  Float_t epsil, stmin, deemax, tmaxfd, stemax;


  // --- Define the various materials for GEANT --- 
  // 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 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 aWater[2]={1.00794,15.9994};
  Float_t zWater[2]={1.,8.};
  Float_t wWater[2]={0.111894,0.888106};
  

  //     Aluminum 
  AliMaterial(9, "Al0$", 26.98, 13., 2.7, 8.9, 37.2);
  AliMaterial(29, "Al1$", 26.98, 13., 2.7, 8.9, 37.2);
  
  //     Stainless Steel 
  AliMixture(19, "STAINLESS STEEL1", asteel, zsteel, 7.88, 4, wsteel);
  AliMixture(39, "STAINLESS STEEL2", asteel, zsteel, 7.88, 4, wsteel);
  AliMixture(59, "STAINLESS STEEL3", asteel, zsteel, 7.88, 4, wsteel);
  //     Iron 
  AliMaterial(10, "Fe0$", 55.85, 26., 7.87, 1.76, 17.1);
  AliMaterial(30, "Fe1$", 55.85, 26., 7.87, 1.76, 17.1);
  
  //     Air 
  AliMixture(15, "AIR0$      ", aAir, zAir, dAir, 4, wAir);
  AliMixture(35, "AIR1$      ", aAir, zAir, dAir, 4, wAir);
  //     Water
  AliMixture(16, "WATER", aWater, zWater, 1., 2, wWater);

  
  // **************** 
  //     Defines tracking media parameters. 
  //     Les valeurs sont commentees pour laisser le defaut 
  //     a GEANT (version 3-21, page CONS200), f.m. 
  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 
  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);
  
  //     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);
  
  //     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);
  //    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);
  //    WATER
  AliMedium(16, "WATER             ", 16, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
}

//_____________________________________________________________________________
void AliMAG::DrawModule() const
{
  //
  // Draw a shaded view of the L3 magnet
  //
}

//_____________________________________________________________________________
void AliMAG::Init()
{
  //
  // Initialise L3 magnet after it has been built
  Int_t i;
  //
  if(AliLog::GetGlobalDebugLevel()>0) {
    printf("\n%s: ",ClassName());
    for(i=0;i<35;i++) printf("*");
    printf(" MAG_INIT ");
    for(i=0;i<35;i++) printf("*");
    printf("\n%s: ",ClassName());
    //
    // Here the MAG initialisation code (if any!)
    for(i=0;i<80;i++) printf("*");
    printf("\n");
  }
}
Commit	Line	Data
4c039060	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
eeacf08b	16	/* $Id$ */
4c039060	17
fe4da5cc	18	///////////////////////////////////////////////////////////////////////////////
	19	// //
	20	// L3 Magnet //
	21	// //
	22	//Begin_Html
	23	/*
1439f98e	24	<img src="picts/AliMAGClass.gif">
fe4da5cc	25	</pre>
	26	<br clear=left>
	27	<font size=+2 color=red>
	28	<p>The responsible person for this module is
	29	<a href="mailto:andreas.morsch@cern.ch">Andreas Morsch</a>.
	30	</font>
	31	<pre>
208c5bb6	32
fe4da5cc	33	*/
	34	//End_Html
	35	// //
	36	///////////////////////////////////////////////////////////////////////////////
	37
88cb7938	38	#include <TVirtualMC.h>
73d5f5aa	39	#include <TGeoMedium.h>
	40	#include <TGeoVolume.h>
	41	#include <TGeoMatrix.h>
	42	#include <TGeoPgon.h>
6b27f3a6	43	#include <TGeoXtru.h>
73d5f5aa	44	#include <TGeoCompositeShape.h>
73d5f5aa	45	#include <TGeoManager.h>
88cb7938	46
fe4da5cc	47	#include "AliMAG.h"
94de3818	48	#include "AliMagF.h"
88cb7938	49	#include "AliRun.h"
fe4da5cc	50
	51	ClassImp(AliMAG)
	52
	53	//_____________________________________________________________________________
	54	AliMAG::AliMAG()
	55	{
	56	//
	57	// Default constructor for L3 magnet
	58	//
	59	}
	60
	61	//_____________________________________________________________________________
	62	AliMAG::AliMAG(const char name, const char title)
b8032157	63	: AliModule(name,title)
fe4da5cc	64	{
	65	//
	66	// Standard constructor for L3 magnet
	67	//
	68	//Begin_Html
	69	/*
1439f98e	70	<img src="picts/aliMAG.gif">
fe4da5cc	71	*/
	72	//End_Html
	73
	74	SetMarkerColor(7);
	75	SetMarkerStyle(2);
	76	SetMarkerSize(0.4);
	77	}
	78
fe4da5cc	79	//_____________________________________________________________________________
	80	void AliMAG::CreateGeometry()
	81	{
	82	//
	83	// Create geometry for L3 magnet
	84	//
	85	//Begin_Html
	86	/*
1439f98e	87	<img src="picts/mag.gif">
fe4da5cc	88	*/
fe4da5cc	89	//End_Html
208c5bb6	90
fe4da5cc	91	//Begin_Html
fe4da5cc	92	/*
1439f98e	93	<img src="picts/tree_mag.gif">
208c5bb6	94	<br> Dimensions taken from drawing: ALIL3___00010
fe4da5cc	95	//End_Html
208c5bb6	96	*/
	97	// Octagon
	98	const Int_t kNSides = 8;
	99	const Float_t kStartAngle = 22.5; // deg
	100	const Float_t kFullAngle = 360.0; // deg
	101	// Mother volume
	102	const Float_t kRBMotherInner = 560.00; // cm
	103	const Float_t kRBMotherOuter = 790.50; // cm
	104	const Float_t kLBMother = 706.00; // cm
	105	// Yoke
	106	const Float_t kRYokeInner = 703.50; // cm
	107	const Float_t kRYokeOuter = 790.50; // cm
	108	const Float_t kLYoke = 620.00; // cm
	109	// Coil
	110	const Float_t kRCoilInner = 593.00; // cm
	111	const Float_t kRCoilOuter = 682.00; // cm
2cadb10d	112	const Float_t kLCoil = 588.00; // cm
	113	// Cooling
	114	const Float_t kRCoolingOuter = 1.70; // cm
	115	const Float_t kRCoolingInner = 1.00; // cm
208c5bb6	116	// Thermal Shield
	117	const Float_t kRThermalShieldInner = 566.00; // cm
	118	const Float_t kRThermalShieldOuter = 571.00; // cm
	119	// Crown
	120	const Float_t kRCrownInner = 560.00; // cm
	121	const Float_t kRCrownOuter = 785.50; // cm
	122	const Float_t kLCrown1 = 605.00; // cm
	123	const Float_t kLCrown2 = 620.00; // cm
	124	const Float_t kLCrown3 = 706.00; // cm
	125	// Door
208c5bb6	126	const Float_t kRDoorOuter = 560.00; // cm
6b27f3a6	127	const Float_t kRPlugInner = 183.50; // cm
208c5bb6	128	const Float_t kLDoor1 = 615.50; // cm
208c5bb6	129	const Float_t kLDoor2 = 714.60; // cm
6b27f3a6	130	//
6b27f3a6	131	const Float_t kDegRad = TMath::Pi()/180.;
208c5bb6	132
208c5bb6	133
73d5f5aa	134	//
	135	// Top volume
	136	TGeoVolume* top = gGeoManager->GetVolume("ALIC");
	137	// Media
2cadb10d	138	TGeoMedium* medAir = gGeoManager->GetMedium("MAG_AIR_C1");
	139	TGeoMedium* medAlu = gGeoManager->GetMedium("MAG_ALU_C1");
	140	TGeoMedium* medAluI = gGeoManager->GetMedium("MAG_ALU_C0");
6b27f3a6	141	TGeoMedium* medSteel = gGeoManager->GetMedium("MAG_ST_C1");
2cadb10d	142	TGeoMedium* medWater = gGeoManager->GetMedium("MAG_WATER");
73d5f5aa	143	//
	144	// Offset between LHC and LEP axis
	145	Float_t os = -30.;
	146
	147	//
208c5bb6	148	// Define Barrel Mother
73d5f5aa	149	//
208c5bb6	150	TGeoPgon* shBMother = new TGeoPgon(kStartAngle, kFullAngle, kNSides, 2);
	151	shBMother->DefineSection(0, -kLBMother, kRBMotherInner, kRBMotherOuter);
	152	shBMother->DefineSection(1, kLBMother, kRBMotherInner, kRBMotherOuter);
	153	//
	154	TGeoVolume* voBMother = new TGeoVolume("L3BM", shBMother, medAir);
	155	//
	156	// Define Thermal Shield
	157	//
2cadb10d	158	// Only one layer
	159	// This can be improved: replace by (protection - shield - insulation) !
	160	//
208c5bb6	161	TGeoPgon* shThermSh = new TGeoPgon(kStartAngle, kFullAngle, kNSides, 2);
	162	shThermSh->DefineSection(0, -kLCoil, kRThermalShieldInner, kRThermalShieldOuter);
	163	shThermSh->DefineSection(1, kLCoil, kRThermalShieldInner, kRThermalShieldOuter);
73d5f5aa	164	//
208c5bb6	165	TGeoVolume* voThermSh = new TGeoVolume("L3TS", shThermSh, medAluI);
208c5bb6	166	voBMother->AddNode(voThermSh, 1, new TGeoTranslation(0., 0., 0.));
73d5f5aa	167	//
2cadb10d	168	// Define Coils and cooling circuits
73d5f5aa	169	//
2cadb10d	170	TGeoPgon* shCoilMother = new TGeoPgon(kStartAngle, kFullAngle, kNSides, 2);
	171	shCoilMother->DefineSection(0, -kLCoil, kRCoilInner - 2. * kRCoolingOuter, kRCoilOuter + 2. * kRCoolingOuter);
	172	shCoilMother->DefineSection(1, kLCoil, kRCoilInner - 2. * kRCoolingOuter, kRCoilOuter + 2. * kRCoolingOuter);
	173	//
	174	// Coils
	175	TGeoVolume* voCoilMother = new TGeoVolume("L3CM", shCoilMother, medAir);
	176	voBMother->AddNode(voCoilMother, 1, new TGeoTranslation(0., 0., 0.));
	177	// Devide into the 168 turns
	178	TGeoVolume* voCoilTurn = voCoilMother->Divide("L3CD", 3, 168, 0., 0.);
208c5bb6	179	TGeoPgon* shCoils = new TGeoPgon(kStartAngle, kFullAngle, kNSides, 2);
2cadb10d	180	shCoils->DefineSection(0, -3., kRCoilInner, kRCoilOuter);
2cadb10d	181	shCoils->DefineSection(1, 3., kRCoilInner, kRCoilOuter);
73d5f5aa	182	//
73d5f5aa	183	TGeoVolume* voCoils = new TGeoVolume("L3C0", shCoils, medAlu);
2cadb10d	184	voCoilTurn->AddNode(voCoils, 1, new TGeoTranslation(0., 0., 0.));
	185	//
	186	// Hexagonal Cooling circuits
	187	//
	188	const Float_t kRCC = kRCoolingOuter;
	189	const Float_t kRCW = kRCoolingInner;
	190	const Float_t kRCL = kRCC * TMath::Tan(30. / 180. * TMath::Pi());
	191	const Float_t kRWL = kRCW * TMath::Tan(30. / 180. * TMath::Pi());
	192	// Outer Circuits
	193	//
	194	// Pipe
	195	TGeoPgon* shCoolingPipeO = new TGeoPgon(kStartAngle, kFullAngle, kNSides, 4);
	196	shCoolingPipeO->DefineSection(0, -kRCC, kRCoilOuter + kRCC, kRCoilOuter + kRCC + 0.01);
	197	shCoolingPipeO->DefineSection(1, -kRCL, kRCoilOuter, kRCoilOuter + 2. * kRCC);
	198	shCoolingPipeO->DefineSection(2, kRCL, kRCoilOuter, kRCoilOuter + 2. * kRCC);
	199	shCoolingPipeO->DefineSection(3, kRCC, kRCoilOuter + kRCC, kRCoilOuter + kRCC + 0.01);
	200	//
	201	TGeoVolume* voCoolingPipeO = new TGeoVolume("L3CCO", shCoolingPipeO, medAlu);
	202	voCoilTurn->AddNode(voCoolingPipeO, 1, new TGeoTranslation(0., 0., 0.));
	203	//
	204	TGeoPgon* shCoolingWaterO = new TGeoPgon(kStartAngle, kFullAngle, kNSides, 4);
	205	shCoolingWaterO->DefineSection(0, -kRCW, kRCoilOuter + kRCC, kRCoilOuter + kRCC + 0.01);
	206	shCoolingWaterO->DefineSection(1, -kRWL, kRCoilOuter + (kRCC - kRCW), kRCoilOuter + kRCC + kRCW);
	207	shCoolingWaterO->DefineSection(2, kRWL, kRCoilOuter + (kRCC - kRCW), kRCoilOuter + kRCC + kRCW);
	208	shCoolingWaterO->DefineSection(3, kRCW, kRCoilOuter + kRCC, kRCoilOuter + kRCC + 0.01);
	209	//
	210	TGeoVolume* voCoolingWaterO = new TGeoVolume("L3CWO", shCoolingWaterO, medWater);
	211	voCoolingPipeO->AddNode(voCoolingWaterO, 1, new TGeoTranslation(0., 0., 0.));
	212
	213	// Inner Circuits
	214	//
	215	// Pipe
	216	TGeoPgon* shCoolingPipeI = new TGeoPgon(kStartAngle, kFullAngle, kNSides, 4);
	217	shCoolingPipeI->DefineSection(0, -kRCC, kRCoilInner - kRCC, kRCoilInner - kRCC + 0.01);
	218	shCoolingPipeI->DefineSection(1, -kRCL, kRCoilInner - 2. * kRCC, kRCoilInner);
	219	shCoolingPipeI->DefineSection(2, kRCL, kRCoilInner - 2. * kRCC, kRCoilInner);
	220	shCoolingPipeI->DefineSection(3, kRCC, kRCoilInner - kRCC, kRCoilInner - kRCC + 0.01);
	221	//
	222	TGeoVolume* voCoolingPipeI = new TGeoVolume("L3CCI", shCoolingPipeI, medAlu);
	223	voCoilTurn->AddNode(voCoolingPipeI, 1, new TGeoTranslation(0., 0., 0.));
	224	//
	225	TGeoPgon* shCoolingWaterI = new TGeoPgon(kStartAngle, kFullAngle, kNSides, 4);
	226	shCoolingWaterI->DefineSection(0, -kRCW, kRCoilInner - kRCC, kRCoilInner - kRCC + 0.01);
	227	shCoolingWaterI->DefineSection(1, -kRWL, kRCoilInner - kRCC - kRCW, kRCoilInner - (kRCC - kRCW));
	228	shCoolingWaterI->DefineSection(2, kRWL, kRCoilInner - kRCC - kRCW, kRCoilInner - (kRCC - kRCW));
	229	shCoolingWaterI->DefineSection(3, kRCW, kRCoilInner - kRCC, kRCoilInner - kRCC + 0.01);
	230	//
	231	TGeoVolume* voCoolingWaterI = new TGeoVolume("L3CWI", shCoolingWaterI, medWater);
	232	voCoolingPipeI->AddNode(voCoolingWaterI, 1, new TGeoTranslation(0., 0., 0.));
	233
73d5f5aa	234	//
208c5bb6	235	// Define Yoke
73d5f5aa	236	//
208c5bb6	237	TGeoPgon* shYoke = new TGeoPgon(kStartAngle, kFullAngle, kNSides, 2);
	238	shYoke->DefineSection(0, -kLYoke, kRYokeInner, kRYokeOuter);
	239	shYoke->DefineSection(1, +kLYoke, kRYokeInner, kRYokeOuter);
73d5f5aa	240	//
6b27f3a6	241	TGeoVolume* voYoke = new TGeoVolume("L3YO", shYoke, medSteel);
208c5bb6	242	voBMother->AddNode(voYoke, 1, new TGeoTranslation(0., 0., 0.));
a00df67b	243
73d5f5aa	244	//
208c5bb6	245	// Define Crown
	246	//
	247	TGeoPgon* shCrown = new TGeoPgon(kStartAngle, kFullAngle, kNSides, 4);
	248	shCrown->DefineSection(0, kLCrown1, kRCrownInner, kRYokeInner);
	249	shCrown->DefineSection(1, kLCrown2, kRCrownInner, kRYokeInner);
	250	shCrown->DefineSection(2, kLCrown2, kRCrownInner, kRCrownOuter);
	251	shCrown->DefineSection(3, kLCrown3, kRCrownInner, kRCrownOuter);
	252	//
6b27f3a6	253	TGeoVolume* voCrown = new TGeoVolume("L3CR", shCrown, medSteel);
6b27f3a6	254
73d5f5aa	255	//
6b27f3a6	256	// Door including "Plug"
73d5f5aa	257	//
6b27f3a6	258	Float_t slo = 2. * kRDoorOuter * TMath::Tan(22.5 * kDegRad);
	259	Float_t sli = 2. * kRPlugInner * TMath::Tan(22.5 * kDegRad);
	260	Double_t xpol1[12], xpol2[12], ypol1[12], ypol2[12];
	261
	262	xpol1[ 0] = 2. ; ypol1[ 0] = kRDoorOuter;
	263	xpol1[ 1] = slo/2. ; ypol1[ 1] = kRDoorOuter;
	264	xpol1[ 2] = kRDoorOuter; ypol1[ 2] = slo/2.;
	265	xpol1[ 3] = kRDoorOuter; ypol1[ 3] = -slo/2.;
	266	xpol1[ 4] = slo/2. ; ypol1[ 4] = -kRDoorOuter;
	267	xpol1[ 5] = 2. ; ypol1[ 5] = -kRDoorOuter;
	268	xpol1[ 6] = 2. ; ypol1[ 6] = -kRPlugInner - os;
	269	xpol1[ 7] = sli/2. ; ypol1[ 7] = -kRPlugInner - os;
	270	xpol1[ 8] = kRPlugInner; ypol1[ 8] = -sli/2. - os;
	271	xpol1[ 9] = kRPlugInner; ypol1[ 9] = sli/2. - os;
	272	xpol1[10] = sli/2. ; ypol1[10] = kRPlugInner - os;
	273	xpol1[11] = 2. ; ypol1[11] = kRPlugInner - os;
	274
	275	TGeoXtru* shL3DoorR = new TGeoXtru(2);
	276	shL3DoorR->DefinePolygon(12, xpol1, ypol1);
	277	shL3DoorR->DefineSection(0, kLDoor1);
	278	shL3DoorR->DefineSection(1, kLDoor2);
	279	TGeoVolume* voL3DoorR = new TGeoVolume("L3DoorR", shL3DoorR, medSteel);
	280
	281	for (Int_t i = 0; i < 12; i++) {
	282	xpol2[i] = - xpol1[11 - i];
	283	ypol2[i] = ypol1[11 - i];
	284	}
	285
	286	TGeoXtru* shL3DoorL = new TGeoXtru(2);
	287	shL3DoorL->DefinePolygon(12, xpol2, ypol2);
	288	shL3DoorL->DefineSection(0, kLDoor1);
	289	shL3DoorL->DefineSection(1, kLDoor2);
	290	TGeoVolume* voL3DoorL = new TGeoVolume("L3DoorL", shL3DoorL, medSteel);
	291	//
	292	// Plug support plate
	293	//
	294	Float_t ro = kRPlugInner + 50.;
	295	slo = 2. * ro * TMath::Tan(22.5 * kDegRad);
	296
	297	xpol1[ 0] = 2. ; ypol1[ 0] = ro - os;
	298	xpol1[ 1] = slo/2. ; ypol1[ 1] = ro - os;
	299	xpol1[ 2] = ro ; ypol1[ 2] = slo/2. - os;
	300	xpol1[ 3] = ro ; ypol1[ 3] = -slo/2.- os;
	301	xpol1[ 4] = slo/2. ; ypol1[ 4] = -ro - os;
	302	xpol1[ 5] = 2. ; ypol1[ 5] = -ro - os;
	303
	304	for (Int_t i = 0; i < 12; i++) {
	305	xpol2[i] = - xpol1[11 - i];
	306	ypol2[i] = ypol1[11 - i];
	307	}
	308
	309
	310	TGeoXtru* shL3PlugSPR = new TGeoXtru(2);
	311	shL3PlugSPR->DefinePolygon(12, xpol1, ypol1);
	312	shL3PlugSPR->DefineSection(0, kLDoor1-10.);
	313	shL3PlugSPR->DefineSection(1, kLDoor1);
	314	TGeoVolume* voL3PlugSPR = new TGeoVolume("L3PlugSPR", shL3PlugSPR, medSteel);
	315
	316	TGeoXtru* shL3PlugSPL = new TGeoXtru(2);
	317	shL3PlugSPL->DefinePolygon(12, xpol2, ypol2);
	318	shL3PlugSPL->DefineSection(0, kLDoor1-10.);
	319	shL3PlugSPL->DefineSection(1, kLDoor1);
	320	TGeoVolume* voL3PlugSPL = new TGeoVolume("L3PlugSPL", shL3PlugSPL, medSteel);
	321
322
208c5bb6	323	// Position crown and door
73d5f5aa	324	TGeoRotation* rotxz = new TGeoRotation("rotxz", 90., 0., 90., 90., 180., 0.);
208c5bb6	325
	326	TGeoVolumeAssembly *l3 = new TGeoVolumeAssembly("L3MO");
	327	voBMother->AddNode(voCrown, 1, new TGeoTranslation(0., 0., 0.));
	328	voBMother->AddNode(voCrown, 2, new TGeoCombiTrans(0., 0., 0., rotxz));
	329	l3->AddNode(voBMother, 1, new TGeoTranslation(0.,0.,0.));
6b27f3a6	330	l3->AddNode(voL3DoorR, 1, new TGeoTranslation(0., 0., 0.));
	331	l3->AddNode(voL3DoorR, 2, new TGeoCombiTrans(0., 0., 0., rotxz));
	332	l3->AddNode(voL3DoorL, 1, new TGeoTranslation(0., 0., 0.));
	333	l3->AddNode(voL3DoorL, 2, new TGeoCombiTrans(0., 0., 0., rotxz));
	334	l3->AddNode(voL3PlugSPR, 1, new TGeoTranslation(0., 0., 0.));
	335	l3->AddNode(voL3PlugSPR, 2, new TGeoCombiTrans(0., 0., 0., rotxz));
	336	l3->AddNode(voL3PlugSPL, 1, new TGeoTranslation(0., 0., 0.));
	337	l3->AddNode(voL3PlugSPL, 2, new TGeoCombiTrans(0., 0., 0., rotxz));
73d5f5aa	338	top->AddNode(l3, 1, new TGeoTranslation(0., os, 0.));
fe4da5cc	339	}
	340
	341	//_____________________________________________________________________________
	342	void AliMAG::CreateMaterials()
	343	{
	344	//
	345	// Create materials for L3 magnet
	346	//
	347
b43eb0dc	348	Int_t isxfld = gAlice->Field()->Integ();
b43eb0dc	349	Float_t sxmgmx = gAlice->Field()->Max();
fe4da5cc	350	Float_t epsil, stmin, deemax, tmaxfd, stemax;
002b3738	351
002b3738	352
fe4da5cc	353	// --- Define the various materials for GEANT ---
6b27f3a6	354	// Steel
	355	Float_t asteel[4] = { 55.847,51.9961,58.6934,28.0855 };
	356	Float_t zsteel[4] = { 26.,24.,28.,14. };
	357	Float_t wsteel[4] = { .715,.18,.1,.005 };
002b3738	358	Float_t aAir[4]={12.0107,14.0067,15.9994,39.948};
	359	Float_t zAir[4]={6.,7.,8.,18.};
	360	Float_t wAir[4]={0.000124,0.755267,0.231781,0.012827};
	361	Float_t dAir = 1.20479E-3;
2cadb10d	362	Float_t aWater[2]={1.00794,15.9994};
	363	Float_t zWater[2]={1.,8.};
	364	Float_t wWater[2]={0.111894,0.888106};
	365
002b3738	366
fe4da5cc	367	// Aluminum
4a9de4af	368	AliMaterial(9, "Al0$", 26.98, 13., 2.7, 8.9, 37.2);
4a9de4af	369	AliMaterial(29, "Al1$", 26.98, 13., 2.7, 8.9, 37.2);
fe4da5cc	370
6b27f3a6	371	// Stainless Steel
	372	AliMixture(19, "STAINLESS STEEL1", asteel, zsteel, 7.88, 4, wsteel);
	373	AliMixture(39, "STAINLESS STEEL2", asteel, zsteel, 7.88, 4, wsteel);
	374	AliMixture(59, "STAINLESS STEEL3", asteel, zsteel, 7.88, 4, wsteel);
fe4da5cc	375	// Iron
4a9de4af	376	AliMaterial(10, "Fe0$", 55.85, 26., 7.87, 1.76, 17.1);
4a9de4af	377	AliMaterial(30, "Fe1$", 55.85, 26., 7.87, 1.76, 17.1);
fe4da5cc	378
fe4da5cc	379	// Air
4a9de4af	380	AliMixture(15, "AIR0$ ", aAir, zAir, dAir, 4, wAir);
4a9de4af	381	AliMixture(35, "AIR1$ ", aAir, zAir, dAir, 4, wAir);
2cadb10d	382	// Water
	383	AliMixture(16, "WATER", aWater, zWater, 1., 2, wWater);
	384
fe4da5cc	385
	386	// ****************
	387	// Defines tracking media parameters.
	388	// Les valeurs sont commentees pour laisser le defaut
	389	// a GEANT (version 3-21, page CONS200), f.m.
	390	epsil = .001; // Tracking precision,
	391	stemax = -1.; // Maximum displacement for multiple scat
	392	tmaxfd = -20.; // Maximum angle due to field deflection
	393	deemax = -.3; // Maximum fractional energy loss, DLS
	394	stmin = -.8;
	395	// ***************
	396
	397	// IRON
	398
b43eb0dc	399	AliMedium(10, "FE_C0 ", 10, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
b43eb0dc	400	AliMedium(30, "FE_C1 ", 30, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
fe4da5cc	401
	402	// ALUMINUM
	403
b43eb0dc	404	AliMedium(9, "ALU_C0 ", 9, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
b43eb0dc	405	AliMedium(29, "ALU_C1 ", 29, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
fe4da5cc	406
	407	// AIR
	408
b43eb0dc	409	AliMedium(15, "AIR_C0 ", 15, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
b43eb0dc	410	AliMedium(35, "AIR_C1 ", 35, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
6b27f3a6	411	// Steel
	412	AliMedium(19, "ST_C0 ", 19, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
	413	AliMedium(39, "ST_C1 ", 39, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
	414	AliMedium(59, "ST_C3 ", 59, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
2cadb10d	415	// WATER
2cadb10d	416	AliMedium(16, "WATER ", 16, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
fe4da5cc	417	}
	418
	419	//_____________________________________________________________________________
ba380e33	420	void AliMAG::DrawModule() const
fe4da5cc	421	{
	422	//
	423	// Draw a shaded view of the L3 magnet
	424	//
	425	}
	426
	427	//_____________________________________________________________________________
	428	void AliMAG::Init()
	429	{
	430	//
	431	// Initialise L3 magnet after it has been built
	432	Int_t i;
	433	//
4951e003	434	if(AliLog::GetGlobalDebugLevel()>0) {
9e1a0ddb	435	printf("\n%s: ",ClassName());
	436	for(i=0;i<35;i++) printf("*");
	437	printf(" MAG_INIT ");
	438	for(i=0;i<35;i++) printf("*");
	439	printf("\n%s: ",ClassName());
	440	//
	441	// Here the MAG initialisation code (if any!)
	442	for(i=0;i<80;i++) printf("*");
	443	printf("\n");
	444	}
fe4da5cc	445	}
fe4da5cc	446